Digital Calibration for Gain, Time Skew, and Bandwidth Mismatch in Under-Sampling Time-Interleaved System

: This paper presents an all-digital background calibration method for gain, time skew, and bandwidth mismatch in M -channel under-sampling time-interleaved analog-to-digital converters (TI-ADCs) systems. Firstly, the characteristics of offset, gain, time skew, and bandwidth mismatch on the TI-ADCs system are analyzed. Secondly, a parameter vector is constructed to correct gain, time skew, and bandwidth mismatch. Then, the constructed parameter vector is calculated with the bandpass fractional delay ﬁlter and least squares (LS) algorithm. Based on the bandpass fractional delay ﬁlter, the proposed technique can work for ultra-high frequency signals. Additionally, the constructed parameter vector has a smaller number of ﬁlter taps than the derivative ﬁlter or Hilbert ﬁlter. Therefore, fewer computing resources are used to correct the input signal after obtaining the proposed parameter vector. Finally, there are matrix inversions in the LS algorithm. Additionally, implementing matrix inversion within FPGA is complex. For this reason, solving a system of linear equations is used to replace matrix inversions. The LS algorithm is affected by quantization error and white Gaussian noise. The simulation results verify the effectiveness of the proposed algorithm when the SNR of sub-ADC is from 30 dB to 100 dB or the ENOB of sub-ADC is from 5-bit to 16-bit. They show that the proposed algorithm is not limited by the ﬁrst sub-ADC Nyquist. Additionally, the measurement results show that the proposed method is effective in the actual time-Interleaved system.


Introduction
In modern communication systems, analog-to-digital converters play an important role.They can convert analog signals to digital signals.Among many kinds of analog-todigital converters (ADCs), the successive-approximation-register (SAR) ADCs are widely accepted due to their low design complexity.Additionally, the newest single-channel SAR analog-to-digital converter (ADC) operates at 12-bit, 1 GS/s [1].However, the center frequency of the 5G signal is above 1 GHz, which requires a higher sampling rate for the SAR ADC.
In 1980, a special framework of ADC, which is named time-interleaved analog-todigital converters (TI-ADCs) has been proposed by W.C. Black [2].Based on this special framework, an M-channel TI-ADCs system can achieve a sampling rate M times faster than a single channel.Nevertheless, there are some mismatches between different sub-ADCs in the special framework.Additionally, the mismatches include offset, gain, time skew, and bandwidth mismatch [3][4][5][6].In the TI-ADCs system, the mismatches lead to a discrepancy between the actual output digital signal and the ideal output digital signal.Therefore, these mismatches should be estimated and corrected to improve the performance of the TI-ADCs system.
Background calibration [7][8][9][10][11][12][13][14][15][16][17][18] and foreground calibration [19][20][21][22] are two of the categories for the calibration techniques of TI-ADCs.Foreground calibration means calibrating mismatches before sampling the target signal in the TI-ADCs system.Additionally, these techniques generally calibrate mismatches by a specific signal [19][20][21].The whole calibration process is completed before the input signal is sampled.If the changes in pressure, volume, and temperature lead to changes in mismatches after the foreground calibrations are completed, the changed mismatches will not be calibrated.Moreover, the performance of foreground calibration is affected by the accuracy of the specific signal.Based on these reasons, this article aims to research a background calibration in the TI-ADCs system.
The background calibration can be divided into the following two steps: mismatches estimation and mismatches correction.The offset and gain can be corrected with digital adders and digital multipliers, respectively [16].The time skew can be corrected by adjusting the sampling clock [8,9] or digitally controlled variable delay lines (VDLs) [15] or the derivative filter [10,11,14,16,17] or the adaptive filter [13].When the time skew is corrected by adjusting the sampling clock, the effect of the correction depends on the accuracy of the adjustment.Therefore, high-precision circuits are required to adjust the sampling clock in the ultra-high-speed TI-ADCs system.The time skew can be corrected with kinds of filters to avoid demanding hardware circuits.However, the derivative filter is applied to correct the time skew, which limits the size of the time skew [13].The bandwidth mismatch can be corrected by adjusting the sample and hold circuit [6] or the cascaded derivative filter or the adaptive filter [4].
Compared with digital correction technologies, using analog circuits to correct time skew and bandwidth mismatch does not involve many digital calculation modules.Therefore, analog correction technologies are suitable for on-chip TI-ADCs.However, they require state-of-the-art manufacturing technologies because of the increasing frequency of the input signals.
Figure 1 is the conceptual graph of an M-channel time-interleaved analog-to-digital conversion system.Here, CLK i is the sample clock of the sub-ADC, M is the number of channels, and Φ i is the initial phase of CLK i .o i , g i , τ i , and H ω i (jω) are offset, gain, time skew [7], and bandwidth mismatch [21], respectively.Firstly, the M sub-ADCs convert the input signal into digital signals by a time-interleaved technique [2].Secondly, the digital signals are transmitted into the field-programmable gate array (FPGA).Then, the digital signals are buffered in the FPGA.Finally, the reorganization of the signals is completed in the FPGA or computer.There are many mathematical operation modules in the FPGA or computer.Additionally, these mathematical operation modules can be used for digital calibration methods.Therefore, it is of practical significance to study digital calibration methods to reduce the accuracy requirements of analog circuits.
Background calibration [7][8][9][10][11][12][13][14][15][16][17][18] and foreground calibration [19][20][21][22] are two of the categories for the calibration techniques of TI-ADCs.Foreground calibration means calibrating mismatches before sampling the target signal in the TI-ADCs system.Additionally, these techniques generally calibrate mismatches by a specific signal [19][20][21].The whole calibration process is completed before the input signal is sampled.If the changes in pressure, volume, and temperature lead to changes in mismatches after the foreground calibrations are completed, the changed mismatches will not be calibrated.Moreover, the performance of foreground calibration is affected by the accuracy of the specific signal.Based on these reasons, this article aims to research a background calibration in the TI-ADCs system.
The background calibration can be divided into the following two steps: mismatches estimation and mismatches correction.The offset and gain can be corrected with digital adders and digital multipliers, respectively [16].The time skew can be corrected by adjusting the sampling clock [8,9] or digitally controlled variable delay lines (VDLs) [15] or the derivative filter [10,11,14,16,17] or the adaptive filter [13].When the time skew is corrected by adjusting the sampling clock, the effect of the correction depends on the accuracy of the adjustment.Therefore, high-precision circuits are required to adjust the sampling clock in the ultra-high-speed TI-ADCs system.The time skew can be corrected with kinds of filters to avoid demanding hardware circuits.However, the derivative filter is applied to correct the time skew, which limits the size of the time skew [13].The bandwidth mismatch can be corrected by adjusting the sample and hold circuit [6] or the cascaded derivative filter or the adaptive filter [4].
Compared with digital correction technologies, using analog circuits to correct time skew and bandwidth mismatch does not involve many digital calculation modules.Therefore, analog correction technologies are suitable for on-chip TI-ADCs.However, they require state-of-the-art manufacturing technologies because of the increasing frequency of the input signals.
Figure 1 is the conceptual graph of an M-channel time-interleaved analog-to-digital conversion system.Here, CLKi is the sample clock of the sub-ADC, M is the number of channels, and Φi is the initial phase of CLKi. , g , τ , and   are offset, gain, time skew [7], and bandwidth mismatch [21], respectively.Firstly, the M sub-ADCs convert the input signal into digital signals by a time-interleaved technique [2].Secondly, the digital signals are transmitted into the field-programmable gate array (FPGA).Then, the digital signals are buffered in the FPGA.Finally, the reorganization of the signals is completed in the FPGA or computer.There are many mathematical operation modules in the FPGA or computer.Additionally, these mathematical operation modules can be used for digital calibration methods.Therefore, it is of practical significance to study digital calibration methods to reduce the accuracy requirements of analog circuits.In the continuous-time wide-sense stationary (WSS) signal, the offset can be estimated using the mean, and the gain can be estimated using the standard deviation.The estimation of time skew and bandwidth mismatch can be divided into three strategies: (a) reference ADC [9]; (b) autocorrelation of the signal [13]; (c) fractional delay filter.There are the following three disadvantages to the strategy of reference ADC: (1) additional reference ADC significantly increases power consumption and hardware overhead [23].(2) The reference channel will generate additional spurs by changing the input impedance of the TI-ADC [24].(3) alternating the reference ADC with a single comparator can result in a slow convergence [25].These disadvantages do not exist in the all-digital calibration strategies of autocorrelation of the signal or fractional delay filter.However, there is high computational complexity, a limitation in the first sub-ADC bandwidth, and slow operational speed for all-digital calibration strategies, which adopt autocorrelation of the signal or fractional delay filter.
A digital calibration technique for timing skew is proposed in the reference [26].However, this technique cannot deal with the gain mismatch.In this paper, a new parameter vector is proposed to correct gain, time skew, and bandwidth mismatch simultaneously.Firstly, compared with the techniques calibrating time skew by the Taylor series [10,13], the proposed approach enables a wide range of time skew mismatch compensation.Secondly, the proposed parameter vector has a smaller number of filter taps than the derivative filter [10,13] or Hilbert filter [27].Therefore, fewer computing resources are used to correct the input signal after obtaining the proposed parameter vector.Then, ref. [16] limits its application to the TI-ADCs with specific channel numbers to ensure there are Hadamard transforms for the given channel number.However, the proposed algorithm is not limited by the number of channels.Finally, the proposed parameter vector is estimated with the bandpass fractional delay filter and least squares (LS) algorithm.Based on the bandpass fractional delay filter, the proposed technique can work for ultra-high frequency signals.
Additionally, there are matrix inversions in the LS algorithm.The implementation of matrix inversion within FPGA is complex.Therefore, solving a system of linear equations is used to replace matrix inversions.
The structure of this paper has the following parts: Section 2 describes the TI-ADCs model, the characteristics of offset, gain, time skew, and bandwidth mismatch, and the correction of these mismatches.Section 3 explains how to apply the LS algorithm to estimate gain, time skew, and bandwidth mismatch.Section 4 explains the bandpass fractional delay filter and the alternative method of the matrix inversion and validates the proposed technique by showing results from simulation and measurement.In the end, Section 5 concludes this article.

System Model and the Correction of Offset, Gain, Time Skew, and Bandwidth Mismatch 2.1. System Model
Figure 1 shows the abstract view of an M-channel time-interleaved analog-to-digital conversion system.Additionally, ADC 0 is assumed as the reference ADC.When there is an offset, gain, and time skew in the TI-ADCs system, the discrete-time Fourier Transform (DTFT) of the i-th channel output is represented as follows [7]: here X i (jω), o i , g i and ∆τ i are respectively the DTFT of the input signal before distortion, the offset, the gain, and the time skew in the i-th channel.
The inverse discrete-time Fourier transform (IDTFT) of e −jω∆τ i can be expressed as follows: When there are bandwidth mismatches in this TI-ADCs system, the DTFT of the output in the i-th channel can be represented as follows: where H ω i (jω) is the frequency response of bandwidth mismatch.Additionally, the IDTFT of H ω i (jω) can be approximated by a single-pole low-pass filter as follows [21]: where ω bi is the cutoff frequency of the sample and hold circuit of the i-th channel.
2.2.Correction of Offset, Gain, Time Skew, and Bandwidth Mismatch ADC 0 is used as the reference ADC.Additionally, offset, gain, time skew, and bandwidth mismatch of the reference ADC are assumed to be zero.Hence, Equation (3) can be rewritten as follows: here ∆o i , ∆g i and ∆τ i are, respectively, the offset, gain, and time skew relative to the reference channel.Additionally, H ∆ω i (jω) can be expressed as follows: where ω br is the cutoff frequency of the sample and the hold circuit of the reference channel.
In the TI-ADCs system, offset and gain can be corrected by the adder and multiplier, respectively.Additionally, Figure 2 is the proposed structure of correction offset and gain.In Figure 2, o i is the value of estimating the offset, g i is the value of estimation of the gain and y icalo,g [n] is the output after calibrating the offset and the gain.When there are bandwidth mismatches in this TI-ADCs system, the DTFT of the output in the i-th channel can be represented as follows: where   is the frequency response of bandwidth mismatch.Additionally, the IDTFT of   can be approximated by a single-pole low-pass filter as follows [21]: where  is the cutoff frequency of the sample and hold circuit of the i-th channel.

Correction of Offset, Gain, Time Skew, and Bandwidth Mismatch
ADC0 is used as the reference ADC.Additionally, offset, gain, time skew, and bandwidth mismatch of the reference ADC are assumed to be zero.Hence, Equation (3) can be rewritten as follows: here ∆o , ∆g and ∆τ are, respectively, the offset, gain, and time skew relative to the reference channel.Additionally,  ∆  can be expressed as follows: where  is the cutoff frequency of the sample and the hold circuit of the reference channel.
In the TI-ADCs system, offset and gain can be corrected by the adder and multiplier, respectively.Additionally, Figure 2 is the proposed structure of correction offset and gain.In Figure 2,  is the value of estimating the offset, g is the value of estimation of the gain and y ,g  is the output after calibrating the offset and the gain.When the filter is used to correct the time skew and the bandwidth mismatch at the same time, the frequency response of the filter can be represented as follows: When the filter is used to correct the time skew and the bandwidth mismatch at the same time, the frequency response of the filter can be represented as follows: The IDTFT of H ∆τ i ,∆ω i (jω) can be expressed as When there are time skews and bandwidth mismatches in the TI-ADCs system, the correction of these mismatches can be expressed as follows: Here * is the convolution and y icalo,g,τ,ω [n] is the output after calibrating the offset, gain, time skew, and bandwidth mismatch.When h ∆τ i ,∆ω i [k] is an N-tap filter, Equation ( 8) can be rewritten as follows: Based on Equation ( 9), Figure 3 shows the proposed structure of correction time skew and bandwidth mismatch.Additionally, N is even in Figure 3. Based on Figures 2 and 3, the offset, gain, time skew, and bandwidth mismatch can be corrected.
gain, time skew, and bandwidth mismatch.When ℎ ∆ ,∆  is an N-tap filter, Equation ( 8) can be rewritten as follows: , when is even.
, when is odd.
Based on Equation ( 9), Figure 3 shows the proposed structure of correction time skew and bandwidth mismatch.Additionally, N is even in Figure 3. Based on Figures 2 and 3, the offset, gain, time skew, and bandwidth mismatch can be corrected.

The Extraction of Offset, Gain, Time Skew, and Bandwidth Mismatch
The input signals in the proposed methods are wide-sense stationary.Hence, the statistical characteristics for each sub-ADC output are the same when there is no mismatch in the TI-ADCs system.

The Extraction of the Offset and Gain
The offset can be extracted from the following formula: where the offset  is approximated by the  samples of the i-th sub-ADC.Similarly, the gain can be extracted from the following formula [28]: here the gain g is approximated by the  g samples of the i-th sub-ADC.
Based on Equations ( 10) and ( 11), the offset and gain can be estimated.However, the process of estimating the gain in Equation ( 11) involves many multiplications.When the offset has been calibrated, Equation ( 5) can be rewritten as follows:

The Extraction of Offset, Gain, Time Skew, and Bandwidth Mismatch
The input signals in the proposed methods are wide-sense stationary.Hence, the statistical characteristics for each sub-ADC output are the same when there is no mismatch in the TI-ADCs system.

The Extraction of the Offset and Gain
The offset can be extracted from the following formula: where the offset o i is approximated by the L o samples of the i-th sub-ADC.Similarly, the gain can be extracted from the following formula [28]: here the gain g i is approximated by the L g samples of the i-th sub-ADC.
Based on Equations ( 10) and ( 11), the offset and gain can be estimated.However, the process of estimating the gain in Equation ( 11) involves many multiplications.When the offset has been calibrated, Equation ( 5) can be rewritten as follows: Based on the Equations ( 8) and ( 12), the correction of the gain, time skew, and bandwidth mismatch can be represented as follows: where h ∆g i ,∆τ i ,∆ω i [k] can be obtained with the following:

The Extraction of the Gain, Time Skew and Bandwidth Mismatch
To extract the gain, time skew, and bandwidth mismatch, the fractional delay filter is introduced to construct the reference.Additionally, the reference of the i-th channel can be represented as follows: Here y 0 [n] is the output of ADC 0 .And the DTFT of h ri [k] is represented as follows: where M is the number of channels.Hence, the error function of estimation gain, time skew, and bandwidth mismatch can be represented as follows: here is the coefficients of an N-tap filter.And N is even.Hence, Equation ( 17) can be represented as follows: The minimum value of ∑(e i [n]) 2 is obtained by the least square method.Hence, h ∆g i ,∆τ i ,∆ω i [k] can be represented as follows: where T is represented as matrix transposition and −1 is represented as matrix inversion.Additionally, y ical can be represented as follows: Here L s is the number of samples.Additionally, y ri can be represented as follows: . .
Hence, the gain, time skew, and bandwidth mismatch can be estimated with Equations ( 19)- (21).Based on Equations ( 10) and ( 19), Figures 2 and 3, the calibration of offset, gain, time skew, and bandwidth mismatch can be shown in Figure 4.

Test Verification
In summary, the offset can be extracted with Equation (10). Figure 4 shows how to calibrate gain, time skew, and bandwidth mismatch.

A Fixed Fractional Delay Filter
The proposed method adopts the fractional delay filters to estimate the gain, time skew, and bandwidth mismatch in each channel.The fractional delay filters are used in each sub-ADC.Therefore, the input signal is limited to the first Nyquist bandwidth of the sub-ADC.To make the proposed algorithm work for the high Nyquist bandwidth of the sub-ADC in the TI-ADCs system, the fractional delay filters should be effective in high sub-ADC Nyquist bandwidth.
When the bandpass (BP) signals are in the  -th Nyquist bandwidth, the frequency can be expressed as follows: Here  1,  , and  are the low and high cutoff frequencies of the signal, respectively.Additionally,  is the sampling frequency of sub-ADC.
The BP signal is constructed by the sum of two complex signals of   and   .
Therefore, the frequency response of the complex signals can be respectively represented as X  and X  .According to the band-pass sampling theory, the original spectrum of the band-pass signal can be folded back to the baseband in the under-sampled TI-ADC converter whose sampling frequency conforms to the Nyquist criterion.In the odd-order Nyquist band, the baseband spectrum of the under-sampling signal has the same shape as the original spectrum.Additionally, the baseband spectrum of even order NB input signal is reversed.In the odd-order Nyquist band, X  can be represented as follows: Here, Y  and Y  is the output of X  and X  , respectively.Furthermore,  can be obtained as follows: where the input signal is the  -th Nyquist bandwidth and ⌊ /2⌋ is the result of rounding down  /2.Therefore, the inverse fast Fourier transform (IFFT) of X  can be expressed as

Test Verification
In summary, the offset can be extracted with Equation (10). Figure 4 shows how to calibrate gain, time skew, and bandwidth mismatch.

A Fixed Fractional Delay Filter
The proposed method adopts the fractional delay filters to estimate the gain, time skew, and bandwidth mismatch in each channel.The fractional delay filters are used in each sub-ADC.Therefore, the input signal is limited to the first Nyquist bandwidth of the sub-ADC.To make the proposed algorithm work for the high Nyquist bandwidth of the sub-ADC in the TI-ADCs system, the fractional delay filters should be effective in high sub-ADC Nyquist bandwidth.
When the bandpass (BP) signals are in the N BP -th Nyquist bandwidth, the frequency can be expressed as follows: Here N BP ≥ 1, f L , and f H are the low and high cutoff frequencies of the signal, respectively.Additionally, f s is the sampling frequency of sub-ADC.
The BP signal is constructed by the sum of two complex signals of x + (t) and x − (t).Therefore, the frequency response of the complex signals can be respectively represented as X + (f ) and X − (f ).According to the band-pass sampling theory, the original spectrum of the band-pass signal can be folded back to the baseband in the under-sampled TI-ADC converter whose sampling frequency conforms to the Nyquist criterion.In the odd-order Nyquist band, the baseband spectrum of the under-sampling signal has the same shape as the original spectrum.Additionally, the baseband spectrum of even order NB input signal is reversed.In the odd-order Nyquist band, X(f ) can be represented as follows: Here, Y + (f ) and Y − (f ) is the output of X + (f ) and X − (f ), respectively.Furthermore, n n can be obtained as follows: where the input signal is the N BP -th Nyquist bandwidth and N BP /2 is the result of rounding down N BP /2.Therefore, the inverse fast Fourier transform (IFFT) of X(f ) can be expressed as According to the Euler formula, x(t) can be rewritten as follows: In the time domain, y + (t) and y − (t) can be obtained as follows: According to Equation ( 27), the reference of the i-th channel can be represented as follows: Therefore, the DFT of Equation ( 28) can be expressed as follows: where, y 0 [n + i/Mf s ] can be represented as follows: And y h0 [n + i/Mf s ] can be represented as follows: In the even Nyquist bands, the DFT of x i (t + i/Mf s ) can be expressed as follows: Figure 5 shows the reference generation of the TI-ADC system with BP signal input.The IDFT of Hilbert filter frequency response can be expressed as follows: Additionally, the causal filters' coefficients are obtained by windowing Equation (33).

Simulation Result
The signal-to-noise distortion ratio (SNDR) is a metric that describes the signal performance.Additionally, SNDR can be calculated by the following: where P Signal is the power of the ideal signal, P Noise is the power of the noise signal, and P Distortion is the spurious power of the output signal.In the TI-ADCs system, the presence of mismatches leads to an increase in P Noise .When and P Noise are fixed, the smaller the calibration error, the higher the SNDR.Therefore, the SNDR is applied to evaluate the effect of error calibration.

Simulation Result
The signal-to-noise distortion ratio (SNDR) is a metric that describes the signal performance.Additionally, SNDR can be calculated by the following: where  is the power of the ideal signal,  is the power of the noise signal, and  is the spurious power of the output signal.In the TI-ADCs system, the presence of mismatches leads to an increase in . When and  are fixed, the smaller the calibration error, the higher the SNDR.Therefore, the SNDR is applied to evaluate the effect of error calibration.
The calibration of offset is simpler than other mismatches.Hence, the offset is set to 0 in the simulation tests.The 12-bit 8 GS/s 8-channel TI-ADCs simulation system is established according to Equation (3) with MATLAB.Additionally, the mismatches of the simulation model are shown in Table 1.Moreover, the parameters of the bandpass fractional filters are shown in Table 2. Additionally, the number of taps of ℎ ∆g ,∆ ,∆  is 73.In the 12-bit 8 GS/s 8-channel TI-ADCs simulation system, the SNDR of the input signal before sampling is 60 dB because of the Gaussian white noise in each channel.

Sub-ADC Number
Offset Gain Time Skew Bandwidth 1 0 1 0 Parameters for the filters.

Parameter Filter Design Method Number of Taps Fractional delay filters
Blackman Window 73 Hilbert transform filters Blackman Window 73 Figure 6 shows the values of SNDR and spurious-free dynamic range (SFDR) for sinusoidal signals with the different frequencies of the first Nyquist bandwidth.The fre- The calibration of offset is simpler than other mismatches.Hence, the offset is set to 0 in the simulation tests.The 12-bit 8 GS/s 8-channel TI-ADCs simulation system is established according to Equation (3) with MATLAB.Additionally, the mismatches of the simulation model are shown in Table 1.Moreover, the parameters of the bandpass fractional filters are shown in Table 2. Additionally, the number of taps of h ∆g i ,∆τ i ,∆ω i [k] is 73.In the 12-bit 8 GS/s 8-channel TI-ADCs simulation system, the SNDR of the input signal before sampling is 60 dB because of the Gaussian white noise in each channel.
Table 1.The mismatches of the TI-ADCs simulation system.

Sub-ADC Number
Offset Gain Time Skew Bandwidth Figure 6 shows the values of SNDR and spurious-free dynamic range (SFDR) for sinusoidal signals with the different frequencies of the first Nyquist bandwidth.The frequency range is limited within 0.05-0.42GHz.Additionally, the frequency spacing between two single-tone sinusoidal signals is 0.01 GHz.Where f c is the sampling frequency of sub-ADC.In Figure 6a, the SNDR is stable above 65 dB in the test frequency range.Additionally, the SFDR is stable above 80 dB in Figure 6b.
between two single-tone sinusoidal signals is 0.01 GHz.Where  is the sampling frequency of sub-ADC.In Figure 6a, the SNDR is stable above 65 dB in the test frequency range.Additionally, the SFDR is stable above 80 dB in Figure 6b.The mismatches shown in Table 3 are set to verify that the proposed methods are still effective when there are gain, time skew, and bandwidth mismatch at the same time.The frequency range is limited within 0.05-0.42GHz.Where  is the frequency of sub-ADC, and  is the frequency of the TI-ADCs system.Additionally, Figure 6 shows the SNDR and SFDR before and after calibration.It can be seen from Figure 7 that the proposed calibration method still has a good effect in the cases of gain, time skew, and bandwidth mismatch.The SNDR before calibration does not exceed 20 dB.Additionally, the SNDR is above 65 dB after calibration.Therefore, the value of SNDR is improved by at least 45 dB after calibration.Moreover, the value of SFDR is above 80 dB after calibration.The mismatches shown in Table 3 are set to verify that the proposed methods are still effective when there are gain, time skew, and bandwidth mismatch at the same time.The frequency range is limited within 0.05-0.42GHz.Where f c is the frequency of sub-ADC, and f s is the frequency of the TI-ADCs system.Additionally, Figure 6 shows the SNDR and SFDR before and after calibration.It can be seen from Figure 7 that the proposed calibration method still has a good effect in the cases of gain, time skew, and bandwidth mismatch.The SNDR before calibration does not exceed 20 dB.Additionally, the SNDR is above 65 dB after calibration.Therefore, the value of SNDR is improved by at least 45 dB after calibration.Moreover, the value of SFDR is above 80 dB after calibration.Figures 8 and 9 show the calibration performance within the second and third Nyquist bandwidths of the sub-ADC, respectively.Additionally, the frequency ranges are limited within 0.55-0.92GHz and 1.05-1.42GHz, respectively.The offset, gain, time skew, and bandwidth mismatch of the TI-ADCs system are shown in Table 3. Figures 8 and 9 show the calibration performance within the second and third Nyquist bandwidths of the sub-ADC, respectively.Additionally, the frequency ranges are limited within 0.55-0.92GHz and 1.05-1.42GHz, respectively.The offset, gain, time skew, and bandwidth mismatch of the TI-ADCs system are shown in Table 3. Figures 8 and 9 show the calibration performance within the second and third Nyquist bandwidths of the sub-ADC, respectively.Additionally, the frequency ranges are limited within 0.55-0.92GHz and 1.05-1.42GHz, respectively.The offset, gain, time skew, and bandwidth mismatch of the TI-ADCs system are shown in Table 3. Figures 8 and 9 show the calibration performance within the second and third Nyquist bandwidths of the sub-ADC, respectively.Additionally, the frequency ranges are limited within 0.55-0.92GHz and 1.05-1.42GHz, respectively.The offset, gain, time skew, and bandwidth mismatch of the TI-ADCs system are shown in Table 3.It can be seen from Figures 7-9 that the SNDR and SFDR of the TI-ADCs system decrease with the increase in the input signal frequency at the same gain, time skew, and bandwidth mismatch.Figures 8 and 9 show that even if the input signal is not located in the first Nyquist bandwidth of the sub-ADC, the proposed methods still have a good calibration effect.Furthermore, Figures 7-9 show that the proposed methods are effective for bandpass signals with different Nyquist bandwidths.Additionally, the SNDR is above 65 dB and the SFDR is above 80 dB in the first, second, and third Nyquist frequency bands.Therefore, the frequency of the input signal has relatively little effect on the proposed technique.
Figure 10 shows the frequency spectrums of the multi-tone signal before and after calibration.Where the frequencies of the multi-tone signal are [2.30GHz, 2.42 GHz, 2.54 GHz, 2.66 GHz, 2.78 GHz]. Figure 10a shows the frequency spectrum of the TI-ADCs system before calibration.Where the offset, gain, time skew, and bandwidth mismatch are shown in Table 3.The performance of the TI-ADCs system shown in Figure 10a is poor.Where the SNDR is −2.77dB and the SFDR is −17.28 dB.Additionally, the performance of the TI-ADCs system is well-improved after calibration.In Figure 10b, the SNDR is 60.19 dB and the SFDR is 75.53 dB.Therefore, the SNDR is improved by 62.96 dB, and the SFDR is improved by 92.81 dB. Figure 10 shows that the proposed calibration algorithm is still effective for multi-tone bandpass signals with high sub-ADC Nyquist bandwidth.
calibration.Where the frequencies of the multi-tone signal are [2.30GHz, 2.42 GHz, 2.54 GHz, 2.66 GHz, 2.78 GHz]. Figure 10a shows the frequency spectrum of the TI-ADCs system before calibration.Where the offset, gain, time skew, and bandwidth mismatch are shown in Table 3.The performance of the TI-ADCs system shown in Figure 10a is poor.Where the SNDR is −2.77dB and the SFDR is −17.28 dB.Additionally, the performance of the TI-ADCs system is well-improved after calibration.In Figure 10b, the SNDR is 60.19 dB and the SFDR is 75.53 dB.Therefore, the SNDR is improved by 62.96 dB, and the SFDR is improved by 92.81 dB. Figure 10 shows that the proposed calibration algorithm is still effective for multi-tone bandpass signals with high sub-ADC Nyquist bandwidth.There are matrix inversions in the LS algorithm.Additionally, the implementation of matrix inversion in FPGA is complex.Therefore, Equation ( 19) is rewritten as follows: Hence, a system of linear equations can be constructed as follows: .
According to Cramer's rule [29], the solution of ℎ ∆g ,∆ ,∆  requires about !  1  1 times of addition and !  1  times of multiplication.In this way, the matrix inversion can be avoided when solving ℎ ∆g ,∆ ,∆  .Additionally, the solution of  requires about   1  times of addition and   times of multiplication.Similarly, the solution of  requires about   1 times of addition and  times of multiplication.Therefore, the estimation in the i-th channel of gain, time skew, and bandwidth mismatch requires about  1   1  1 !  1  1  times of addition and   1  !  1  times of multiplication.Where L is the number of bandpass fractional delay filter taps.There are matrix inversions in the LS algorithm.Additionally, the implementation of matrix inversion in FPGA is complex.Therefore, Equation ( 19) is rewritten as follows: Hence, a system of linear equations can be constructed as follows: According to Cramer's rule [29], the solution of h ∆g i ,∆τ i ,∆ω i [k] requires about N!(N − 1)(N + 1) times of addition and N!(N + 1)N times of multiplication.In this way, the matrix inversion can be avoided when solving h ∆g i ,∆τ i ,∆ω i [k].Additionally, the solution of A N×N requires about N(L s − 1)N times of addition and NL s N times of multiplication.Similarly, the solution of b N×1 requires about N(L s − 1) times of addition and NL s times of multiplication.Therefore, the estimation in the i-th channel of gain, time skew, and bandwidth mismatch requires about (N + 1)(N(L s − 1) + (N − 1)N!) + (M − 1)(L − 1)L s times of addition and N(N + 1)(L s + N!) + (M − 1)LL s times of multiplication.Where L is the number of bandpass fractional delay filter taps.
The comparison results of the proposed technique with the prior state-of-the-art are summarized in Table 4.The bandwidth mismatch can be calibrated with the proposed technique.In the reference [13], the calibration of time skew is based on the correlation between channels.However, Equation (36) of reference [13] is established under the necessary condition that the time skew of the adjacent channels is much smaller than the sampling period of the sub-ADC.To summarize, the calibration techniques of reference [13] must meet the necessary condition that the time skew between adjacent channels or compared with the reference channel is much smaller than the sampling period of the sub-ADC.Similarly, there are restrictions on the size of the time skew in the references [30,31].References [16,27] and this work approach enables a wide range of time skew mismatch compensation.Therefore, references [16,27], and this work have a good effect on the improvement of SNDR and SFDR.Work [13,16,27], and this work needs 1045.5 K, 277.5 K, 720 K, and 128.1 K multiplication modules to estimate time skew, respectively.Additionally, the multiplication modules occupy most of the computing resources in the techniques of time skew estimation or correction.Therefore, it can be considered that the proposed estimation technique in this work has the lowest power consumption.Work [13], work [16], work [27], work [30], work [31], and this work needs 31, 33, 12, 25, 21, and 5 multiplication modules to correct time skew respectively.Therefore, the power consumption of the proposed technique of correction time skew is far less than techniques [13,16,27,30,31].ENOB represents the final performance of the TI-ADC system after calibration.Compared with work [16,27], and [30], the proposed technique has the best performance in ENOB.The quantization error is ignored in the work [31].Therefore, the ENOB after calibration is not compared between work [31] and others.  T is time skew. 4B is bandwidth mismatch.M is the channel number.L is filter length.L s is the number of samples.N = 5 is the length of the calibration coefficient.
In Equation (19), quantization noise and Gaussian white noise are not considered.Therefore, the calibration of time skew is affected by quantization noise and Gaussian white noise.Figure 11 shows the performances of calibration time skew in a 2-channel 6.4 GS/s TI-ADCS under different white noise levels.
O is offset.G is gain.T is time skew.B is bandwidth mismatch.M is the channel number.L is filter length.s L is the number of samples.

5
N  is the length of the calibration coefficient.
In Equation (19), quantization noise and Gaussian white noise are not considered.Therefore, the calibration of time skew is affected by quantization noise and Gaussian white noise.Figure 11 shows the performances of calibration time skew in a 2-channel 6.4 GS/s TI-ADCS under different white noise levels.As shown in Figure 11, there is a positive correlation basically between the SNDR of TI-ADCS and the SNR of sub-ADC when the SNR of sub-ADC is in the range of 0~100 dB.Therefore, the proposed technique can effectively calibrate the gain, time skew, and bandwidth mismatch at the SNR of 0~100 dB.In Figure 12, there is a positive correlation between the SNDR of TI-ADCS and the ENOB of sub-ADC when the ENOB of sub-ADC is in the range of 1 bit~16 bit.It can be seen from Figures 11 and 12 that the proposed technique is effective when the SNR of sub-ADC is in the range of 0~100 dB or the ENOB of sub-ADC is in the range of 1~16 bits.As shown in Figure 11, there is a positive correlation basically between the SNDR of TI-ADCS and the SNR of sub-ADC when the SNR of sub-ADC is in the range of 0~100 dB.Therefore, the proposed technique can effectively calibrate the gain, time skew, and bandwidth mismatch at the SNR of 0~100 dB.In Figure 12, there is a positive correlation between the SNDR of TI-ADCS and the ENOB of sub-ADC when the ENOB of sub-ADC is in the range of 1 bit~16 bit.It can be seen from Figures 11 and 12 that the proposed technique is effective when the SNR of sub-ADC is in the range of 0~100 dB or the ENOB of sub-ADC is in the range of 1~16 bits.The minimum value of ∑   is obtained by the least square method.Hence, the cumulative error ∑   is affected by quantization error and white Gaussian noise.In Figures 13 and 14, the cumulative errors of the proposed technique are respectively shown under different quantization errors and Gaussian white noise of sub-ADC.It can be seen from Figures 13 and 14 that the cumulative error is stable and approaches 0 when the SNR is greater than 30 dB and the ENOB is greater than 5 bits.It can be concluded from Figures 11, 12, 13 and 14 that the proposed technique will be affected by the SNR and ENOB.However, the proposed technique has a good effect in the 30~100 dB of SNR or the 5~16 bit of ENOB.The minimum value of ∑(e i [n]) 2 is obtained by the least square method.Hence, the cumulative error ∑(e i [n]) 2 is affected by quantization error and white Gaussian noise.In Figures 13 and 14, the cumulative errors of the proposed technique are respectively shown under different quantization errors and Gaussian white noise of sub-ADC.It can be seen from Figures 13 and 14 that the cumulative error is stable and approaches 0 when the SNR is greater than 30 dB and the ENOB is greater than 5 bits.It can be concluded from Figures 11-14 that the proposed technique will be affected by the SNR and ENOB.However, the proposed technique has a good effect in the 30~100 dB of SNR or the 5~16 bit of ENOB.  the cumulative error is stable and approaches 0 when the SNR is greater than 30 dB and the ENOB is greater than 5 bits.It can be concluded from Figures 11, 12, 13 and 14 that the proposed technique will be affected by the SNR and ENOB.However, the proposed technique has a good effect in the 30~100 dB of SNR or the 5~16 bit of ENOB.

Measurement Result
A commercial 12-bit 800 Msps 2-channel TI-ADC was used to validate the propo calibration algorithm under real conditions.The measurement setup is illustrated in ure 15, and it consists of a signal generator (SDG6032X), a commercial 2-channel TI-A (ADC12DL3200EVM), and a field-programmable gate array (TSWDL3200) for data tr mission.

Measurement Result
A commercial 12-bit 800 Msps 2-channel TI-ADC was used to validate the proposed calibration algorithm under real conditions.The measurement setup is illustrated in Figure 15, and it consists of a signal generator (SDG6032X), a commercial 2-channel TI-ADCs (ADC12DL3200EVM), and a field-programmable gate array (TSWDL3200) for data transmission.
The features of experimental equipment can be shown in Table 5. Figure 16 shows the frequency spectrum before and after calibration.Additionally, the input is a single sinusoidal signal.The frequency of the input signal is 70 MHz.As shown in Figure 16, the SNDR and SFDR are, respectively, 17.15 dB and 17.17 dB.Additionally, the SNDR and SFDR are respectively 43.19 dB and 65.19 dB after calibration by the proposed algorithm.Additionally, it can also be seen from Figure 16 that the spurs of time skew drop below the noise floor.Therefore, the measurement result demonstrates that the proposed method is effective against time skew.

Measurement Result
A commercial 12-bit 800 Msps 2-channel TI-ADC was used to validate the proposed calibration algorithm under real conditions.The measurement setup is illustrated in Figure 15, and it consists of a signal generator (SDG6032X), a commercial 2-channel TI-ADCs (ADC12DL3200EVM), and a field-programmable gate array (TSWDL3200) for data transmission.The features of experimental equipment can be shown in Table 5. Figure 16 shows the frequency spectrum before and after calibration.Additionally, the input is a single sinusoidal signal.The frequency of the input signal is 70 MHz.As shown in Figure 16, the SNDR and SFDR are, respectively, 17.15 dB and 17.17 dB.Additionally, the SNDR and SFDR are respectively 43.19 dB and 65.19 dB after calibration by the proposed algorithm.Additionally, it can also be seen from Figure 16 that the spurs of time skew drop below the noise floor.Therefore, the measurement result demonstrates that the proposed method is effective against time skew.

Parameter Test Conditions Typical Values Unit
The In Table 6, the calibration effects of [13,16,27,30,31], and this paper are shown.After calibration, there is the highest ENOB in the proposed technique compared with [13,16,27,30,31].Figure 16 shows the frequency spectrum of the commercial 2-channel TI- In Table 6, the calibration effects of [13,16,27,30,31], and this paper are shown.After calibration, there is the highest ENOB in the proposed technique compared with [13,16,27,30,31].Figure 16 shows the frequency spectrum of the commercial 2-channel TI-ADCs before and after calibration for the two-tone sinusoidal signals.The frequencies of the input signal are 70 MHz and 170 MHz.As shown in Figure 17, the SNDR and SFDR are respectively 11.83 dB and 10.03 dB.Additionally, the SNDR and SFDR are, respectively, 42.38 dB and 56.35 dB after calibration by the proposed algorithm.Therefore, the proposed algorithm also has a good calibration effect for actual two-tone signals.

Figure 1 .
Figure 1.Conceptual graph of an M-channel time-interleaved analog-to-digital conversion system.Figure 1. Conceptual graph of an M-channel time-interleaved analog-to-digital conversion system.

Figure 1 .
Figure 1.Conceptual graph of an M-channel time-interleaved analog-to-digital conversion system.Figure 1. Conceptual graph of an M-channel time-interleaved analog-to-digital conversion system.

Figure 2 .
Figure 2. Proposed structure of correction offset and gain.

Figure 2 .
Figure 2. Proposed structure of correction offset and gain.

Figure 3 .
Figure 3. Proposed structure of correction time skew and bandwidth mismatch.

Figure 3 .
Figure 3. Proposed structure of correction time skew and bandwidth mismatch.

Figure 4 .
Figure 4. Proposed structure of correction of the offset, gain, time skew, and bandwidth mismatch.

Figure 4 .
Figure 4. Proposed structure of correction of the offset, gain, time skew, and bandwidth mismatch.

Figure 6 .
Figure 6.Performances of bandwidth mismatch calibration with the proposed method.(a) The SNDR of TI-ADCs with different frequencies; (b) the SFDR of TI-ADC with different frequencies.

Figure 6 .
Figure 6.Performances of bandwidth mismatch calibration with the proposed method.(a) The SNDR of TI-ADCs with different frequencies; (b) the SFDR of TI-ADC with different frequencies.

Figure 7 .
Figure 7. Performances of the proposed method.(a) The SNDR of TI-ADCs with different frequencies; (b) the SFDR of TI-ADC with different frequencies.

Figure 7 .
Figure 7. Performances of the proposed method.(a) The SNDR of TI-ADCs with different frequencies; (b) the SFDR of TI-ADC with different frequencies.

Figure 7 .
Figure 7. Performances of the proposed method.(a) The SNDR of TI-ADCs with different frequencies; (b) the SFDR of TI-ADC with different frequencies.

Figure 8 .Figure 9 .Figure 8 .Figure 7 .
Figure 8. Performances of the proposed method in the second Nyquist bandwidth of the sub-ADC.(a) The SNDR of TI-ADCs with different frequencies; (b) the SFDR of TI-ADC with different frequencies.

Figure 8 .Figure 9 .Figure 9 .
Figure 8. Performances of the proposed method in the second Nyquist bandwidth of the sub-ADC.(a) The SNDR of TI-ADCs with different frequencies; (b) the SFDR of TI-ADC with different frequencies.

Figure 10 .
Figure 10.Performances of mismatches calibration for a multi-tone signal.(a) Before calibration; (b) after calibration.

Figure 10 .
Figure 10.Performances of mismatches calibration for a multi-tone signal.(a) Before calibration; (b) after calibration.

Figure 11 .
Figure 11.The performance of the proposed technique with different SNR of the input signal.

Figure 11 .
Figure 11.The performance of the proposed technique with different SNR of the input signal.

20 Figure 12 .
Figure 12.The performance of the proposed technique with different ENOB of the sub-ADC.

Figure 12 .
Figure 12.The performance of the proposed technique with different ENOB of the sub-ADC.

Figures 13 and 14 ,
Figures13 and 14, the cumulative errors of the proposed technique are respectively shown under different quantization errors and Gaussian white noise of sub-ADC.It can be seen from Figures13 and 14that the cumulative error is stable and approaches 0 when the SNR is greater than 30 dB and the ENOB is greater than 5 bits.It can be concluded from Figures11, 12, 13 and 14 that the proposed technique will be affected by the SNR and ENOB.However, the proposed technique has a good effect in the 30~100 dB of SNR or the 5~16 bit of ENOB.
Figures13 and 14, the cumulative errors of the proposed technique are respectively shown under different quantization errors and Gaussian white noise of sub-ADC.It can be seen from Figures13 and 14that the cumulative error is stable and approaches 0 when the SNR is greater than 30 dB and the ENOB is greater than 5 bits.It can be concluded from Figures11, 12, 13 and 14 that the proposed technique will be affected by the SNR and ENOB.However, the proposed technique has a good effect in the 30~100 dB of SNR or the 5~16 bit of ENOB.

Figure 13 .
Figure 13.The cumulative errors under different SNR of the sub-ADC in the proposed technique.Figure 13.The cumulative errors under different SNR of the sub-ADC in the proposed technique.

Figure 13 . 16 Figure 14 .
Figure 13.The cumulative errors under different SNR of the sub-ADC in the proposed technique.Figure 13.The cumulative errors under different SNR of the sub-ADC in the proposed technique.Appl.Sci.2022, 12, x FOR PEER REVIEW 16

Figure 14 .
Figure 14.The cumulative errors under different ENOB of the sub-ADC in the proposed technique.

Figure 15 .
Figure 15.Experiment setup which is aimed to validate the proposed method.

Figure 15 .
Figure 15.Experiment setup which is aimed to validate the proposed method.

Figure 16 .
Figure 16.The frequency spectrum of the commercial 2-channel TI-ADCs before and after calibration.

Figure 16 .
Figure 16.The frequency spectrum of the commercial 2-channel TI-ADCs before and after calibration.

Figure 17 .
Figure 17.The frequency spectrum of the two-tone signals in the commercial 2-channel TI-ADCs before and after calibration.

Table 1 .
The mismatches of the TI-ADCs simulation system.

Table 2 .
Parameters for the filters.

Table 3 .
The mismatches of the TI-ADCs system.

Table 3 .
The mismatches of the TI-ADCs system.

Table 4 .
Compared with other techniques.

Table 5 .
The features of experimental equipment.
Analog differential input full-scale range Default 800 mV Appl.Sci.2022, 12, x FOR PEER REVIEW 17 of 20

Table 5 .
The features of experimental equipment

Table 6 .
Comparison of calibration effects in the ADC12DL3200EVM.

Table 6 .
Comparison of calibration effects in the ADC12DL3200EVM