A method, apparatus, equipment, medium, and product for determining channel error values ​​in a data acquisition system.

Abstract

This invention discloses a method, apparatus, device, medium, and product for determining the error value of a channel in a data acquisition system. The method includes: determining the parameters of a first median frequency and a second median frequency within the frequency range corresponding to each channel; if the parameters of the first median frequency and the second median frequency corresponding to the channel do not satisfy the iteration stopping condition, updating the frequency range corresponding to the channel based on the parameters of the first median frequency and the second median frequency corresponding to the channel; and, based on the updated frequency range corresponding to the channel, returning to execute an operation based on a three-parameter sine fitting algorithm to determine the parameters of the first median frequency and the second median frequency within the frequency range corresponding to the channel, until the parameters of the first median frequency and the second median frequency satisfy the iteration stopping condition; and determining the error value of the channel based on the parameters of the first median frequency and the second median frequency within the frequency range corresponding to the channel.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of data processing technology, and in particular to a method, apparatus, device, medium and product for determining the error value of a channel in a data acquisition system. Background Technology

[0002] Data acquisition systems can increase the speed of analog-to-digital converters (ADCs) without reducing the data resolution of individual ADCs; however, errors exist between channels in the data acquisition system. These inter-channel errors can severely impact the overall performance of the data acquisition system.

[0003] How to accurately estimate the error value of the channel in a data acquisition system is a technical problem that urgently needs to be solved by those skilled in the art. Summary of the Invention

[0004] This invention provides a method, apparatus, device, medium, and product for determining the error value of a channel in a data acquisition system, so as to achieve accurate estimation of the error value of a channel in the data acquisition system.

[0005] According to one aspect of the present invention, a method for determining the error value of a channel in a data acquisition system is provided, comprising:

[0006] Based on the initial frequency of the target sine wave and the frequency offset value of each channel in the data acquisition system, the frequency range interval corresponding to each channel is determined, wherein the target sine wave includes low-frequency sine waves and / or high-frequency sine waves.

[0007] Based on the three-parameter sine fitting algorithm, the parameters of the first median frequency and the second median frequency in the frequency range corresponding to each channel are determined.

[0008] For any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, then the frequency range interval corresponding to the channel is updated according to the parameters of the first median frequency and the second median frequency corresponding to the channel. Based on the updated frequency range interval corresponding to the channel, the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel is returned to be executed until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition. Based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, the error value of the channel is determined.

[0009] According to another aspect of the present invention, an error value determination device for a channel in a data acquisition system is provided, the error value determination device for a channel in the data acquisition system comprising:

[0010] The frequency range interval determination module is used to determine the frequency range interval corresponding to each channel based on the initial frequency of the target sine wave and the frequency offset value corresponding to each channel in the data acquisition system, wherein the target sine wave includes a low-frequency sine wave and / or a high-frequency sine wave.

[0011] The first median frequency parameter and second median frequency parameter determination module is used to determine the first median frequency parameter and second median frequency parameter in the frequency range interval corresponding to each channel based on a three-parameter sine fitting algorithm.

[0012] The channel error value determination module is used to, for any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, update the frequency range interval corresponding to the channel according to the parameters of the first median frequency and the second median frequency corresponding to the channel, and return to execute the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel based on the updated frequency range interval corresponding to the channel, until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition, and determine the error value of the channel according to the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel.

[0013] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0014] At least one processor; and

[0015] A memory communicatively connected to the at least one processor; wherein,

[0016] The memory stores a computer program that can be executed by the at least one processor, which enables the at least one processor to perform the channel error value determination method in the data acquisition system according to any embodiment of the present invention.

[0017] According to another aspect of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium storing computer instructions for causing a processor to execute and implement the method for determining the error value of a channel in a data acquisition system according to any embodiment of the present invention.

[0018] According to another aspect of the present invention, a computer program product is provided, which, when executed by a processor, implements a method for determining the error value of a channel in a data acquisition system as described in any of the embodiments of the present invention.

[0019] This invention, in its embodiments, determines the frequency range interval for each channel based on the initial frequency of the target sine wave and the frequency offset value corresponding to each channel in the data acquisition system. Using a three-parameter sine fitting algorithm, it determines the parameters of the first median frequency and the second median frequency within the frequency range interval for each channel. For any channel, if the parameters of the first median frequency and the second median frequency do not meet the iteration stopping condition, the frequency range interval for the channel is updated based on these parameters. The process then returns to executing the three-parameter sine fitting algorithm to determine the parameters of the first median frequency and the second median frequency within the frequency range interval for the channel, continuing until these parameters meet the iteration stopping condition. Finally, the error value of the channel is determined based on these parameters, enabling accurate estimation of the error values ​​for each channel in the data acquisition system.

[0020] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description

[0021] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as a limitation on the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0022] Figure 1 This is a flowchart of a method for determining the error value of a channel in a data acquisition system according to an embodiment of the present invention;

[0023] Figure 2 This is a flowchart of another method for determining the error value of a channel in a data acquisition system according to an embodiment of the present invention;

[0024] Figure 3 This is a schematic diagram of the structure of a channel error value determination device in a data acquisition system according to an embodiment of the present invention;

[0025] Figure 4 This is a schematic diagram of the structure of an electronic device according to an embodiment of the present invention. Detailed Implementation

[0026] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0027] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0028] It is understood that before using the technical solutions disclosed in the various embodiments of this disclosure, users should be informed of the types, scope of use, and usage scenarios of the personal information involved in this disclosure in an appropriate manner in accordance with relevant laws and regulations, and user authorization should be obtained.

[0029] It should be noted that the data acquisition system can be an analog-to-digital converter (TIADC). In existing technologies, there are two main methods for obtaining error values ​​in data acquisition systems.

[0030] The first method involves applying a specific test signal to the data acquisition system and obtaining the output with sampling time mismatch for that channel from the sub-channel analog-to-digital converter (ADC). The sampling time mismatch value of the sub-channel ADC can then be obtained by comparing the ideal output value with the actual value. However, this method requires high accuracy of the test signal; low accuracy will result in low accuracy of the calculated sampling time mismatch value.

[0031] The second approach involves a least-squares sine curve fitting algorithm to estimate parameters such as channel delay, data acquisition system gain, offset, and phase. This algorithm requires the sine wave frequency to be known in order to estimate the gain, offset, and phase parameters. However, since the frequency of a sine wave is usually unknown or inaccurate, the three parameters obtained using this algorithm will inevitably contain estimation errors. This results in low accuracy of the estimated error values ​​for the data acquisition system.

[0032] Example 1

[0033] Figure 1 This is a flowchart illustrating a method for determining the error value of a channel in a data acquisition system according to an embodiment of the present invention. This embodiment is applicable to situations where the error value of a channel in a data acquisition system is determined. This method can be executed by a channel error value determination device in the data acquisition system according to this embodiment of the present invention. This device can be implemented in software and / or hardware, such as... Figure 1 As shown, the method specifically includes the following steps:

[0034] S110, based on the initial frequency of the target sine wave and the frequency offset value of each channel in the data acquisition system, determine the frequency range interval corresponding to each channel.

[0035] In this embodiment, the target sine wave includes a low-frequency sine wave and / or a high-frequency sine wave. The frequency range interval is the frequency range interval corresponding to each channel in the data acquisition system.

[0036] In a specific example, if the initial frequency of the target sine wave is f0, the frequency offset value corresponding to channel 1 in the data acquisition system is Δf. max Then the frequency range corresponding to channel 1 in the data acquisition system is [f0-Δf max ,f0+Δf max ].

[0037] Optionally, before determining the frequency range for each channel based on the initial frequency of the target sine wave and the frequency offset values ​​corresponding to each channel in the data acquisition system, the method further includes:

[0038] The initial frequency of the target sine wave is obtained by performing a Fourier transform on the target sine wave.

[0039] In this embodiment, a Fourier transform is performed on the target sine wave to find the center frequency point, and the frequency corresponding to the center frequency point is determined as the initial frequency.

[0040] The target sine wave is input into the data acquisition system to obtain the total number of sampling points and the total sampling time for each channel.

[0041] In this embodiment, the target sine wave is input into the data acquisition system. Based on the sampled data output from each channel in the data acquisition system, the total number of sampling points and the total sampling time for each channel are determined. In a specific example, the sampled data output from channel 1 in the data acquisition system includes: known times t1, t2, ..., t n The corresponding sampled data y1, y2, ..., y n .

[0042] The frequency offset value for each channel is determined based on the total number of sampling points and the total sampling time for each channel.

[0043] In this embodiment, the frequency offset value corresponding to each channel can be determined based on the total number of sampling points and the total sampling time of each channel as follows: the sampling rate corresponding to each channel is determined based on the total number of sampling points and the total sampling time of each channel, and the ratio of the sampling rate corresponding to each channel to the total number of sampling points is determined as the frequency offset value corresponding to each channel.

[0044] In a specific example, based on the formula The sampling rate corresponding to channel 1 is calculated, where v is the sampling rate of channel 1, n is the total number of sampling points for channel 1, and t is the total sampling time for channel 1. Based on the formula Δf... max =v / n, calculate the frequency offset value corresponding to channel 1, Δf max This is the frequency offset value corresponding to channel 1.

[0045] S120, based on the three-parameter sine fitting algorithm, determines the parameters of the first median frequency and the second median frequency in the frequency range corresponding to each channel.

[0046] In this embodiment, the parameters include: gain, phase, DC offset, and RMS value of the fitted residual.

[0047] In this embodiment, given the frequency, gain, DC offset, and phase, a sine function can be uniquely determined. Let the sine function model be:

[0048]

[0049] Among them, f in Let be a sine wave with a known input frequency, 'a' be the gain of the sine wave, and 'C' be the DC offset of the sine wave. The phase of the sine wave. It is the standard sine wave expression, Asin(2πf in t)+Bcos(2πf in t)+C is a three-parameter expression, for A, B, a, The relationship can be represented as:

[0050]

[0051] Given the total number of sampling points and the total sampling time for each channel, for example, channel 1, times t1, t2, ..., t n n sampling points y1, y2, ..., y2 were obtained. n It should be noted that the frequency of the input target sine wave is a known parameter, and A, B, and C are used to minimize the sum of squared residuals E to complete the fitting.

[0052]

[0053] It should be noted that since the above formula is linearly closed, the parameters A, B, and C can be uniquely determined. Construct the following matrix:

[0054]

[0055]

[0056] The above formula, expressed in matrix form, can be represented as:

[0057] E = (y - Ds) T (y-Ds);

[0058] Among them, (y-Ds) T Taking the transpose of matrix (y-Ds), we can obtain the final estimated three parameters as follows:

[0059]

[0060] Based on the above formulas, the gain, phase, and DC offset of the sinusoidal wave with known frequency corresponding to each channel are calculated. The effective value of the fitting residual is then calculated using the following formula:

[0061]

[0062] Where E is the sum of squared residuals, n is the total number of sampling points, and ρ is the effective value of the fitted residuals.

[0063] S130, for any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, then update the frequency range interval corresponding to the channel according to the parameters of the first median frequency and the second median frequency corresponding to the channel. Based on the updated frequency range interval corresponding to the channel, return to execute the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel using the three-parameter sine fitting algorithm, until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition. Based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, determine the error value of the channel.

[0064] In this embodiment, for any channel, if the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency do not meet the iteration stopping condition, then the frequency range interval corresponding to the channel is updated according to the parameters of the first median frequency and the second median frequency. Based on the updated frequency range interval corresponding to the channel, the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel is returned to be executed, until the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency meet the iteration stopping condition. Based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, the error value of the channel is determined.

[0065] In this embodiment, the parameter of the first median frequency can be the effective value of the fitting residual of the first median frequency, and the parameter of the second median frequency can be the effective value of the fitting residual of the second median frequency.

[0066] In this embodiment, for any channel, if the effective values ​​of the fitting residuals of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, then a frequency range interval update rule is determined based on the parameters of the first median frequency and the second median frequency corresponding to the channel. The frequency range interval corresponding to the channel is updated according to the frequency range interval update rule. Based on the updated frequency range interval corresponding to the channel, the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel is returned to be executed, until the effective values ​​of the fitting residuals of the first median frequency and the second median frequency meet the iteration stopping condition. Based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, the error value of the channel is determined.

[0067] It should be noted that the frequency range interval update rule can be determined based on the parameters of the first median frequency and the second median frequency corresponding to the channel as follows: If the effective value of the fitting residual of the first median frequency is less than the effective value of the fitting residual of the second median frequency, then the second median frequency of the frequency range interval before the update is determined as the left boundary frequency of the frequency range interval after the update, and the right boundary frequency of the frequency range interval before the update is determined as the right boundary frequency of the frequency range interval after the update; if the effective value of the fitting residual of the first median frequency is greater than the effective value of the fitting residual of the second median frequency, then the first median frequency of the frequency range interval before the update is determined as the right boundary frequency of the frequency range interval after the update, and the left boundary frequency of the frequency range interval before the update is determined as the left boundary frequency of the frequency range interval after the update.

[0068] Optionally, updating the frequency range corresponding to the channel based on the parameters of the first median frequency and the second median frequency includes:

[0069] If the effective value of the fitting residual of the first median frequency is less than the effective value of the fitting residual of the second median frequency, then the second median frequency of the frequency range before the update is determined as the left boundary frequency of the frequency range after the update, and the right boundary frequency of the frequency range before the update is determined as the right boundary frequency of the frequency range after the update.

[0070] If the effective value of the fitting residual of the first median frequency is greater than the effective value of the fitting residual of the second median frequency, then the first median frequency of the frequency range interval before the update is determined as the right boundary frequency of the frequency range interval after the update, and the left boundary frequency of the frequency range interval before the update is determined as the left boundary frequency of the frequency range interval after the update.

[0071] In this embodiment, the left boundary frequency of the frequency range interval is the minimum frequency in the frequency range interval, and the right boundary frequency of the frequency range interval is the maximum frequency in the frequency range interval.

[0072] In this embodiment, for any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, the frequency range interval corresponding to the channel is updated according to the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency. The three-parameter sine fitting algorithm is executed again to obtain the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition, and the iteration stops.

[0073] In this embodiment, the method for updating the frequency range interval corresponding to the channel based on the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency can be as follows: If the effective value of the fitting residual of the first median frequency is less than the effective value of the fitting residual of the second median frequency, then the second median frequency of the frequency range interval before the update is determined as the left boundary frequency of the updated frequency range interval, and the right boundary frequency of the frequency range interval before the update is determined as the right boundary frequency of the updated frequency range interval. For example, the updated frequency range interval can be [f T ,f R ], where f T f is the left boundary frequency (i.e., the minimum frequency) of the updated frequency range interval, which is the second median frequency of the original frequency range interval. R The right boundary frequency (i.e., the maximum frequency) of the original frequency range interval is used. The first median frequency of the updated frequency range interval is determined based on the following formula:

[0074] f M =f T +0.618(f R -f T );

[0075] Among them, f M f is the first median frequency of the updated frequency range. R f is the right boundary frequency of the frequency range before the update. T This is the left boundary frequency of the updated frequency range.

[0076] It should be noted that the second median frequency of the updated frequency range is equal to the first median frequency of the original frequency range.

[0077] If the effective value of the fitting residual of the first median frequency is greater than the effective value of the fitting residual of the second median frequency, then the first median frequency of the frequency range interval before the update is determined as the right boundary frequency of the updated frequency range interval, and the left boundary frequency of the frequency range interval before the update is determined as the left boundary frequency of the updated frequency range interval. For example, the updated frequency range interval could be [f L ,f M ], where f M f is the right boundary frequency of the updated frequency range, which is also the first median frequency of the original frequency range. L The left boundary frequency of the frequency range interval before the update is given. The second median frequency of the updated frequency range interval is determined based on the following formula:

[0078] f T =f M +0.618(fR -f M );

[0079] Among them, f T f is the second median frequency of the updated frequency range. R f is the right boundary frequency of the frequency range before the update. M This is the right boundary frequency of the updated frequency range.

[0080] It should be noted that the first median frequency of the updated frequency range interval is equal to the second median frequency of the original frequency range interval. In this embodiment, the error value of the channel is determined based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, as follows: the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel are determined based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel; the error value of the channel is then determined based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel.

[0081] In this embodiment, the error value includes at least one of the following: mismatch error value, gain mismatch error value, sampling time mismatch error value, and bandwidth mismatch error value.

[0082] In practical data sampling systems, the frequency of the input target sine wave is often unknown or inaccurate during channel error estimation. This unknown or inaccurate frequency severely affects the accuracy of the fitted gain, DC offset, and phase. Consequently, the sampling time mismatch and bandwidth mismatch calculated from the frequency, gain, DC offset, and phase will be severely distorted, impacting the performance of the data acquisition system. This invention, based on the three-parameter sine fitting algorithm, adds a frequency iterative search algorithm. The final target frequency is infinitely close to the true frequency of the input target sine wave, thus obtaining the optimal target gain, target DC offset, and target phase. Then, based on the target frequency, target gain, target DC offset, and target phase of the target sine wave corresponding to each channel, the error value of each channel is determined.

[0083] Optionally, the parameters include: the effective value of the fitted residual, and the iteration stopping condition includes:

[0084] The absolute value of the difference between the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency, and the ratio of the effective value of the fitting residual of the second median frequency, are greater than or equal to a set ratio threshold.

[0085] In this embodiment, the iteration stopping condition is:

[0086] |(ρM -ρ T ) / ρ T | <h e ;

[0087] Where, ρ M ρ is the effective value of the fitted residual for the first median frequency. T h is the effective value of the fitted residual for the second median frequency. e To set a ratio threshold.

[0088] Optionally, the parameters may also include: gain, phase, and DC offset;

[0089] The error value of the channel is determined based on the parameters corresponding to the first median frequency and the second median frequency within the frequency range corresponding to the channel, including:

[0090] Based on the parameters of the first median frequency and the second median frequency in the frequency range corresponding to the channel, the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel are determined.

[0091] In this embodiment, the method for determining the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel based on the parameters of the first median frequency and the second median frequency within the frequency range corresponding to the channel can be as follows: If the effective value of the fitting residual of the first median frequency corresponding to the channel is less than the effective value of the fitting residual of the second median frequency, then the first median frequency corresponding to the channel is determined as the target frequency of the target sine wave corresponding to the channel, the gain of the first median frequency corresponding to the channel is determined as the target gain of the target sine wave corresponding to the channel, and the phase of the first median frequency corresponding to the channel is determined as the target phase of the target sine wave corresponding to the channel. The DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the target sine wave corresponding to the channel; if the effective value of the fitting residual of the first median frequency corresponding to the channel is greater than or equal to the effective value of the fitting residual of the second median frequency, then the second median frequency corresponding to the channel is determined as the target frequency of the target sine wave corresponding to the channel, the gain of the second median frequency corresponding to the channel is determined as the target gain of the target sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the target sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the target sine wave corresponding to the channel.

[0092] It should be noted that, based on the initial frequency of the target sine wave and the frequency offset value corresponding to each channel in the data acquisition system, the frequency range interval corresponding to each channel is determined. Based on the three-parameter sine fitting algorithm, the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel is iteratively executed until the parameters of the first median frequency and the second median frequency satisfy the iteration stopping condition. Based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel are determined.

[0093] The error value of the channel is determined based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel.

[0094] In this embodiment, the error value of a channel can be determined based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel. The mismatch value of the channel is determined based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel. The difference between the mismatch value of the channel and the mismatch value of the reference channel is determined as the error value of the channel.

[0095] Optionally, based on the parameters of the first median frequency and the second median frequency within the frequency range corresponding to the channel, the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel are determined, including:

[0096] If the effective value of the fitting residual of the first median frequency corresponding to the channel is less than the effective value of the fitting residual of the second median frequency, then the first median frequency corresponding to the channel is determined as the target frequency of the target sine wave corresponding to the channel, the gain of the first median frequency corresponding to the channel is determined as the target gain of the target sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the target sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the target sine wave corresponding to the channel.

[0097] If the effective value of the fitting residual of the first median frequency corresponding to the channel is greater than or equal to the effective value of the fitting residual of the second median frequency, then the second median frequency corresponding to the channel is determined as the target frequency of the target sine wave corresponding to the channel, the gain of the second median frequency corresponding to the channel is determined as the target gain of the target sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the target sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the target sine wave corresponding to the channel.

[0098] In a specific example, within the frequency range [f0-Δf] corresponding to the channel max ,f0+Δf max Within [the framework], the absolute convergence of the four-parameter sine fitting process is achieved using a three-parameter sine fitting algorithm.

[0099] The four-parameter fitting process is as follows:

[0100] (1) Set the iteration stopping condition as |(ρ M -ρ T ) / ρ T | <h e ;

[0101] (2) Perform a Fourier transform on the target sine wave to obtain its initial frequency f0; input the target sine wave into the data acquisition system, and determine the total number of sampling points n and the total sampling time t for each channel based on the sampling data output from each channel in the data acquisition system; determine the frequency offset value of the target sine wave corresponding to each channel based on the total number of sampling points and the total sampling time for each channel. Δf max =v / n,

[0102] (3) Determine the frequency range of each channel based on the initial frequency and frequency offset value of the target sine wave corresponding to each channel. Left boundary frequency: Right boundary frequency: Median frequency: f M =f L +0.618(f R -f L );f T =f R -0.618(f R -f L ).

[0103] (4) In f M Perform a three-parameter sine curve fitting operation to obtain: A M θ M C M ρ M ; in f T Perform a three-parameter sine curve fitting to obtain A. T θ T C T ρ T .

[0104] (5) If ρ M <ρ T Then ρ = ρ M , there exists f0∈[f T ,fR ], f T The second median frequency of the frequency range before the update; f R The frequency range before the update is the right boundary frequency; the frequency range after the update is f. M =f T +0.618(f R -f T The second median frequency of the updated frequency range is equal to the first median frequency of the original frequency range. If ρ M >ρ T Then ρ = ρ T , there exists f0∈[f L ,f M The updated frequency range interval f T =f M -0.618(f R -f M The first median frequency of the updated frequency range interval is equal to the second median frequency of the original frequency range interval.

[0105] (6) Determine whether |(ρ) M -v T ) / ρ T | <h e If yes, then stop the iteration. If no, then repeat steps (3) to (5). If ρ M <ρ T Then ρ = ρ M The parameters of the four-parameter fitted sine curve are obtained as A = A M θ=θ M C=C M f = f M ρ M The fitting process ends. If ρ M >ρ T Then ρ = ρ T The parameters of the four-parameter fitted sine curve are obtained as A = A T θ=θ T C=C T f = f T ρ T The fitting process is now complete.

[0106] Optionally, the data acquisition system includes: a reference channel;

[0107] Based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel, determine the channel's error value, including:

[0108] The mismatch value of the channel is determined based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel.

[0109] In this embodiment, the mismatch value includes at least one of the following: gain mismatch value, offset mismatch value, bandwidth mismatch value, and sampling time mismatch value. It should be noted that the sampling time mismatch value is the time delay caused by sampling time mismatch. The bandwidth mismatch value is the time delay caused by bandwidth mismatch.

[0110] In data acquisition systems, channel mismatch mainly includes the following types of mismatch, and each type of mismatch affects different performance indicators as follows:

[0111] (1) Gain mismatch

[0112] 1) Gain differences can lead to amplitude errors, thereby reducing overall measurement accuracy;

[0113] 2) Gain mismatch may increase spurious components, thereby reducing the effective dynamic range;

[0114] 3) It affects the linear response of the system, causing nonlinear distortion.

[0115] (2) Mismatch (DC offset mismatch)

[0116] 1) Offset mismatch introduces a constant DC error;

[0117] 2) Increased noise levels in the low-frequency range;

[0118] 3) Because the signal is not properly zeroed, it may affect the smallest detectable signal change.

[0119] (3) Sampling time offset

[0120] 1) High-frequency signals are more severely affected by time mismatch, which leads to spectral broadening;

[0121] 2) It causes signal distortion or blurring in the time domain, affecting the quality of signal reconstruction.

[0122] (4) Phase mismatch

[0123] 1) Affects the phase measurement of the signal, especially in multi-channel analysis;

[0124] 2) It leads to an increase in harmonic distortion, thereby affecting the signal's spectrum analysis;

[0125] 3) It has a negative impact on applications that require precise phase information (such as communication demodulation).

[0126] In this embodiment, the method for determining the channel mismatch value based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel can be as follows: the target gain of the target sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the target sine wave corresponding to the channel is determined as the offset mismatch value of the channel; the bandwidth mismatch value of the channel is determined based on the target frequency and target gain of the target sine wave corresponding to the channel; and the difference between the target phase of the target sine wave corresponding to the channel and the bandwidth mismatch value is determined as the sampling time mismatch value of the channel.

[0127] The difference between the mismatch value of the channel and the mismatch value of the reference channel is determined as the channel error value.

[0128] In this embodiment, the difference between the channel mismatch value and the reference channel mismatch value can be used to determine the channel error value as follows: the difference between the channel gain mismatch value and the reference channel gain mismatch value can be used to determine the channel gain mismatch error value; the difference between the channel offset mismatch value and the reference channel offset mismatch value can be used to determine the channel offset mismatch error value; the difference between the channel bandwidth mismatch value and the reference channel bandwidth mismatch value can be used to determine the channel bandwidth mismatch error value; and the difference between the channel sampling time mismatch value and the reference channel sampling time mismatch value can be used to determine the channel sampling time mismatch error value.

[0129] In this embodiment, the error value includes at least one of the following: mismatch error value, gain mismatch error value, sampling time mismatch error value, and bandwidth mismatch error value.

[0130] In a specific example, using a four-channel oscilloscope, channel 1 is used as the reference channel, and channels 2, 3, and 4 are used as measurement channels. A four-parameter sine fitting algorithm is used to estimate the target gain, target phase, target DC offset, and target frequency of the target sine wave for each channel. Then, based on the target gain, target phase, target DC offset, and target frequency of the target sine wave for each channel, the mismatch value for each channel is calculated, resulting in the offset mismatch value, gain mismatch value, sampling time mismatch value, and bandwidth mismatch value for each channel. The offset mismatch error value of the measurement channel is determined by the difference between the offset mismatch value of the measurement channel and the offset mismatch value of the reference channel; the gain mismatch error value of the measurement channel is determined by the difference between the gain mismatch value of the measurement channel and the gain mismatch value of the reference channel; the sampling time mismatch error value of the measurement channel is determined by the difference between the sampling time mismatch value of the measurement channel and the sampling time mismatch value of the reference channel; and the bandwidth mismatch error value of the measurement channel is determined by the difference between the bandwidth mismatch value of the measurement channel and the bandwidth mismatch value of the reference channel.

[0131] Optionally, the channel mismatch value is determined based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel, including:

[0132] The target gain of the target sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the target sine wave corresponding to the channel is determined as the mismatch value of the channel.

[0133] In this embodiment, if the target sine wave is a low-frequency sine wave, the target gain of the low-frequency sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the low-frequency sine wave corresponding to the channel is determined as the offset mismatch value of the channel. If the target sine wave is a high-frequency sine wave, the target gain of the high-frequency sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the high-frequency sine wave corresponding to the channel is determined as the offset mismatch value of the channel. If the target sine wave includes both low-frequency and high-frequency sine waves, the target gain of the low-frequency sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the low-frequency sine wave corresponding to the channel is determined as the offset mismatch value of the channel.

[0134] The bandwidth mismatch value of the channel is determined based on the target frequency and target gain of the target sine wave corresponding to the channel.

[0135] In this embodiment, if the target sine wave is a high-frequency sine wave, the bandwidth mismatch value of the channel is determined based on the target frequency and target gain of the high-frequency sine wave corresponding to the channel. If the target sine wave is a low-frequency sine wave, the bandwidth mismatch value of the channel is determined based on the target frequency and target gain of the low-frequency sine wave corresponding to the channel. If the target sine wave includes both low-frequency and high-frequency sine waves, the bandwidth mismatch value of the channel is determined based on the target frequency and target gain of the high-frequency sine wave corresponding to the channel.

[0136] The difference between the target phase of the target sine wave corresponding to the channel and the bandwidth mismatch value is determined as the sampling time mismatch value of the channel.

[0137] In this embodiment, if the target sine wave is a high-frequency sine wave, the difference between the target phase of the high-frequency sine wave corresponding to the channel and the bandwidth mismatch value is determined as the sampling time mismatch value of the channel. If the target sine wave is a low-frequency sine wave, the difference between the target phase of the low-frequency sine wave corresponding to the channel and the bandwidth mismatch value is determined as the sampling time mismatch value of the channel. If the target sine wave includes both low-frequency and high-frequency sine waves, the difference between the target phase of the high-frequency sine wave corresponding to the channel and the bandwidth mismatch value is determined as the sampling time mismatch value of the channel.

[0138] Optionally, the bandwidth mismatch value of the channel is determined based on the target frequency and target gain of the target sine wave corresponding to the channel, including:

[0139] The bandwidth of the sample-and-hold circuit in the channel is determined based on the target gain of the target sine wave corresponding to the channel and the angular frequency of the target frequency of the target sine wave.

[0140] In this embodiment, the bandwidth of the sample-and-hold circuit in the channel can be determined based on the target gain of the target sine wave corresponding to the channel and the angular frequency corresponding to the target frequency of the target sine wave, using the formula: ω in =2πf / v, obtain the angular frequency corresponding to the target frequency of the target sine wave corresponding to the channel, and determine the bandwidth of the sample-and-hold circuit in the channel based on the target gain of the target sine wave corresponding to the channel and the angular frequency corresponding to the target frequency of the target sine wave.

[0141] The bandwidth mismatch value of the channel is determined based on the angular frequency of the target sine wave corresponding to the channel and the bandwidth of the sample-and-hold circuit in the channel.

[0142] In this embodiment, the arctangent of the ratio of the angular frequency of the target frequency of the target sine wave corresponding to the channel to the bandwidth of the sample-and-hold circuit in the channel is determined as the bandwidth mismatch value of the channel.

[0143] Optionally, the bandwidth of the sample-and-hold circuit in the channel is determined based on the target gain of the target sine wave corresponding to the channel and the angular frequency corresponding to the target frequency of the target sine wave, including:

[0144] Based on the target gain of the target sine wave corresponding to the channel, the angular frequency of the target sine wave, and a first formula, the bandwidth of the sample-and-hold circuit in the channel is determined, wherein the first formula is:

[0145]

[0146] Where, ω in ω is the angular frequency of the target frequency corresponding to the channel. b For the bandwidth of the sample-and-hold circuit in the channel, g i The target gain is the target sine wave corresponding to the channel.

[0147] Optionally, the bandwidth mismatch value of the channel is determined based on the angular frequency corresponding to the target frequency of the target sine wave and the bandwidth of the sample-and-hold circuit in the channel, including:

[0148] Based on the angular frequency of the target sine wave corresponding to the channel, the bandwidth of the sample-and-hold circuit in the channel, and a second formula, the bandwidth mismatch value of the channel is determined, wherein the second formula is:

[0149]

[0150] Where, θ fBW This represents the bandwidth mismatch value.

[0151] In this embodiment, the arctangent of the ratio of the angular frequency of the target frequency of the target sine wave corresponding to the channel to the bandwidth of the sample-and-hold circuit in the channel is determined as the bandwidth mismatch value of the channel.

[0152] In a specific example, taking a two-channel data acquisition system as an example, channel 0 is used as the reference channel and channel 1 is used as the measurement channel.

[0153] After a four-parameter fitting and iterative process, the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to channel 0 are obtained:

[0154] The target gain g0 = A0 of the target sine wave corresponding to channel 0, the target phase θ0 = θ0 of the target sine wave corresponding to channel 0, the target DC offset o0 = C0 of the target sine wave corresponding to channel 0, and the target frequency f0 = f of the target sine wave corresponding to channel 0.

[0155] Target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to channel 1:

[0156] The target gain g1 = A1 of the target sine wave corresponding to channel 1, the target phase θ1 = θ1 of the target sine wave corresponding to channel 1, and the target DC offset o1 = C1 of the target sine wave corresponding to channel 1. Since the input sine wave remains unchanged, the target frequency of the target sine wave corresponding to channel 1 is still f1 = f.

[0157] With channel 0 as the reference channel, the offset mismatch error value, gain mismatch error value, sampling time mismatch error value, and bandwidth mismatch error value of channel 1 can all be determined by the difference between the gain mismatch value, offset mismatch value, bandwidth mismatch value, and sampling time mismatch value of channel 1 and the corresponding gain mismatch value, offset mismatch value, bandwidth mismatch value, and sampling time mismatch value of channel 0.

[0158] In this embodiment, the gain mismatch and offset mismatch values ​​can be determined by inputting a low-frequency sine wave f into the data acquisition system. lowIterative fitting estimation is used to obtain the results because gain mismatch and offset mismatch are low-frequency types of mismatch, and the mismatch value does not change with the frequency of the input sine wave. For sampling time mismatch, at the same sampling clock, the higher the frequency of the input sine wave, the greater the error caused by sampling time mismatch. Therefore, for high-frequency sine signals f... high Estimating the sampling time and bandwidth mismatch under the given input conditions can better improve the performance of the data acquisition system.

[0159] Under low-frequency sinusoidal signal input conditions, the gain mismatch error value of the measurement channel is obtained based on the difference between the gain mismatch value of the reference channel and the gain mismatch value of the measurement channel; the offset mismatch error value of the measurement channel is obtained based on the difference between the offset mismatch value of the reference channel and the offset mismatch value of the measurement channel. Distinguishing between the time delay of sampling time mismatch and the time delay of bandwidth mismatch at high frequencies is possible because bandwidth mismatch is related to the input signal frequency; the higher the input frequency, the greater the bandwidth mismatch. Therefore, at high frequencies, it is easier to distinguish the time delay in sampling time mismatch and bandwidth mismatch. In this embodiment of the invention, a determined high-frequency sine wave is input to the data acquisition system. Through a four-parameter fitting iterative process, the target gain, target phase, and target frequency of the high-frequency sine wave corresponding to each channel are obtained. The sampling time mismatch value and bandwidth mismatch value of each channel are determined. Then, the sampling time mismatch error value of the measurement channel is obtained based on the difference between the sampling time mismatch value of the reference channel and the sampling time mismatch value of the measurement channel; the bandwidth mismatch error value of the measurement channel is obtained based on the difference between the bandwidth mismatch value of the reference channel and the bandwidth mismatch value of the measurement channel.

[0160] In a data acquisition system, the transfer function of the sample-and-hold circuit in the channel can be:

[0161]

[0162] Where A0 is the target gain of the target sine wave, ω b This determines the bandwidth of the sample-and-hold circuit in the channel.

[0163] The bandwidth ω of the sample-and-hold circuit in the channel b Represented as:

[0164]

[0165] Among them, g i ω represents the target gain value of the target sine wave corresponding to the channel. in The angular frequency is the target frequency of the target sine wave corresponding to the channel.

[0166] According to ω b The bandwidth mismatch value of the channel is calculated.

[0167]

[0168] By using four-parameter fitting, the target phase of the target sine wave corresponding to the channel is estimated to include sampling time mismatch and bandwidth mismatch values. Therefore, we can obtain:

[0169]

[0170] Where, θ s This represents the sampling time mismatch value for the channel.

[0171] The technical solution provided in this embodiment can estimate the error value of the channel in the data acquisition system, thereby improving the consistency and accuracy of the multi-channel data acquisition system.

[0172] It should be noted that channel mismatch calibration is crucial for ensuring high-precision and high-performance data acquisition, mainly in the following aspects:

[0173] High-frequency digital signal analysis: Any channel mismatch will lead to incorrect signal demodulation or spectrum analysis.

[0174] Wideband signal acquisition: For applications that require the acquisition of signals over a wide frequency range, such as ultra-wideband radio frequency signals, channel mismatches must be corrected for accurate analysis in the frequency domain.

[0175] High-precision measurement: The equipment requires extremely high measurement accuracy, and even slight errors can have a significant impact on the results.

[0176] Multi-sensor measurement: In systems such as radar and sonar that require multiple sensors to work together, any mismatch between channels will affect the accuracy of signal fusion.

[0177] Semiconductor test equipment: It needs to process a large number of parallel signals, and inter-channel calibration can ensure the reliability and consistency of measurements.

[0178] The technical solution of this embodiment determines the frequency range interval corresponding to each channel based on the initial frequency of the target sine wave and the frequency offset value corresponding to each channel in the data acquisition system; based on a three-parameter sine fitting algorithm, it determines the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to each channel; for any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, the frequency range interval corresponding to the channel is updated according to the parameters of the first median frequency and the second median frequency corresponding to the channel, and the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel is returned based on the updated frequency range interval corresponding to the channel, until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition, and the error value of the channel is determined according to the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel. This approach addresses the issue that the frequency of the input target sine wave is often unknown or inaccurate. Unknown or inaccurate frequencies severely affect the accuracy of the fitted gain, DC offset, and phase, leading to significant distortion in the sampling time and bandwidth mismatches calculated from these parameters, thus impacting the performance of the data acquisition system. The method uses an iterative algorithm to obtain the target frequency, target gain, target DC offset, and target phase of the target sine wave for each channel. Based on these parameters, the error value for each channel is determined, thereby improving the performance of the data acquisition system.

[0179] Example 2

[0180] Figure 2 A flowchart of another method for determining the channel error value in a data acquisition system provided by an embodiment of the present invention is shown below. Figure 2 As shown, the method specifically includes the following steps:

[0181] S201 performs a Fourier transform on the low-frequency sine wave to obtain the initial frequency of the low-frequency sine wave.

[0182] S202 inputs a low-frequency sine wave into the data acquisition system to obtain the total number of sampling points and the total sampling time for each channel.

[0183] S203, determine the frequency offset value corresponding to each channel based on the total number of sampling points and the total sampling time for each channel.

[0184] S204. Based on the initial frequency of the low-frequency sine wave and the frequency offset value corresponding to each channel, determine the frequency range interval corresponding to each channel.

[0185] S205, based on a three-parameter sine fitting algorithm, determines the gain, phase, DC offset, and effective value of the fitting residual at the first median frequency in the frequency range corresponding to each channel, and the gain, phase, DC offset, and effective value of the fitting residual at the second median frequency.

[0186] S206, for any channel, determine whether the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency satisfy the iteration stopping condition. If they satisfy the condition, proceed to S207. If they do not satisfy the condition, update the frequency range interval corresponding to the channel according to the parameters of the first median frequency and the second median frequency. Based on the updated frequency range interval corresponding to the channel, return to S205 until the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency satisfy the iteration stopping condition.

[0187] S207, determine whether the effective value of the fitting residual of the first median frequency is less than the effective value of the fitting residual of the second median frequency. If yes, execute S208; otherwise, execute S209.

[0188] S208, the first median frequency corresponding to the channel is determined as the target frequency of the low-frequency sine wave corresponding to the channel, the gain of the first median frequency corresponding to the channel is determined as the target gain of the low-frequency sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the low-frequency sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the low-frequency sine wave corresponding to the channel.

[0189] S209, the second median frequency corresponding to the channel is determined as the target frequency of the low-frequency sine wave corresponding to the channel, the gain of the second median frequency corresponding to the channel is determined as the target gain of the low-frequency sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the low-frequency sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the low-frequency sine wave corresponding to the channel.

[0190] S210, the target gain of the low-frequency sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the low-frequency sine wave corresponding to the channel is determined as the mismatch value of the channel.

[0191] S211 performs a Fourier transform on the high-frequency sine wave to obtain the initial frequency of the high-frequency sine wave.

[0192] S212 inputs a high-frequency sine wave into the data acquisition system to obtain the total number of sampling points and the total sampling time for each channel.

[0193] S213, determine the frequency offset value corresponding to each channel based on the total number of sampling points and the total sampling time for each channel.

[0194] S214. Based on the initial frequency of the high-frequency sine wave and the frequency offset value corresponding to each channel, determine the frequency range interval corresponding to each channel.

[0195] S215, based on a three-parameter sine fitting algorithm, determines the effective values ​​of gain, phase, DC offset, and fitting residual corresponding to the first median frequency in the frequency range interval corresponding to each channel, and the effective values ​​of gain, phase, DC offset, and fitting residual corresponding to the second median frequency.

[0196] S216, for any channel, determine whether the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency satisfy the iteration stopping condition. If they satisfy the condition, proceed to S217. If they do not satisfy the condition, update the frequency range interval corresponding to the channel according to the parameters of the first median frequency and the second median frequency. Based on the updated frequency range interval corresponding to the channel, return to S215 until the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency satisfy the iteration stopping condition.

[0197] In this embodiment, the method of updating the frequency range interval corresponding to the channel based on the parameters of the first median frequency and the second median frequency can be as follows: If the effective value of the fitting residual of the first median frequency is less than the effective value of the fitting residual of the second median frequency, then the second median frequency of the frequency range interval before the update is determined as the left boundary frequency of the updated frequency range interval, and the right boundary frequency of the frequency range interval before the update is determined as the right boundary frequency of the updated frequency range interval; if the effective value of the fitting residual of the first median frequency is greater than the effective value of the fitting residual of the second median frequency, then the first median frequency of the frequency range interval before the update is determined as the right boundary frequency of the updated frequency range interval, and the left boundary frequency of the frequency range interval before the update is determined as the left boundary frequency of the updated frequency range interval.

[0198] It should be noted that if the effective value of the fitted residual of the first median frequency is less than the effective value of the fitted residual of the second median frequency, then the first median frequency of the updated frequency range interval is determined based on the following formula:

[0199] f M =f T +0.618(f R -f T );

[0200] Among them, f M f is the first median frequency of the updated frequency range. Rf is the right boundary frequency of the frequency range before the update. T This is the left boundary frequency of the updated frequency range.

[0201] It should be noted that the second median frequency of the updated frequency range is equal to the first median frequency of the original frequency range.

[0202] If the effective value of the fitted residual for the first median frequency is greater than the effective value of the fitted residual for the second median frequency, then the second median frequency of the updated frequency range interval is determined based on the following formula:

[0203] f T =f M +0.618(f R -f M );

[0204] Among them, f T f is the second median frequency of the updated frequency range. R f is the right boundary frequency of the frequency range before the update. M This is the right boundary frequency of the updated frequency range.

[0205] It should be noted that the first median frequency of the updated frequency range is equal to the second median frequency of the original frequency range.

[0206] S217, determine whether the effective value of the fitting residual of the first median frequency is less than the effective value of the fitting residual of the second median frequency. If yes, execute S218; otherwise, execute S219.

[0207] S218, the first median frequency corresponding to the channel is determined as the target frequency of the high-frequency sine wave corresponding to the channel, the gain of the first median frequency corresponding to the channel is determined as the target gain of the high-frequency sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the high-frequency sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the high-frequency sine wave corresponding to the channel.

[0208] S219, the second median frequency corresponding to the channel is determined as the target frequency of the high-frequency sine wave corresponding to the channel, the gain of the second median frequency corresponding to the channel is determined as the target gain of the high-frequency sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the high-frequency sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the high-frequency sine wave corresponding to the channel.

[0209] S220 determines the bandwidth mismatch value of the channel based on the target frequency and target gain of the high-frequency sine wave corresponding to the channel.

[0210] In this embodiment, the bandwidth mismatch value of the channel can be determined based on the target frequency and target gain of the high-frequency sine wave corresponding to the channel as follows: the bandwidth of the sample-and-hold circuit in the channel is determined based on the target gain and the angular frequency of the target frequency of the high-frequency sine wave corresponding to the channel; the bandwidth mismatch value of the channel is determined based on the angular frequency of the target frequency of the high-frequency sine wave corresponding to the channel and the bandwidth of the sample-and-hold circuit in the channel.

[0211] In this embodiment, the bandwidth of the sample-and-hold circuit in the channel is determined based on the target gain and angular frequency of the target frequency of the high-frequency sine wave corresponding to the channel, including:

[0212] Based on the target gain of the high-frequency sine wave corresponding to the channel, the angular frequency of the target frequency of the high-frequency sine wave, and the first formula, the bandwidth of the sample-and-hold circuit in the channel is determined, wherein the first formula is:

[0213]

[0214] Where, ω in ω is the angular frequency corresponding to the target frequency of the high-frequency sine wave. b For the bandwidth of the sample-and-hold circuit in the channel, g i The target gain for a high-frequency sine wave.

[0215] In this embodiment, the bandwidth mismatch value of the channel is determined based on the angular frequency of the target frequency of the high-frequency sine wave corresponding to the channel and the bandwidth of the sample-and-hold circuit in the channel, including:

[0216] Based on the angular frequency of the target frequency of the high-frequency sine wave corresponding to the channel, the bandwidth of the sample-and-hold circuit in the channel, and the second formula, the bandwidth mismatch value of the channel is determined, wherein the second formula is:

[0217]

[0218] in, This represents the bandwidth mismatch value of the channel.

[0219] S221, the difference between the target phase of the high-frequency sine wave corresponding to the channel and the bandwidth mismatch value is determined as the sampling time mismatch value of the channel.

[0220] S222, the difference between the gain mismatch value of the channel and the gain mismatch value of the reference channel is determined as the gain mismatch error value of the channel; the difference between the offset mismatch value of the channel and the offset mismatch value of the reference channel is determined as the offset mismatch error value of the channel; the difference between the bandwidth mismatch value of the channel and the bandwidth mismatch value of the reference channel is determined as the bandwidth mismatch error value of the channel; the difference between the sampling time mismatch value of the channel and the sampling time mismatch value of the reference channel is determined as the sampling time mismatch error value of the channel.

[0221] In existing technologies, the calibration of data acquisition system outputs mainly includes two methods. One method directly measures the error caused by clock delay mismatch at each channel output point of the data acquisition system using test signals, and then uses this error to calibrate the output sampling points of the data acquisition system. This direct measurement method requires multiple measurement calibrations for different mismatch parameters, and different test signals need to be input for different mismatches. Furthermore, the input test signals need to be of high quality; otherwise, the accuracy of the measured mismatch value will be limited. Direct measurement of channel mismatch also has some significant engineering drawbacks, such as:

[0222] (1) Additional precision measurement equipment and reference signal sources are required to ensure accurate measurement of each ADC channel. This increases the complexity of the overall system.

[0223] (2) Each channel must be measured and adjusted individually, which can be very time-consuming in a multi-channel system.

[0224] (3) The measurement process is easily affected by environmental factors such as temperature changes and electromagnetic interference, which can lead to errors; the accuracy and stability of the measuring equipment itself will directly affect the accuracy of the measurement results.

[0225] (4) Direct measurements are usually performed under static or low-frequency conditions, which cannot fully characterize the characteristics and requirements of high-speed dynamic signal processing.

[0226] (5) High-precision measuring instruments and standard signal sources are required, which increases the cost of equipment investment.

[0227] (6) It requires specialized technical personnel to operate and analyze, resulting in high labor costs.

[0228] Another approach is to estimate the mismatch between channels by fitting a sine curve. Traditional sine curve fitting algorithms use a three-parameter fitting algorithm for amplitude, phase, and DC offset. This method requires the frequency of the sine wave input to the data acquisition system to be known, which places higher demands on the sine wave in practical applications. It is difficult to guarantee the ideal accuracy of the sine wave frequency, which ultimately leads to low accuracy of the estimated mismatch value.

[0229] In this embodiment of the invention, to address the above-mentioned problems, a frequency parameter is added to the original three-parameter sine fitting estimation algorithm. The target frequency that is closest to the true frequency of the target sine wave is obtained through iteration, along with the target gain, target DC offset, and target phase. Then, the error value of each channel is determined based on the target gain, target DC offset, target frequency, and target phase corresponding to each channel.

[0230] Compared with the three-parameter sinusoidal fitting calibration algorithm, the four-parameter calibration algorithm has the following advantages:

[0231] 1. The four-parameter fitting includes gain, phase, DC offset, and frequency, which can adapt to more diverse signal characteristics;

[0232] 2. An additional frequency parameter has been added, allowing direct estimation of the true frequency of a sine wave without relying on an external, precise signal source. This improves the tolerance for frequency errors.

[0233] 3. Four-parameter fitting can more accurately describe sine waves, and even if a sine wave with ideal accuracy cannot be provided, it can still obtain accurate mismatch values. This simplifies the calibration process and improves the level of automation.

[0234] 4. It does not rely on the absolute accuracy of sine waves, reducing the dependence on high-precision signal generation equipment in the production process, and is more robust to changes in environmental factors such as temperature and electromagnetic interference.

[0235] 5. Four-parameter fitting reduces the requirement for high-precision signal generators, thereby reducing equipment procurement and maintenance costs. Although the initial implementation complexity may increase slightly, it provides long-term stability and calibration accuracy, reducing subsequent maintenance expenses.

[0236] The method for determining the channel error value in the data acquisition system proposed in this invention has the following advantages: structurally, it does not rely on a high-precision signal source; functionally, it simplifies the entire channel mismatch estimation process and improves the level of automation; in terms of manufacturing, it is more robust to changes in factors such as temperature and electromagnetic interference; and in terms of cost, it requires less manual intervention and adjustment.

[0237] It should be noted that, to facilitate the analysis of the dynamic calibration process, the parameter values ​​of each iteration are stored, and a four-parameter sine fitting dynamic calibration graph is plotted. For example, the effective values ​​of the fitting residuals can be displayed. In practice, the real-time display of the effective values ​​of the fitting residuals clearly indicates the stopping position of the sine fitting iteration and the entire fitting process.

[0238] The embodiments of the present invention can be applied to the accurate estimation of channel mismatch in single-channel and multi-channel data acquisition systems, and based on this, channel mismatch correction can be performed on the data acquisition system to improve the performance of the data acquisition system.

[0239] Example 3

[0240] Figure 3 This is a schematic diagram of a channel error value determination device in a data acquisition system according to an embodiment of the present invention. This embodiment is applicable to the determination of channel error values ​​in a data acquisition system. The device can be implemented using software and / or hardware, and can be integrated into any device that provides channel error value determination functionality in a data acquisition system, such as… Figure 3 As shown, the channel error value determination device in the data acquisition system specifically includes: a frequency range interval determination module 310, a parameter determination module for the first median frequency and a parameter determination module for the second median frequency 320, and a channel error value determination module 330.

[0241] The frequency range interval determination module is used to determine the frequency range interval corresponding to each channel based on the initial frequency of the target sine wave and the frequency offset value corresponding to each channel in the data acquisition system. The target sine wave includes low-frequency sine waves and / or high-frequency sine waves.

[0242] The first median frequency parameter and second median frequency parameter determination module is used to determine the first median frequency parameter and second median frequency parameter in the frequency range interval corresponding to each channel based on a three-parameter sine fitting algorithm.

[0243] The channel error value determination module is used to, for any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, update the frequency range interval corresponding to the channel according to the parameters of the first median frequency and the second median frequency corresponding to the channel, and return to execute the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel based on the updated frequency range interval corresponding to the channel, until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition, and determine the error value of the channel according to the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel.

[0244] The above-described products can perform the methods provided in any embodiment of the present invention, and have the corresponding functional modules and beneficial effects for performing the methods.

[0245] Example 4

[0246] Figure 4A schematic diagram of an electronic device 10 that can be used to implement embodiments of the present invention is shown. The electronic device is 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. The electronic device can also represent various forms of mobile devices, such as personal digital processors, cellular phones, smartphones, wearable devices (e.g., helmets, glasses, watches, etc.), 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 invention described and / or claimed herein.

[0247] like Figure 4 As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 or a random access memory (RAM) 13, communicatively connected to the at least one processor 11. The memory stores computer programs executable by the at least one processor. The processor 11 can perform various appropriate actions and processes based on the computer program stored in the ROM 12 or loaded from storage unit 18 into the RAM 13. The RAM 13 may also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are interconnected via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0248] Multiple components in electronic device 10 are connected to I / O interface 15, including: input unit 16, such as keyboard, mouse, etc.; output unit 17, such as various types of displays, speakers, etc.; storage unit 18, such as disk, optical disk, etc.; and communication unit 19, such as network card, modem, wireless transceiver, etc. Communication unit 19 allows electronic device 10 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0249] Processor 11 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 11 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 processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 11 performs the various methods and processes described above, such as the error value determination method for a channel in a data acquisition system.

[0250] In some embodiments, the channel error value determination method in the data acquisition system can be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the channel error value determination method in the data acquisition system described above can be performed. Alternatively, in other embodiments, processor 11 can be configured to perform the channel error value determination method in the data acquisition system by any other suitable means (e.g., by means of firmware).

[0251] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), payload-programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.

[0252] Computer programs used to implement the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be performed. The computer programs may be executed entirely on a machine, partially on a machine, or as a standalone software package, partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

[0253] In the context of this invention, a computer-readable storage medium can be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, apparatus, or device. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0254] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0255] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or computing systems that include middleware components (e.g., application servers), or computing systems that include frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0256] A computing system can include clients and servers. Clients and servers are generally located far apart and typically interact through communication networks. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system to address the shortcomings of traditional physical hosts and VPS services, such as high management difficulty and weak business scalability.

[0257] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.

[0258] This invention also provides a computer program product, including a computer program that, when executed by a processor, implements the method for determining the error value of a channel in a data acquisition system according to any embodiment of the invention.

[0259] In implementing the computer program product, computer program code for performing the operations of this invention can be written in one or more programming languages ​​or a combination thereof. Programming languages ​​include object-oriented programming languages ​​such as Java, Smalltalk, and C++, as well as conventional procedural programming languages ​​such as C or similar languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0260] The specific embodiments described above do not constitute a limitation on the scope of protection of this invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this invention should be included within the scope of protection of this invention.

Claims

1. A method for determining the error value of a channel in a data acquisition system, characterized in that, include: Based on the initial frequency of the target sine wave and the frequency offset value of each channel in the data acquisition system, the frequency range interval corresponding to each channel is determined, wherein the target sine wave includes low-frequency sine waves and / or high-frequency sine waves. Based on the three-parameter sine fitting algorithm, the parameters of the first median frequency and the second median frequency in the frequency range corresponding to each channel are determined. For any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, then the frequency range interval corresponding to the channel is updated according to the parameters of the first median frequency and the second median frequency corresponding to the channel. Based on the updated frequency range interval corresponding to the channel, the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel is returned to be executed until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition. Based on the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel, the error value of the channel is determined.

2. The method according to claim 1, characterized in that, Before determining the frequency range for each channel based on the initial frequency of the target sine wave and the frequency offset values ​​corresponding to each channel in the data acquisition system, the following steps are also included: Perform a Fourier transform on the target sine wave to obtain the initial frequency of the target sine wave; The target sine wave is input into the data acquisition system to obtain the total number of sampling points and the total sampling time for each channel. The frequency offset value for each channel is determined based on the total number of sampling points and the total sampling time for each channel.

3. The method according to claim 1, characterized in that, The parameters include: the effective value of the fitted residual; the iteration stopping condition includes: The absolute value of the difference between the effective value of the fitting residual of the first median frequency and the effective value of the fitting residual of the second median frequency, and the ratio of the effective value of the fitting residual of the second median frequency, are greater than or equal to a set ratio threshold.

4. The method according to claim 3, characterized in that, The parameters also include: gain, phase, and DC offset; The error value of the channel is determined based on the parameters corresponding to the first median frequency and the second median frequency within the frequency range corresponding to the channel, including: Based on the parameters of the first median frequency and the second median frequency in the frequency range corresponding to the channel, determine the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel; The error value of the channel is determined based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel.

5. The method according to claim 4, characterized in that, Based on the parameters of the first median frequency and the second median frequency within the frequency range corresponding to the channel, the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel are determined, including: If the effective value of the fitting residual of the first median frequency corresponding to the channel is less than the effective value of the fitting residual of the second median frequency, then the first median frequency corresponding to the channel is determined as the target frequency of the target sine wave corresponding to the channel, the gain of the first median frequency corresponding to the channel is determined as the target gain of the target sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the target sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the target sine wave corresponding to the channel. If the effective value of the fitting residual of the first median frequency corresponding to the channel is greater than or equal to the effective value of the fitting residual of the second median frequency, then the second median frequency corresponding to the channel is determined as the target frequency of the target sine wave corresponding to the channel, the gain of the second median frequency corresponding to the channel is determined as the target gain of the target sine wave corresponding to the channel, the phase of the first median frequency corresponding to the channel is determined as the target phase of the target sine wave corresponding to the channel, and the DC offset of the first median frequency corresponding to the channel is determined as the target DC offset of the target sine wave corresponding to the channel.

6. The method according to claim 4, characterized in that, The data acquisition system includes: a reference channel; Based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel, determine the channel's error value, including: Based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel, the mismatch value of the channel is determined, wherein the mismatch value includes at least one of the following: gain mismatch value, offset mismatch value, bandwidth mismatch value, and sampling time mismatch value; The difference between the mismatch value of the channel and the mismatch value of the reference channel is determined as the channel error value, wherein the error value includes at least one of the following: mismatch error value, gain mismatch error value, sampling time mismatch error value, and bandwidth mismatch error value.

7. The method according to claim 6, characterized in that, Based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel, determine the channel mismatch value, including: The target gain of the target sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the target sine wave corresponding to the channel is determined as the mismatch value of the channel. The bandwidth mismatch value of the channel is determined based on the target frequency and target gain of the target sine wave corresponding to the channel. The difference between the target phase of the target sine wave corresponding to the channel and the bandwidth mismatch value is determined as the sampling time mismatch value of the channel.

8. The method according to claim 7, characterized in that, Based on the target frequency and target gain of the target sine wave corresponding to the channel, determine the bandwidth mismatch value of the channel, including: The bandwidth of the sample-and-hold circuit in the channel is determined based on the target gain of the target sine wave corresponding to the channel and the angular frequency of the target sine wave. The bandwidth mismatch value of the channel is determined based on the angular frequency of the target sine wave corresponding to the channel and the bandwidth of the sample-and-hold circuit in the channel.

9. The method according to claim 8, characterized in that, Based on the target gain of the target sine wave corresponding to the channel and the angular frequency corresponding to the target frequency of the target sine wave, the bandwidth of the sample-and-hold circuit in the channel is determined, including: Based on the target gain of the target sine wave corresponding to the channel, the angular frequency of the target sine wave, and a first formula, the bandwidth of the sample-and-hold circuit in the channel is determined, wherein the first formula is: Where, ω in ω is the angular frequency corresponding to the target frequency. b For the bandwidth of the sample-and-hold circuit in the channel, g i The target gain is the target sine wave.

10. The method according to claim 8, characterized in that, Based on the angular frequency of the target sine wave corresponding to the channel and the bandwidth of the sample-and-hold circuit in the channel, the bandwidth mismatch value of the channel is determined, including: Based on the angular frequency of the target sine wave corresponding to the channel, the bandwidth of the sample-and-hold circuit in the channel, and a second formula, the bandwidth mismatch value of the channel is determined, wherein the second formula is: in, This represents the bandwidth mismatch value of the channel.

11. The method according to claim 6, characterized in that, The target sine wave includes low-frequency sine waves and high-frequency sine waves; Accordingly, based on the target gain, target DC offset, target frequency, and target phase of the target sine wave corresponding to the channel, the channel mismatch value is determined, including: The target gain of the low-frequency sine wave corresponding to the channel is determined as the gain mismatch value of the channel, and the target DC offset of the low-frequency sine wave corresponding to the channel is determined as the mismatch value of the channel. Determine the bandwidth mismatch value of the channel based on the target frequency and target gain of the high-frequency sine wave corresponding to the channel; The difference between the target phase of the high-frequency sine wave corresponding to the channel and the bandwidth mismatch value of the channel is determined as the sampling time mismatch value of the channel.

12. The method according to claim 3, characterized in that, The frequency range interval corresponding to the channel is updated based on the parameters of the first median frequency and the second median frequency, including: If the effective value of the fitting residual of the first median frequency is less than the effective value of the fitting residual of the second median frequency, then the second median frequency of the frequency range before the update is determined as the left boundary frequency of the frequency range after the update, and the right boundary frequency of the frequency range before the update is determined as the right boundary frequency of the frequency range after the update. If the effective value of the fitting residual of the first median frequency is greater than the effective value of the fitting residual of the second median frequency, then the first median frequency of the frequency range interval before the update is determined as the right boundary frequency of the frequency range interval after the update, and the left boundary frequency of the frequency range interval before the update is determined as the left boundary frequency of the frequency range interval after the update.

13. A device for determining the error value of a channel in a data acquisition system, characterized in that, include: The frequency range interval determination module is used to determine the frequency range interval corresponding to each channel based on the initial frequency of the target sine wave and the frequency offset value corresponding to each channel in the data acquisition system, wherein the target sine wave includes a low-frequency sine wave and / or a high-frequency sine wave. The first median frequency parameter and second median frequency parameter determination module is used to determine the first median frequency parameter and second median frequency parameter in the frequency range interval corresponding to each channel based on a three-parameter sine fitting algorithm. The channel error value determination module is used to, for any channel, if the parameters of the first median frequency and the second median frequency corresponding to the channel do not meet the iteration stopping condition, update the frequency range interval corresponding to the channel according to the parameters of the first median frequency and the second median frequency corresponding to the channel, and return to execute the operation of determining the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel based on the updated frequency range interval corresponding to the channel, until the parameters of the first median frequency and the second median frequency meet the iteration stopping condition, and determine the error value of the channel according to the parameters of the first median frequency and the second median frequency in the frequency range interval corresponding to the channel.

14. An electronic device, characterized in that, The electronic device includes: At least one processor; and A memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, the computer program being executed by the at least one processor to enable the at least one processor to perform the method for determining the error value of a channel in the data acquisition system according to any one of claims 1-12.

15. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions that, when executed by a processor, implement the method for determining the error value of a channel in the data acquisition system according to any one of claims 1-12.

16. A computer program product, characterized in that, The computer program product includes a computer program that, when executed by a processor, implements the method for determining the error value of a channel in a data acquisition system according to any one of claims 1-12.

Images