专利名称：一种基于拉格朗日插值方法的时间交替模拟数字转换(tiadc)系统时间失配实时补偿算法的制作方法
技术领域：
本发明涉及一种基于拉格朗日插值方法的时间交替模拟数字转换CriADC,
Time-Interleaved Analog-to-digital Converter)系统时间失配实时补偿算法，属于高速高精度模拟 数字转换技术领域。
背景技术：
高速高精度模数转换是雷达、医疗仪器、通信系统等现代电子系统的重要部件，而传统 的ADC (如流水线结构ADC)由于器件工艺的限制很难同时具有高精度和高速度两种性能。 一种新的解决方案是使用并行交替ADC结构(TIADC)。这种并行交替ADC由M个独立的并 行子通道(子ADC)构成，各个子通道以/VM的采样频率对相同的输入信号进行分时交替采 样,然后M个子通道的输出重组成一组数字输出信号。这样，整个并行交替ADC的采样频率 为乂。理想情况下，M个子通道为相同线性电路，并且具有相同的电路特性。例如，所有子 通道具有相同的增益，相同的采样周期(1《)等。但由于实际制造工艺的限制，使得各个子通 道之间不可避免地产生以下几种失配增益失配(Gain mismatch)、直流偏置失配(Offset mismatch)和时间失配(Timing mismatch)。这些通道失配使得TIADC的动态无伪谱范围(SFDR) 大大降低，严重影响TIADC的性能。其中，增益失配和直流偏置失配较易解决，只需在各个 子通道输出信号通路上增加至多一个加法器和乘法器。而时间失配的补偿则要困难许多。
国际上针对时间失配的研究很多，提出了许多时间失配的补偿技术。如Jin和Lee提出 了基于Neville算法的插值补偿算法(H. Jin and E. Lee， "A digital-background calibration technique for minimizing timing-error effects in time-interleaved ADC，s，" IEEE Trans. Cicuits Syst" vol.47,no. 7, pp. 603~613， Jul.2000)， Prendergast， Levy, Hurst提出了一种新的采用多率滤 波器组的方法实现TIADC输出序列重构，从而实现时间失配的补偿(R. S. Prendergast, B. C. Levy, and P. J. Hurst, "Reconstruction of Band-Limited Periodic Nonuniformly Sampled Signals Through Multirate Filter Banks," IEEE Trans. Circuits and Systems I， vol. CAS-I-51, no. 8, pp.1612-1622, Aug 2004)。研究表明，这些已发表的补偿算法能有效地实现TIADC系统的时间 失配补偿，但都具有计算复杂度大，难以实现硬件实时补偿。因此，研究一种新的补偿性能好、 计算复杂度低，易于硬件实现的TIADC系统时间失配实时补偿算法具有重要现实意义。
本发明的目的是提出一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系统 时间失配实时补偿算法，该补偿算法具有高速、高精度、计算复杂度低的特点，易于硬件实现，能对TIADC系统的输出进行实时时间失配补偿。
一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系统时间失配实时补偿算 法其思想是，采用TIADC系统实际非均匀采样输出值，利用拉格朗日插值的数学方法实现 TIADC系统理想的均匀采样输出值的估计，并采用多通道补偿滤波器实现实时补偿。具体原理 表述如下。对于《+l个TIADC系统实际非均匀采样输出数据点(x。j。),...其中a 为釆样时间，》为采样时间^时的采样输出值，x,^c//#)，在时刻f的W点拉格朗日插值估计 由下式计算<formula>formula see original document page 6</formula>
其中，W为拉格朗日插值阶数，f为欲求理想均匀采样点的时间坐标值，/2,仍为对应于时 刻Z的第Z个拉格朗日插值系数，/A^)为对应于/时刻的拉格朗日插值(输出估计值)。利用
TIADC系统各个子通道的周期均匀采样特性，可以证明公式(2)求得的^=[/^,...,/^]对于 同一通道所有理想均匀采样点的拉格朗日插值系数为一常数矢量。因此，基于拉格朗日插值 方法的TIADC时间失配插值补偿算法可采用多通道补偿滤波器实现。对于TIADC系统的第 附个子通道，其多通道补偿滤波器系数即为^ 6'=0，/,2,...，Aa其值由TIADC系统的采样周 期r,、各通道时间失配参数zl/m (m=0,7,...M-/)、以及所用拉格朗日插值多项式的阶数iV共 同决定。因而，对于处于一定工作环境下的特定TIADC系统，所有M个多通道补偿滤波器 系数均为N点常数矢量，依据公式(2)，采用离线计算方法进行确定。
一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系统时间失配实时补偿算 法，其步骤是
(a) 确定采用的拉格朗日插值多项式的阶数W,其值等于所用插值点数减l。这里假定为 偶数，即采用的插值点数为奇数。
(b) 求前^+^M+iV/2个TIADC系统实际非均匀采样点的时间坐标私栌O丄...， (7+^)M+iW2-7 。其中，/=ceil((iV+l)/2/W)， M为TIADC的通道数。
根据以下公式求a:
<formula>formula see original document page 6</formula>
其中，7;是TIADC的采样周期，水h。dM为通道(A:M^MJ的时间失配参数，该参数可以离线精确测得，这里假定为已知。
(c) 求理想均匀采样点的时间坐标^/=/*似/*^+/，...,(/+/力1/-/。其中，/=ceil((iV+l)/2A0。 根据以下公式求t-
,' = 0,/'2,... (4)
(d) 求用于补偿各个通道时间失配的多通道补偿滤波器的W+7点拉格朗日插值系数向量 I w ,..,,M-/。
求好m的公式如下
好m = [Vo 、2…附=。,7' 2,…'(5)
M/+_/ +—W/2 z一x
、,=ft 一丄 /=0'/,2，...，W w="7，2，.."M-7 (6)
乂'藩+m-JV/2+Z
/ = ce,7(^~(7) 、:2A/ 、
其中，m为通道标志。理想均匀采样时刻表示为^=/*r,， /=0乂2，...，而/"z,)为在 理想均匀采样时刻^的w阶拉格朗日插值输出估计值。
(e) 利用多通道补偿滤波器实现理想均匀采样时刻的拉格朗日插值输出估计值，即求 細。
计算/WO按照以下公式进行 /w(^+ ,) = [; :T"m P=ft人2"." / =0'7,2"."M-7 "=pM+m^V/2 (8)
-W/2少"-W/2+l…凡…少w+W/2-l少"+W/2] (9) 公式(8)由公式(1)推导得来。其中，/pA/^为第W通道所属第p个输出时的理想
均匀采样时间点。[义:T为多通道补偿滤波器的输入向量，即由W+7个TIADC系统 实际输出采样值形成的向量。依次求出的/;^i)， /= iV/2+l， JV/2+2，...，即为补偿后的 TIADC的输出。因为求/MW时要用到其前面的个输出值0,-w/2, ...,力-》，所
以真正的有效补偿输出是从第^/2+/点开始的。
值得一提的是，步骤(a)、 (b)、 (c)、 (d)都不涉及TIADC系统的采样输出值，只与TIADC 系统的参数K、 ^U和所采用的拉格朗日多项式的阶数iV有关。因此，步骤(a)、 (b)、 (c)、 (d)可以在实现TIADC系统时间失配实时补偿前执行，即离线执行。步骤(e)在TIADC系统进行 时间失配实时补偿时执行。因该算法实现一个TIADC系统输出补偿值只需执行个乘法 和W个加法，无迭代运算，计算复杂度为0(JV+7),故而步骤(e)可以实现TIADC系统的时间 失配实时补偿。
本发明的有益效果在于本发明所述的算法能对TIADC系统的时间失配进行实时补偿， 算法由多通道滤波器实现，计算复杂度低，易于硬件实现。对于4通道12bit 320MHz釆样率 的TIADC硬件系统，采用XILINXXUPVirtexIIProFPGA实现了本发明提出的算法(N-7), 并实现了 TIADC系统的时间失配实时补偿。测试结果表明，本发明提出的算法能够提升该 TIADC硬件系统的SFDR (动态无伪谱范围)达30dB。
图1为时间交替模拟数字转换(TIADC)系统的示意图。 图2为TIADC系统的非均匀采样示意图。
图3为本发明所述算法的结构框图。图中的多通道滤波器有多个输入和一个输出，输入 分别来自TIADC系统的M个通道ADC的输出采样值。M个多通道滤波器的输出通过一个多 路选择器最后组成一个输出。
图4为实现四通道TIADC系统的时间失配3点插值补偿时，力(W、力(b)、力(G)、力(")的 计算结构图。其中
具体实施例方式
下面将结合附图和实施例对本发明作进一步的说明。
实施例1:以对四通道(^/=2)采样周期为r,通道失配参数为水,,,(/ =0乂2，3)的TIADC进行 3点插值补偿(A^2)为例。
(a) 确定采用的拉格朗日阶数^=2;
(b) 求前9个TIADC系统实际非均匀采样输出的时间坐标，此时 /=ceil((iV+l)/2M)=ceil((2+l)/2*4)-l。
根据公式(3)可求得A的值如下<formula>formula see original document page 9</formula>(c) 求4个理想均匀采样时间坐标/= 3,《5，6。 根据公式(4)求得f,值如下<formula>formula see original document page 9</formula>
(d) 求拉格朗日插值系数向量双 ，附=0,7,2,3。
根据公式(5)和(6)求得四个通道各自的拉格朗日插值系数向量(多通道滤波器 系数)如下<formula>formula see original document page 9</formula> (9)
附<formula>formula see original document page 9</formula> (10)
(e) 对TIADC系统的输出进行实时补偿。 补偿按照以下公式进行
<formula>formula see original document page 9</formula>(ii)<formula>formula see original document page 9</formula> (12)
如对TIADC系统的第1、 2、 3、 4个输出点进行补偿时，计算式如下(见图4):
<formula>formula see original document page 9</formula>
依次求出的力(f,)， /=1， 2，...，即为TIADC系统时间失配补偿后的输出。
权利要求
1.一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系统时间失配实时补偿算法，其特征在于，其思想是，采用TIADC系统实际非均匀采样输出值，利用拉格朗日插值的数学方法实现TIADC系统理想的均匀采样输出值的估计，并采用多通道补偿滤波器实现实时补偿。具体原理表述如下对于K+1个TIADC系统实际非均匀采样输出数据点(x0，y0)，...，(xk，yk)，其中xk为采样时间，yk为采样时间xk时的采样输出值，当xi≠j(i≠j)，在时刻t的N点拉格朗日插值估计由下式计算<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>f</mi> <mi>N</mi></msub><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>0</mn> </mrow> <mi>N</mi></munderover><msub> <mi>h</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><msub> <mi>y</mi> <mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo></mrow> </mrow>]]></math></maths><maths id="math0002" num="0002" ><math><![CDATA[ <mrow><msub> <mi>h</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><mo>=</mo><munderover> <munder><mi>Π</mi><mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn></mrow> </munder> <mrow><mi>j</mi><mo>≠</mo><mi>i</mi> </mrow> <mi>N</mi></munderover><mfrac> <mrow><mi>t</mi><mo>-</mo><msub> <mi>x</mi> <mi>j</mi></msub> </mrow> <mrow><msub> <mi>x</mi> <mi>i</mi></msub><mo>-</mo><msub> <mi>x</mi> <mi>j</mi></msub> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo></mrow> </mrow>]]></math></maths>其中，N为拉格朗日插值阶数，t为欲求理想均匀采样点的时间坐标值，hi(t)为对应于时刻t的第i个拉格朗日插值系数，fN(t)为对应于t时刻的拉格朗日插值(输出估计值)。利用TIADC系统各个子通道的周期均匀采样特性，可以证明公式(2)求得的H＝[h0，...，hN]对于同一通道所有理想均匀采样点的拉格朗日插值系数为一常数矢量。因此，基于拉格朗日插值方法的TIADC时间失配插值补偿算法可采用多通道补偿滤波器实现。对于TIADC系统的第m个子通道，其多通道补偿滤波器系数即为hi(i＝0，1，2，...，N)，其值由TIADC系统的采样周期Ts、各通道时间失配参数Δtm(m＝0，1，...M-1)、以及所用拉格朗日插值多项式的阶数N共同决定。因而，对于处于一定工作环境下的特定TIADC系统，所有M个多通道补偿滤波器系数均为N点常数矢量，依据公式(2)，采用离线计算方法进行确定。
2. 根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系 统时间失配实时补偿算法，其特征在于，其步骤是(a) 确定采用的拉格朗日插值多项式的阶数W，其值等于所用插值点数减l。这里假定为 偶数，即采用的插值点数为奇数。(b) 求前^+/> 7^+^/2个TIADC系统实际非均匀采样点的时间坐标a， hO丄...， (7+ /JM+iV/2-7 。其中，/=ceil((AT+l)/2M), M为TIADC的通道数。根据以下公式求;^<formula>formula see original document page 2</formula>(3)其中，r,是TIADC的采样周期，Ah。dM为通道(AM^姊的时间失配参数，该参数可以离线精确测得，这里假定为已知。(c) 求理想均匀采样点的时间坐标^1'=/*械/*^+7，...,(/+^^-7。其中，/=ceil((W+l)/2JW)。 根据以下公式求^<formula>formula see original document page 3</formula> (4)(d) 求用于补偿各个通道时间失配的多通道补偿滤波器的点拉格朗日插值系数向量求及m的公式如下<formula>formula see original document page 3</formula>其中，m为通道标志。理想均匀采样时刻表示为/产/*7;， /=0，1，2,...,而/w(A)为在 理想均匀采样时刻A的iV阶拉格朗日插值输出估计值。(e) 利用多通道补偿滤波器实现理想均匀采样时刻的拉格朗日插值输出估计值，即求 細。计算//^,)按照以下公式进行<formula>formula see original document page 3</formula>公式(8)由公式(1)推导得来。其中，Z;^+附为第W通道所属第p个输出时的理想均匀采样时间点。[y f为多通道补偿滤波器的输入向量，即由JV+J个TIADC系统实 际输出采样值形成的向量。
3.根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系 统时间失配实时补偿算法，其特征在于，其步骤(a)在TIADC正常工作前离线执行。
4. 根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换CTIADC)系统时间失配实时补偿算法，其特征在于，其步骤(b)在TIADC正常工作前离线执行。
5. 根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)时间失配实时补偿算法，其特征在于，其步骤(c)在TIADC正常工作前离线执行。
6. 根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)时间失配实时补偿算法，其特征在于，其步骤(d)在TIADC正常工作前离线执行。
7. 根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)时间失配实时补偿算法，其特征在于，其步骤(e)在TIADC正常工作时实时执行。
8. 根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)时间失配实时补偿算法，其特征在于，步骤(d)中的系数向量及m是一次性求出，其值仅由TIADC系统的采样周期r,、各通道时间失配参数zl^、所用拉格朗日插值多项式的阶数W决定。当r,、 Am、 W不变时，Wm保持不变，不需另求。
9. 根据权利要求1所述的一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系统时间失配实时补偿算法，其特征在于，因为求/"/,)时要用到其前面的个输出值O,.w","-w/2+7,...j^)，所以有效时间失配补偿输出是从第iV/2+7点开始。即有效补偿后的TIADC系统的输出为/^,)， …。
全文摘要
本发明涉及一种基于拉格朗日插值方法的时间交替模拟数字转换(TIADC)系统时间失配实时补偿算法。所述算法主要利用TIADC系统各个子通道分别对输入信号进行周期采样的特点推导得出，算法执行前需要事先获取TIADC系统各通道的时间失配信息。所述算法采用易于硬件实现的多通道补偿数字滤波器实现TIADC系统的时间失配补偿。多通道滤波器的系数由TIADC系统采样频率、各通道时间失配信息和所用拉格朗日插值的点数共同决定。因此，在某一工作环境下对于特定的TIADC系统来说，多通道滤波器的系数是常数，只需要离线计算一次。本发明所述算法计算复杂度为O(N+1)，硬件实现时，占用资源少，执行速度快，能够对TIADC系统的输出进行实时补偿，非常适用于TIADC系统芯片设计。
文档编号H03M1/12GK101656538SQ20091010948
公开日2010年2月24日 申请日期2009年8月21日 优先权日2009年8月21日
发明者张尚良, 邹月娴 申请人:北京大学深圳研究生院