Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Research on Multi-Fault Diagnosis Method Based on Time Domain Features of Vibration Signals
Next Article in Special Issue
Operational State Recognition of a DC Motor Using Edge Artificial Intelligence
Previous Article in Journal
Graph Convolutional Network for 3D Object Pose Estimation in a Point Cloud
Previous Article in Special Issue
Static Load Characteristics of Hydrostatic Journal Bearings: Measurements and Predictions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
National Key Laboratory for Electronic Measurement Technology, North University of China, Taiyuan 030051, China
2
School of Electronic Information Engineering, Taiyuan University of Science and Technology, Taiyuan 030024, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(21), 8170; https://doi.org/10.3390/s22218170
Submission received: 22 September 2022 / Revised: 15 October 2022 / Accepted: 21 October 2022 / Published: 25 October 2022
(This article belongs to the Collection Instrument and Measurement)

Abstract

:
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.

1. 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].
To suppress spectral leakage due to non-coherent sampling, researchers have proposed several approaches, the most commonly used of which are the windowing technique [8], the interpolated discrete Fourier transform method [9], the fundamental identification and replacement method [10], and the four-parameter sine-fitting method [11].
The windowing technique involves multiplying the digital output of the ADC with a suitable window function to mitigate the leakage of fundamental power. Common window functions include the Hanning, Blackman, Nuttall, and maximum attenuation (MSD) windows [12,13,14,15]. The spectrum of the window function includes a main lobe and multiple side lobes. When using the windowing technique, it is critical to select an appropriate window function for the specific situation, as the use of unsuitable window functions often causes difficulties in testing. A good side lobe performance of the window function means more concentrated power, low harmonics influence, and accurate power calculation, but it also means a wide main lobe and low-frequency resolution. This means that testers need to have a great deal of a priori knowledge of window functions to select the right type to use. At the same time, the industry currently imposes extremely strict restrictions on the selection criteria of the window function, requiring the first side lobe to be lower than the noise floor, which means that most common window functions cannot be applied to high-precision ADCs, while the special windows that meet the requirements are often very complex expressions, which are not conducive to the calculation.
The interpolated discrete Fourier transform method, which uses the fast Fourier transform for parameter estimation by discrete spectral analysis, can effectively solve the spectral leakage problem in some situations [16,17]. Interpolation of spectral lines by adding complementary zeros (zero-filling) after the effective signal in the time domain can make the spectral lines denser and avoid the limitation of frequency resolution, successfully reducing the fence effect to a certain extent. However, the amount of data increases after zero-filling, which leads to longer processing times and poor real-time performance [18].
Researchers have also proposed many solutions to address the difficulties of implementing coherent sampling in high-resolution ADCs during measurements [19,20,21]. Since non-coherent sampling is difficult to avoid, the non-coherent sampling problem is usually studied and analyzed when considering how to relax other strict standard test conditions [22]. For example, when studying input test signals with harmonic distortion, the non-coherent sampling problem is often considered at the same time. The FIRE algorithm proposed by [23,24,25] accomplishes the testing of ADC dynamic parameters under non-coherent sampling conditions by identifying and replacing the fundamental wave of the output signal under non-coherent sampling conditions. However, the Newton iteration part of the FIRE algorithm consumes a large amount of processing time, which is unfavorable for the mass production of ADC chips and will increase the cost of testing.
Another widely used method is the four-parameter sine-fitting method with accurate estimates of spectral parameters [26,27], such as the total harmonic distortion (THD), signal-to-noise ratio (SNR), signal-to-noise and distortion (SINAD), spurious free dynamic range (SFDR), and the effective number of bits (ENOB). However, this approach is limited by its computational inefficiency, the need for pre-estimation of the four fitted parameters, and the need for sufficiently high estimation accuracy [28,29,30]. In addition, the convergence of the four-parameter iterative search method has been problematic. There is no clear conclusion or criterion for judging when it converges, and there is no clear interval of convergence. The closer the initial value is to the optimal value, the more the convergence of the iterative process is guaranteed; it will otherwise be more likely to lead to a divergent iterative process.
All the above methods suffer from one or more problems, such as a large amount of data, long computation time, inaccurate frequency estimation, a lack of robustness in the Nyquist range, a dependence of the results on the selected window function type, or an inability to perform a full-spectrum test. Therefore, there is a need to develop a method that can solve all the above problems and perform spectral tests accurately without the need for coherent sampling.
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.

2. 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:
x ( t ) = A cos ( 2 π f s i g t + φ ) + h = 2 H A h cos ( 2 π h f s i g t + φ h ) + w ( t )
where A is the amplitude of the signal x ( t ) , f s i g 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:
x [ n ] = A cos ( 2 π J M n + φ ) + h = 2 H A h cos ( 2 π h J M n + φ h ) + w [ n ]
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:
X k = 1 M n = 0 M 1 x [ n ] e j 2 π k M n , k = 0 , 1 , , M 1
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:
f s i g f s a m p = J M
where f s a m p 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.

3. 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 Figure 2 and Figure 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.

3.1. 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.

3.1.1. 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:
x ( t ) = A 0 cos ( 2 π f s i g t + φ ) + c 0 = a 0 cos ( 2 π f s i g t ) + b 0 sin ( 2 π f s i g t ) + c 0
where a 0 = A 0 cos ( φ ) , b 0 = A 0 sin ( φ ) , 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 s i g / v , then the signal can be expressed in the following discrete form:
x ( i ) = A 0 cos ( ω i + φ 0 ) + c 0 = a 0 cos ( ω i ) + b 0 sin ( ω i ) + c 0
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 :
ε = i = 1 n [ x i a cos ( ω i ) b sin ( ω i ) c ] 2
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:
x ^ ( i ) = acos ( ω i ) + bsin ( ω i ) + c
That is:
x ^ ( i ) = A cos ( ω i + θ ) + c
A = a 2 + b 2
θ = { arctan ( b a ) , a 0 arctan ( b a ) + π , a < 0
The fitted residuals RMS value can be written:
ρ = ε n ,   ε = i = 1 n ( x i x ^ ( i ) ) 2
As this is a closed algorithm, there is no convergence problem.

3.1.2. 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:
z ( i ) = A 0 cos ( ω i + φ 0 )     i = 1 , 2 , , n
If we let the error of the measurement point x i be γ i , and q = 2 cos ω , then:
x i = x ( i ) + γ i = z ( i ) + c + γ i
z ( i ) + z ( i 2 ) = ( 2 cos ω ) z ( i 1 ) = q z ( i 1 )
x ( i ) c + x ( i 2 ) c = q ( x ( i 1 ) c )
x i c γ i + x i 2 c γ i 2 = q ( x i 1 c γ i 1 )
x i + 1 c γ i + 1 + x i 1 c γ i 1 = q ( x i c γ i )
which we can rearrange into:
λ k , i = x i + k x i
η k , i = γ i + k γ i
x i + k + 1 c γ i + k + 1 + x i + k 1 c γ i + k 1 = q k ( x i + k c γ i + k )
λ k , i + 1 + λ k , i 1 q k λ k , i = η k , i + 1 + η k , i 1 q k η k , i
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:
ρ k = i = 2 n k 1 ( λ k , i + 1 + λ k , i 1 q k λ k , i ) 2 = m i n
Let d ρ k d q k = 0 , then we obtain:
q k = λ k , n k 1 λ k , n k + λ k , 1 λ k , 2 + 2 i = 2 n k 2 λ k , i λ k , i + 1 i = 2 n k 2 λ k , i 2
ω k = arccos ( q k 2 )
ω 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 three-parameter 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:
x t = A cos ( 2 π f t + φ ) + c

3.2. 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:
x ( t ) = A cos ( 2 π J + δ M f s a m p t + φ ) + c
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:
f = J + δ M f s a m p
From Equation (28), we obtain J + δ as:
J + δ = M f f s a m p
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 n e w , the signal frequency of the coherent test condition is obtained from J n e w :
f r = J n e w M f s a m p
Replacing the frequency f in Equation (26) with f r yield:
x r = A cos ( 2 π f r t + φ ) + c
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 n e w [ n ] :
x n e w [ n ] = x [ n ] x t [ n ] + x r [ n ]
DFT of the test data provides the frequency domain information x new k , 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-curve-fitting algorithm.

4. 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.

4.1. Functionality

Firstly, the output signal of a 16-bit ADC is simulated under the conditions of non-coherent 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.

4.2. Robustness

To verify the reliability of the proposed algorithm under different degrees of non-coherence, 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.

4.3. 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.

5. 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.
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.

5.1. AD976 Test

The proposed algorithm was first tested using a 16-bit high-precision AD976. The procedure for the test was as follows:
  • 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 i n = 44,055.15234375 Hz with the sampling frequency set to f s a m p = 200 kHz.
  • Power up the device and mode control, and connect the device’s digital output to the digital acquisition terminal through the high-speed interface.
  • 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.
  • Use the logic analyzer/evaluation board to set the AD976 for dynamic conversion and acquisition of the digital output signals of the device.
  • 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.

5.2. 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.
  • Apply a 1.8 V power supply.
  • Apply a voltage of −1 dBFS amplitude to the ADC under test, then a sine wave source at f i n = 10,002,700 Hz, and filter the sine wave input to remove distortion and random noise from the input signal.
  • Apply a sample clock with the specified sampling frequency f s a m p = 75,366,400 Hz to the ADC under test.
  • After the ADC has stabilized, collect 32,768 output conversion data points.
  • 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.
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.

6. 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.

Author Contributions

Conceptualization, X.P. and J.L.; methodology, D.Z.; software, J.J.; validation, C.H.; formal analysis, N.S.; investigation, J.L.; resources, D.Z.; data curation, J.J.; writing—original draft preparation, C.H.; writing—review and editing, X.P.; visualization, D.Z.; supervision, J.L.; project administration, C.H.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No: 61973280), the National Natural Science Foundation of China (No: 62003316), China Postdoctoral Science Foundation (No: 2019M661069) and Shanxi Province Key R&D Program (No: 202003D111007).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Varier, V.; Sun, N. High-Precision ADC Testing with Relaxed Reference Voltage Stationarity. IEEE Trans. Instrum. Meas. 2021, 70, 2001809. [Google Scholar] [CrossRef]
  2. Islam, R.; Maktoomi, M.H.; Ren, H.; Arigong, B. Spectrum Aggregation Dual-Band Real-Time RF/Microwave Analog Signal Processing from Microstrip Line High-Frequency Hilbert Transformer. IEEE Trans. Microw. Theory Tech. 2021, 69, 4647–4657. [Google Scholar] [CrossRef]
  3. Poczekajlo, P.; Suszynski, R. Modern computing methods for digital signal processing engineering systems. In Proceedings of the 25th KES International Conference on Knowledge-Based and Intelligent Information and Engineering Systems, KES 2021, Szczecin, Poland, 8–10 September 2021; pp. 3534–3541. [Google Scholar]
  4. IEEE Std 1241TM-2010; IEEE Standard for Terminology and Test Methods for Analog-to-Digital Converters. IEEE: New York, NY, USA, 2011; pp. 1–139.
  5. 1658-2011; IEEE Standard for Terminology and Test Methods of Digital-to-Analog Converter Devices. IEEE: New York, NY, USA, 2012.
  6. Kuwalek, P.; Otomanski, P. The effect of the phenomenon of ‘spectrum leakage’ on the measurement of power quality parameters. In Proceedings of the 12th International Conference on Measurement, Smolenice, Slovakia, 27–29 May 2019; pp. 70–73. [Google Scholar]
  7. Xiao, Q.Y.; Li, J.; Sun, J.D.; Feng, H.; Jin, S.J. Natural-gas pipeline leak location using variational mode decomposition analysis and cross-time-frequency spectrum. Measurement 2018, 124, 163–172. [Google Scholar] [CrossRef]
  8. Jwo, D.-J.; Chang, W.-Y.; Wu, I.H. Windowing Techniques, the Welch Method for Improvement of Power Spectrum Estimation, Computers. Mater. Contin. 2021, 67, 3983–4003. [Google Scholar]
  9. Yang, C.; Li, B.; Hu, X.; Fu, Z. Accurate Frequency Estimation Based on Analytical Interpolated Discrete Time Fourier Transform. Diangong Jishu Xuebao/Trans. China Electrotech. Soc. 2022, 37, 1338–1350. [Google Scholar]
  10. Sudani, S.; Wu, M.; Chen, D. A novel robust and accurate spectral testing method for non-coherent sampling. In Proceedings of the International Test Conference 2011, Anaheim, CA, USA, 18–23 September 2011. [Google Scholar]
  11. Bilau, T.Z.; Megyeri, T.; Sarhegyi, A.; Markus, J.; Kollar, I. Four-Parameter Fitting of Sine Wave Testing Result: Iteration and Convergence; Elsevier: Amsterdam, The Netherlands, 2004; pp. 51–56. [Google Scholar]
  12. Agarwal, P.; Kumar, S.; Singh, S.P. Closed form solutions of various window functions in fractional fourier transform domain. In Proceedings of the 6th International Conference on Computing for Sustainable Global Development, INDIACom 2019, New Delhi, India, 13–15 March 2019; pp. 64–68. [Google Scholar]
  13. Albrecht, H.-H. Tailoring of minimum sidelobe cosine-sum windows for high-resolution measurements. Open Signal Process. J. 2010, 3, 20–29. [Google Scholar] [CrossRef]
  14. Wang, X.; Lei, K.; Yang, X.; Li, M.; Wang, X. Harmonic analysis based on Blackman Harris self-multiplication window. In Proceedings of the 5th Asia Conference on Power and Electrical Engineering, ACPEE 2020, Chengdu, China, 4–7 June 2020; pp. 2165–2169. [Google Scholar]
  15. Sulistyaningsih, S.; Putranto, P.; Qurrachman, T.; Desvasari, W.; Daud, P.; Wijayanto, Y.N.; Mahmudin, D.; Kurniadi, D.P.; Rahman, A.N.; Hardiati, S.; et al. Performance comparison of Blackman, bartlett, Hanning, and kaiser window for radar digital signal processing. In Proceedings of the 4th International Conference on Information Technology, Information Systems and Electrical Engineering, ICITISEE 2019, Yogyakarta, Indonesia, 20–21 November 2019; pp. 391–394. [Google Scholar]
  16. Belega, D.; Petri, D.; Dallet, D. Iterative sine-wave frequency estimation by generalized fourier interpolation algorithms. In Proceedings of the 11th International Symposium on Electronics and Telecommunications, ISETC 2014, Timisoara, Romania, 14–15 November 2014. [Google Scholar]
  17. Borkowski, J.; Matusiak, A. Dynamic characteristics of Interpolated-DFT-based frequency estimation method with Generalized Maximum Sidelobes Decay Windows. In Proceedings of the 19th IEEE International Conference on Environment and Electrical Engineering and 2019 IEEE Industrial and Commercial Power Systems Europe, EEEIC/I and CPS Europe 2019, Genoa, Italy, 11–14 June 2019. [Google Scholar]
  18. Yang, C.; Li, B.; Hu, X.; Fu, Z. High-precision Frequency Estimation Based on Iterative DTFT Interpolation. Power Syst. Technol. 2021, 45, 3955–3963. [Google Scholar]
  19. Sudani, S.; Chen, D.; Geiger, R. A method for accurate full spectrum testing without requiring coherency. In Proceedings of the 2012 IEEE 55th International Midwest Symposium on Circuits and Systems, MWSCAS 2012, Boise, ID, USA, 5–8 August 2012; pp. 346–349. [Google Scholar]
  20. Sudani, S.K.; Xu, L.; Chen, D. A comparative study of state-of-The-Art high-performance spectral test methods. IEEE Des. Test 2015, 32, 26–35. [Google Scholar] [CrossRef]
  21. Wang, X.; Lei, K.; Yang, X.; Wang, X. An effective harmonic parameters estimation algorithm based on transformed FFT sequence. In Proceedings of the 2nd International Conference on Information Technologies and Electrical Engineering, ICITEE 2019, Zhuzhou, China, 6–7 December 2019. [Google Scholar]
  22. Zhuang, Y.; Chen, D. Accurate and Robust Spectral Testing with Relaxed Instrumentation Requirements; Springer International Publishing: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  23. Zhuang, Y.; Chen, D. Algorithms for Accurate Spectral Analysis in the Presence of Arbitrary Noncoherency and Large Distortion. IEEE Trans. Instrum. Meas. 2017, 66, 2556–2565. [Google Scholar] [CrossRef]
  24. Sudani, S.K.; Chen, D. FIRE: A fundamental identification and replacement method for accurate spectral test without requiring coherency. IEEE Trans. Instrum. Meas. 2013, 62, 3015–3025. [Google Scholar] [CrossRef]
  25. Sudani, S.; Xu, L.; Chen, D. Accurate full spectrum test robust to simultaneous non-coherent sampling and amplitude clipping. In Proceedings of the 44th IEEE International Test Conference, ITC 2013, Anaheim, CA, USA, 10–12 September 2013. [Google Scholar]
  26. Simon, G.; Pintelon, R.; Sujbert, L.; Schoukens, J. An efficient nonlinear least square multisine fitting algorithm. IEEE Trans. Instrum. Meas. 2002, 51, 750–755. [Google Scholar] [CrossRef]
  27. Kollar, I.; Blair, J.J. Improved determination of the best fitting sine wave in ADC testing. IEEE Trans. Instrum. Meas. 2005, 54, 1978–1983. [Google Scholar] [CrossRef]
  28. Da Silva, M.F.; Ramos, P.M.; Serra, A. A new four parameter sine fitting technique. Measurement 2004, 35, 131–137. [Google Scholar] [CrossRef]
  29. Ramos, P.M.; Radil, T.; Serra, A.C. Four-parameter sine-fitting algorithm for detection and classification of transients and waveform distortions. Metrol. Meas. Syst. 2008, 15, 441–456. [Google Scholar]
  30. Liang, Z.; Ren, D.; Sun, J.; Zhu, Z. Fitting algorithm of sine wave with partial period waveforms and non-uniform sampling based on least-square method. In Proceedings of the 13th International Conference on Vibration Measurements by Laser and Noncontact Techniques, Ancona, Italy, 20–22 June 2018. [Google Scholar]
  31. Wu, M.; Ban, C.; Xu, J.; Geng, L.; Chen, D. MIRE: A Multitone Identification and Replacement Method for Multitone Spectral Test Without Requiring Coherent Sampling. IEEE Trans. Instrum. Meas. 2019, 69, 4578–4591. [Google Scholar] [CrossRef]
Figure 1. Coherent (red) and non-coherent (blue) sampling spectra.
Figure 1. Coherent (red) and non-coherent (blue) sampling spectra.
Sensors 22 08170 g001
Figure 2. Spectral leakage in the DFT of the data for non-coherent sampling.
Figure 2. Spectral leakage in the DFT of the data for non-coherent sampling.
Sensors 22 08170 g002
Figure 3. Residual spectrum after removing the fundamental frequency.
Figure 3. Residual spectrum after removing the fundamental frequency.
Sensors 22 08170 g003
Figure 4. Algorithm process flow.
Figure 4. Algorithm process flow.
Sensors 22 08170 g004
Figure 5. The simulated spectrum of a 16-bit ADC under three sets of conditions (original, coherently sampled, and after processing with the proposed algorithm).
Figure 5. The simulated spectrum of a 16-bit ADC under three sets of conditions (original, coherently sampled, and after processing with the proposed algorithm).
Sensors 22 08170 g005
Figure 6. Algorithm robustness verification.
Figure 6. Algorithm robustness verification.
Sensors 22 08170 g006
Figure 7. Overall scheme of the ADC spectrum testing setup.
Figure 7. Overall scheme of the ADC spectrum testing setup.
Sensors 22 08170 g007
Figure 8. Test evaluation board.
Figure 8. Test evaluation board.
Sensors 22 08170 g008
Figure 9. Spectrum tests using different test methods: (a) full spectrum; and (b) amplified fundamental spectrum [13].
Figure 9. Spectrum tests using different test methods: (a) full spectrum; and (b) amplified fundamental spectrum [13].
Sensors 22 08170 g009
Figure 10. Test equipment for the AD9230 test.
Figure 10. Test equipment for the AD9230 test.
Sensors 22 08170 g010
Figure 11. Spectrogram of AD9230 test data.
Figure 11. Spectrogram of AD9230 test data.
Sensors 22 08170 g011
Table 1. Computation time and functionality comparison. ( J = 3413.33, M = 32,768, 16-b ADC).
Table 1. Computation time and functionality comparison. ( J = 3413.33, M = 32,768, 16-b ADC).
MethodTime (ms)Functionality
Direct DFT39Inaccurate
Hanning40Inaccurate
Blackman-Nuttall (4-term)42Inaccurate
Window in [13]87Accurate
FIRE80Accurate
Proposed Method72Accurate
Table 2. Dynamic test results.
Table 2. Dynamic test results.
Date SetMethodSNR (dB)SINAD (dB)THD (dB)SFDR (dB)
Test 1Coherent75.06575.145−90.24881.142
Non-coherent + Algorithm74.11273.373−88.56679.523
Non-coherent + Window in [13]73.55372.469−87.93577.946
Test 2Coherent75.21375.314−90.56880.997
Non-coherent + Algorithm74.24873.674−89.02379.653
Non-coherent + Window in [13]73.97371.466−88.57578.795
Test 3Coherent74.96575.665−89.84481.616
Non-coherent + Algorithm73.56874.472−88.26579.945
Non-coherent + Window in [13]72.76172.429−86.96378.149
Table 3. Results of the AD9320 test.
Table 3. Results of the AD9320 test.
Dynamic IndicatorsStandard Values in DatasheetTest Results
SNR63.8 dB (Min); 64.6 dB (Typ);62.329 dB
SINAD63.7 dB (Min); 64.5 dB (Typ);63.097 dB
ENOB10.6 bits (Typ)9.198 bits
WORST HARMONIC (Second OR Third)−82 dB (Typ); −78 dB (Max);−78.280 dB
WORST OTHER (SFDR Excluding Second and Third)−89 dB (Typ); −84 dB (Max);−85.247 dB
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peng, X.; Li, J.; Zhang, D.; Hu, C.; Sun, N.; Jiang, J. High-Precision ADC Spectrum Testing under Non-Coherent Sampling Conditions. Sensors 2022, 22, 8170. https://doi.org/10.3390/s22218170

AMA Style

Peng X, Li J, Zhang D, Hu C, Sun N, Jiang J. High-Precision ADC Spectrum Testing under Non-Coherent Sampling Conditions. Sensors. 2022; 22(21):8170. https://doi.org/10.3390/s22218170

Chicago/Turabian Style

Peng, Xiaofei, Jie Li, Debiao Zhang, Chenjun Hu, Ning Sun, and Jie Jiang. 2022. "High-Precision ADC Spectrum Testing under Non-Coherent Sampling Conditions" Sensors 22, no. 21: 8170. https://doi.org/10.3390/s22218170

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop