# A Self-Adaptive Frequency Response Compensation Method for a TIADC System

Authors: Peng Ye; Huiqing Pan; Hao Zeng; Min Li; Wuhuang Huang
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Time interleaving is one of the most efficient techniques employed in high-speed sampling systems. However, the frequency response mismatch among different channels will create distortion tones that degrade the system performance. In this paper, a self-adaptive frequency response mismatch compensation method is presented, where the design of compensation filter is optimized with a self-adapting strategy. This digital postprocessing technique realizes the compensation of frequency response effectively and also the increase of the digital bandwidth of the acquisition system. MATLAB-based simulation and an actual two-channel acquisition system test verify the effectiveness of the algorithm.


---

## Body

## 1. Introduction

High-speed digitizing of signals plays an important part in broadband communication and electronic measurement instruments, in particular, digital storage oscilloscopes (DSO). Increasing requirements are put forward for these systems. Because of the limits of existing electronic components and manufacturing technologies, the time-interleaved analog-to-digital converter (TIADC) parallel sampling technique is the most effective solution to high sampling rate and high sampling resolution needs.Unfortunately, the performance of this TIADC system is degraded by following factors:(1) the mismatching of amplitude among various channels and  (2) the different phase delay of sampling circuits. Generally, these mismatches are related to frequency and are also called frequency response mismatch. The frequency response mismatch will lead to additional spurious components in the spectrum and degrade the system’s signal-to-noise ratio (SNR) [1, 2].For data acquisition systems, as channel bandwidth becomes wider, both the acquisition bandwidth of the analog signals and the signal fidelity of the system increase. However, increasing the channel bandwidth is always difficult. The low integration due to manufacturing restrictions and the different input bandwidth of multiple electronic components easily lead to the deterioration of the system’s bandwidth [3]. On the other hand, the inconsistency of frequency response of each sampling channel in TIADC structure also exacerbates the unflatness of the system’s frequency response, which results in the decrease in measurement bandwidth of the whole system.To address this problem, many approaches are proposed. Among these approaches, digital postprocessing is the most popular technique, as analog matching techniques are either too imprecise or too expensive [4]. In [5], the channel of the TIADC system is fitted as a first-order RC circuit and it formulates the impact of bandwidth mismatching on SNR. Tsai, Satarzadeh et al. [6, 7] present corresponding calibration algorithms to such first-order modeled circuit. However, due to the diversity of the ADCs and other electronic components, the actual frequency response of acquisition channel is too complicated to be simply modeled as a first-order system. Based on WLS algorithm, Seo et al. [8] developed a bandwidth compensation filter. By employing LS algorithm Johansson and Löwenborg [9] suggested a technique for the filter design. However, whether in the WLS or LS, with the increase of test frequency points, the design of filter would be computationally burdensome, and the algorithm will become more complicated.In this paper, we develop a digital postprocessing algorithm to compensate the frequency response in the TIADC system. Based on frequency domain sampling, an optimized design technique is proposed by constructing the error function based on a self-adapting strategy. With the proposed algorithm, it is feasible to compensate frequency response and increase the bandwidth of the sampling system, and the fidelity of data acquisition system could also be ensured.In the remainder of this paper, the frequency response mismatch is analyzed, and the compensation filter is presented in Section2. The implementation details of the algorithm are given in Section 3. Several experiments using the proposed algorithm are reported in Section 4. Further discussions are summarized in Section 5.


## 2. Compensation of Frequency Response Mismatch

### 2.1. Frequency Response of Sampling System

Due to the limitation of front-end electronic elements, the sampling channel trends to be a low-pass filter. The passing capacity for different frequency components of the input signal is inconsistent; that is, the amplitude attenuation rate of high-frequency components is larger than that of low-frequency components. Meanwhile, the manufacturing processes also cause the inconsistent response of various parts of the acquisition channel, resulting in the unflatness of channel frequency response and the decrease of measurement bandwidth of the overall system. These bring the new challenge to the signal reconstruction of the TIADC system.Figure1 depicts a TIADC sampling system with M channels. It consists of M parallel ADCs, which is called “channel” and works at the same sampling rate of fs/M driven by a master sampling clock. After reorganization, the M ADCs work equivalently as a single ADC operating at an M times higher sampling rate. The impulse response of signal conditioning components of each acquisition channel is defined as si(t). Unfortunately, due to the mismatching of si(t), the actual output of TIADC y[n] can be expressed as:
(1)y[n]=∑m=0M-1x(t)*sm(t)∑n=0N-1δ(t-(nM+m)Ts)=yd[n]+e[n],
where δ(t) is the dirac delta function, Ts is the sampling period, and * denotes the convolution operation. yd[n] represents the desired output, and e[n] is an error term introduced by frequency response mismatch among channels.Figure 1
Block diagram of a typical parallel TIADC system.e
[
n
] is a frequency-associated error introduced by different responses to different frequency components of input signals through the acquisition channels. As the phase-frequency response may not change till very high sampling frequency, we only consider the amplitude-frequency response mismatch in this paper. If it is possible to estimate the error e^[n] by certain construction function, the compensation of the frequency response mismatch can be realized.


### 2.2. Frequency Response Compensation

In order to compensate the errore[n] introduced by the frequency response mismatch of each channel, a compensation filter h[n] may be built for postprocessing the digital signals collected by the  TIADC. The fundamental principle is described as below.Equation (1) may be expressed as
(2)y[n]=∑n=0N-1∑m=0M-1x(nMTs+mTs)*s(nMTs+mTs)·δ(t-nMTs-mTs)
The Fourier transform of (1) can be written as
(3)Y(ω)=1MTs∑k=-∞∞∑m=0M-1Sm(ω-k2πMTs)X(ω-k2πMTs)×e-jk(2π/M)m.If the frequency response of various channels is consistent, the ideal TIADC output frequency response will be(4)Yd(ω)=1MTs∑k=-∞∞∑m=0M-1X(ω-k2πMTs)e-jk(2π/M)mThen it is possible to obtain the estimation of output’s frequency responseY^(ω) from the TIADC system after the frequency response mismatch is calibrated by compensation filter h[n]:
(5)Y^(ω)=H(ω)Y(ω)≈Yd(ω),
where  H(ω) refers to the frequency response of compensation filter. Hence, it can be obtained according to (3), (4), and (5); it can be rewritten as
(6)H(ω)=∑k=-∞∞1∑m=0M-1Sm(ω-k(2π/MTs)).The frequency response errore[n] may be compensated by constructing the calibration filter H(ω).In the TIADC system, we could make full use of parallel features, and then signal processing could be decomposed to single channel, and Figure2 can be converted into Figure 3.Figure 2
Fundamental principle of frequency response mismatch compensation.Figure 3
Fundamental principle of frequency response error compensation.Obviously, in order to design the compensation filter, in (6) the frequency response function Sm(ω) of corresponding channel should be determined. We noted that the response of each channel is related to not only the parameters of the actual devices but also the attenuation of channel. Meanwhile, the design of the digital compensation filter is related to real-time sampling rate of the channel. For the acquisition systems which has different sampling rate fs, the frequency response compensation filter needs to be redesigned. On the other hand, a different response mode will provide the system with different features. Therefore, a suitable reference response should be selected according to the testing requirements, such as the flat response, bessel response or gaussian response.The fundamental principle that was used to construct the compensation system is shown in Figure4. Based on frequency domain compensation techniques, the correction can be made to the amplitude frequency mismatch in the parallel acquisition system.Figure 4
Fundamental principle of frequency response compensation filter design.


### 2.3. Design of Frequency Response Compensation Filter

We adopt the adaptive method to design the compensation filter. Considering the advantage of parallel structure in the TIADC acquisition system, we design the compensation algorithm by parallel computation.We use the frequency responseHmd(ω) that meets the design objectives as the reference for the mth channel. The objective of the design is to minimize the mean square error (MSE) of the performance function
(7)J=E{em(ω)em(ω)T}=1N∑ω(Hm(ω)-Hmd(ω))2,
where  Hm(ω) refers to the frequency response of the filter to be designed.TheHmd(ω) can be calculated as
(8)Hmd(ω)=Hmr(ω)Sm(ω),
where Hmr(ω) refers to the reference frequency response mode.The minimax algorithm could be used for filter design that is, the design objective is achieved by the constrained conditions that the maximum mean squares error (MMSE) is minimal. Base on an iterative method, the filter that satisfies the minimize requirements of (7) could be obtained.For reducing of multiplying operations in the filtering process, the symmetric or antisymmetric characteristic is used as constraint for subfilters, that is, the implementation complex can be reduced approximately by a factor of two.In summary, the design criteria of the frequency response compensation filter could be expressed as;(9)Δa=maxω∈ωp||Hm(ω)|-|Hmd(ω)||≤δa,(10)hm(n)=±hm(N-1-n),
where δa refers to the amplitude error limits N refers to the length of each compensation filter. If both (9) and (10) are satisfied, the frequency response compensation filter may fulfill the compensation within the bandwidth ωp.


## 2.1. Frequency Response of Sampling System

Due to the limitation of front-end electronic elements, the sampling channel trends to be a low-pass filter. The passing capacity for different frequency components of the input signal is inconsistent; that is, the amplitude attenuation rate of high-frequency components is larger than that of low-frequency components. Meanwhile, the manufacturing processes also cause the inconsistent response of various parts of the acquisition channel, resulting in the unflatness of channel frequency response and the decrease of measurement bandwidth of the overall system. These bring the new challenge to the signal reconstruction of the TIADC system.Figure1 depicts a TIADC sampling system with M channels. It consists of M parallel ADCs, which is called “channel” and works at the same sampling rate of fs/M driven by a master sampling clock. After reorganization, the M ADCs work equivalently as a single ADC operating at an M times higher sampling rate. The impulse response of signal conditioning components of each acquisition channel is defined as si(t). Unfortunately, due to the mismatching of si(t), the actual output of TIADC y[n] can be expressed as:
(1)y[n]=∑m=0M-1x(t)*sm(t)∑n=0N-1δ(t-(nM+m)Ts)=yd[n]+e[n],
where δ(t) is the dirac delta function, Ts is the sampling period, and * denotes the convolution operation. yd[n] represents the desired output, and e[n] is an error term introduced by frequency response mismatch among channels.Figure 1
Block diagram of a typical parallel TIADC system.e
[
n
] is a frequency-associated error introduced by different responses to different frequency components of input signals through the acquisition channels. As the phase-frequency response may not change till very high sampling frequency, we only consider the amplitude-frequency response mismatch in this paper. If it is possible to estimate the error e^[n] by certain construction function, the compensation of the frequency response mismatch can be realized.


## 2.2. Frequency Response Compensation

In order to compensate the errore[n] introduced by the frequency response mismatch of each channel, a compensation filter h[n] may be built for postprocessing the digital signals collected by the  TIADC. The fundamental principle is described as below.Equation (1) may be expressed as
(2)y[n]=∑n=0N-1∑m=0M-1x(nMTs+mTs)*s(nMTs+mTs)·δ(t-nMTs-mTs)
The Fourier transform of (1) can be written as
(3)Y(ω)=1MTs∑k=-∞∞∑m=0M-1Sm(ω-k2πMTs)X(ω-k2πMTs)×e-jk(2π/M)m.If the frequency response of various channels is consistent, the ideal TIADC output frequency response will be(4)Yd(ω)=1MTs∑k=-∞∞∑m=0M-1X(ω-k2πMTs)e-jk(2π/M)mThen it is possible to obtain the estimation of output’s frequency responseY^(ω) from the TIADC system after the frequency response mismatch is calibrated by compensation filter h[n]:
(5)Y^(ω)=H(ω)Y(ω)≈Yd(ω),
where  H(ω) refers to the frequency response of compensation filter. Hence, it can be obtained according to (3), (4), and (5); it can be rewritten as
(6)H(ω)=∑k=-∞∞1∑m=0M-1Sm(ω-k(2π/MTs)).The frequency response errore[n] may be compensated by constructing the calibration filter H(ω).In the TIADC system, we could make full use of parallel features, and then signal processing could be decomposed to single channel, and Figure2 can be converted into Figure 3.Figure 2
Fundamental principle of frequency response mismatch compensation.Figure 3
Fundamental principle of frequency response error compensation.Obviously, in order to design the compensation filter, in (6) the frequency response function Sm(ω) of corresponding channel should be determined. We noted that the response of each channel is related to not only the parameters of the actual devices but also the attenuation of channel. Meanwhile, the design of the digital compensation filter is related to real-time sampling rate of the channel. For the acquisition systems which has different sampling rate fs, the frequency response compensation filter needs to be redesigned. On the other hand, a different response mode will provide the system with different features. Therefore, a suitable reference response should be selected according to the testing requirements, such as the flat response, bessel response or gaussian response.The fundamental principle that was used to construct the compensation system is shown in Figure4. Based on frequency domain compensation techniques, the correction can be made to the amplitude frequency mismatch in the parallel acquisition system.Figure 4
Fundamental principle of frequency response compensation filter design.


## 2.3. Design of Frequency Response Compensation Filter

We adopt the adaptive method to design the compensation filter. Considering the advantage of parallel structure in the TIADC acquisition system, we design the compensation algorithm by parallel computation.We use the frequency responseHmd(ω) that meets the design objectives as the reference for the mth channel. The objective of the design is to minimize the mean square error (MSE) of the performance function
(7)J=E{em(ω)em(ω)T}=1N∑ω(Hm(ω)-Hmd(ω))2,
where  Hm(ω) refers to the frequency response of the filter to be designed.TheHmd(ω) can be calculated as
(8)Hmd(ω)=Hmr(ω)Sm(ω),
where Hmr(ω) refers to the reference frequency response mode.The minimax algorithm could be used for filter design that is, the design objective is achieved by the constrained conditions that the maximum mean squares error (MMSE) is minimal. Base on an iterative method, the filter that satisfies the minimize requirements of (7) could be obtained.For reducing of multiplying operations in the filtering process, the symmetric or antisymmetric characteristic is used as constraint for subfilters, that is, the implementation complex can be reduced approximately by a factor of two.In summary, the design criteria of the frequency response compensation filter could be expressed as;(9)Δa=maxω∈ωp||Hm(ω)|-|Hmd(ω)||≤δa,(10)hm(n)=±hm(N-1-n),
where δa refers to the amplitude error limits N refers to the length of each compensation filter. If both (9) and (10) are satisfied, the frequency response compensation filter may fulfill the compensation within the bandwidth ωp.


## 3. Implementation

Considering the computational capacity of the existing computers, we adopt an optimization technique to design a compensation filter for each channel. The purpose is to minimize the design deviation under the condition to ensure the accuracy of channel frequency response correction with less resource consumption.(i)
Step 1: obtain the frequency response of channelsSm(ω), m=0,1,…,M-1 and select the reference frequency response mode Hmr(ω).(ii)
Step 2: The optimal design of compensation filter, which meets the constraints of (9) and (10), could be designed by combination with the frequency sampling filter design method and Parks-McClellan optimization algorithm [10, 11].(iii)
Step 2(a): a proper filter should be firstly searched for as initial filter used for the iterative algorithm. As the amplitude-frequency response and the phase-frequency response have been known, the initial filterc0(n) could be simply designed. The amplitude-frequency response of the filter within the given pass band |ω|<ωp will be |Hmd(ω)| and within the stop band will be zero. And the phase response is linear a corroding to constraint (10). With N0 frequency response sampling data, it is very quick to obtain the initial filter coefficient c0(n),  n=0,…,N0-1 by the frequency sampling design method and the IFFT transformation.(iv)
Step 2(b): compute actual frequency response of initial filterc0(n). According to the Parks-McClellan optimization algorithm, if the actual frequency response of filter c0(n) with order N0 does not satisfy the constraint condition (9) and (10), then the filter order N1=N0+1 will be updated. The minimax algorithm is used to design the new filter c1(n) with order N1, until the constrained conditions are satisfied. And then, the N-order filter is the frequency response compensation filter hm(n) for the mth channel.(v)
Step 3: repeat Steps 2(a) and 2(b), until the frequency response compensation filtershm(n) for all channels (m=0,1,…,M-1) are designed. Then, the design of frequency response compensation filter is completed.In the design process, in order to utilize the filter design criteria simply, the pass band|ω|<ωp is divided into I equivalent portions at equal spaces, ωi∈[0,ωp],i=0, 1,…,I-1. Therefore, the constraint of (9) may be expressed as
(11)Δa=max0≤i≤I-1||Hm(ωi)|-|Hmd(ωi)||≤δa.


## 4. Experimental Results

In this section, we performed several simulations results to demonstrate the feasibility and effectiveness of the proposed algorithm.Experiment
(
1
)
. The feasibility of proposed calibration algorithm has been evaluated in a two-channel TIADC system with equivalent sampling rate at 1GSPS supported by two 8 Bit-500MSPS ADCs (Ideal_8_Bit.adc) [12]. Assume that a certain bandwidth mismatch exists among channels. Based on assumption, the frequency response of the mth channel can be defined as [4]
(12)Hm(ω)=11+jωτm
Equation (12) contains both magnitude and phase contributions, where τm is the sampling time constant of mth channel, and τ1 = 0.8 fs, τ2=0.9 fs. Figure 5 shows the frequency response of each channel in the TIADC system.Figure 5
Frequency response of two channels.A sine wave with frequencyf0=10MHz is used as a test signal. The frequency spectra of the original output sequence before and after compensation are depicted in Figures 6(a) and 6(b), respectively.Frequency spectrum of 10 MHz sine wave captured by two-channel TIADC before and after compensation.
(a)
(b)The experimental results indicate that the error spectra caused by frequency response mismatches are with very weak power after compensation. Now the calculated value of SNR is approximately 41.95 dB leading to an improvement of 9.6 dB, the ENOB (effective number of bits) is increased by 6.67 Bit–5.09 Bit, and the SFDR (spurious free dynamic range) is enhanced from 32.83 dB to 54.29 dB. The additional spurious components in the spectrum have been reduced efficiently. The rest of the additional spurious components may be caused by phase-frequency response mismatch that is not calibrated in this paper. This performance of this TIADC system verifies the effectiveness of the proposed algorithm.Experiment
(
2
)
. The frequency response calibration is already used in an actual underdesigned digital storage oscilloscope (DSO) (shown in Figure 7). This system is a double-channel acquisition system with a whole sampling frequency of fs=6GSPS and designed bandwidth of fc=1GHz.Figure 7
Under designed DSO(double-channel acquisition system) test platform.The amplitude-frequency response of each channel is characterized by a sinusoid sweeping test. Figure8 shows the experimental result.Figure 8
Amplitude-frequency response of actual double-channel DSO before compensation.In Figure8, the “Ideal” curve gives the reference amplitude-frequency response of the Gaussian response mode, curve with “·” and curve with “*” represent the amplitude-frequency response of two actual sampling channels respectively. Clearly, the gains of each channel with respect to different frequency components of input signals are changing, and it results in the unflatness of the actual frequency response. Meanwhile, the insufficiency of designed bandwidth also degrades the actual system’s performance.According to the proposed design procedure, the amplitude-frequency response of compensation filtersh0(n) and h1(n) are shown in Figure 9.Figure 9
Amplitude-frequency response of compensation filters.After applying the compensation, the amplitude-frequency response of the whole system is also evaluated using a sinusoid sweeping test, and the result is shown in Figure10.Figure 10
Amplitude-frequency response of system after compensation.Obviously, it can be seen from Figure10 that the bandwidth of the system after compensation is fc=1 GHz (ωc=0.33π), and, in the pass band, the flatness of the amplitude-frequency response also meets the requirements.


## 5. Discussions and Conclusions

In the wideband TIADC acquisition system, the frequency response mismatch at the high-frequency part is not avoidable. The high fidelity of signals sampling could be ensured by employing the digital postprocessing technique. With the proposed algorithm, it is possible to improve the measurement bandwidth and compensate the frequency response mismatch. The algorithm is realized by an adaptive optimal filter that could meet the design requirements. By the digital signal postprocessing technique, frequency response can be compensated online. And it breaks the limitation of the analog calibration and improves the system’s measurement bandwidth so as to satisfy the testing requirements of wideband signals.In this paper, we only consider the amplitude-frequency response mismatch of the TIADC system. As the phase-frequency response will change at a higher sampling frequency, our future work may focus on both amplitude-frequency response and phase-frequency response mismatches for high sampling system.



---

*Source: 613025-2013-04-11.xml*