High-Precision ADC Spectrum Testing under Non-Coherent Sampling Conditions

Realizing coherent sampling is one of the major bottlenecks in high-precision ADC spectrum testing. In spectrum analysis, if coherent sampling is not implemented, spectral leakage will result, which in turn leads to inaccurate test results. In this paper, a combined four-parameter sine-curve-fitting algorithm is proposed incorporating non-coherent sampling, with the amplitude, initial phase, and frequency parameters of the sine wave being obtained by fitting. The corresponding coherent sine wave is then calculated and replaced according to the obtained sine wave to reconstruct the new test data, eliminating the requirement of coherent sampling. Numerous simulations demonstrated the functionality and robustness of the algorithm, which was then used to process and analyze the measured data of two commercial high-precision ADCs. The results show that our algorithm can achieve accurate testing of ADC parameters under relaxed test conditions, which verifies the effectiveness and superiority of the scheme.


Introduction
With the rapid improvement of digital signal processing and the speed of computing, the accompanying increased requirement for system sensitivity has led, in turn, to an increased requirement for the indicators of high-speed and high-precision ADC (Analog to Digital Converter) testing [1,2]. In recent years, ADC spectrum testing has become an area of intense research interest in the mixed-signal IC test field. Of the dynamic parameter test methods available, spectrum analysis is the most commonly used [3]. This method calculates the dynamic parameters of the device to be measured by obtaining the power of the fundamental frequency and all harmonics in the spectrum. However, the challenge with this approach is that the IEEE Standard for Terminology and Test Methods for Analog-to-Digital Converters [4] and IEEE Standard for Terminology and Test Methods of Digital-to-Analog Converter Devices [5] impose several strict requirements, one of which is the implementation of coherent sampling.
The realization of this requirement for coherent sampling in spectrum analysis means that, for example, the input sine wave must be pure enough so that its distortion is much lower than that of the ADC under test. This is a particular challenge when the test environment (such as temperature), the aging of the equipment, and the inability to perform accurate period interception of the signal at the output (i.e., the number of measurement periods is not an integer), lead to skirt effects appearing on the output signal spectrum, as the quantization index corresponding to the fundamental wave will no longer be the unique index. Spectral leakage occurs, and the formula for calculating the spectrum parameters under standard conditions is no longer applicable, so accurate parameters cannot be obtained [6,7].
In this paper, a combined four-parameter sine-fitting algorithm is proposed for the theoretical derivation of the spectral leakage caused by non-coherent sampling. First, a three-parameter sine fit is performed using the input signal frequency as the known frequency to obtain the four parameters of amplitude, frequency, phase, and DC component. Then, using this as the initial value, the proposed four-parameter sine fit algorithm is used without iteration and absolute convergence to obtain the real test signal. Finally, the corresponding coherent sine wave is calculated and replaced according to the obtained test signal to reconstruct the new test data, eliminating the influence of spectral leakage of the fundamental on the spectrum test. The proposed algorithm is verified by simulation and experiment, and, finally, the accurate testing of high-precision ADC dynamic parameters is achieved under non-coherent sampling test conditions. The rest of the paper is divided into the following sections. Section 2 discusses the fundamentals of the ADC spectrum testing method. Section 3 presents the proposed algorithm. Section 4 describes the simulation and robustness verification of the algorithm. Section 5 verifies the applicability of the algorithm to determining ADC performance using a test board. Section 6 concludes the paper.

ADC Spectrum Testing
The basic principle of spectrum analysis is to input a high-purity sine signal to the ADC under test and then acquire the output quantized digital signal and perform a fast Fourier transform on this data to obtain the spectrum of the ADC to calculate its relevant dynamic parameters.
If we let x(t) be the time-domain expression of the simulated input signal, modeled with harmonic distortion, then the time-domain expression is: where A is the amplitude of the signal x(t), f sig is the input signal frequency, ϕ is the initial phase of the signal, H is the total number of harmonics, A h and ϕ h are the amplitude and initial phase of the harmonics, A h << A for all values of h taken 2 ≤ h ≤ H, whose initial phase ϕ h is distributed in the interval [0, 2π], and w is the noise in the signal. If we let x[n] be the analog interpretation of the digital output obtained from an ADC whose gain error and offset error have been calibrated, it can be expressed as: where n = 0, 1, 2, . . . , M − 1, M is the total number of sampling points, generally taken to the power of 2, J is the number of periods of the intercepted sampled signal, also known as the number of measurement periods, and w[n] is the noise at the nth sample point. Harmonics in the ADC output signal represent a nonlinear distortion of the ADC [21]. Discrete Fourier transform (DFT) with M sampling points is performed on the sampling sequence to obtain the exact spectral parameters. The expression for the Fourier transform of x[n] is: where k is the index of the frequency code value. Coherent sampling requires that the sampling time be an integer multiple of the input signal period [22], which has the advantage of avoiding spectral leakage. The international ADC test standard defines coherent sampling as: where f samp is the sampling frequency when the number of measurement periods J is an integer and the total number of sampling points M is coprime with J. To avoid repeated sampling, the same point is taken at different periods to achieve coherent sampling. Ideally, information on the fundamental frequency and all harmonics can be obtained in the spectrum of the ADC to calculate various dynamic parameters such as SNR, SFDR, or THD [24]. However, the sampling process of the input signal requires the input signal frequency to be accurate to 8-12 decimal places, and this condition is difficult to achieve, resulting in the coherent sampling condition not being met and spectral leakage, which will lead to the leakage of fundamental power to nearby frequencies, affecting or even drowning the harmonics, thus, preventing the accurate calculation of the fundamental frequency and all harmonics, and, ultimately, a failure to obtain the correct parameter test results. Figure 1 shows the spectrum obtained under coherent and non-coherent sampling conditions. where samp f is the sampling frequency when the number of measurement periods J is an integer and the total number of sampling points M is coprime with J . To avoid repeated sampling, the same point is taken at different periods to achieve coherent sampling. Ideally, information on the fundamental frequency and all harmonics can be obtained in the spectrum of the ADC to calculate various dynamic parameters such as SNR, SFDR, or THD [24]. However, the sampling process of the input signal requires the input signal frequency to be accurate to 8-12 decimal places, and this condition is difficult to achieve, resulting in the coherent sampling condition not being met and spectral leakage, which will lead to the leakage of fundamental power to nearby frequencies, affecting or even drowning the harmonics, thus, preventing the accurate calculation of the fundamental frequency and all harmonics, and, ultimately, a failure to obtain the correct parameter test results. Figure 1 shows the spectrum obtained under coherent and non-coherent sampling conditions.

Fundamental Waveform Estimation and Reconstruction Algorithm
Before describing our new approach to spectral analysis under non-coherent sampling conditions in detail, we will briefly describe the fundamental frequency estimation and reconstruction methods [22,23].
When DFT is applied to test data under non-coherent conditions, the leakage in the spectrum is mainly due to leakage of the fundamental frequency. For high-resolution Spectrum Power in(dB)

Fundamental Waveform Estimation and Reconstruction Algorithm
Before describing our new approach to spectral analysis under non-coherent sampling conditions in detail, we will briefly describe the fundamental frequency estimation and reconstruction methods [22,23].
When DFT is applied to test data under non-coherent conditions, the leakage in the spectrum is mainly due to leakage of the fundamental frequency. For high-resolution ADC tests, the leakage at all frequencies other than the fundamental frequency is significantly lower than the total noise power of the ADC. This effect is illustrated in Figures 2 and 3. Figure 2 shows the DFT spectrum of the non-coherent test data. There is significant spectral leakage around the fundamental. However, if the non-coherent fundamental frequency in this data is estimated and removed, accurate harmonic and noise information can be obtained from the remaining spectrum ( Figure 3). This demonstrates that an accurate estimation of the non-coherent fundamental frequency is required to obtain the correct spectral results. ADC tests, the leakage at all frequencies other than the fundamental frequency is significantly lower than the total noise power of the ADC. This effect is illustrated in    Figure 2 shows the DFT spectrum of the non-coherent test data. There is significant spectral leakage around the fundamental. However, if the non-coherent fundamental frequency in this data is estimated and removed, accurate harmonic and noise information can be obtained from the remaining spectrum ( Figure 3). This demonstrates that an accurate estimation of the non-coherent fundamental frequency is required to obtain the correct spectral results. ADC tests, the leakage at all frequencies other than the fundamental frequency is significantly lower than the total noise power of the ADC. This effect is illustrated in Figures 2 and 3.   Figure 2 shows the DFT spectrum of the non-coherent test data. There is significant spectral leakage around the fundamental. However, if the non-coherent fundamental frequency in this data is estimated and removed, accurate harmonic and noise information can be obtained from the remaining spectrum ( Figure 3). This demonstrates that an accurate estimation of the non-coherent fundamental frequency is required to obtain the correct spectral results.

Fitting of Sine Test Signals
There are several methods for estimating non-coherent fundamental waves, such as the interpolated discrete Fourier transform method [16], which first adds a window to the non-coherent sampled data and then interpolates to accurately estimate the fundamental and the previously-discussed FIRE algorithm [24], with its time-intensive Newton iteration. However, we propose a combined four-parameter sine-fitting algorithm, which first uses the frequency of the input signal as the known frequency as the basis of three-parameter sine fitting to obtain the amplitude, initial phase, DC component, and fitting residuals of the fundamental frequency and then combines frequency estimation and three-parameter sine fitting to implement a four-parameter least-squares fitting algorithm to finally obtain an accurate test signal.

Three-Parameter Sine-Fitting Method
The input signal is modeled in the same way as Equation (1), and the input signal expression is obtained by splitting Equation (1) using the sum and difference product formula in trigonometric functions as follows: , and c 0 are the DC components. The equal interval uniform data acquisition sequence is x 1 , x 2 , · · · , x n . If the sampling interval is taken as ∆t, the acquisition rate v, the acquisition time t i = i × ∆t = i/v of the sampling point x i , i = 1, 2, · · · , n, and the digital angular frequency ω = 2π f sig /v, then the signal can be expressed in the following discrete form: The three-parameter sine-curve-fitting process, using the frequency of the input signal as the initial value, minimizes the sum of squares of the residuals ε described in Equation (6) by searching for a, b, and c: Then a, b, and c are the least squares fits of a 0 , b 0 , and c 0 . The fitted function is obtained when the value ε reaches its minimum: That is: The fitted residuals RMS value can be written: As this is a closed algorithm, there is no convergence problem.

Four-Parameter Combination Estimation Algorithm
The input signal frequency estimation is combined with the above three-parameter sine-curve-fitting method to obtain the estimation results of the four parameters of the sine waveform. However, when the sampling interval is sufficiently small, the difference between the amplitude values of adjacent sampling points will be negligible. When the difference is less than a minimum amplitude quantization step, adjacent analog value points may have the same quantization step code value. Therefore, the difference between the cumulative formula of the difference between adjacent sampling points and the theoretical value becomes too large, leading to the error of the obtained digital angular frequency ω becoming too large, which affects the accuracy of other parameters obtained by the three-parameter fitting algorithm based on ω. Here, we propose to replace the subtraction of adjacent sampling points with the subtraction of two adjacent sampling points of k, to avoid the influence of too many adjacent sampling points with the same quantization code value and, thus, form the following improved algorithm.
The sampling sequence of the sine signal waveform x 1 , x 2 , · · · , x n , whose function relationship is shown in Equation (6) can be written: If we let the error of the measurement point x i be γ i , and q = 2 cos ω, then: which we can rearrange into: where k is a natural number ranging from one to N, and N is the number of samples in each sine wave period. The accuracy of the estimated parameters changes periodically as the value of k increases, and the change period is the number of sampling points per cycle N. The best estimation results are obtained when k ∼ = N/2. In this case, the accuracy of the estimated values of each parameter is the highest. Here, q k = 2 cos ω k and η k,i is the random error, which is obtained by selecting q k such that: Let dρ k dq k = 0, then we obtain: ω k is the least square estimation of the digital angular frequency of sine wave, and the real input signal frequency f can be obtained from ω k = 2π f /v. After using the threeparameter estimation method, the estimated values of the other three parameters can be obtained as A for the amplitude, ϕ for the initial phase, and c for the DC component, thus, completing the four-parameter estimation. The resulting test signal obtained by fitting can then be:

Reconstructing Coherent Sine Signals
To obtain coherent test data, the test signal obtained by fitting is replaced in the test data by the reconstructed coherent test signal, and then the reconstructed test data is subjected to DFT to obtain accurate spectral test results [20]. The test signal expression, which is derived from Equation (4), can be written: where δ = 0 for coherent sampling, and δ = 0 for non-coherent sampling. As Equation (26) is the test signal obtained by fitting, the frequencies in Equations (26) and (27) are equal, and it is obtained that: From Equation (28), we obtain J + δ as: where f is the true frequency obtained from the test data. To ensure that all samples are valid and there are no redundant samples, the total number of sample points M is mutually prime with the number of periods J, and the length of the test data at the time of testing M is an integer power of two. When this is the case, any odd value of J satisfies the mutual prime. Reassigning the nearest odd value of J + δ to J new , the signal frequency of the coherent test condition is obtained from J new : Replacing the frequency f in Equation (26) with f r yield: where x r is the coherent test signal, and the non-coherent sampling fundamental signal is removed from the original test data and replaced with a coherent sampling fundamental signal by Equations (26) and (31) to obtain the new test data x new [n]: DFT of the test data provides the frequency domain information x newk , which in turn yields accurate ADC dynamic parameters [31]. The algorithm can perform spectrum testing without windowing and requires no iterations, absolute convergence, or high computational speed, which saves time and reduces testing costs. The flowchart in Figure 4 shows the process flow of spectrum testing using a combined four-parameter sine-curvefitting algorithm.
yields accurate ADC dynamic parameters [31]. The algorithm can perfo ing without windowing and requires no iterations, absolute convergen tational speed, which saves time and reduces testing costs. The flow shows the process flow of spectrum testing using a combined four-par fitting algorithm.

Numerical Simulation
In this section, effective simulation experiments are performed to v ality and robustness of the proposed algorithm. The computation tim is also compared with other non-coherent sampling algorithms in the l

Functionality
Firstly, the output signal of a 16-bit ADC is simulated under the coherent sampling and harmonic distortion of the input signal. The freq signal is 50 kHz, the sampling rate is 480 kS/s, and the total number of 32768. The spectrum of the original output signal, the coherently samp signal after processing by the proposed algorithm are shown in Figur

Numerical Simulation
In this section, effective simulation experiments are performed to verify the functionality and robustness of the proposed algorithm. The computation time of this algorithm is also compared with other non-coherent sampling algorithms in the literature.

Functionality
Firstly, the output signal of a 16-bit ADC is simulated under the conditions of noncoherent sampling and harmonic distortion of the input signal. The frequency of the input signal is 50 kHz, the sampling rate is 480 kS/s, and the total number of sampling points is 32,768. The spectrum of the original output signal, the coherently sampled signal, and the signal after processing by the proposed algorithm are shown in Figure 5. The test signal is a single frequency signal source. According to the actual test situation, the amplitude of harmonics and noise has been large enough. In the real test situation, the harmonics and noise generated by the power supply and the peripheral circuit of the system are far smaller than those set in the text. The simulation conditions have been set enough to cover the actual application test.
For the above signal, the number of measurement periods calculated by Equation (4) is 3413.33, i.e., a non-integer number; therefore, the original output signal is obtained under non-coherent sampling conditions; thus, spectral leakage will occur in the spectrum, which cannot be directly used to calculate the dynamic parameters by the standard calculation formula, as shown in the blue curve in Figure 5. Processing the non-coherent sampling data represented by the blue curve with the proposed algorithm results in the green curve. A comparison of the green and blue curves shows that the use of the algorithm eliminates the spectral leakage caused by non-coherent sampling completely, and the clean spectrum has the same effect as the output signal spectrum of the coherent test signal (red curve) under ideal test conditions, which verifies the effectiveness of the algorithm.

Robustness
To verify the reliability of the proposed algorithm under different degrees of noncoherence, we performed many random simulations. As a reference, the parameters obtained from a set of coherently sampled data were used as the real values. The fractional part of the sampling period delta was randomly generated in the range of -0.5 to 0.5, and all other simulation settings were the same as in the functional verification. Further, 1000 simulations were performed. The non-coherent data from each simulation were processed using our algorithm. The error between the obtained results and the real values for the 1000 simulations is shown in Figure 6. After processing with our algorithm, the maximum error of each parameter is less than 2 dB in absolute value, and the algorithm can still obtain accurate results even with severe non-coherent sampling (when the delta is very close to 0.5). This test confirms that, for any degree of non-coherent sampling, the errors of each parameter are within the acceptable range, and the algorithm is reliable. For the above signal, the number of measurement periods calculated by Equation (4) is 3413.33, i.e., a non-integer number; therefore, the original output signal is obtained under non-coherent sampling conditions; thus, spectral leakage will occur in the spectrum, which cannot be directly used to calculate the dynamic parameters by the standard calculation formula, as shown in the blue curve in Figure 5.
Processing the non-coherent sampling data represented by the blue curve with the proposed algorithm results in the green curve. A comparison of the green and blue curves shows that the use of the algorithm eliminates the spectral leakage caused by non-coherent sampling completely, and the clean spectrum has the same effect as the output signal spectrum of the coherent test signal (red curve) under ideal test conditions, which verifies the effectiveness of the algorithm.

Robustness
To verify the reliability of the proposed algorithm under different degrees of noncoherence, we performed many random simulations. As a reference, the parameters obtained from a set of coherently sampled data were used as the real values. The fractional part of the sampling period delta was randomly generated in the range of -0.5 to 0.5, and all other simulation settings were the same as in the functional verification. Further, 1000 simulations were performed. The non-coherent data from each simulation were processed using our algorithm. The error between the obtained results and the real values for the 1000 simulations is shown in Figure 6. After processing with our algorithm, the maximum error of each parameter is less than 2 dB in absolute value, and the algorithm can still obtain accurate results even with severe non-coherent sampling (when the delta is very close to 0.5). This test confirms that, for any degree of non-coherent sampling, the errors of each parameter are within the acceptable range, and the algorithm is reliable.

Computation Time
In Table 1, we list the computation time of the direct DFT using different window functions and the FIRE algorithm. Among all the methods listed, the algorithm proposed in this paper provides accurate test results with the least cost in terms of computation time. It can also be seen from the table that only one window function can accurately test the 16-bit ADC. This demonstrates the degree to which the use of the window function for spectrum testing depends on prior knowledge of the ADC resolution, while the special windows that meet the requirements are often very complex expressions, which are not conducive to computation. The FIRE algorithm provides accurate test results but takes a long time to perform the Newton iteration, which is not conducive to mass production testing of ADC chips.

Experimental Verification
To further demonstrate the utility of the proposed algorithm, we conducted a spectral test of non-coherent sampling by using the proposed algorithm for two commercial high-precision ADCs. A Kintex-7-FPGA based A/D converter test system was designed and built for testing the dynamic parameters of the ADCs, with the overall scheme shown in Figure 7.

Computation Time
In Table 1, we list the computation time of the direct DFT using different window functions and the FIRE algorithm. Among all the methods listed, the algorithm proposed in this paper provides accurate test results with the least cost in terms of computation time. It can also be seen from the table that only one window function can accurately test the 16-bit ADC. This demonstrates the degree to which the use of the window function for spectrum testing depends on prior knowledge of the ADC resolution, while the special windows that meet the requirements are often very complex expressions, which are not conducive to computation. The FIRE algorithm provides accurate test results but takes a long time to perform the Newton iteration, which is not conducive to mass production testing of ADC chips.

Experimental Verification
To further demonstrate the utility of the proposed algorithm, we conducted a spectral test of non-coherent sampling by using the proposed algorithm for two commercial highprecision ADCs. A Kintex-7-FPGA based A/D converter test system was designed and built for testing the dynamic parameters of the ADCs, with the overall scheme shown in Figure 7.  The test evaluation board used to test the performance of the ADC is shown in Figure  8 below. It uses a high-speed Kintex-7 FPGA as the core device, which carries out the sampling of the output data of the high-performance A/D converter. It also allows the caching of data and the processing of data across clock domains. DDR3 SDRAM is used as the memory for the sampling data during A/D conversion, with a data throughput rate of up to 20 Gbits/s, which meets the requirements of high bandwidth, high storage depth, and high-speed real-time data reading and writing for high-performance testing. The SDRAM uses EMMC to realize high speed and large capacity storage, with a total effective storage capacity of 2TB, four pieces of storage speed up to 250MB/s*4 = 1GB/s. The board is completed by ethernet and fiber optic interface protocol implementation and data transmission. In addition, the JESD204B standard transceiver interface has been integrated into the FPGA, which provides great convenience for data sampling and the sending of highperformance A/D converters. The test evaluation board used to test the performance of the ADC is shown in Figure 8 below. It uses a high-speed Kintex-7 FPGA as the core device, which carries out the sampling of the output data of the high-performance A/D converter. It also allows the caching of data and the processing of data across clock domains. DDR3 SDRAM is used as the memory for the sampling data during A/D conversion, with a data throughput rate of up to 20 Gbits/s, which meets the requirements of high bandwidth, high storage depth, and high-speed real-time data reading and writing for high-performance testing. The SDRAM uses EMMC to realize high speed and large capacity storage, with a total effective storage capacity of 2TB, four pieces of storage speed up to 250MB/s*4 = 1GB/s. The board is completed by ethernet and fiber optic interface protocol implementation and data transmission. In addition, the JESD204B standard transceiver interface has been integrated into the FPGA, which provides great convenience for data sampling and the sending of high-performance A/D converters.

AD976 Test
The proposed algorithm was first tested using a 16-bit high-precision AD976. The procedure for the test was as follows:

AD976 Test
The proposed algorithm was first tested using a 16-bit high-precision AD976. The procedure for the test was as follows: 1.
Connect the AD976 board to be tested to the DC power supply and the evaluation board in the specified environment; input a high-precision sine wave f in = 44, 055.15234375 Hz with the sampling frequency set to f samp = 200 kHz.

2.
Power up the device and mode control, and connect the device's digital output to the digital acquisition terminal through the high-speed interface.

3.
Provide a frequency-specific analog input signal through a high-performance RF signal source, and connect a fixed-frequency filter to the AD976 analog input. 4.
Use the logic analyzer/evaluation board to set the AD976 for dynamic conversion and acquisition of the digital output signals of the device. 5.
DFT to the obtained data to obtain the frequency domain information.
As expected, there was significant spectral leakage in the non-coherently sampled data. However, applying our algorithm, clean and accurate spectra were obtained, as is shown by the red spectrum in Figure 9, which matches the blue spectrum (coherent) and provides accurate spectral results for the high-resolution ADC. To show the effect of the use of a window function, the black plot in Figure 9 shows the spectrum obtained when using the Window in [13] function for the same non-coherent sampled data. Table 2 compares the SNR, SINAD, THD, and SFDR values of the ADC using our algorithm and the windowing technique with the values obtained using a coherent sampling method. The algorithm accurately estimates the parameters and suppresses the spectral leakage somewhat better than the Window in [13] function, and the higher the number of ADC bits, the better the suppression of leakage.
x FOR PEER REVIEW 14 of 17

AD9230 Test
To further verify the proposed algorithm, another commercial high-precision ADC was also tested ( Figure 10). The procedure, in this case, was as follows. The obtained spectrum is shown in Figure 11. Our algorithm the spectral leakage. Table 3 shows the test results, which demon proposed algorithm, achieving accurate estimation of the dy SINAD, and ENOB.

1.
Apply a 1.8 V power supply.

2.
Apply a voltage of −1 dBFS amplitude to the ADC under test, then a sine wave source at f in = 10,002,700 Hz, and filter the sine wave input to remove distortion and random noise from the input signal.

3.
Apply a sample clock with the specified sampling frequency f samp = 75,366,400 Hz to the ADC under test.

4.
After the ADC has stabilized, collect 32,768 output conversion data points. 5.
Use the proposed algorithm and DFT to obtain the spectrum.
The obtained spectrum is shown in Figure 11. Our algorithm successfully eliminated the spectral leakage. Table 3 shows the test results, which demonstrate the accuracy of the proposed algorithm, achieving accurate estimation of the dynamic parameters SNR, SINAD, and ENOB.
The obtained spectrum is shown in Figure 11. Our algorithm successfully eliminated the spectral leakage. Table 3 shows the test results, which demonstrate the accuracy of the proposed algorithm, achieving accurate estimation of the dynamic parameters SNR, SINAD, and ENOB. Figure 11. Spectrogram of AD9230 test data. In summary, our proposed algorithm can accurately identify and reconstruct the fundamental frequency of the output signal under non-coherent sampling and input signal with harmonic distortion test conditions and achieve accurate testing of the ADC dynamic parameters. Non-coherent + Algorithm Coherent + DFT Figure 11. Spectrogram of AD9230 test data. In summary, our proposed algorithm can accurately identify and reconstruct the fundamental frequency of the output signal under non-coherent sampling and input signal with harmonic distortion test conditions and achieve accurate testing of the ADC dynamic parameters.

Conclusions
Under the test conditions of non-coherent sampling and an input signal with harmonic distortion, the ADC dynamic parameters cannot be calculated directly due to spectral leakage. In this paper, a combined four-parameter sine-curve-fitting algorithm is proposed to obtain the amplitude, initial phase, and frequency parameters of the sine wave by fitting. This is used to construct a coherent sampling signal of the fundamental frequency that replaces the non-coherent sampling signal of the fundamental frequency in the test data, and the new signal obtained greatly improved the spectral leakage. Unlike windowing, this algorithm can be used to test any ADC output without prior knowledge of the ADC resolution. The accuracy and robustness of the proposed algorithm were verified by numerical simulation, and a simulation-based computation time comparison verifies the operational efficiency of the algorithm. Finally, two commercial high-precision ADCs, the AD976 and the AD9230, were evaluated, and the results were SNR: 74.321 dB, SINAD: 73.654 dB, SFDR: 79 dB, and THD: 89.024 dB for the AD976, and SNR: 62.329 dB, SINAD: 63.097 dB, and ENOB: 9.198 bits for the AD9230. These results show the practical effectiveness of the proposed algorithm. The algorithm can achieve accurate testing of ADC dynamic parameters under arbitrary non-coherent sampling conditions without complex calculations or tedious manual calibration, greatly reducing test time and test cost.