用于估计时间交错a/d转换器系统的时间误差的方法和设备的制作方法

专利名称：用于估计时间交错a/d转换器系统的时间误差的方法和设备的制作方法
技术领域：
现代的信号处理应用需要A/D转换器具有越来越高的采样率。增加采样率的一种方法是使用多个并行的、时间交错的A/D转换器。该解决方法的一个问题在于由于并行结构而引入时间失配误差。本发明涉及一种用于估计时间交错A/D转换器或ADC系统的时间误差的方法和设备。与其他的方法和设备相比，该方法和设备在后台运行并且不需要任何特殊的校准信号。所述方法和设备产生的估计也将渐进无偏。
发明的技术背景诸如无线电基站或VDSL调制解调器那样的许多数字信号处理应用需要具有非常高的采样率和非常高的精度的A/D转换器。为了获得足够高的采样率，能够使用时间交错的M个A/D转换器或ADC的阵列。参照

图1，每一个ADC应以所期望的采样率[7]的第1/M进行工作。在图1中，示出了并行的ADC系统，其由多个单个或独立的ADC 1(也称为ADC单元或ADC通道)所组成，这些ADC在周期处理中转换连续采样的值，每一单元中的转换与其他单元中的转换并行执行。每一单元中的转换处理在用于连续采样的模拟值的连续时间上启动。模拟信号u(t)被输入到线路2上的设备。采样时钟信号被输入到线路3上的ADC单元的控制端子并且通过连接到时钟输入线路上的延迟元件4来连续地延迟。被提供给控制端子的信号控制相应ADC单元中转换程序的启动。来自ADC单元的M个输出信号y1至yM被连接到MUX 5上。MUX将信号y1至yM放在一起以形成一个输出信号y。ADC单元的延迟元件和输入电路引起时间偏移误差。
由于交错结构而引入了三种失配误差时间误差或静态抖动不同A/D转换器的时钟延迟时间是不相等的。这意味着信号将被周期性地、但不均匀地采样。
幅度偏移误差在不同A/D转换器中地电平可能稍微不同。这意味着每一A/D转换器中存在恒定的幅度偏移。
增益误差对于不同的A/D转换器来说，从模拟输入到数字输出的增益可能是不同的。
以上所列出的误差被假设是静态的，以致于从一个周期到下一周期相同A/D转换器中的误差是相同的。由于例如热噪声和量化，所以也存在随机误差，该随机误差从一个采样到下一采样是不同的。这些误差与A/D转换器的并行结构没有关系并且由于它们的随机特性，它们是不能够被估计的。然而，对于研究估计算法的稳健性并且对于计算估计精度的下限来说，所述随机误差是重要的。
A/D转换器中的主要的随机误差是量化如果已知输入信号，那么这是确定性的误差。但是对于大多数信号来说，其能够被视为与输入信号和均匀分布[11]无关联的附加的白噪声。
随机抖动由于时钟信号中的噪声，在采样情况下存在随机误差。这些误差能够被视为在采样情况下的高斯白噪声。
发明概述在一种用于计算并行A/D转换器中的单元的时间偏移值的方法中，输入信号基本上是不知道的。所以，在没有任何校准信号的情况下运行算法。时间偏移估计适用于频带被限制于尼奎斯特频率的信号。本发明涉及一种用于并行A/D转换器系统的时间误差估计和信号重构的方法和设备。这是对[3，2，5，4]中所介绍并且所评价的估计方法的改进。在此处所介绍的估计方法中，只要输入信号是频带受限的，就渐进地消除偏差估计误差。这意味着仅取决于数据数量，估计精度就能够是任意地好。在例如[1]和[8]中也提供了用于估计时间误差的方法，但是这些方法需要已知的校准信号。A/D转换器的校准是耗费时间的并且昂贵的。所以，如果能够自动地估计并且在偏移时补偿ADC中的误差，那么能够节省很多成本。
附图简述图1是受控于同一主时钟的M个并行ADC的框图。
图2是具有根据本发明的时间误差估计和校正的时间交错ADC系统的框图。
图3是示出在时间误差估计和补偿之前和之后的正弦信号的频谱的图。
图4是示出估计精度如何随着数据数量的增加而提高的图。
优选实施例的详细描述没有时间误差的额定采样间隔被表示为Ts。M表示时间交错阵列中A/D转换器的数量。时间误差参数被表示为Δti，i＝0，...，M-1。这些误差的估计被表示为ti，i＝0，...，M-1，并且真误差被表示为△ti0，i＝，...，M-1。下列符号被用于有关以下内容的信号u(t)是模拟输入信号，u[k]表示在无时间误差的情况下采样的人工信号，ui[k]表示u[k]和yi[k]的M个子序列是在具有时间误差的情况下采样的M个A/D转换器的输出子序列。
yi[k]=u((kM+i)Ts+Δti0)---(1)]]>y[m]是所有A/D转换器的多路复用输出信号。
z(Δi)[m]表示利用误差参数Δti重构的输出信号y[m]。
zi(△ti)是z(Δi)[m]的子序列。假设u(t)频带被限制于奈奎斯特频率。
在随后的部分中，建立一些稍后将使用的定义。[准稳态信号[9]]如果m=limN→∞1N∫0NE(u(t))dt]]>Ru(τ)=limN→∞1N∫0NE(u(t)u(t+τ))dt---(3)]]>存在，其中在可能的随机部分上取期望值，那么信号u(t)是准稳态的。稳态随机处理是准稳态的，其中m和气Ru(τ)分别表示平均值和协方差函数。[准稳态频谱[9]假设u(t)是准稳态的。于是s(t)的功率谱被定义为Φu(ω)=∫-∞∞Ru(τ)e-jωτdτ---(4)]]> 假设对于函数g(·，·，)来说存在g‾yi1,yi2,=limN→∞1NΣi=1Ng(yi1[t],yi2[t],...)]]>i1，i2，...＝0，...，M-1 (5)那么如果g‾i1,i2,=g‾(i1+l)modM,(i2+l)modM,)]]>l∈{...，-1，0，1，...} (6)则y相对于g是模M准稳态的。模M准稳态特性保证在时间交错阵列中的所有A/D转换器中信号具有相同的统计特性。
补偿后的输出信号的平均值被定义为m‾zi(Δti)=limN→∞1NΣk=lNE{zi(Δti)[k]}---(7)]]>补偿后的输出信号的均方值被定义为(σ‾zi2)(Δti)=limN→∞1NΣk=1NE{(zi(Δti)[k])2}---(8)]]>第i个A/D转换器与第j个A/D转换器的补偿后的输出之间的均方差被定义为 其中 表示整数部分。
如果输入信号的频带被限制于奈奎斯特频率 ，并且时间误差参数是已知的，那么根据不规则的采样能够理想地重构输入信号。根据[10]，随后的结果将是1.一般信号重构首先为了符号的简明性而引入符号αi=-M-12+i+Δti.]]>为Hi(ω，t)而求解等式系统
Σi=0M-1ejαiωHi(ω,t)=1]]>Σi=0M-1ejαi(ω+2πMTs)Hi(ω,t)=ej2πMTst]]>Σi=0M-1ejαi(ω+(M-1)2πMTs)Hi(ω,t)=ej(M-1)2πMTst---(10)]]>然后可以如下在任意时间情况下计算输入信号u(t)=Σk=-∞∞Σi=0M-1yi[k]hi(t-kMTs)]]>wherehi(t)=MTs2π∫-π/Ts-π/Ts+2π/(MTs)Hi(ω,t)ejωtdω---(11)]]>然而，该信号的重构在计算上是费力的，如果仅在额定的采样情况下需要该信号，那么其能够被简化。
2.采样信号的重构假设在时间情况t＝(kM+l)Tsl＝0，...，M-1，k＝...，-1，0，1，...(12)下重构信号。于是(10)的右手边不依赖于k。此外，所述右手边能够被因式分解为一个取决于ω的对角矩阵和一个不依赖于ω的矩阵AE(ω)H(i)(ω)＝BiA=1···1ejα02πMTs···ejαM-12πMTs···ejα0(M-1)2πMTs···ejαM-1(M-1)2πMTs---(14)]]>E(ω)=ejα0ω0···00ejα1ω···0···00···ejαM-1ω---(15)]]>Bl=[1ej2πl/M···ej2π(M-1)l/M]T---(16)]]>
由于仅E(ω)取决于ω并且去除了时间相依性，所以能够容易地计算系数hi(l)(k)=hi((kM+l)Ts)]]>h(l)(k)=MTs2π∫γ1γ2E-1(ω)ejω(kM+l)TsdωA-1Bl]]>γ1＝-π/Ts，γ2＝-π/Ts+2π/(MTs) (17)从这里开始假设M是偶数，M是奇数则产生相似的计算。计算子序列h(l)k的TDFT产生H(l)(ejωMTs)=MTsΣk=-∞∞h(l)(k)]]>=(MTs)22π∫-π/Ts-π/Ts+2π/(MTs)E-1(γ)ejγlTsΣk=-∞∞ejγkMTse-jωKMTsdγA-1Bl]]>=MTsΣr=-∞∞∫-π/Ts-π/Ts+2π/(MTs)E-1(γ)ejγlTsδ(γ-ω+r2πMTs)dγA-1Bl]]>=MTsE-1(ω-πTs)ejωlTs(-1)lA-1Bl]]>0≤ω<2πMTs---(18)]]>然后可以如下计算子序列 其中YT(ejωMTs)=[Y0(ejωMTs)···YM-1(ejωMTs)]]]>然后可以根据其子序列[6]计算经时间误差补偿的信号的TDFT、即Z(Δt0)(ejωTs):]]>Z(Δt0)(ejωTs)=Σl=0M-1Zl(Δi0)(ej(ωMTs2π))e-jlωTs---(20)]]>在接下来的部分中，介绍一种用于估计时间交错A/D转换器系统的时间误差的算法。当使用A/D转换器时所述估计算法在后台运行，也就是说不需要特殊的校准信号。将假设在第一A/D转换器中误差将为零，Δt00=0,]]>也就是所有的其他误差都与第一ADC的时间误差有关系。这不是限制，因为在所有采样之间仅需要相同的间隔，并且绝对采样时间并不重要。参数估计适合于随机梯度算法对时间误差估计的要求1.该算法除了应针对整个ADC系统将输入信号频带限制于尼奎斯特频率之外不需要有关输入信号的知识。
2.真误差参数必须小于采样间隔的一半，即Δti0<Ts/2.]]>这基本上仅意味着必须按正确的顺序进行采样。
初始化1.从每一A/D转换器收集一批N个采样，yi[k]，i＝0，...，M-1。
2.初始化随机梯度算法的步长μ。
3.初始化参数估计，Δ^ti(0)=0,i=0,...,M-1.]]>4.初始化迭代计数器n＝0。
5.定义在估计损失函数中应使用多少时间滞后L。
6.定义损失函数权重参数β。
调整1.计算利用Δ^ti(n)=0,i=0,...,M-1]]>校正的输出信号ziΔ^ti(n)=IDFT{ZiΔ^ti(n)---(21)]]>其中根据部分3由y[k]计算 2.计算相邻ADC之间的均方差 i＝0，...，M-1，l＝0，...，L (22)3.计算每一个补偿后的ADC输出的均方(σ‾zi2)(Δti(n))=limN→∞1NΣk=1NE{(zi(Δti(n))[k])2}]]>i＝0，...，M-1(23)4.计算损失函数V((i))V(Δ^(i))=Σl=0LΣi,j(Rzi,zi-1(Δt(0))[l]-Rzi,zi-1(Δt(n))[l])2+βΣi,j((σ‾zi2)(Δti(n))-(σ‾zj2)(Δtj(n)))2---(24)]]>5.计算损失函数的梯度V((n))。
6.沿负梯度方向更新参数估计
(n+1)＝(n)-μV((n)) (25)7.计算V((n+1))，如果V((n+1))＞V((n))，那么更新μ＝μ/2，并且返回到点6。
8.收集新的一批数据并且从点2继续估计更新。
图2中描述了具有时间误差估计和校正的时间交错ADC系统。所述系统包括多个单个或独立的ADC 6(也称为ADC单元或ADC通道)ADC1至ADCM，这些ADC在周期处理中转换连续采样的值，每一单元中的转换与其他单元中的转换并行执行。每一单元中的转换处理在用于被连续采样的模拟值的连续时间上启动。每一ADC单元以MTs的采样间隔而工作。整个系统的采样率是Ts。模拟信号u(t)被输入到线路7上的设备。采样时钟被输入到线路8上的ADC单元的控制端子。在ADC之间的时钟信号上使用延迟为Ts的延迟元件9。被提供给控制端子的信号控制相应ADC单元中转换程序的启动。输入信号u被连接到所有的ADC上。必须针对整个ADC系统将输入信号的频带限制于尼奎斯特频率。M个输出信号y1至yM被连接到校正设备10上，该校正设备10利用所估计的时间误差delta t、即Δt来校正信号，并且产生M个信号z1至zM。在估计算法设备11中使用信号z1至zM，其中时间误差delta t、即Δt被估计。然后将所估计的时间误差delta t、即Δt反向连接到校正设备上。信号z1至zM也被连接到MUX 12上，其中所述信号被多路复用在一起以形成一个信号z。根据时间误差来校正信号z并且以采样间隔Ts来采样所述信号z。
已经利用具有三种不同输入信号正弦信号、多正弦信号以及频带受限的白噪声的仿真来评价了该估计算法。也已经利用输入上的加性白噪声和时钟上的随机抖动来测试了该算法。所述仿真示出该估计算法适用于宽范围的信号并且其对噪声和随机抖动来说也是稳健的。在图3中，单个正弦信号被用作为ADC系统的输入信号u(t)。上面的绘图示出了在时间误差估计和校正之前ADC系统的输出信号y的频谱。下面的绘图示出了在时间误差估计和补偿之后校正设备的输出信号z。在图4中，针对不同的数据数量示出了时间误差估计精度。在此，输入信号u(t)是单个正弦曲线。时间误差估计精度在该图中也与利用那些数据数量而能够获得的最好的Cramer-Rao界限进行了比较。
在许多应用中需要高速A/D转换器。提高采样速度的一种方法是使用多个并行的、时间交错的A/D转换器。于是每一A/D转换器能够以较低速度运行，同时整个系统实现高采样率。时间交错结构的问题在于引入了失配误差。在本说明书中已经介绍了用于时间误差的估计和补偿的方法和设备。当使用A/D转换器并且A/D转换器不需要特殊的校准信号时，能够在线进行所述估计。估计算法也已利用仿真被验证。所述仿真已说明该估计方法很好地适用于许多不同种类的信号，唯一的要求是信号频带被限制于尼奎斯特频率。
参考文献[1]J.J.Corcoran.Timing and amplitude error estimation fortimeinterleaved analog-to-digital converters.US专利号5294926，1992年10月。
J.Elbornsson.Equalization of Distortion in A/DConverters.Lic.thesis 883，Department of ElectricalEngineering，Linkoping University，Linkoping，Sweden，2001年4月。
J.Elbornsson and J.-E.Eklund.Blind estimation of timingerrors in interleaved AD converters.In Proc.ICASSP 2001，卷6，第3913-3916页，IEEE，2001年。
J.Elbornsson，K.Folkesson，and J.-E.Eklund.Measurementverification of estimation method for time errors in atime-interleaved A/D converter system.in Proc.ISCAS2002.IEEE，2002年。
J.Elbornsson，F.Gustafsson，and J.-E.Eklund.Amplitudeand gain error inuence on time error estimation algorithm fortime interleaved A/D converter system.in Proc.ICASSP2002.IEEE，2002年。
F.Gustafsson，L.Ljung，and M.Millnert.DigitalSignalbehandling.Studentlitteratur，2001年。(瑞典语)[7]Y-C Jenq.Digital spectra of nonuniformly sampledsignalsA robust sampling time offset estimation algorithm forultra high-speed waveform digitizers using interleaving.IEEETransactions on Instrumentation and Measurement.39(1)71{75，1990年2月。
H.Jin and E.K.Lee.A digital-background calibrationtechnique for minimizing timing-error effects in time-interleaved ADC’s.IEEE Transactions on Cicuits and Systems，47(7)603{613，2000年7月。
L.Ljung.System Identification，Theory for the user.Prentice-Hall，第2版，1999年。
A.Papoulis.Signal Analysis.McGraw-Hill，1977年。
B.Widrow，I.Kollar，and M.-C.Liu.Statistical theory ofquantization.IEEE Transactions on Instrumentation andMeasurement，45(2)353{361，1996年4月。
权利要求
1.一种用于估计时间交错A/D转换器系统的时间误差的方法，其特征在于计算利用参数估计校正的输出信号，计算相邻ADC之间的均方差，计算每一补偿后的ADC输出的均方，计算损失函数，计算损失函数的梯度，以及更新参数估计。
2.一种用于估计时间交错A/D转换器系统的时间误差的设备，其特征在于输出信号y1至yM被连接到校正设备(10)上，所述校正设备(10)被配置用于利用所估计的时间误差来校正所述信号并且被配置用于产生M个信号z1至zM，这些信号被提供用于在估计算法设备中用来估计时间误差，以及所估计的时间误差被提供用于连接到所述校正设备上。
3.依据权利要求2所述的设备，其特征在于所述估计算法设备被布置用于利用根据新数据的估计来自适应地更新时间误差估计。
全文摘要
本发明涉及一种用于估计时间交错A/D转换器系统的时间误差的设备。输出信号(y1，y2，…，yM)被连接到校正设备(5)上，该校正设备被配置用于利用所估计的时间误差来校正所述信号并且被配置用于产生M个信号(z1，z2，…，zM)。这些信号被提供用于在估计算法设备中用来估计时间误差并且所估计的时间误差被提供用于连接到校正设备(5)上。
文档编号H03M1/12GK1739241SQ200480002215
公开日2006年2月22日 申请日期2004年1月15日 优先权日2003年1月15日
发明者J·埃尔博恩松 申请人:因芬尼昂技术股份公司