Mixed-Domain Adaptive Blind Correction of High-Resolution Time-Interleaved ADCs
Mixed-Domain Adaptive Blind Correction of High-Resolution Time-Interleaved ADCs
ETRI Journal. 2014. Oct, 36(6): 894-904
Copyright © 2014, Electronics and Telecommunications Research Institute(ETRI)
  • Received : January 24, 2014
  • Accepted : July 18, 2014
  • Published : October 01, 2014
Export by style
Cited by
About the Authors
Munkyo Seo
Eunsoo Nam
Mark Rodwell

Blind mismatch correction of time-interleaved analog-to-digital converters (TI-ADC) is a challenging task. We present a practical blind calibration technique for low-computation, low-complexity, and high-resolution applications. Its key features are: dramatically reduced computation; simple hardware; guaranteed parameter convergence with an arbitrary number of TI-ADC channels and most real-life input signals, with no bandwidth limitation; multiple Nyquist zone operation; and mixed-domain error correction. The proposed technique is experimentally verified by an M = 4 400 MSPS TI-ADC system. In a single-tone test, the proposed practical blind calibration technique suppressed mismatch spurs by 70 dB to 90 dB below the signal tone across the first two Nyquist zones (10 MHz to 390 MHz). A wideband signal test also confirms the proposed technique.
I. Introduction
A time-interleaved analog-to-digital converter (TI-ADC) is a scalable architecture for very high sampling rates. A number of subconverters cyclically sample the input signal, and multiple outputs are combined to yield a single digital stream ( Fig. 1 (a) ). The collective sampling rate is, therefore, proportional to the number of subconverters. It has been well known, however, that the spectral performance of a TI-ADC is limited by aliasing spectra due to mismatches in subconverter gain, sampling time, and so on [1] [6] .
Currently known mismatch correction techniques can be categorized into training (foreground) [2] , [6] and blind (background) methods [7] [26] . Training methods are suitable for high-resolution application in general, since they are capable of correcting general linear mismatches [6] , but at the cost of suspension of data acquisition during each calibration. They are also subject to post-calibration detuning due to temperature variation, aging, and so on [6] . Blind methods, on the other hand, use normal input signals for calibration purposes; therefore, they do not require a dedicated calibration period. Errors that slowly vary as a function of time may also be tracked. There have been proposed a variety of blind methods with different accomplishments and limitations. One group of techniques performs error detection and error correction entirely in the digital domain [7] [17] , as illustrated in Fig. 1(b) . Another class of techniques use both the analog and digital domain [18] [26] .
Previous mixed-domain methods typically involve special analog signal processing (for example, adding a known signal to the input, as in Fig. 1 (c) ), to facilitate mismatch estimation, which may potentially compromise input signal integrity. Purely digital techniques keep the analog signal path intact, but their computational cost is highly demanding. This is partly due to the complex parameter search algorithm, but mostly due to digital correction of sampling time delay, since matrix inversion or transformation between the time and frequency domains is usually required for correction filter calculation. Further, the resulting correction filter is long (sinc( t ) decays only as 1/ t ), increasing hardware cost and complexity. Depending on blind correction algorithms, the estimation of timing error may also be computationally complex. For both classes of techniques, special assumptions are necessary (for example, reduced input bandwidth, limited number of TI-ADC channels, or wide-sense stationarity) to make the problem of blind estimation solvable. Full-digital techniques tend to require stronger assumptions; thus, they are more restrictive than mixed-domain methods. Using an extra ADC may provide a convenient calibration reference, but the analog input path is subject to switching, which raises concerns about signal integrity.
The above discussion naturally leaves us one interesting option in the mixed-domain: timing error correction in the analog domain by directly tuning sampling clocks ( Fig. 1 (d) ). This obviates long digital filters as well as their online calculation for timing correction. In addition, timing correction can now be perfect over multiple Nyquist zones. In contrast, a single set of digital filters can only provide approximate timing correction over a single Nyquist zone, due to limited usable bandwidth and in-band ripples, among others. We must still develop a reliable blind estimation method to close the feedback loop, preferably with weaker working assumptions so as to allow application to the widest range of signals.
PPT Slide
Lager Image
(a) TI-ADC system, (b) TI-ADC with full digital-domain blind correction, (c) conventional mixed-domain method with analog input pre-processing, and (d) proposed mixed-domain method with sampling clock tuning.
In accordance with the above discussion, the authors previously reported a new mixed-domain blind method with analog tuning of sampling clocks [26] . This paper presents a full-fledged discussion with theory and new experimental results. Dramatically reduced computational complexity and exceptionally wide applicability are among the contributions of the paper. The proposed blind method is based on the assumption that the input signal is wide-sense stationary (WSS). Under the input WSS assumption, the mismatch estimates are guaranteed to converge to true parameters. No further restriction on the input signal is necessary. Specifically, the TI-ADC can have an arbitrary number of channels, with no need of an additional subconverter for calibration reference. The input spectrum can cover the full Nyquist bandwidth. In fact, the proposed blind method works in any Nyquist zone.
Similar mixed-domain approaches have been proposed [23] [25] , with their own assumptions and restrictions: it is unclear whether the correction method in [23] applies to wideband signals. The proposed method in [24] only applies to random data; for example, digitally modulated signals. An additional subconverter channel is necessary in [25] for calibration reference.
Section II describes our system model. Section III introduces the proposed error detection method, and Section IV develops the proposed error detection method’s adaptive implementation for parameter estimation. Section V discusses the experimental results, and Section VI concludes the paper.
II. System Model
Figure 2 shows a block diagram of the M = 4 TI-ADC with the proposed mismatch correction scheme. Each of the four subconverters successively samples the input signal, x ( t ), every 4 Ts such that the overall sampling rate is fs (= 1/ Ts ). Listed below are pertinent assumptions and clarifications relating to the proposed mismatch correction scheme.
  • ▪ The inputx(t) is WSS and bandlimited from dc tofs/2. No further information aboutx(t) is known.
  • ▪ The mismatch in subconverter dc offsets is independently corrected by first measuring and then subsequent subtracting.
  • ▪ Thekth channel subconverter hasintrinsicgainG*kand sampling time error Δt*k, both of which are unknown.
  • ▪ Theestimateof the intrinsic gain and sampling time errors isG˜kandΔt˜k, respectively.
  • ▪ Correction of sampling time mismatches is achieved by tuning individual sampling clocks to the estimateΔt˜k.
  • ▪ Correction of gain error is performed by digitally dividing the subconverter output by the gain estimateG˜k.
  • ▪Residualmismatch error is defined as the difference between intrinsic and estimated parameters.
  • ▪ The magnitude of intrinsic mismatches is small. The precise interpretation will be made clear in the context.
The output of the k th subconverter, yk [ n ], is
y k [ n ]= G k * x( ( Mn+k ) T s +Δ t k * Δ t ˜ k ),
where M is the number of TI-ADC channels. After gain correction, we have
z k [ n ]= y k [ n ] G ˜ k .
Therefore, the subsequence zk [ n ] is a scaled, time-shifted, and undersampled version of x ( t ) and is given as
z k [ n ]=[ G k * G ˜ k ][ x( ( Mn+k ) T s +Δ t k * Δ t ˜ k ) ]  for  k=0,  ...  ,M1.
PPT Slide
Lager Image
M = 4 TI-ADC system with proposed mismatch correction scheme.
The overall TI-ADC output, z [ n ], is obtained by taking zk [ n ]’s in a cyclic fashion as follows:
z[ n ]= z n mod M [ floor( n M ) ],
where floor( x ) is the greatest integer less than or equal to x .
- 1. TI-ADC Mismatch Model
If there is no residual mismatch, then it follows from (2)–(3) that z [ n ] = x ( nTs ) — a perfect reconstruction of the input. Otherwise, z [ n ] is modulated by residual gain or sampling time errors. In the frequency domain, this modulation manifests itself as frequency-shifted input spectra, known as aliasing error. Aliasing effects of gain and sampling time error have been extensively studied [1] , [3] [5] . In this section, we briefly review the generalized mismatch model [6] . Let Hk ( f ) be the frequency response of the k th channel (which includes static gain and time delay as a special case). The spectrum of z [ n ] consists of linearly weighted frequency-shifted input spectra [6] and is as follows:
Z( e j2πf )= m=0 M1 c m ( f m M f s )X( f m M f s ) ,
where Z (e j2πf ) and X ( f ) are the Fourier transform of z [ n ] and x ( t ), respectively. The conversion gains, cm ( f )’s, are obtained by taking the discrete Fourier transform (DFT) of the Hk ( f )’s with respect to k . If there is no mismatch, then the Hk ( f )’s are all equal. The only nonzero DFT coefficient in this case is c 0 ( f ) — the average of the Hk ( f )’s.
- 2. Adjacent-Channel Timing Offset
It will prove to be useful to use the adjacent-channel timing offset , δk , between two cyclically neighboring channels; thus, we have
δ k Δ t k+1 Δ t k          for  k=0,  ...  ,M2.
The δk ’s do not retain the common timing offset in sampling instances, but the timing mismatch information is still preserved. To retrieve the Δ tk ’s from the δk ’s (to drive the tunable sampling clock), we need an additional constraint. Two reasonable choices are as follows:
Δ t k = n=0 k1 δ n    assuming Δ t 0 =0
Δ t k = n=0 k1 δ n 1 M m=0 M2 ( M1m ) δ m    assuming  k=0 M1 Δ t k =0.
Timing conversion in this work follows (7), the one based on the center of offset, since it evenly distributes the timing error across M tunable clocks, requiring a smaller delay tuning range than (6), where the first channel is chosen as reference.
III. Stationarity-Based Blind Method
In general, looking at the TI-ADC output alone does not uniquely determine converter mismatches, since there are many input-mismatch combinations that will yield the same TI-ADC output. It is necessary, therefore, to constrain the permissible input signal to a proper subset of all Nyquist-bandlimited signals. We want these constraints to be weak enough to admit as large a class of signals as possible, yet strong enough to enable blind mismatch detection. How to constrain the input signal is an important question and one worthy of attention, since this determines the practicality and complexity of the blind algorithm. One of the previously proposed approaches is to constrain the input signal bandwidth in the frequency domain, inspired by the frequency-translation action of mismatch in (4). The unoccupied portion of the input spectrum plays the role of aliasing (thus mismatch) detector. Some recent techniques require only a small fraction of extra bandwidth, minimizing input spectrum loss. It remains unclear, however, as to the effect of out-of-band interferers (for example, signal harmonics, adjacent-channel residual power, thermal noise, and so on) in real-life signals, especially when the anti-aliasing filter is not perfect. In this paper, we take a time-domain view of TI-ADC artifacts; that is, we focus on a periodic modulation of the input by converter mismatches.
Under the assumption of WSS TI-ADC input, output autocorrelation plays the role of a mismatch indicator. Thus, the proposed technique can be best introduced by first examining the properties of TI-ADC output autocorrelation.
Since the input x ( t ) is a WSS, its autocorrelation is shift-independent and thus depends only on the time lag between two samples. Consequently, the autocorrelation of x ( t ) is
R x ( τ )=E[ x( t+τ )x( t ) ]     for all t.
With nonzero residual mismatches, the TI-ADC output z [ n ] no longer satisfies WSS properties, and its autocorrelation is shift-dependent. Specifically, we focus on the subset of TI-ADC output autocorrelations with zero and unit lag, Rz,k [0] and Rz,k [1], referenced to each channel as follows:
R z,k [ 0 ]=E[ z k 2 [ n ] ],
R z,k [ 1 ]={ E[ z k [ n ] z k+1 [ n ] ]          for  k=0,  ...  ,M2, E[ z k [ n ] z 0 [ n+1 ] ]       for  k=M1.
These can be rewritten in terms of the input autocorrelation, Rx ( τ ) , by using (2), as follows:
R z,k [ 0 ]= ( G k * G ˜ k ) 2 R x ( 0 ),
R z,k [ 1 ]= G k * G ( k+1 ) mod M * G ˜ k G ˜ ( k+1 ) mod M ( R x ( T s )+( δ k * δ ˜ k ) d R x dτ | τ= T s ).
The first-order approximation in (11) is valid if δ * k and
δ ˜ k
are much smaller than Ts . Note that Rz,k [0] is a function of only gain mismatch, but Rz,k [1] depends on both gain and timing errors. The following important observations are made from (10)–(11):
  • 1)No residual mismatch condition. If there are no residual mismatches, thenRz,0[m] =Rz,1[m] … =Rz,M−1[m] form= 0. In other words, theRz,k[m]’s are equalized across channels, and shift-independence has been attained.
  • 2)Equalized autocorrelation condition. If all theRz,k[m]’s are equalized (m= 0, 1), then gain and timing estimates are equal to their intrinsic namesakes up to a common scale factor and time delay, respectively; that is,G˜0/G0*=G˜1/G1*= ... G˜M−1/GM−1*, and Δt˜0−Δt0*=Δt˜1−Δt1*= ... Δt˜M−1−ΔtM−1*.
  • 3) Therefore, if we disregard common time delay and scaling, 1) and 2) establish the following: the attainment of equalization ofRz,k[0]’s andRz,k[1]’s is necessary and sufficient for perfect mismatch correction.
The equivalence condition in 3) is a key result that the proposed method is based on. Next, the actual adaptive algorithm used to achieve the output correlation equalization will be discussed.
IV. Adaptive Estimation Algorithm
For adaptation to time-varying mismatch errors, the algorithm will be iterative in nature. Starting from an initial estimate, the calibration loop will gradually refine parameter estimates until the zero- and unit-lag output correlation coefficients are all equalized. The equivalence result in Section III then guarantees that the estimates are equal to true parameters up to a common scale and time delay.
- 1. Empirical Output Autocorrelation
The calibration cycle starts with the procurement of the output autocorrelation coefficients ( Rz,k [0]’s and Rz,k [1]). Given a batch of subconverter outputs, yk [ n ], the gain-corrected stream zk [ n ] is calculated from (1). The output correlation coefficients are then empirically obtained by
R z,k ( i ) [ 0 ]= 1 N n=0 N1 z k ( i ) [ n ] 2   ,
R z,k ( i ) [ 1 ]={ 1 N n=0 N2 z k ( i ) [ n ] z k+1 ( i ) [ n ]           for     k=0,  ...  ,  M2, 1 N n=0 N2 z k ( i ) [ n ] z 0 ( i ) [ n+1 ]      for     k=M1,
where the superscript indicates that the calculation is based on the i th iteration batch data. For simpler notation, this superscript will be dropped afterwards unless necessary for clarity.
- 2. Equalization Reference
Given the two sets of empirical correlation coefficients from (12)–(13), we want to choose an appropriate equalization reference for each set against which empirical coefficients are compared. One of the sensible choices is the following average coefficients:
R z,ref [ 0 ]= 1 M k=0 M1 G ˜ k 2 R z,k [ 0 ]  ,
R z,ref [ 1 ]= 1 M k=0 M1 R z,k [ 1 ]  ,
where Rz,ref [0] and Rz,ref [1] is the equalization reference for the Rz,k [0]’s and Rz,k [1]’s, respectively. The reference coefficients in (14)–(15) are basically an average across channels, except for the weighting by the gain estimate in (14). In fact, Rz,ref [0] is equal to the average correlation observed before digital gain correction; that is, at yk [ n ]. However, using (14) is more efficient than separately observing yk [ n ] and calculating its autocorrelation coefficients. After plugging (10)–(11) into (14)–(15) and neglecting the common scaling factor, Rz,ref [0] and Rz,ref [1] can be rewritten as a function of the TI-ADC input autocorrelation; that is,
R z,ref [ 0 ]= R x ( 0 ),
R z,ref [ 1 ]= R x ( T s )       if no residual gain error.
- 3. Parameter Recursion
By subtracting (16)–(17) from (10)–(11), we have
R z,k [ 0 ] R z,ref [ 0 ]= R x ( 0 )( ( G k * G ˜ k ) 2 1 ),
R z,k [ 1 ] R z,ref [ 1 ]= d R x dτ | τ= T s ( δ k * δ ˜ k ).
Here, Rx (0) is the input signal power; thus, it is strictly positive for nonzero signals. It can be shown that the derivative dRx / at τ = Ts is strictly negative for signals in the odd-ordered Nyquist zones and strictly positive in the even-ordered zones. Therefore, the sign of the left-hand side of (18) and (19) uniquely determines if the current gain or timing estimate is greater or smaller than its intrinsic namesake. This leads to the following parameter update rule:
G ˜ k ( i+1 ) = G ˜ k ( i ) + β g ( R z,k ( i ) [ 0 ] R z,ref ( i ) [ 0 ] ),
δ ˜ k ( i+1 ) = δ ˜ k ( i ) β t ( R z,k ( i ) [ 1 ] R z,ref ( i ) [ 1 ] ),
  • βg> 0
  • βt> 0 ifx(t) is in odd-order Nyquist zones,
  • βt< 0 otherwise.
G ˜ k ( i )
δ ˜ k ( i )
are the estimate of k th channel gain and adjacent-channel timing error at the i th iteration, respectively. Once the Rz,k [0]’s and Rz,k [1]’s are all equalized, then the driving term in the parenthesis in (20) and (21) is zero; hence, convergence is achieved. The stability and speed of convergence is controlled by βg and βt , which will be referred to as convergence parameters. Note that the use of adjacent-channel timing parameters decouples the parameter updates for each channel.
The calculation of 2 M correlation coefficients is practically all that is necessary for a single parameter update. This is significantly more efficient than previous blind techniques. The recursion rule in (20) and (21) also yields faster convergence than finite-difference methods or general search algorithms, because parameter adjustment is made with a priori known direction of decreasing residual error.
- 4. Convergence Analysis
Let’s define γk (i) and εk (i) as residual gain and adjacent-channel timing error at the i th iteration, respectively. Then, we have
γ k ( i ) = G ˜ k ( i ) G k * ,
ε k ( i ) = δ ˜ k ( i ) δ k * .
Then, from (18)–(21), it can be shown that γk (i) and εk (i) each follow a uniquely different geometric series under a small-mismatch regime.
γ k ( i ) = [ 12 β g R x ( 0 ) ] i γ k ( 0 ) ,
ε k ( i ) = [ 1 β t | d R x dτ | τ= T s ] i ε k ( 0 ) .
The effect of residual gain error is neglected in (25), for simplicity. The magnitude of the geometric ratios in (24)–(25) should be less than one so as to ensure convergence. For monotonic convergence, which is usually preferred over an oscillatory one, βg and βt must lie in the following range:
0< β g < 1 2 R x ( 0 ) ,
0< β t < 1 | d R x dτ | τ= T s .
It is seen from the above that small values of βg and βt will, in general, guarantee monotonic convergence. Large values of βg and βt will speed up convergence as long as (26) and (27) are satisfied. This will, however, necessarily amplify noise from the driving term (the one in the parenthesis in (20) and (21)), making the parameter estimates also noisy. With βg and βt fixed, the gain and timing estimate will converge faster with higher input power and fast-changing input signal, respectively, since the residual errors in (24) and (25) will diminish faster.
- 5. Other Considerations
For simplicity, previous analysis assumed no gain mismatch when discussing the timing estimate convergence. Equation (11) suggests that nonzero residual gain error may bias timing error estimates. However, as long as the gain calibration loop is in action, the timing estimates will eventually converge to a true parameter (see Section V for experimental results).
There are several factors that will affect parameter estimation: ADC quantization noise, sampling clock random jitter, autocorrelation estimation error due to finite observation, finite-resolution sampling clock tuning, and so on. Using small values of βg and βt , or increasing the batch size N , in general, decreases the contribution of these noise sources. However, the sampling clock quantization error can only be reduced by increasing its resolution. The minimum tuning resolution should be commensurate with the target signal-to-noise ratio (SNR) or Spurious-Free Dynamic Range (SFDR) level. For example, 80 dB of target SNR requires 0.00003 Ts of timing control resolution (assuming the input sinusoid is at fs /2). If the delay control range is 0.03 Ts , for instance, then approximately 1,000 quantization levels (that is, 10-bit level) will be necessary. Other than the minimum resolution requirement, the proposed calibration scheme is tolerant to analog imperfections in clock tuning circuitry. For example, the precise tuning curve need not be known. It can even change between calibration cycles since the feedback action of the blind calibration will eventually track such variations.
V. Experimental Results
Figure 3 shows the M = 4 TI-ADC experimental setup. Four 14-bit 100 MSPS commercial ADC’s (AD6645 from Analog Devices, Inc.) are used to achieve 400 MSPS of overall sampling rate. The logic analyzer performs both data acquisition and digital signal processing. The four-phase sampling clock is derived from a single 100 MHz reference clock, followed by a voltage-controlled delay line. Each tunable line consists of a single 50 Ω T-section ( L-C-L ) with varactor diodes (MV104 from ON Semiconductor) for delay control. The delay line provides 0.2 Ts ( Ts = 2.5 ns) of delay tuning range across 0 V to 3 V of tuning voltage.
PPT Slide
Lager Image
Experimental setup for M = 4 400-MSPS TI-ADC with proposed blind adaptive calibration loop. Sampling clock is fine-tuned by a single L-C-L section of varactor-loaded delay line. 6 dB attenuators minimize impedance variation with varactor tuning, thus preventing undesirable cross-line tuning effects.
- 1. Narrowband Input Test
Sinusoids are used as a representative narrowband input signal. Convergence parameters are chosen to be one tenth of the stability limit given by (26)–(27) ( βg , βt ≈ 0.2). The batch size for all narrowband tests is N = 4,096.
First, a 171.567 MHz sinusoid is applied at the TI-ADC input, and the blind calibration loop is initiated. Figure 4 shows the measured convergence curve for gain and sampling time error estimate. It is seen that the measurement (solid lines) closely follows the prediction (dotted lines). The small discrepancy in timing convergence plot is due to the slight curvature in the delay tuning characteristic (the prediction assumes linear tuning). Figure 5 shows the TI-ADC output spectrum before and after 300 calibration iterations. We define the following performance-evaluation metric, disregarding mismatch-irrelevant spurs:
PPT Slide
Lager Image
Measured convergence plot of (a) gain and (b) timing error estimates with a 171.567 MHz sinusoidal input. Solid lines represent measurement. Dotted lines denote predicted curves by (22)–(25) with intrinsic gain [0.99865, 1.0038, 1.0005, 0.99695], and sampling timing errors [0.0224Ts, 0.0012Ts, 0.0058Ts, −0.0294Ts] characterized at 172.8 MHz by a training method [6].
PPT Slide
Lager Image
Measured TI-ADC output spectrum with a 171.567 MHz input (marked with ‘1’): (a) before and (b) after 300 iterations (N = 4,096). Gain and sampling time mismatch spurs are labeled with ‘X.’ The offset spurs are represented by ‘O,’ and the input signal harmonics up to the 7th order are also shown as numbers.
Mismatch-limited SFDR =  signal power max (mismatch spur power) .
Thus, mismatch-limited SFDR (ML-SFDR) is the maximum achievable SFDR in the absence of other spurious tones. The proposed calibration suppressed mismatch spurs by more than 40 dB, achieving a value of more than 80 dB of ML-SFDR ( Fig. 6 ). Note that the initial linear increase (in dB scale) comes from exponential parameter convergence. The equalization process of the output correlation coefficients is shown in Fig. 7 .
PPT Slide
Lager Image
Improvement of mismatch-limited SFDR during 300 iterations with a 171.567 MHz sinusoid input (N = 4,096).
PPT Slide
Lager Image
Measured equalization plot of (a) unit-lag, (b) zero-lag correlation coefficients, and (c) their cross-channel variance. Input signal is a sinusoid at 171.567 MHz.
Next, the frequency of the input sinusoid is swept across the first two Nyquist zones (10 MHz to 390 MHz), with the same batch size of N = 4,096 ( Fig. 8 ). Up to 70 dB to 90 dB of SFDR is obtained.
PPT Slide
Lager Image
Single-tone test result across the first two Nyquist zones (N = 4,096). The sign of the correlation derivative and convergence parameters is also shown (see (26), (27)). The input anti-aliasing filter in Fig. 3 (cutoff = 187 MHz) was removed for this test.
There are three frequencies in the first Nyquist zone where the stationary input assumption fails for M = 4: 50 MHz, 100 MHz, and 150 MHz. As the input sinusoid comes closer to one of these frequencies, a longer batch size is required to maintain calibration performance. Calibration performance versus batch size near to 150 MHz is investigated in Fig. 9 . No performance loss is observed down to 150.5 MHz. If the input is at 150.05 MHz, however, then the ML-SFDR drops down to 50 dB (still 20 dB of calibration gain, though), and at least N = 10 5 is necessary to regain 80 dB of the ML-SFDR level. In summary, the batch size N = 4,096 can meet 80 dB of single-tone ML-SFDR requirement over 98.5% of the entire Nyquist band. If a longer observation of up to N = 10 5 is allowed, then the spectral performance is maintained over 99.85% of the bandwidth.
PPT Slide
Lager Image
Single-tone calibration performance near 150 MHz input, which is one of the failure frequencies (50 MHz, 100 MHz, and 150 MHz). The closer the input frequency approaches to 150 MHz, the greater the required batch size. Dotted line is uncalibrated performance.
- 2. Wideband Input Test
An independent, identically distributed sequence is first generated by an arbitrary waveform generator and then filtered by a 10-tap finite impulse response (FIR) filter. Its occupied bandwidth is approximately 180 MHz. Unlike sinusoidal inputs, wideband input signals mostly overlap with their own aliasing products. For the purpose of alias identification, a deep notch is created at 80 MHz. Before calibration, there exist significant channel mismatches, and the notch at 80 MHz is partially filled with frequency-shifted input spectra, as given by (4) ( Fig. 10(a) ). After 1,200 iterations, the calibration suppresses any aliasing products, and the deep notch is restored, as seen in Fig. 10(b) . Parameter convergence plots are shown in Fig. 11 . Note that the batch size N = 131,072 is much longer than the narrowband case. This is because empirical correlation coefficients for wideband signals are noisier than narrowband ones, in general.
PPT Slide
Lager Image
Measured TI-ADC output spectrum with a dc to 180 MHz wideband input signal: (a) before and (b) after 1,200 iterations (N = 131,072). The original deep notch at 80 MHz is masked by aliasing products in (a), but restored after calibration in (b), implying that the mismatch is corrected.
PPT Slide
Lager Image
Measured convergence plot of (a) gain and (b) timing error estimates with the wideband input signal in Fig. 10 (N = 131,072).
- 3. Effects of Residual Gain Error on Timing Curve
As a final test, the timing calibration loop is first initiated with a 171.567 MHz sine input, with gain correction turned off ( Fig. 12 ). In Fig. 12(b) , timing estimates are seen to converge with bias in the presence of residual gain error. After 200 iterations, a gain correction loop begins ( Fig. 12(a) ), when timing parameters also start readjusting themselves. After another 100 iterations, both gain and timing estimates achieved convergence to a true parameter. Under a small-mismatch regime, in general, residual gain error does not significantly affect the convergence of timing parameters, as experimentally verified by Figs. 12(b) and (c) .
PPT Slide
Lager Image
Measured convergence plot with a 171.567 MHz sine input: (a) gain and (b) timing error estimate. Gain adaptation is intentionally delayed by 200 iterations to see its effect on timing error convergence.
VI. Conclusion
We have demonstrated a new adaptive blind technique for multi-channel TI-ADCs. The analog-domain correction of timing mismatches, combined with autocorrelation-based error detection, dramatically reduces hardware and computational complexity. Specifically, empirical calculation of 2 M output autocorrelation coefficients is practically enough for a single parameter update. After blind mismatch correction, 70 dB to 90 dB of SFDR was experimentally achieved across the first two Nyquist zones. Proof of parameter convergence is given under the WSS input assumption. There is no restriction in the number of TI-ADC channels or input signal bandwidth.
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea Government (MEST) (No. 2013-056424).
Corresponding Author
Munkyo Seo received his BS and MS degrees in electronics engineering from Seoul National University, Rep. of Korea, in 1994 and 1996, respectively. He received his PhD degree in electrical engineering from the University of California, Santa Barbara (UCSB), CA, USA, in 2007. From 1997 to 2002, he was an RF engineer with LG Electronics Inc., Anyang, Rep. of Korea, designing microwave subsystems for wireless communication. He was an assistant project scientist with UCSB from 2008 to 2009. In 2009, he joined Teledyne Scientific Company, Thousand Oaks, CA, USA, where he worked on the design of various millimeter-wave, terahertz, and high-speed mixed-signal circuits. Since 2013, he has been with Sungkyunkwan University, Suwon, Rep. of Korea, as an assistant professor.
Eunsoo Nam received his BS degree in physics from Kyungpook National University, Daegu, Rep. of Korea, in 1983 and his MS and PhD degrees in physics from the State University of New York, Buffalo, NY, USA, in 1992 and 1994, respectively. In 1994, he joined the Electronics and Telecommunications Research Institute (ETRI), Daejeon, Rep. of Korea, where his research involved compound semiconductor devices, including GaAs-and GaN-based high-frequency power device platforms, long wavelength InP semiconductor photonic devices, optoelectronic integrated circuits for microwave circuits, and SDN network devices. In 2006, he was a visiting scholar at the Division of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA. His current research interest is to combine radio and optics technologies to develop millimeter-wave photonic components and integrated functions for a number of applications, including broadband communications, radar, security, and instrumentation. He is now a senior vice president managing the materials and devices R&D department at ETRI.
Mark Rodwell received his PhD degree from Stanford University, CA, USA, in 1988. He holds the Doluca Family Endowed Chair in Electrical and Computer Engineering at the University of California, Santa Barbara (UCSB), CA, USA. He directs the UCSB node of the NSF Nanofabrication Infrastructure Network, and the SRC Nonclassical CMOS Research Center. His research group works to extend the operation of integrated circuits to the highest feasible frequencies. He received the 2010 IEEE Sarnoff Award for the development of InP-based bipolar IC technology for mm-wave and sub-mm-wave applications. His group’s work on GaAs Schottky-diode ICs for subpicosecond/mm-wave instrumentation was awarded the 1997 IEEE Microwave Prize.
Black W.C. , Hodges D.A. 1980 “Time Interleaved Converter Arrays,” IEEE J. Solid-State Circuits 15 (6) 1022 - 1029    DOI : 10.1109/JSSC.1980.1051512
Poulton K. 2003 “A 20 GS/s 8 b ADC with a 1 MB Memory in 0.18 μm CMOS,” IEEE Int. Solid-State Circuits Conf. San Francisco, CA, USA Feb. 13 318 - 496    DOI : 10.1109/ISSCC.2003.1234315
Petraglia A. , Mitra S.K. 1991 “Analysis of Mismatch Effects among A/D Converters in a Time-Interleaved Waveform Digitizer,” IEEE Trans. Instrum. Meas. 40 (5) 831 - 835    DOI : 10.1109/19.106306
Kurosawa N. “Explicit Analysis of Channel Mismatch Effects in Time-Interleaved ADC Systems,” IEEE Trans. Circuits Syst. I, Fundam. Theory Appl. 48 (3) 261 - 271    DOI : 10.1109/81.915383
Vogel C. , Johansson H. “Time-Interleaved Analog-to-Digital Converters: Status and Future Directions,” IEEE Int. Symp. Circuits Syst. Kos Island, Greece May 21–24, 2006 3386 - 3389    DOI : 10.1109/ISCAS.2006.1693352
Seo M. , Rodwell M. , Madhow U. 2005 “Comprehensive Digital Correction of Mismatch Errors for a 400-Msamples/s, 80-db SFDR Time-Interleaved Analog-to-Digital Converter,” IEEE Trans. Microw. Theory Techn. 53 (3) 1072 - 1082    DOI : 10.1109/TMTT.2005.843487
Jamal S.M. 2002 “A 10-b 120-Msample/s Time-Interleaved Analog-to-Digital Converter with Digital Background Calibration,” IEEE J. Solid-State Circuits 37 (12) 1618 - 1627    DOI : 10.1109/JSSC.2002.804327
Jamal S.M. 2004 “Calibration of Sample-Time Error in a Two-Channel Time-Interleaved Analog-to-Digital Converter,” IEEE Trans. Circuits Syst. I, Reg. Papers 51 (1) 130 - 139
Matsuno J. 2013 “All-Digital Background Calibration Technique for Time-Interleaved ADC Using Pseudo Aliasing Signal,” IEEE Trans. Circuits Syst. I, Reg. Papers 60 (5) 1113 - 1121    DOI : 10.1109/TCSI.2013.2249176
Divi V. , Wornell G. 2009 “Blind Calibration of Timing Skew in Time-Interleaved Analog-to-Digital Converters,” IEEE J. Sel. Topics Signal Process. 3 (3) 509 - 522    DOI : 10.1109/JSTSP.2009.2020269
Elbornsson J. , Gustafsson F. , Eklund J.-E. 2004 “Blind Adaptive Equalization of Mismatch Errors in a Time-Interleaved A/D Converter System,” IEEE Trans. Circuits Syst. I, Reg. Papers 51 (1) 151 - 158    DOI : 10.1109/TCSI.2003.821300
Saleem S. , Vogel C. 2011 “Adaptive Blind Background Calibration of Polynomial-Represented Frequency Response Mismatches in a Two-Channel Time-Interleaved ADC,” IEEE Trans. Circuits Syst. I, Reg. Papers 58 (6) 1300 - 1310    DOI : 10.1109/TCSI.2010.2094330
Law C. , Hurst P. , Lewis S. 2010 “A Four-Channel Time-Interleaved ADC with Digital Calibration of Interchannel Timing and Memory Errors,” IEEE J. Solid-State Circuits 45 (10) 2091 - 2103    DOI : 10.1109/JSSC.2010.2061630
Huang S. , Levy B.C. 2006 “Adaptive Blind Calibration of Timing Offset and Gain Mismatch for Two-Channel Time-Interleaved ADCs Converters,” IEEE Trans. Circuits Syst. I, Reg. Papers 53 (6) 1278 - 1288    DOI : 10.1109/TCSI.2006.875180
Huang S. , Levy B.C. 2007 “Blind Calibration of Timing Offsets for Four-Channel Time-Interleaved ADCs,” IEEE Trans. Circuits Syst. I, Reg. Papers 544 (4) 863 - 876    DOI : 10.1109/TCSI.2006.888770
M. Seo , M. Rodwell , U. Madhow 2005 “Blind Correction of Gain and Timing Mismatches for a Two-Channel Time-Interleaved Analog-to-Digital Converter,” Asilomar Conf. Signals, Syst. Comput. Pacific Grove, CA, USA 1121 - 1125
Seo M. , Rodwell M. , Madhow U. “Blind Correction of Gain and Timing Mismatches for a Two-Channel Time-Interleaved Analog-to-Digital Converter: Experimental Verification,” IEEE Int. Symp. Circuits Syst. Kos Island, Greece May 21–24, 2006 3394 - 3397    DOI : 10.1109/ISCAS.2006.1693354
Fu D. 1998 “A Digital Background Calibration Technique for Time-Interleaved Analog-to-Digital Converters,” IEEE J. Solid-State Circuits 33 (12) 1904 - 1911    DOI : 10.1109/4.735530
Dyer K.C. 1998 “An Analog Background Calibration Technique for Time-Interleaved Analog-to-Digital Converters,” IEEE J. Solid-State Circuits 33 (12) 1912 - 1919    DOI : 10.1109/4.735531
Jin H. , Lee E.K.F. 2000 “A Digital-Background Calibration Technique for Minimizing Timing-Error Effects in Time-Interleaved ADCs,” IEEE Trans. Circuits Syst. II, Analog Digit. Signal Process. 47 (7) 603 - 613    DOI : 10.1109/82.850419
Iroaga E. , Murmann B. , Nathawad L. “A Background Correction Technique for Timing Errors in Time-Interleaved Analog-to-Digital Converters,” IEEE Int. Symp. Circuits Syst. Kobe, Japan May 23–26, 2005 6 5557 - 5560    DOI : 10.1109/ISCAS.2005.1465896
Huang C.-C. , Wang C.-Y. , Wu J.-T. 2011 “A CMOS 6-Bit 16-GS/s Time-Interleaved ADC Using Digital Background Calibration Techniques,” IEEE J. Solid-State Circuits 46 (4) 848 - 858    DOI : 10.1109/JSSC.2011.2109511
Camarero D. 2008 “Mixed-Signal Clock-Skew Calibration Technique for Time-Interleaved ADCs,” IEEE Trans. Circuits Syst. I, Reg. Papers 55 (11) 3676 - 3687
Haftbaradaran A. , Martin K. 2008 “A Background Sample-Time Error Calibration Technique Using Random Data for Wide-Band High-Resolution Time-Interleaved ADCs,” IEEE Trans. Circuits Syst. II, Exp. Briefs 55 (3) 234 - 238    DOI : 10.1109/TCSII.2008.918970
El-Chammas M. , Murmann B. 2011 “A 12-GS/s 81-mW 5-Bit Time-Interleaved Flash ADC with Background Timing Skew Calibration,” IEEE J. Solid-State Circuits 46 (4) 838 - 847
Seo M. , Rodwell M. , Madhow U. “A Low Computation Adaptive Blind Correction for Time-Interleaved ADCs,” IEEE Int. Midwest Symp. Circuits Syst. San Juan, Puerto Rico Aug. 6–9, 2006 292 - 296