Ion cyclotron resonance transmission line signal processing method and device, electronic equipment and medium

Abstract

The application discloses an ion cyclotron transmission line signal processing method and device, electronic equipment and medium, and relates to the technical field of ion cyclotron radio frequency systems and the like. The ion cyclotron transmission line signal processing method comprises the following steps: synchronously collecting a measurement point signal based on a voltage probe, correcting the measurement point signal based on a current calibration coefficient to obtain a corrected measurement point voltage signal; performing preliminary signal fitting through a forward wave and a reflected wave superposition model based on the measurement point voltage signal to obtain a measurement point fitting signal; calculating a measurement point difference index based on the measurement point fitting signal and the measurement point voltage signal, and generating a measurement point health score according to the measurement point difference index; determining a measurement point weight of each measurement point based on the measurement point health score, and performing standing wave reconstruction based on the measurement point weight to obtain a target standing wave voltage. According to the health score, the application reduces the weight or shields the measurement point, reduces the pollution of abnormal measurement points on the forward wave, the reflected wave and the standing wave distribution reconstruction result, and improves the accuracy of the target standing wave voltage.

Description

Technical Field

[0001] This application relates to the technical fields of ion cyclotron radio frequency systems, and in particular to a signal processing method, apparatus, electronic device and medium for an ion cyclotron transmission line. Background Technology

[0002] In an ion cyclotron radio frequency system, radio frequency energy is delivered to the antenna load via a transmitter, transmission line, and matching network. The standing wave voltage amplitude, antinode position, node position, and standing wave ratio on the transmission line directly affect matching status determination, insulation margin assessment, and overvoltage protection operation. In engineering practice, multiple voltage probes are typically arranged along the transmission line, and the standing wave voltage distribution is reconstructed by measuring the voltage amplitude, phase, or envelope signal at each measurement point.

[0003] Under conditions of high-power radio frequency, long pulse, strong electromagnetic interference and dynamic load changes, voltage probes and acquisition links may experience gain drift, phase deviation, zero drift, saturation, increased noise, loose connectors, channel abnormalities or temperature drift of the front-end conditioning circuit, which may cause the judgment of VSWR peak value, antinode position and VSWR to deviate from the true state.

[0004] Many related technologies rely on periodic offline calibration, single-point threshold alarms, or fixed-weight fusion, which makes it difficult to distinguish between real standing wave changes and measurement link anomalies in a timely manner, and also makes it difficult to maintain output reliability when some measurement points are abnormal. Summary of the Invention

[0005] Therefore, the purpose of this application is to propose a signal processing method, apparatus, electronic device and medium for ion cyclotron transmission lines, which performs health scoring on the actual acquired signals at each measurement point and dynamically adjusts the weight of each measurement point according to the health scoring results, so as to reduce or shield the measurement points, reduce the contamination of the forward wave, reflected wave and standing wave distribution reconstruction results by abnormal measurement points, and improve the accuracy of the target standing wave voltage.

[0006] This application provides a signal processing method for an ion cyclotron transmission line. The method includes: synchronously acquiring measurement point signals based on a voltage probe, and correcting the measurement point signals based on a current calibration coefficient to obtain a corrected measurement point voltage signal; performing preliminary signal fitting based on the measurement point voltage signal using a forward wave and reflected wave superposition model to obtain a measurement point fitted signal; calculating a measurement point difference index based on the measurement point fitted signal and the measurement point voltage signal, and generating a measurement point health score based on the measurement point difference index; determining the measurement point weight of each measurement point based on the measurement point health score, and performing standing wave reconstruction based on the measurement point weight to obtain a target standing wave voltage.

[0007] For example, the measurement point difference index includes at least two of the following: amplitude residual index, phase residual index, noise index, saturation index, drift index, and spatial consistency index. The step of generating a measurement point health score based on the measurement point difference index includes: weighting and normalizing at least two of the amplitude residual index, phase residual index, noise index, saturation index, drift index, and spatial consistency index based on index weighting coefficients to obtain the measurement point health score.

[0008] For example, determining the weight of each measuring point based on the measuring point health score includes: if the consecutive over-limit count of the current measuring point is greater than or equal to the consecutive over-limit number threshold, or if the measuring point health score of the current measuring point is less than the shielding threshold, determining the measuring point weight of the current measuring point to be an integer 0; otherwise, calculating the sum of the preset weight sensitivity coefficient powers of the measuring point health scores of each unshielded measuring point, and determining the ratio of the power of the measuring point health score of the current measuring point to the summation result as the measuring point weight of the current measuring point.

[0009] For example, the method further includes updating the current calibration coefficient when the reference signal is valid, the sampling saturation index is less than a preset threshold, the power change rate is less than a power threshold, the model residual is less than a residual threshold, and the health score of the measurement point is greater than a score threshold.

[0010] For example, the current calibration coefficient includes at least one of a recalibration coefficient, an amplitude calibration coefficient, a phase compensation amount, and an offset compensation amount. Updating the current calibration coefficient includes at least one of the following: constructing an instantaneous recalibration estimate based on the current window standing wave model fitting value, the measurement point voltage signal, and the current offset compensation amount; updating the amplitude calibration coefficient based on the amplitude of the instantaneous recalibration estimate; updating the phase compensation amount based on the phase of the instantaneous recalibration estimate; obtaining an update step size based on the product of the base step size and the measurement point health score; calculating a first product result of the update step size and the instantaneous recalibration estimate; calculating a second product result of the difference between a preset integer and the update step size and the current recalibration coefficient; updating the recalibration coefficient based on the amplitude limiting result obtained by summing the first and second product results; calculating a third product result of the offset update step size and the sample mean within a preset window; calculating a fourth product result of the difference between the preset integer and the offset update step size and the current offset compensation amount; and updating the offset compensation amount based on the sum of the third and fourth product results.

[0011] For example, the method further includes: obtaining the number of trusted measurement points, the number of model conditions, and the safety judgment voltage; and performing alarm protection based on the number of trusted measurement points, the number of model conditions, and the safety judgment voltage.

[0012] For example, the alarm protection based on the number of trusted measurement points, the number of model conditions, and the safety judgment voltage includes: triggering a data acquisition link health alarm when the number of trusted measurement points is less than a quantity threshold or the number of model conditions is greater than a condition number threshold; triggering an overvoltage alarm or interlocking protection when the number of trusted measurement points is greater than or equal to the quantity threshold and the number of model conditions is less than or equal to the condition number threshold, and the safety judgment voltage is greater than or equal to a protection threshold; and outputting the target standing wave voltage normally when the number of trusted measurement points is greater than or equal to the quantity threshold and the number of model conditions is less than or equal to the condition number threshold, and the safety judgment voltage is less than the protection threshold.

[0013] For example, the amplitude residual index is calculated based on the difference between the amplitude of the fitted signal at the measurement point and the amplitude of the voltage signal at the measurement point; the phase residual index is calculated based on the difference between the phase of the fitted signal at the measurement point and the phase of the voltage signal at the measurement point; the noise index is calculated based on the root mean square of the residual noise after removing the fundamental frequency; the saturation index is calculated based on the proportion of the number of points at the upper limit of sampling; the drift index is calculated based on the difference between the recalibration coefficients of the current window and the previous window; and the spatial consistency index is calculated based on the difference between the predicted value obtained by interpolating the voltages of adjacent measurement points and the amplitude of the voltage signal at the measurement point.

[0014] For example, before correcting the measurement point signal based on the current calibration coefficient, the method further includes performing at least one of complex phase extraction or envelope extraction on the measurement point signal.

[0015] Another embodiment of this application provides a signal processing device for an ion cyclotron transmission line. The device includes: a calibration module, used to synchronously acquire measurement point signals based on a voltage probe and correct the measurement point signals based on a current calibration coefficient to obtain a corrected measurement point voltage signal; a fitting module, used to perform preliminary signal fitting based on the measurement point voltage signal using a forward wave and reflected wave superposition model to obtain a measurement point fitted signal; a generation module, used to calculate a measurement point difference index based on the measurement point fitted signal and the measurement point voltage signal, and generate a measurement point health score based on the measurement point difference index; and a reconstruction module, used to determine the measurement point weight of each measurement point based on the measurement point health score, and perform standing wave reconstruction based on the measurement point weight to obtain a target standing wave voltage.

[0016] Another embodiment of this application provides an electronic device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of the above-described signal processing method for ion cyclotron transport lines.

[0017] Another embodiment of this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the signal processing method for the ion cyclotron transport line described above.

[0018] In the above embodiments, the signal processing method of the ion cyclotron transmission line includes: synchronously acquiring measurement point signals based on voltage probes, and correcting the measurement point signals based on the current calibration coefficients to obtain corrected measurement point voltage signals; performing preliminary signal fitting based on the measurement point voltage signals using a forward wave and reflected wave superposition model to obtain measurement point fitted signals; calculating measurement point difference indices based on the measurement point fitted signals and measurement point voltage signals, and generating measurement point health scores based on the measurement point difference indices; determining the measurement point weights of each measurement point based on the measurement point health scores, and performing standing wave reconstruction based on the measurement point weights to obtain the target standing wave voltage; performing health scores on the actual acquired signals of each measurement point, and dynamically adjusting the weights of each measurement point based on the health score results to achieve weight reduction or shielding of measurement points, reduce the contamination of forward wave, reflected wave and standing wave distribution reconstruction results by abnormal measurement points, and improve the accuracy of the target standing wave voltage. Attached Figure Description

[0019] Figure 1 A flowchart of a signal processing method for an ion cyclotron transmission line provided in an embodiment of this application; Figure 2 A schematic diagram of the structure of the ion cyclotron transmission line standing wave voltage acquisition and diagnosis system provided in this application embodiment; Figure 3 A schematic diagram of standing wave voltage reconstruction under multi-point voltage sampling provided for an embodiment of this application; Figure 4 A schematic diagram of a health scoring, anomaly shielding, and protection alarm state machine provided for the implementation of this application; Figure 5 A schematic diagram of a signal processing apparatus for an ion cyclotron transmission line provided in an embodiment of this application; Figure 6 A block diagram of an electronic device provided in an embodiment of this application. Detailed Implementation

[0020] The embodiments of this application are described in detail below. Examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.

[0021] Figure 1 This is a flowchart of a signal processing method for an ion cyclotron transmission line according to an embodiment of this application.

[0022] As an example, such as Figure 1 As shown, the signal processing method for ion cyclotron transmission lines includes: S101 synchronously acquires the measurement point signal based on the voltage probe, and corrects the measurement point signal based on the current calibration coefficient to obtain the corrected measurement point voltage signal.

[0023] S102, based on the voltage signal at the measurement point, performs preliminary signal fitting using a forward wave and reflected wave superposition model to obtain the fitted signal at the measurement point.

[0024] S103 calculates the measurement point difference index based on the measurement point fitting signal and the measurement point voltage signal, and generates the measurement point health score based on the measurement point difference index.

[0025] S104: Determine the weight of each measuring point based on the health score of the measuring point, and perform standing wave reconstruction based on the measuring point weight to obtain the target standing wave voltage.

[0026] For example, multiple voltage probes can be arranged on an ion cyclotron transmission line, with the probe positions designated as measurement points. Measurement point signals are synchronously acquired using these voltage probes. For instance, voltage probe signals for each measurement point can be acquired within a preset acquisition window. Of course, the measurement point signals also include the reference voltage, known calibration signals, and other signals. Considering the influence of gain, phase, bias, and channel delay on the voltage probes and the front-end acquisition link, the measurement point signals are corrected based on the current calibration coefficients. These calibration coefficients are used to eliminate the influence of gain, phase, bias, and channel delay on the voltage probes and the front-end acquisition link. Based on the corrected measurement point voltage signals, preliminary signal fitting is performed using a forward wave and reflected wave superposition model to obtain the measurement point fitted signal. It can be understood that, based on the forward wave and reflected wave superposition model, the forward wave and reflected wave parameters for the current window are inferred from the measured values ​​of multiple measurement points. Then, the forward wave and reflected wave parameters are substituted back into the transmission line model to obtain the theoretically corresponding signal at each measurement point location. This theoretical signal is the measurement point fitted signal. The measurement point difference index is calculated based on the fitted signal and voltage signal of the measurement point, and a measurement point health score is generated based on the measurement point difference index. In essence, the health level of a measurement point is measured by the difference between the theoretical voltage signal and the actual acquired voltage signal, thus generating a measurement point health score. Different measurement point health scores correspond to different measurement point weights. The measurement point weight represents the degree of participation of that measurement point. For example, if the anomaly of a measurement point is significant, it can be directly blocked. Alternatively, if the anomaly of a measurement point is slight, its corresponding measurement point weight can be reduced. Signal reconstruction is performed based on the measurement point weight of each measurement point, ultimately outputting the target standing wave voltage. The target standing wave voltage is not a simple voltage value, but rather the spatial distribution of the complex phasor of the standing wave voltage reconstructed along the transmission line position coordinates, including multiple features such as amplitude and phase.

[0027] The signal processing method for ion cyclotron transmission lines in this application introduces an online measurement point health scoring mechanism. By mapping the measurement point health score to the measurement point weight, normal or reliable measurement points occupy a higher weight in the standing wave reconstruction, while abnormal measurement points are automatically downweighted or shielded, reducing the contamination of the forward wave, reflected wave and standing wave distribution reconstruction results by abnormal measurement points.

[0028] Figure 2 This is a schematic diagram of the structure of an ion cyclotron transmission line standing wave voltage acquisition and diagnosis system according to an embodiment of this application.

[0029] like Figure 2 As shown, the system includes multiple voltage probes installed at different locations on the ion cyclotron transmission line. It is configured with a front-end conditioning circuit, a synchronous data acquisition module, a reference voltage verification module, a complex phasor / envelope extraction module, a standing wave voltage reconstruction module, an online calibration module, a health diagnostic module, a weighted fusion module, a protection alarm interface, and a log recording module. The position parameters of each voltage probe are entered into the system during installation or commissioning. The transmission line phase constant and characteristic impedance can be obtained from design parameters, offline calibration, or by looking up tables at the operating frequency. The front-end conditioning circuit is used to couple the high-voltage radio frequency signal into a low-voltage signal that the acquisition module can receive, and may include attenuation, isolation, filtering, detection, or synchronous phase acquisition units. The synchronous data acquisition module acquires the voltage probe signals, reference voltage, known verification signal, transmission line position parameters, and sampling channel status at each measurement point within a preset acquisition window, calibration window, stable operating range, or pulse gap. The reference voltage verification module can inject a known verification signal into the acquisition link during pulse gaps or calibration windows, or select historical calibration data from a verified reference measurement point or under a standard load as a reference. The complex phasor / envelope extraction module filters, synchronizes, compensates for delays, and extracts the fundamental frequency component from the voltage sampling signal. The standing wave voltage reconstruction module converts multi-point measurements into forward, reverse, and standing wave distributions based on the measurement point locations and the transmission line traveling wave superposition model, and mitigates the impact of abnormal measurement points on the reconstruction results under health weight constraints. The online calibration module estimates and updates the gain, phase, and bias compensation coefficients. The health diagnosis module generates a health score based on residuals, noise, saturation, drift, and spatial consistency. The weight fusion module maps the health score to reconstruction weights. The protection alarm interface outputs maintenance prompts, insufficient confidence alarms, or protection interlock signals based on the reliable standing wave voltage, measurement point reliability, model condition number, and safety margin.

[0030] As an example, before correcting the measurement point signal based on the current calibration coefficient, the signal processing method for the ion cyclotron transmission line also includes: performing at least one of complex phase extraction or envelope extraction on the measurement point signal.

[0031] For example, the extraction of complex phasors from the measurement point signal is used as an illustration. Within the acquisition window, the output signal of the i-th voltage probe is subjected to synchronous demodulation, digital down-conversion, discrete Fourier transform, or equivalent envelope extraction. Let the measurement point position be zi, and the sampled signal be... [n], sampling time is The radio frequency angular frequency is The voltage complex quantity is as follows:

[0032] in, Let N be the voltage complex phase of the i-th measurement point within the current acquisition window, and N be the number of sampling points in the window. Let be the radio frequency angular frequency, j be the imaginary unit, and e be the natural constant, where the exponent term e is a complex exponential function with base e.

[0033] If the acquisition link only provides the envelope signal, then take =| As a characteristic of amplitude, Let be the voltage amplitude at the i-th measurement point; if a synchronous phase is provided, then ∠ is retained. .

[0034] Of course, complex phasor extraction can be achieved using synchronous demodulation, phase-locked detection, digital downconversion, discrete Fourier transform, Hilbert envelope extraction, or other equivalent fundamental frequency or envelope extraction methods.

[0035] After extracting the complex phase signal from the measurement point, the measurement point signal is corrected based on the current calibration coefficient to obtain the corrected measurement point voltage signal.

[0036] As an example, considering the effects of gain, phase, bias, and channel delay on the voltage probe and the front-end acquisition link, a measurement point correction model is established:

[0037] in, The corrected measurement point voltage signal, also known as the measurement point voltage complex phasor, For recalibration coefficients, For amplitude calibration coefficient, This is the phase compensation amount. This is the offset compensation amount. The channel delay compensation amount is calculated as follows: The recalibration coefficient, amplitude calibration coefficient, phase compensation amount, and offset compensation amount are estimated and updated by the online calibration module. The channel delay compensation amount is calibrated during the commissioning phase using a synchronization reference pulse, a known verification signal, or the phase difference between channels, and can be used as a fixed compensation amount during operation. The amplitude calibration coefficient is determined by the magnitude of the recalibration coefficient, and the phase compensation amount is determined by the phase angle of the recalibration coefficient. In another embodiment, the amplitude component and the phase component can also be updated separately.

[0038] After correcting the measurement point signal, preliminary signal fitting is performed using a forward wave and reflected wave superposition model based on the corrected measurement point voltage signal to obtain the measurement point fitted signal.

[0039] As an example, to avoid the problem that single-point threshold judgment cannot reflect the voltage distribution of the entire transmission line, this application represents the transmission line voltage as the superposition of forward and reverse waves. Let the complex amplitude of the forward voltage wave be... The complex amplitude of the reflected wave is The transmission line phase constant is Then the voltage complex phasor at any position z is:

[0040]

[0041] in, Let be the measurement residual at the i-th measurement point. The voltage complex phasors of the above M measurement points can be written in matrix form as follows:

[0042] in, The observation vector is composed of the voltage complex phasors from multiple measurement points. The transmission line observation matrix is ​​constructed based on the measurement point location zi and the transmission line phase constant β. Let ε be the vector of complex amplitude parameters of the forward and reflected waves to be estimated, and let ε be the measurement residual vector. Substituting the voltage signal at the measurement point obtained from the actual acquired measurement point correction model into the superposition model of the forward and reflected waves, we obtain the estimated complex amplitude parameters of the forward and reflected waves. Then, based on the estimated complex amplitude parameters of the forward and reflected waves, we perform preliminary signal fitting to obtain the fitted signal at the measurement point.

[0043] The measurement point difference index is calculated based on the measurement point fitting signal and the measurement point voltage signal, and a measurement point health score is generated based on the measurement point difference index.

[0044] As an example, the measurement point difference index includes at least two of the following: amplitude residual index, phase residual index, noise index, saturation index, drift index, and spatial consistency index. A measurement point health score is generated based on the measurement point difference index, including: The health score of the measurement point is obtained by weighting and normalizing at least two of the amplitude residual index, phase residual index, noise index, saturation index, drift index and spatial consistency index based on the index weighting coefficient.

[0045] For example, the health diagnosis module uses the fitted signal at the measurement point as a reference to calculate a difference index, including at least two of the following: amplitude residual index, phase residual index, noise index, saturation index, drift index, and spatial consistency index. This application illustrates the calculation of the measurement point health score using a combination of six indices, but of course, any two or more difference indices can be used to calculate the measurement point health score.

[0046] For example, let the initially fitted measurement point fitting signal be denoted as The amplitude residual index is The amplitude residual index is calculated based on the difference between the amplitude of the fitted signal at the measuring point and the amplitude of the voltage signal at the measuring point. The specific formula is as follows:

[0047] in, The amplitude of the voltage signal at the measuring point. The amplitude of the fitted signal at the measurement point. To prevent small positive numbers with a denominator of zero, the following formula is similar.

[0048] For example, the phase residual index is denoted as The phase residual index is calculated based on the difference between the phase of the fitted signal at the measuring point and the phase of the voltage signal at the measuring point. The specific formula is as follows:

[0049] in, The phase of the voltage signal at the measuring point. The phase of the fitted signal at the measurement point, The phase residual normalization threshold can be a positive number, and its unit is the same as the phase difference. It can be determined based on the phase measurement uncertainty of the synchronous acquisition system, the allowable error of channel delay, the offline calibration results, or the allowable phase deviation at the corresponding operating frequency. wrap(·) means to limit the phase difference to the interval [-π, π].

[0050] For example, the noise index is denoted as The noise index is calculated based on the root mean square of the residual noise after removing the fundamental frequency, as shown in the following formula:

[0051] in, To remove the root mean square of residual noise after removing the fundamental frequency, The amplitude of the voltage signal at the measuring point. To prevent small positive numbers with a denominator of zero.

[0052] For example, the saturation index is denoted as The saturation index is calculated based on the proportion of points at the upper limit of sampling, and the specific formula is as follows:

[0053] in, The number of points is the upper limit of sampling, which can be understood as the number of points that reach or approach the upper limit of sampling. N is the number of sampling points in the window.

[0054] For example, the drift index is denoted as The drift index is calculated based on the difference between the recalibration coefficients of the current window and the previous window, and the specific formula is shown below:

[0055] in, The recalibration coefficient for the current window. These are the recalibration coefficients from the previous window. To prevent small positive numbers with a denominator of zero.

[0056] For example, let the spatial consistency index be denoted as The spatial consistency index is calculated based on the difference between the predicted value obtained by interpolating the voltages of adjacent measuring points and the amplitude of the voltage signal at the measuring point. The specific calculation formula is as follows:

[0057] in, The predicted value of the i-th measuring point is obtained by interpolation of adjacent measuring points or model fitting. To prevent small positive numbers with a denominator of zero.

[0058] The health score for each measurement point is obtained by weighted fusion of the above multiple differential indicators.

[0059] As an example, the health score at the measurement point is... The health score at each testing point is obtained by weighting the above six differential indicators, and the specific calculation formula is shown below:

[0060] in, The weighting coefficients for the amplitude residual index are... The weighting coefficients for the phase residual index are... These are the weighting coefficients for the noise index. The weighting coefficient for the saturation index is... The weighting coefficients for the drift index are... The weighting coefficients are used for spatial consistency indicators. The weighting coefficients for each indicator are set based on engineering experience or calibration results. The exponential processing aims to limit the health score to the range [0,1]. Of course, other types of normalization processing are also acceptable. A higher health score indicates a more reliable measurement point. This application sets a health score threshold; measurement points with health scores greater than this threshold are considered reliable measurement points. Of course, for processing exponential health scores, piecewise linear, fuzzy rule-based, expert rule-based, or interpretable machine learning models can also be used.

[0061] As an example, in addition to calculating health scores, the system can also identify anomaly types based on combinations of the aforementioned indicators. For instance, if the saturation indicator continuously increases and a high proportion of sampling points reach the upper limit, it is determined to be channel saturation. Another example is if the noise indicator increases significantly and the amplitude and phase residuals show no synchronization pattern, indicating noise enhancement or grounding interference. Yet another example is if the phase residual changes abruptly while the amplitude residual is small, indicating channel delay or phase reference anomaly. Furthermore, if the spatial consistency residual is significantly higher than adjacent measurement points while other measurement points remain consistent, it is determined to be a loose connection, local probe anomaly, or front-end conditioning anomaly.

[0062] The above-mentioned anomaly identification does not rely on a single fixed threshold, but uses a combination of transmission line model residuals, spatial consistency, health scores, and continuous over-limit counts for judgment, which makes it easy to set thresholds based on power level, frequency, and load status in the engineering field.

[0063] This application uses a multi-indicator health score to comprehensively identify anomalous types such as amplitude anomalies, phase anomalies, noise enhancement, saturation, drift, and spatial consistency disruption.

[0064] As an example, the weight of each measurement point is determined based on the measurement point health score, including: If the consecutive over-limit count of the current measuring point is greater than or equal to the consecutive over-limit number threshold, or if the measuring point health score of the current measuring point is less than the shielding threshold, the measuring point weight of the current measuring point is determined to be an integer 0. Otherwise, the sum of the preset weight sensitivity coefficients of the health scores of each unmasked measurement point is calculated, and the ratio of the power of the current measurement point's health score to the sum is determined as the measurement point weight of the current measurement point.

[0065] For example, this application introduces a shielding threshold for the measuring point, which can be denoted as: If the current measurement point's health score Less than the shielding threshold If the signal collected by the i-th measurement point is unreliable, the measurement point weight is set to an integer 0, meaning the measurement point is blocked. A threshold for the number of consecutive exceedances is also introduced. This means that the measurement point signal has corresponding upper limits, such as voltage upper limits, phase upper limits, saturation upper limits, etc. If the i-th measurement point continuously exceeds the limit... Greater than or equal to the threshold for consecutive over-limit times If so, the measurement point weight is determined to be an integer 0, that is, the measurement point is masked.

[0066] For example, for other measurement points, the measurement point weight is determined based on the specific measurement point health score. The sum of the preset weight sensitivity coefficients raised to the power of the measurement point health scores of each unmasked measurement point is calculated. The ratio of the power of the current measurement point's measurement point health score to the sum is determined as the measurement point weight of the current measurement point. The preset weight sensitivity coefficient is denoted as γ. The formula for the measurement point weight of the current measurement point can be expressed as follows:

[0067] like or ; Other situations in, Let M be the weight of the i-th measurement point, and M be the total number of measurement points. To preset the weight sensitivity coefficient, This is the activation flag for measurement point i. To measure the health score, To prevent small positive numbers with a denominator of zero, For the consecutive exceedance count at the i-th measurement point, This is the threshold for the number of consecutive exceedances. The preset weight sensitivity coefficient, masking threshold, and consecutive exceedance threshold can be set according to the operating conditions.

[0068] After obtaining the measurement point weights for each measurement point, standing wave reconstruction is performed to obtain the target standing wave voltage.

[0069] As an example, when the health status of each measuring point is different, this invention maps the health score to a weight matrix and uses weighted least squares inversion to solve for the forward and reflected waves. Measuring points with higher health scores have a larger weight in the solution, while measuring points with lower or abnormal health scores are reduced in weight or masked, thereby reducing the impact of abnormal measuring points on the estimation results of the forward, reflected, and standing wave distributions. The complex amplitude parameter vectors of the forward and reverse waves to be estimated are shown below:

[0070] in, This represents the vector of complex amplitude parameters for the forward and reverse waves to be estimated. The flag indicating whether the measuring point is shielded is as described above. In matrix form, A W represents the conjugate transpose of the transmission line observation matrix A; W = diag( () is a diagonal matrix composed of the health weights of each measurement point; The positive regularization coefficient is non-negative. This method is used to improve solution stability under conditions of degraded measurement point layout or noise. After obtaining estimates of the complex amplitude parameters of the forward and reverse waves based on the above formulas, multiple characteristics of the target standing wave voltage are calculated, including the magnitude of the reflection coefficient. Standing Population Maximum value of standing wave voltage and minimum value The calculation formulas are as follows:

[0071]

[0072] in, This is an estimate of the complex amplitude of the forward voltage wave. This is an estimate of the complex amplitude of the reflected wave. Additionally, the positions of antinodes and nodes can be determined based on... and The relative phase and transmission line phase constant Calculated and All of these are small positive numbers to prevent the denominator from being zero.

[0073] Of course, standing wave voltage reconstruction can be achieved using multi-point fitting, table lookup correction, transmission line model inversion, Kalman filtering, recursive least squares, robust state estimation, or other data-driven models.

[0074] Figure 3 This is a schematic diagram of standing wave voltage reconstruction under multi-point voltage sampling according to an embodiment of this application.

[0075] like Figure 3 As shown, preliminary fitting observations are performed based on the actual acquired voltage signals. The observation method can be a continuous observation model or a multi-point sampling observation model to obtain the measurement point fitting signal. The measurement point fitting signal is obtained by fitting the corrected voltage, measurement point position, and transmission line phase constant of multiple measurement points within the current acquisition window, and is not predicted from the previous acquisition window to the next window. Preliminary fitting can employ robust fitting, the previous window's reliable weights, or other reliable measurement points after removing the measurement points to be evaluated. Based on the difference between the fitted measurement point signal and the actual measurement point signal, a health score is calculated for each measurement point. The health score is mapped to the measurement point's participation weight, and standing wave reconstruction is performed to obtain the final target standing wave voltage. The standing wave reconstruction output includes the peak value, valley value, and standing wave ratio. Information such as antinode and node positions, and abnormal measurement point status is collected. After the initial health score, the data is reconstructed again based on the updated weights until the preset convergence condition is met or the preset number of iterations is reached.

[0076] This application uses a transmission line standing wave model and weighted least squares inversion to determine the maximum and minimum values ​​of the standing wave voltage, the location of the antinodes, and the standing wave ratio by multiple measurements, thereby reducing the impact of single-point deviations on the overall judgment.

[0077] As an example, signal processing methods for ion cyclotron transmission lines also include: The current calibration coefficients are updated when the reference signal is valid, the sampling saturation index is less than the preset threshold, the power change rate is less than the power threshold, the model residual is less than the residual threshold, and the measurement point health score is greater than the score threshold.

[0078] For example, to avoid misinterpreting actual load changes as probe drift, this application sets a stability gating within the calibration window. The current calibration coefficients are updated only when the reference signal is valid, the sampling saturation index is less than a preset threshold, the power change rate is less than a power threshold, the model residual is less than a residual threshold, and the measurement point health score is greater than a score threshold. This can be understood as allowing the calibration coefficients to be updated when the reference signal is valid, sampling is not saturated, power changes are stable, the model residual is normal, and the measurement point health score meets the requirements.

[0079] For example, it can be expressed as a formula, namely:

[0080] in, This indicates a calibration permission flag; for example, a value of 1 allows updates, while a value of 0 prevents updates. As a valid reference signal, As a saturation index, The preset threshold for the sampling saturation index, Power threshold For the overall model residuals, The residual threshold, This is the scoring threshold.

[0081] As an example, if the current calibration coefficients include at least one of the following: recalibration coefficients, amplitude calibration coefficients, phase compensation, and bias compensation, then updating the current calibration coefficients includes at least one of the following: Construct an instantaneous recalibration estimate based on the current window standing wave model fitting value, the voltage signal at the measurement point, and the current bias compensation amount; update the amplitude calibration coefficient based on the amplitude of the instantaneous recalibration estimate. The phase compensation amount is updated based on the phase estimated by instantaneous recalibration; The updated step size is obtained by multiplying the base step size and the health score of the measurement point. The first product result of the updated step size and the instantaneous recalibration estimate is calculated. The second product result of the difference between the preset integer and the updated step size and the current recalibration coefficient is calculated. The recalibration coefficient is updated based on the amplitude limiting result after summing the first product result and the second product result. Calculate the third product of the bias update step size and the sample mean within the preset window, the fourth product of the difference between the preset integer and the bias update step size and the current bias compensation amount, and update the bias compensation amount based on the sum of the third and fourth product results.

[0082] For example, when calibration permits, the online calibration module constructs an instantaneous recalibration estimate using the reference voltage, a known verification signal, or the current window standing wave model fitting value reconstructed from a reliable measurement point, as shown below:

[0083] in, This is the instantaneous recalibration estimate obtained for the i-th measurement point within the current calibration window. Y is the reference voltage complex phasor (current window VSWR model fitted value) given by the reference voltage, known verification signal, or reliable model prediction. i b represents the original complex phase signal (measurement point voltage signal) extracted in the current window. i This is the bias compensation value before this update (current bias compensation amount). To prevent small positive numbers with a denominator of zero.

[0084] The amplitude of the instantaneous recalibration estimate is determined to be the updated amplitude calibration coefficient, and the phase of the instantaneous recalibration estimate is determined to be the phase compensation amount.

[0085] For example, to suppress parameter jumps caused by single noise or short-term disturbances, a limited exponential recursive update is used, as shown in the following equation:

[0086]

[0087] Where k is the discrete calibration window number, These are the recalibration coefficients before the update. These are the updated recalibration coefficients; Let i be the update step size for the i-th measurement point. Base step size; Calibration permission flag; {·} is the limiting operator, used to limit the amplitude of a single gain or phase change; Update the step size for the bias; mean{·} represents the sample mean within the specified window. The default integer is 1.

[0088] It is understandable that the above recursive update is equivalent to a first-order low-pass filter or exponential smoothing, and the new estimate is only based on... The proportion of historical calibration values ​​is included in the updated results. Retained when the health score of the measuring point decreases or the calibration gate is not met. Decrease or zero, freeze or slowly change the calibration coefficients, and the limiting operator can further constrain single-step variables, thereby suppressing parameter abrupt changes caused by noise spikes and short-term disturbances. This application uses the above recursive update to suppress calibration coefficient jumps caused by single noise or short-term disturbances. The bias compensation is updated only within the no-RF or low-power reference window, and is frozen or limited during high-power operation.

[0089] This application updates the calibration coefficients by using calibration gating and recursive updates when the reference signal is valid, the sampling is not saturated, and the operating state is stable, thereby reducing the impact of gain drift, phase deviation, and zero drift on voltage measurement during long pulse operation.

[0090] As an example, signal processing methods for ion cyclotron transmission lines also include: Obtain the number of reliable measurement points, the number of model conditions, and the safety judgment voltage; Alarm protection is implemented based on the number of reliable measurement points, the number of model conditions, and the safety judgment voltage.

[0091] For example, the protection alarm interface determines whether to allow the output of a reliable standing wave voltage based on the number of reliable measurement points, the number of model conditions, the reconstruction residual, and the voltage margin. The number of reliable measurement points and the number of model conditions are:

[0092] in, For the number of reliable measurement points, The condition number is the model condition number. cond(·) represents the condition number operation, which can be understood as the ratio of the largest singular value to the smallest singular value of the matrix. The larger the model condition number, the more sensitive the inversion is to noise and measurement point errors.

[0093] For example, to reduce the risk of overvoltage failure due to reconstruction error, this application superimposes an uncertainty margin on the maximum standing wave voltage to obtain a safe judgment voltage, as shown below:

[0094] in, The safety determination voltage can be understood as a conservative voltage value used for protection determination. To obtain the maximum standing wave voltage through reconstruction, The voltage uncertainty is estimated from the fitting residuals and parameter covariance. The confidence margin coefficient is... For engineering safety margin.

[0095] This application can output reliable overvoltage alarms or data acquisition link health alarms based on the number of reliable measurement points, the number of model conditions, and the safety judgment voltage, thereby reducing the risk of false alarms and missed alarms.

[0096] As an example, alarm protection based on the number of trusted measurement points, the number of model conditions, and the safety judgment voltage includes: If the number of trusted measurement points is less than the quantity threshold, or the number of model conditions is greater than the condition number threshold, a health alarm for the data acquisition link will be triggered. If the number of reliable measurement points is greater than or equal to the number threshold and the number of model conditions is less than or equal to the number of conditions threshold, and the safety judgment voltage is greater than or equal to the protection threshold, then an overvoltage alarm or interlock protection will be triggered. If the number of reliable measurement points is greater than or equal to the quantity threshold and the number of model conditions is less than or equal to the condition number threshold, the target standing wave voltage will be output normally if the safety judgment voltage is less than the protection threshold.

[0097] For example, if the number of reliable measurement points is less than the quantity threshold, or the number of model conditions is greater than the condition number threshold, the acquired signal can be considered unreliable, triggering a health alarm for the acquisition link. Conversely, if the number of reliable measurement points is greater than or equal to the quantity threshold and the number of model conditions is less than or equal to the condition number threshold, the acquired signal is considered reliable. Under this premise, further judgment is made based on the safety judgment voltage. If the safety judgment voltage is greater than or equal to the protection threshold, an overvoltage alarm or interlocking protection is triggered; if the safety judgment voltage is less than the protection threshold, the target standing wave voltage is output normally.

[0098] For example, the protection threshold is denoted as The quantity threshold is The condition number threshold is When the safety judgment voltage is greater than or equal to the protection threshold, the number of reliable measurement points is greater than or equal to the quantity threshold, and the number of model conditions is less than or equal to the condition number threshold, an overvoltage alarm or interlock protection is triggered, as shown in the following formula: ,like and and

[0099] in, This is the overvoltage alarm or interlock protection flag. When ≥ When the number of reliable measurement points and the number of model conditions meet the requirements, a reliable overvoltage alarm or interlock signal is output. If the number of reliable measurement points is insufficient, the number of model conditions exceeds the limit, or the number of measurement points with low health scores exceeds the set threshold, the system outputs a health alarm for the acquisition link to avoid using unreliable measurement results directly as the basis for protection actions.

[0100] The alarm protection in this application can distinguish between actual standing wave changes and measurement link anomalies.

[0101] Figure 4 This is a schematic diagram of a health scoring, anomaly shielding, and protection alarm state machine according to an embodiment of this application.

[0102] like Figure 4 As shown, a normalized weighted fusion of multiple indicators, including amplitude residual, phase residual, noise, saturation, drift, and spatial consistency, is used to generate a health score and corresponding weights for the measurement points. The shielding threshold for the aforementioned measurement points is then considered. Based on this, a warning threshold for slightly lower measuring points was also introduced. Warning threshold Greater than the shielding threshold Health scores can be obtained based on the measurement points. The testing points are categorized, and the health score of the testing points is... ≥ When it is a normal measuring point, when ≤ < At that time, the measuring point was found to have a slight anomaly and its weight was reduced. < At that time, the measuring point was severely abnormal and shielded, and the count of consecutive exceedances occurred at the i-th measuring point. ≥ When the number of consecutive over-limits reaches a certain threshold, the measurement point is also considered severely abnormal and is shielded. The abnormal measurement point status in this application is determined by a health score, shielding flag, consecutive over-limit count, and amplitude residual, phase residual, noise, saturation, drift, and spatial consistency indices. Additionally, when the safety judgment voltage is greater than or equal to the protection threshold, the number of reliable measurement points is greater than or equal to the quantity threshold, and the number of model conditions is less than or equal to the condition number threshold, a reliable overvoltage alarm or interlock signal is output. When the number of reliable measurement points is less than the quantity threshold, or the number of model conditions is greater than the condition number threshold, a data acquisition link health alarm is output, and a log is recorded. The log recording module records the original sampling summary, complex phasor / envelope, reference signal status, calibration parameters, health score, reconstruction result, abnormality type, shielding status, protection alarm action, and manual confirmation mark according to the data acquisition window, for use in false alarm review, threshold verification, and operation and maintenance.

[0103] As an example, the overall process of this application can be described as follows: 1) Enter the acquisition window, calibration window, stable operation range, or pulse gap.

[0104] 2) Synchronously collect voltage probe signals, reference voltage, verification signals, sampling channel status, and transmission line position parameters at each measurement point.

[0105] 3) Filter, synchronize, delay compensate, and extract fundamental frequency complex phase / envelope from the signals at each measurement point.

[0106] 4) Calculate the corrected measurement point voltage based on the current calibration coefficient, bias, and channel delay.

[0107] 5) Reconstruct the forward and reflected waves based on the transmission line model, and calculate the peak value, valley value, standing wave ratio, and antinode / node position of the standing wave.

[0108] 6) Calculate residual, noise, saturation, drift and spatial consistency indices to generate a health score for the measurement point.

[0109] 7) Generate collection weights based on health scores, reduce the weight of slightly abnormal measurement points, and block severely abnormal or continuously exceeding measurement limits measurement points.

[0110] 8) Update the calibration coefficients recursively when the calibration gate is satisfied, and freeze the calibration coefficients when it is not satisfied.

[0111] 9) Output a reliable standing wave voltage or protection alarm status based on the number of reliable measurement points, the number of model conditions, and the safety judgment voltage.

[0112] 10) Record the acquisition window, original sampling volume, calibration coefficient, reconstruction result, health score, acquisition weight, anomaly type, shielding status, and protection action.

[0113] This application sets calibration gating for reference signal validity, saturation state, power change rate, model residual, and health score within the calibration window or stable operating range. It updates the recalibration coefficients through adaptive step size and limits single gain and phase changes. It recursively updates bias compensation during low power, no RF, or calibration windows. It freezes or limits bias updates during high power operation. It reduces the weight of measurement points when there are minor anomalies. It shields measurement points when there are serious anomalies or continuous over-limits. It outputs an alarm for insufficient reliability of the acquisition link when there are insufficient reliable measurement points or the model is ill-conditioned. It outputs a reliable overvoltage alarm or interlock signal when the safety judgment voltage exceeds the protection threshold. It records the anomaly type, shielding status, calibration parameters, and protection actions in a correlated manner.

[0114] This application also proposes a signal processing device for an ion cyclotron transmission line.

[0115] As an example, such as Figure 5As shown, the signal processing device for the ion cyclotron transmission line includes: a correction module 501, used to synchronously acquire measurement point signals based on voltage probes and correct the measurement point signals based on the current calibration coefficients to obtain the corrected measurement point voltage signals; a fitting module 502, used to perform preliminary signal fitting based on the measurement point voltage signals using a forward wave and reflected wave superposition model to obtain the measurement point fitted signal; a generation module 503, used to calculate the measurement point difference index based on the measurement point fitted signal and the measurement point voltage signal, and generate a measurement point health score based on the measurement point difference index; and a reconstruction module 504, used to determine the measurement point weights of each measurement point based on the measurement point health score, and perform standing wave reconstruction based on the measurement point weights to obtain the target standing wave voltage.

[0116] This application also proposes a computer-readable storage medium.

[0117] In this embodiment, a computer program is stored on a computer-readable storage medium, and when the computer program is executed by a processor, it implements the steps of the signal processing method for the ion cyclotron transport line described above.

[0118] Figure 6 A block diagram of an electronic device provided in an embodiment of this application.

[0119] This application provides an electronic device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the above-described signal processing method for ion cyclotron transport lines.

[0120] like Figure 6 As shown, for ease of understanding, embodiments of this application illustrate a specific electronic device.

[0121] Electronic devices are intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. Electronic devices can also represent various forms of mobile devices, such as personal digital processors, cellular phones, smartphones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the present disclosure described and / or claimed herein.

[0122] like Figure 6As shown, the device includes a computing unit 601, which can perform various appropriate actions and processes based on a computer program stored in a read-only memory (ROM) 602 or a computer program loaded into a random access memory (RAM) 603 from a storage unit 608. The RAM 603 may also store various programs and data required for the operation of the electronic device. The computing unit 601, ROM 602, and RAM 603 are interconnected via a bus 604. An input / output (I / O) interface 605 is also connected to the bus 604.

[0123] Multiple components in the electronic device are connected to the I / O interface 605. These components include: an input unit 606, such as a keyboard or mouse; an output unit 607, such as various types of displays or speakers; a storage unit 608, such as a disk or optical disk; and a communication unit 609, such as a network interface card (NIC), a modem, or a wireless transceiver. The communication unit 609 allows the electronic device to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0124] The computing unit 601 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the computing unit 601 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various computing units running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The computing unit 601 performs the various methods described above, such as the signal processing method for ion cyclotron transmission lines. For example, in some embodiments, the signal processing method for ion cyclotron transmission lines can be implemented as a computer software program tangibly contained in a machine-readable medium, such as storage unit 608. In some embodiments, part or all of the computer program can be loaded and / or installed on an electronic device via ROM 602 and / or communication unit 609. When the computer program is loaded into RAM 603 and executed by the computing unit 601, the signal processing method for ion cyclotron transmission lines described above can be performed. Alternatively, in other embodiments, the computing unit 601 can be configured to perform the signal processing method for ion cyclotron transmission lines by any other suitable means (e.g., by means of firmware).

[0125] It should be noted that the logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be specifically implemented in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this application, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: electrical connections (electronic devices) having one or more wires, portable computer disk drives (magnetic devices), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Furthermore, computer-readable media can even be paper or other suitable media on which programs can be printed, because programs can be obtained electronically, for example, by optically scanning the paper or other media, followed by editing, interpreting, or otherwise processing as necessary, and then stored in computer memory.

[0126] It should be understood that various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0127] In the description of this application, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this application, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples.

[0128] In the description of this application, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential", etc., indicating the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings, are only for the convenience of describing this application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of this application.

[0129] Furthermore, the terms "first," "second," etc., used in the embodiments of this application are for descriptive purposes only and should not be construed as indicating or implying relative importance, or implicitly specifying the number of technical features indicated in this embodiment. Therefore, features defined with terms such as "first" and "second" in the embodiments of this application can explicitly or implicitly indicate that the embodiment includes at least one of those features. In the description of this application, the word "multiple" means at least two or more, such as two, three, four, etc., unless otherwise explicitly and specifically defined in the embodiments.

[0130] In this application, unless otherwise explicitly specified or limited in the embodiments, the terms "installation," "connection," "joining," and "fixing" appearing in the embodiments should be interpreted broadly. For example, a connection can be a fixed connection, a detachable connection, or an integral part; it can also be a mechanical connection, an electrical connection, etc. Of course, it can also be a direct connection, or an indirect connection through an intermediate medium, or it can be the internal communication between two components, or the interaction between two components. Those skilled in the art can understand the specific meaning of the above terms in this application based on the specific implementation.

[0131] In this application, unless otherwise expressly specified and limited, "above" or "below" the second feature can mean that the first feature is in direct contact with the second feature, or that the first feature is in indirect contact with the second feature through an intermediate medium. Furthermore, "above," "on top of," and "over" the second feature can mean that the first feature is directly above or diagonally above the second feature, or simply that the first feature is at a higher horizontal level than the second feature. "Below," "below," and "under" the second feature can mean that the first feature is directly below or diagonally below the second feature, or simply that the first feature is at a lower horizontal level than the second feature.

[0132] Although embodiments of this application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting this application. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of this application.

Claims

1. A signal processing method for an ion cyclotron transmission line, characterized in that, The method includes: The voltage probe synchronously acquires the measurement point signal, and the measurement point signal is corrected based on the current calibration coefficient to obtain the corrected measurement point voltage signal; Based on the voltage signal at the measurement point, preliminary signal fitting is performed using a forward wave and reflected wave superposition model to obtain the fitted signal at the measurement point. The measurement point difference index is calculated based on the measurement point fitting signal and the measurement point voltage signal, and a measurement point health score is generated based on the measurement point difference index. The weight of each measuring point is determined based on the health score of the measuring point, and the standing wave is reconstructed based on the measuring point weight to obtain the target standing wave voltage.

2. The signal processing method for the ion cyclotron transmission line according to claim 1, characterized in that, The measurement point difference index includes at least two of the following: amplitude residual index, phase residual index, noise index, saturation index, drift index, and spatial consistency index. Generating a measurement point health score based on the measurement point difference index includes: Based on the index weighting coefficients, at least two of the amplitude residual index, the phase residual index, the noise index, the saturation index, the drift index, and the spatial consistency index are weighted and normalized to obtain the measurement point health score.

3. The signal processing method for the ion cyclotron transmission line according to claim 1, characterized in that, The determination of the weight of each measuring point based on the health score of the measuring points includes: If the consecutive over-limit count of the current measuring point is greater than or equal to the consecutive over-limit number threshold, or if the measuring point health score of the current measuring point is less than the shielding threshold, the measuring point weight of the current measuring point is determined to be an integer 0. Otherwise, the sum of the preset weight sensitivity coefficients of the health scores of each unmasked measurement point is calculated, and the ratio of the power of the current measurement point's health score to the sum is determined as the measurement point weight of the current measurement point.

4. The signal processing method for the ion cyclotron transmission line according to claim 1, characterized in that, The method further includes: The current calibration coefficients are updated when the reference signal is valid, the sampling saturation index is less than a preset threshold, the power change rate is less than a power threshold, the model residual is less than a residual threshold, and the health score of the measurement point is greater than a score threshold.

5. The signal processing method for the ion cyclotron transmission line according to claim 4, characterized in that, The current calibration coefficients include at least one of the following: recalibration coefficients, amplitude calibration coefficients, phase compensation amount, and bias compensation amount. Updating the current calibration coefficients includes at least one of the following: An instantaneous recalibration estimate is constructed based on the current window standing wave model fitting value, the voltage signal at the measurement point, and the current bias compensation amount. The amplitude calibration coefficient is then updated based on the amplitude of the instantaneous recalibration estimate. The phase compensation amount is updated based on the phase estimated by the instantaneous recalibration; The updated step size is obtained by multiplying the base step size and the health score of the measurement point. The first product result of the updated step size and the instantaneous recalibration estimate is calculated. The second product result of the difference between the preset integer and the updated step size and the current recalibration coefficient is calculated. The recalibration coefficient is updated based on the amplitude limiting result after summing the first product result and the second product result. Calculate the third product of the bias update step size and the sample mean within the preset window, and the fourth product of the difference between the preset integer and the bias update step size and the current bias compensation amount. Update the bias compensation amount based on the sum of the third product and the fourth product.

6. The signal processing method for the ion cyclotron transmission line according to claim 1, characterized in that, The method further includes: Obtain the number of reliable measurement points, the number of model conditions, and the safety judgment voltage; Alarm protection is provided based on the number of reliable measurement points, the number of model conditions, and the safety judgment voltage.

7. The signal processing method for the ion cyclotron transmission line according to claim 6, characterized in that, The alarm protection based on the number of trusted measurement points, the number of model conditions, and the safety judgment voltage includes: If the number of trusted measurement points is less than the quantity threshold, or if the number of model conditions is greater than the condition number threshold, a data acquisition link health alarm is triggered. If the number of reliable measurement points is greater than or equal to the number threshold and the number of model conditions is less than or equal to the number threshold, and the safety judgment voltage is greater than or equal to the protection threshold, then an overvoltage alarm or interlock protection is triggered. If the number of reliable measurement points is greater than or equal to the number threshold and the number of model conditions is less than or equal to the condition number threshold, and the safety judgment voltage is less than the protection threshold, then the target standing wave voltage is output normally.

8. The signal processing method for the ion cyclotron transmission line according to claim 2, characterized in that, The amplitude residual index is calculated based on the degree of difference between the amplitude of the fitted signal at the measuring point and the amplitude of the voltage signal at the measuring point. The phase residual index is calculated based on the degree of difference between the phase of the fitted signal at the measurement point and the phase of the voltage signal at the measurement point. The noise index is calculated based on the root mean square of the residual noise after removing the fundamental frequency. The saturation index is calculated based on the percentage of points at the upper limit of sampling. The drift index is calculated based on the degree of difference between the recalibration coefficients of the current window and the previous window; The spatial consistency index is calculated based on the degree of difference between the predicted value obtained by interpolating the voltages of adjacent measurement points and the amplitude of the voltage signal at the measurement point.

9. The signal processing method for the ion cyclotron transmission line according to claim 1, characterized in that, Before correcting the measurement point signal based on the current calibration coefficient, the method further includes: The measurement point signal is subjected to at least one of complex phase extraction or envelope extraction.

10. A signal processing device for an ion cyclotron transmission line, characterized in that, The device includes: The calibration module is used to synchronously acquire the measurement point signal based on the voltage probe, and to correct the measurement point signal based on the current calibration coefficient to obtain the corrected measurement point voltage signal; The fitting module is used to perform preliminary signal fitting based on the voltage signal at the measurement point using a forward wave and reflected wave superposition model to obtain the fitted signal at the measurement point. The generation module is used to calculate the measurement point difference index based on the measurement point fitting signal and the measurement point voltage signal, and generate a measurement point health score based on the measurement point difference index. The reconstruction module is used to determine the measurement point weight of each measurement point based on the measurement point health score, and to perform standing wave reconstruction based on the measurement point weight to obtain the target standing wave voltage.

11. An electronic device, characterized in that, It includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the signal processing method for the ion cyclotron transmission line according to any one of claims 1-9.

12. A computer-readable storage medium, characterized in that, It stores a computer program, which, when executed by a processor, implements the steps of the signal processing method for the ion cyclotron transmission line according to any one of claims 1-9.

Images