CN105720983A - 用于时间交织模数转换***的误差估计方法和装置 - Google Patents
用于时间交织模数转换***的误差估计方法和装置 Download PDFInfo
- Publication number
- CN105720983A CN105720983A CN201610044047.7A CN201610044047A CN105720983A CN 105720983 A CN105720983 A CN 105720983A CN 201610044047 A CN201610044047 A CN 201610044047A CN 105720983 A CN105720983 A CN 105720983A
- Authority
- CN
- China
- Prior art keywords
- error
- function
- frequency
- time
- fourier transformation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000012360 testing method Methods 0.000 claims abstract description 26
- 230000009466 transformation Effects 0.000 claims description 73
- 238000005070 sampling Methods 0.000 claims description 70
- 238000012545 processing Methods 0.000 claims description 34
- 230000001143 conditioned effect Effects 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 abstract description 10
- 230000006870 function Effects 0.000 description 91
- 238000012937 correction Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000009966 trimming Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/10—Calibration or testing
- H03M1/1071—Measuring or testing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Analogue/Digital Conversion (AREA)
Abstract
本发明涉及一种用于时间交织模数转换***的误差估计方法和装置,通过在多通道TIADC***中输入测试信号获得误差信号,然后对误差信号进行傅里叶变换得到其频谱幅值数据,在时间误差存在的情况下对各个子通道的增益误差与偏置误差进行估计,然后对增益误差与偏置误差进行校正,消除增益误差和偏置误差,最后对时间误差进行估计。本发明可以在线同时估计三种失配误差;当外界的环境改变时,无需调整滤波器***或者重新设计硬件电路,只需重新输入测试信号即可得到三种失配误差;对输入信号的频率没有限制;估计的三种失配误差精度高。
Description
技术领域
本发明涉及信号采样技术领域,特别是涉及一种用于时间交织模数转换***的误差估计方法和装置。
背景技术
随着集成电路技术的不断发展,数字化技术的推广,对模数转换器件ADC的采样速率以及采样精度的要求越来越高,不仅要求数据采集***有高的采样率,还要有高的采样精度。在实际的运用中,对实时采样速率以及采样精度有极高的依赖性。然而ADC的最大采样速率受限于它的分辨率,分辨率与采样速率之间是一对矛盾体,高采样速率要求较短的转换时间,而高分辨率则要求较长的转换时间。根据目前的IC设计(IntegratedCircuitDesign,集成电路设计)工艺,要实现更高速的采样速率,我们需要探索一种基于新结构和新方法的ADC。一种实现超高速采样的重要方法就是利用多通道时间交织模数转换***TIADC(Time-interleavedAnalog-to-DigitalConverter)。
如图1所示,TIADC是利用M片有着相同采样率fs的单个ADC,采用并行的结构,每片ADC以相隔1/(M*fs)的时间间隔进行采样,以达到采样率为M*fs(总采样率f=M*fs)的效果。理论上,这种对于M通道的并行交替采样的ADC结构能够使得整个***采样率达到单个DAC的M倍。但是由于制造工艺本身固有的缺点,不可能使得每一片ADC完全一模一样,所以必然会使得各个通道ADC之间存在失配误差,从而严重降低了整个ADC***的信噪比。
主要的失配误差为时间误差(Time-skewerror)、增益误差(Gainerror)和偏置误差(offseterror)。偏置误差相当于一个加性的误差,它会引起频谱上单个通道采样率整数倍的频点上出现伪谱线。增益误差相当于一个乘性的误差,各个通道之间的增益误差可以引起采样信号的幅度调制。时间误差相对于采样周期来说是一个很小的量,但是它会导致信号存在较大的杂散分量,它会在频谱上导致原始信号频谱在单个通道采样率的整数倍频点上的一个频谱偏移。
国内外早期基于多通道时间交织ADC***的失配修正一般是利用对前端电路的修调,通过精心布局的线路来减少失配误差的影响。这种方法的缺点就是当随着时间的推移,温度的变化,电器元件的老化会使得电路的修正效果失效。为了克服这种前端修正的方法,可以利用后端处理的方法。
目前基于多通道时间交织ADC***的失配误差及其数字后端处理的修正算法是未来发展的关键。国内外的研究机构针对时间误差、增益误差和偏置误差的修正算法有很多,但是传统的修正算法都是时间误差、增益误差和偏置误差分开单独校正的,即使可以实现多个误差的同时校正,也需要输入多种频率的测试信号,用数据库的方式记录误差,或者需对***输入一个固定频率的测试信号,这类方法不仅需要额外的存储空间而且对输入的测试信号要求严格,所以研究一种同时对三种误差同时估计,对输入测试信号限制少的用于多通道时间交织ADC***的失配修正方法非常有实际意义。
发明内容
基于此,有必要针对上述问题,提供一种用于时间交织模数转换***的误差估计方法和装置,能够实现适配误差的同时估计并且对测试信号没有限制。
为了达到上述目的,本发明采取的技术方案如下:
一种用于时间交织模数转换***的误差估计方法,包括步骤:
在多通道时间交织模数转换***中输入测试信号,获取输出的误差信号;
将所述误差信号进行傅里叶变换,获得傅里叶变换后增益误差出现的第一频点和偏置误差出现的第二频点;
根据所述误差信号的傅里叶变换在第一频点的值,获得增益误差的傅里叶变换;根据所述误差信号的傅里叶变换在第二频点的值,获得偏置误差的傅里叶变换;
根据增益误差的傅里叶变换和偏置误差的傅里叶变换,获得各个子通道的增益误差和偏置误差;
根据原采样函数、增益误差和偏置误差得到各个子通道的新的采样函数;
将各个子通道的新的采样函数进行傅里叶变换,得到第一函数;
将各个子通道的新的采样函数合并后进行傅里叶变换,得到第二函数;
根据预设条件对所述第一函数和所述第二函数进行处理;
根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道;
根据第三频点和第三频点对应的通道,获得各个子通道的时间误差。
一种用于时间交织模数转换***的误差估计装置,包括:
误差信号获取模块,用于在多通道时间交织模数转换***中输入测试信号,获取输出的误差信号;
频点确定模块,用于将所述误差信号进行傅里叶变换,获得傅里叶变换后增益误差出现的第一频点和偏置误差出现的第二频点;
傅里叶变换模块,用于根据所述误差信号的傅里叶变换在第一频点的值,获得增益误差的傅里叶变换;根据所述误差信号的傅里叶变换在第二频点的值,获得偏置误差的傅里叶变换;
增益和偏置误差确定模块,用于根据增益误差的傅里叶变换和偏置误差的傅里叶变换,获得各个子通道的增益误差和偏置误差;
采样函数更新模块,用于根据原采样函数、增益误差和偏置误差得到各个子通道的新的采样函数;
第一函数确定模块,用于将各个子通道的新的采样函数进行傅里叶变换,得到第一函数;
第二函数确定模块,用于将各个子通道的新的采样函数合并后进行傅里叶变换,得到第二函数;
函数处理模块,用于根据预设条件对所述第一函数和所述第二函数进行处理;
频点和通道确定模块,用于根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道;
时间误差确定模块,用于根据第三频点和第三频点对应的通道,获得各个子通道的时间误差。
本发明用于时间交织模数转换***的误差估计方法和装置,通过在多通道TIADC***中输入测试信号获得误差信号,然后对误差信号进行傅里叶变换得到其频谱幅值数据,在时间误差存在的情况下对各个子通道的增益误差与偏置误差进行估计,然后对增益误差与偏置误差进行校正,消除增益误差和偏置误差,最后对时间误差进行估计。本发明可以在线同时估计三种失配误差;当外界的环境改变时,无需调整滤波器***或者重新设计硬件电路,只需重新输入测试信号即可得到三种失配误差;对输入信号的频率没有限制;估计的三种失配误差精度高。
附图说明
图1为多通道时间交织模数转换***实施例一的结构示意图;
图2为本发明用于时间交织模数转换***的误差估计方法实施例的流程示意图;
图3为多通道时间交织模数转换***实施例二的结构示意图;
图4为单个通道中三种失配误差示意图;
图5为输入信号与存在三种失配误差的输出信号(时域)的示意图;
图6为输入信号与存在三种失配误差的输出信号(频域)的示意图;
图7为本发明校正增益误差与偏置误差后的信号与误差信号(时域)的示意图;
图8为本发明输入信号与校正增益误差与偏置误差的校正后信号(频域)的示意图;
图9为本发明RMSE(dB)与时间误差的示意图;
图10为本发明RMSE(dB)与增益误差的示意图;
图11为本发明RMSE(dB)与偏置误差的示意图;
图12为本发明用于时间交织模数转换***的误差估计装置实施例的结构示意图;
图13为本发明函数处理模块实施例的结构示意图;
图14为本发明频点和通道确定模块实施例的结构示意图;
图15为本发明时间误差确定模块实施例的结构示意图。
具体实施方式
为更进一步阐述本发明所采取的技术手段及取得的效果,下面结合附图及较佳实施例,对本发明的技术方案,进行清楚和完整的描述。
如图1所示,一种用于时间交织模数转换***的误差估计方法,包括步骤:
S101、在多通道时间交织模数转换***中输入测试信号,获取输出的误差信号;
S102、将所述误差信号进行傅里叶变换,获得傅里叶变换后增益误差出现的第一频点和偏置误差出现的第二频点;
S103、根据所述误差信号的傅里叶变换在第一频点的值,获得增益误差的傅里叶变换;根据所述误差信号的傅里叶变换在第二频点的值,获得偏置误差的傅里叶变换;
S104、根据增益误差的傅里叶变换和偏置误差的傅里叶变换,获得各个子通道的增益误差和偏置误差;
S105、根据原采样函数、增益误差和偏置误差得到各个子通道的新的采样函数;
S106、将各个子通道的新的采样函数进行傅里叶变换,得到第一函数;
S107、将各个子通道的新的采样函数合并后进行傅里叶变换,得到第二函数;
S108、根据预设条件对所述第一函数和所述第二函数进行处理;
S109、根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道;
S110、根据第三频点和第三频点对应的通道,获得各个子通道的时间误差。
在步骤S101中,输入的测试信号可以为余弦信号Acos(2πfint+θ),其中A为余弦信号的幅值,fin为频率,θ为初始相位,t为时间。输入的测试信号还可以为其它信号,例如正弦信号等。
当测试信号为余弦信号时,输出的误差信号为:
其中,ei(t)为误差信号,gi表示第i通道的增益误差,ai表示第i通道的偏置误差,Δti表示第i通道的时间误差,M代表总的通道数,i=(0,1.....M-1),δ()代表狄拉克δ函数,Ts=1/fs,fs代表每一个通道即单片ADC的采样频率,n表示范围,为模数转换中数字信号的表达方式。
需要说明的是,上述误差信号仅仅为了对测试信号为余弦信号的情况进行示例,并不对具体的误差信号做出限定,当测试信号为其它信号时,误差信号做相应的改变。后续公式的示意是在测试信号为余弦信号的基础上得到的,并不对具体的形式加以限定,后面不再一一说明。
在步骤S102中,对误差信号ei(t)进行傅里叶变换,例如进行快速傅里叶变换FFT。需要说明的是,后续公式中仅对FFT进行了示意,但是并不对傅里叶变换的具体形式加以限定,后续不再一一说明。
在误差信号的傅里叶变换上选取增益误差和偏置误差出现的频点:
增益误差gi在频域上出现在第一频点 偏置误差ai在频域上出现在第二频点
在步骤S103中,利用误差函数ei(t)的傅里叶变换在第一频点 和第二频点的值,分别获得:
增益误差gi的傅里叶变换为:
偏置误差ai的傅里叶变换为:
在步骤S104中,根据以下公式估计各个子通道的增益误差和偏置误差:
增益误差
偏置误差
其中,IFFT为FFT的逆变换。
在步骤S105中,根据以下表达式得到各个子通道的新的采样函数:
其中,i=(0.....M-1),ai为偏置误差,gi为增益误差,为各个子通道的新的采样函数,yADCi为原采样函数。
将原来的采样函数更新为新的采样函数,就可以消除各个子通道的增益误差和偏置误差。
在步骤S106中,对新的采样函数进行傅里叶变换,例如快速傅里叶变换,则第一函数
在步骤S107中,将新的采样函数合并成函数再对函数进行快速傅里叶变换,则第二函数
在步骤S108中,在一个实施例中,根据预设条件对所述第一函数和所述第二函数进行处理的步骤可以包括:
将第一函数Yi中小于ε×MAX(Yi)的值设为0,其余值不变,其中ε为预设常数,例如ε可以设置为0.06等,Yi为第一函数中各通道的函数,i=(0.....M-1),M为总的通道数,MAX为最大值;
将第二函数Y中小于ε×MAX(Y)的值设为0,其余值不变。
在步骤S109中,在一个实施例中,所述测试信号为余弦信号,根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道的步骤可以包括:
获取Y0在(0,π)中使得幅值|Y0(ejω)|不为零的频点ωk1,其中Y0为处理后的第一函数中第一通道的函数,ω为角频率,e为自然指数,j为虚数单位;
从k=1,...,M中选取使不为零的k,将选取的k作为kp1,其中M为总的通道数,Y为处理后的第二函数;
获取Y0在(π,2π)中使得幅值|Y0(ejω)|不为零的频点ωk2;
从k=1,...,M中选取使不为零的k,将选取的k作为kp2。
寻找通道0中Y0(即Yi中i=0)中没有混频的欠采样频点。由于Y0为零散时间的傅里叶变换,周期为2π,所以分为(0,π)和(π,2π)两部分寻找Y0中没有混频的欠采样频点。
搜索Y0在(0,π)中使得幅值|Y0(ejω)|不会为零的频点ωk1。由于输入函数为余弦信号,有且只有一个ωk对应于使得Y的幅值函数
搜索Y0在(π,2π)中的非混频的欠采样频点,其中寻找方法同理,此时Y的幅值函数变为同理找到欠采样频点ωk2和对应的kp2。
在步骤S110中,在一个实施例中,根据第三频点和第三频点对应的通道,获得各个子通道的时间误差的步骤可以包括:
根据ωk1和kp1,以及表达式:得到Δt1,其中i=(0.....M-1),Yi为处理后的第一函数中各通道的函数, ω<2π>=((ω+π)mod2π)-π,mod为求余运算,T=1/f,f为总采样率,即f=fs×M;
根据ωk2和kp2,以及表达式:得到Δt2;
根据Δt1和Δt2的平均值获得各个子通道的时间误差。
上述Δt1和Δt2的计算过程为:将ωk1替换表达式中的ωk,kp1替换表达式中的Kp,得到Δt1;将ωk2替换表达式中的ωk,kp2替换表达式中的Kp,得到Δt1。
为了更好的理解本发明的实施方式,下面结合一个M=4通道的TIADC的具体实施例进行详细说明。
如图3所示,为本发明的多通道TIADC时间交织***示意图,输入信号输入M通道的TIADC,每条通道以相同的采样率但不同的采样时刻(相邻通道相差T时刻)对高速输入信号采样,最终合并出输出信号,以此实现高速采样的模数转化。
如图4所示,为各个子通道中的失配误差情况。由于现实中布线等原因会产生时间失配误差Δti,同时会产生加性的失配偏置误差ai和乘性失配增益误差gi。
输入测试信号为余弦函数Acos(2πfint+θ),设置A=1,f为总采样率,m为正整数,nfft为2的正整数次幂,nfft就是总采样数。这里设m=10027,nfft=216,f=16×109。每条子通道的偏置误差a=[15,20,25,30],增益误差g=[6,5,6,5],时间误差
在没有采取本发明方案时,仿真结果如图5所示,当输入的余弦信号时,4通道中含有三种失配误差(偏置失配误差,增益失配误差,时间失配误差),输出信号严重影响了原信号的波形。由图6可以清楚地看出,在频域上原信号的频谱出现了偏移并出现了很多伪谱线。
本发明的主要思想是:先估计增益误差与偏置误差,然后再对信号校正,消除增益误差与偏置误差,最后估计时间误差。
首先对误差信号 进行65536个点的FFT快速傅里叶变换,在固定的频点上得到与由公式 就可以估计出的每条子通道的偏置误差ak=[15.00,20.00,25.00,30.00]与增益误差gk=[6.00,5.00,6.00,5.00]。然后每条子通道经过公式可以消除每条子通道内的增益误差和偏置误差。
如图7所示,为校正增益误差与偏置误差后的信号在时域内的表示与原信号与校正之后的信号之差(即误差信号)在时域内的表示。可以看到误差信号已经控制在±2%之内。
如图8所示,为输入信号在频域上的表示与校正增益误差与偏置误差的校正后信号在频域上的表示。相对于原来同时存在三种失配误差信号频谱上的误差,可以看到消除了大部分的伪谱线和频线搬移,还剩下由于时间失配误差产生的误差。
然后再对每条子通道进行16384个点的FFT快速傅里叶变换,得到Yi,i=(0,1,2,3),再把4条通道的时域信号合并成一个信号,再进行65536个点的FFT快速傅里叶变化,得到Y,然后对Yi和Y预处理,其中ε可以设置为0.06。寻找非混频的频点ωk1和相应的Kp1,ωk2和相应的Kp2。根据相应的ωk和Kp,可以得到Δt1,Δt2。再求Δt1,Δt2的均值就可以得到每条子通道的时间误差估计。最后仿真得到的时间误差估计为得到的时间误差估计可以通过后续的方法进行消除,以达到同时消除三种失配误差的目的。本发明不涉及时间误差后续的消除校正方法。
可以对时间误差估计的效果进行评估,考虑到RMSE(root-mean-squareerror)作为标准,估计实际的时间误差与估计误差之间的误差效果。与之前的实例数据一致,时间误差以每个子通道0.01的步长增加,观察其RMSE的变化。可以从图9中看到,总体的误差RMSE低于-114dB。而且估计的精度随着时间失配误差的增大而增大,证明了本发明的有效性。
同理也可以对估计的增益误差与偏置误差的效果进行评估。在评估增益误差的估计效果时,每条子通道的偏置误差a=[15,20,25,30],时间误差Δt/T=[0,0.03,0.01,-0.02],增益误差g=[6,5,6,5],每个子通道的增益误差以步长为5而增加。其结果如图10所示,随着增益误差的增加,其RMSE也增加,但是总体低于-250dB。在评估偏置误差的估计效果时,每条子通道的增益误差g=[6,5,6,5],时间误差Δt/T=[0,0.03,0.01,-0.02],偏置误差a=[15,20,25,30],每个子通道的偏置误差以步长为5而增加。其结果如图11所示,总体的RMSE低于-276dB。
从以上的实验结果可以看出本发明对多通道时间交织模数转换***的增益误差、偏置误差和时间误差的估计方法的有效性与高精度。
基于同一发明构思,本发明还提供一种用于时间交织模数转换***的误差估计装置,下面对本发明装置的具体实施方式做详细描述。
如图12所示,一种用于时间交织模数转换***的误差估计装置,包括:
误差信号获取模块101,用于在多通道时间交织模数转换***中输入测试信号,获取输出的误差信号;
频点确定模块102,用于将所述误差信号进行傅里叶变换,获得傅里叶变换后增益误差出现的第一频点和偏置误差出现的第二频点;
傅里叶变换模块103,用于根据所述误差信号的傅里叶变换在第一频点的值,获得增益误差的傅里叶变换;根据所述误差信号的傅里叶变换在第二频点的值,获得偏置误差的傅里叶变换;
增益和偏置误差确定模块104,用于根据增益误差的傅里叶变换和偏置误差的傅里叶变换,获得各个子通道的增益误差和偏置误差;
采样函数更新模块105,用于根据原采样函数、增益误差和偏置误差得到各个子通道的新的采样函数;
第一函数确定模块106,用于将各个子通道的新的采样函数进行傅里叶变换,得到第一函数;
第二函数确定模块107,用于将各个子通道的新的采样函数合并后进行傅里叶变换,得到第二函数;
函数处理模块108,用于根据预设条件对所述第一函数和所述第二函数进行处理;
频点和通道确定模块109,用于根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道;
时间误差确定模块110,用于根据第三频点和第三频点对应的通道,获得各个子通道的时间误差。
在一个实施例中,所述采样函数更新模块105可以根据以下表达式得到各个子通道的新的采样函数:
其中,i=(0.....M-1),ai为偏置误差,gi为增益误差,为各个子通道的新的采样函数,yADCi为原采样函数。
在一个实施例中,如图13所示,所述函数处理模块108包括:
第一函数处理单元1081,用于将第一函数中小于ε×MAX(Yi)的值设为0,其中ε为预设常数,Yi为第一函数中各通道的函数,i=(0.....M-1),M为总的通道数,MAX为最大值;
第二函数处理单元1082,用于将第二函数中小于ε×MAX(Y)的值设为0,其中Y为第二函数。
在一个实施例中,所述测试信号为余弦信号,如图14所示,所述频点和通道确定模块109可以包括:
第一获取单元1091,用于获取Y0在(0,π)中使得幅值|Y0(ejω)|不为零的频点ωk1,其中Y0为处理后的第一函数中第一通道的函数,ω为角频率,e为自然指数,j为虚数单位;
第一通道确定单元1092,用于从k=1,...,M中选取使不为零的k,将选取的k作为kp1,其中M为总的通道数,Y为处理后的第二函数;
第二获取单元1093,用于获取Y0在(π,2π)中使得幅值|Y0(ejω)|不为零的频点ωk2;
第二通道确定单元1094,用于从k=1,...,M中选取使不为零的k,将选取的k作为kp2。
在一个实施例中,如图15所示,所述时间误差确定模块110可以包括:
第一时间误差确定单元1101,用于根据ωk1和kp1,以及表达式: 得到Δt1,其中i=(0.....M-1),Yi为处理后的第一函数中各通道的函数,ω<2π>=((ω+π)mod2π)-π,mod为求余运算,T=1/f,f为总采样率;
第二时间误差确定单元1102,用于根据ωk2和kp2,以及表达式: 得到Δt2;
第三时间误差确定单元1103,用于根据Δt1和Δt2的平均值获得各个子通道的时间误差。
本发明装置其它技术特征与本发明方法相同,在此不予赘述。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种用于时间交织模数转换***的误差估计方法,其特征在于,包括步骤:
在多通道时间交织模数转换***中输入测试信号,获取输出的误差信号;
将所述误差信号进行傅里叶变换,获得傅里叶变换后增益误差出现的第一频点和偏置误差出现的第二频点;
根据所述误差信号的傅里叶变换在第一频点的值,获得增益误差的傅里叶变换;根据所述误差信号的傅里叶变换在第二频点的值,获得偏置误差的傅里叶变换;
根据增益误差的傅里叶变换和偏置误差的傅里叶变换,获得各个子通道的增益误差和偏置误差;
根据原采样函数、增益误差和偏置误差得到各个子通道的新的采样函数;
将各个子通道的新的采样函数进行傅里叶变换,得到第一函数;
将各个子通道的新的采样函数合并后进行傅里叶变换,得到第二函数;
根据预设条件对所述第一函数和所述第二函数进行处理;
根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道;
根据第三频点和第三频点对应的通道,获得各个子通道的时间误差。
2.根据权利要求1所述的用于时间交织模数转换***的误差估计方法,其特征在于,根据预设条件对所述第一函数和所述第二函数进行处理的步骤包括:
将第一函数中小于ε×MAX(Yi)的值设为0,其中ε为预设常数,Yi为第一函数中各通道的函数,i=(0…..M-1),M为总的通道数,MAX为最大值;
将第二函数中小于ε×MAX(Y)的值设为0,其中Y为第二函数。
3.根据权利要求1所述的用于时间交织模数转换***的误差估计方法,其特征在于,所述测试信号为余弦信号,根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道的步骤包括:
获取Y0在(0,π)中使得幅值|Y0(ejω)|不为零的频点ωk1,其中Y0为处理后的第一函数中第一通道的函数,ω为角频率,e为自然指数,j为虚数单位;
从k=1,…,M中选取使不为零的k,将选取的k作为kp1,其中M为总的通道数,Y为处理后的第二函数;
获取Y0在(π,2π)中使得幅值|Y0(ejω)|不为零的频点ωk2;
从k=1,…,M中选取使不为零的k,将选取的k作为kp2。
4.根据权利要求3所述的用于时间交织模数转换***的误差估计方法,其特征在于,根据第三频点和第三频点对应的通道,获得各个子通道的时间误差的步骤包括:
根据ωk1和kp1,以及表达式:得到Δt1,其中i=(0…..M-1),Yi为处理后的第一函数中各通道的函数, ω<2π>=((ω+π)mod2π)-π,mod为求余运算,Y=1/f,f为总采样率;
根据ωk2和kp2,以及表达式:得到Δt2;
根据Δt1和Δt2的平均值获得各个子通道的时间误差。
5.根据权利要求1至4任意一项所述的用于时间交织模数转换***的误差估计方法,其特征在于,根据表达式得到各个子通道的新的采样函数,其中i=(0…..M-1),ai为偏置误差,gi为增益误差,为各个子通道的新的采样函数,yADCi为原采样函数。
6.一种用于时间交织模数转换***的误差估计装置,其特征在于,包括:
误差信号获取模块,用于在多通道时间交织模数转换***中输入测试信号,获取输出的误差信号;
频点确定模块,用于将所述误差信号进行傅里叶变换,获得傅里叶变换后增益误差出现的第一频点和偏置误差出现的第二频点;
傅里叶变换模块,用于根据所述误差信号的傅里叶变换在第一频点的值,获得增益误差的傅里叶变换;根据所述误差信号的傅里叶变换在第二频点的值,获得偏置误差的傅里叶变换;
增益和偏置误差确定模块,用于根据增益误差的傅里叶变换和偏置误差的傅里叶变换,获得各个子通道的增益误差和偏置误差;
采样函数更新模块,用于根据原采样函数、增益误差和偏置误差得到各个子通道的新的采样函数;
第一函数确定模块,用于将各个子通道的新的采样函数进行傅里叶变换,得到第一函数;
第二函数确定模块,用于将各个子通道的新的采样函数合并后进行傅里叶变换,得到第二函数;
函数处理模块,用于根据预设条件对所述第一函数和所述第二函数进行处理;
频点和通道确定模块,用于根据处理后的第一函数和第二函数,获得第一通道中没有混频的第三频点,以及与第三频点对应的通道;
时间误差确定模块,用于根据第三频点和第三频点对应的通道,获得各个子通道的时间误差。
7.根据权利要求6所述的用于时间交织模数转换***的误差估计装置,其特征在于,所述函数处理模块包括:
第一函数处理单元,用于将第一函数中小于ε×MAX(Yi)的值设为0,其中ε为预设常数,Yi为第一函数中各通道的函数,i=(0…..M-1),M为总的通道数,MAX为最大值;
第二函数处理单元,用于将第二函数中小于ε×MAX(Y)的值设为0,其中Y为第二函数。
8.根据权利要求6所述的用于时间交织模数转换***的误差估计装置,其特征在于,所述测试信号为余弦信号,所述频点和通道确定模块包括:
第一获取单元,用于获取Y0在(0,π)中使得幅值|Y0(ejω)|不为零的频点ωk1,其中Y0为处理后的第一函数中第一通道的函数,ω为角频率,e为自然指数,j为虚数单位;
第一通道确定单元,用于从k=1,…,M中选取使不为零的k,将选取的k作为kp1,其中M为总的通道数,Y为处理后的第二函数;
第二获取单元,用于获取Y0在(π,2π)中使得幅值|Y0(ejω)|不为零的频点ωk2;
第二通道确定单元,用于从k=1,…,M中选取使不为零的k,将选取的k作为kp2。
9.根据权利要求8所述的用于时间交织模数转换***的误差估计装置,其特征在于,所述时间误差确定模块包括:
第一时间误差确定单元,用于根据ωk1和kp1,以及表达式: 得到Δt1,其中i=(0…..M-1),Yi为处理后的第一函数中各通道的函数,ω<2π>=((ω+π)mod2π)-π,mod为求余运算,T=1/f,f为总采样率;
第二时间误差确定单元,用于根据ωk2和kp2,以及表达式: 得到Δt2;
第三时间误差确定单元,用于根据Δt1和Δt2的平均值获得各个子通道的时间误差。
10.根据权利要求6至9任意一项所述的用于时间交织模数转换***的误差估计装置,其特征在于,所述采样函数更新模块根据表达式得到各个子通道的新的采样函数,其中i=(0…..M-1),ai为偏置误差,gi为增益误差,为各个子通道的新的采样函数,yADCi为原采样函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044047.7A CN105720983A (zh) | 2016-01-22 | 2016-01-22 | 用于时间交织模数转换***的误差估计方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044047.7A CN105720983A (zh) | 2016-01-22 | 2016-01-22 | 用于时间交织模数转换***的误差估计方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105720983A true CN105720983A (zh) | 2016-06-29 |
Family
ID=56153720
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610044047.7A Pending CN105720983A (zh) | 2016-01-22 | 2016-01-22 | 用于时间交织模数转换***的误差估计方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105720983A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105846822A (zh) * | 2016-03-24 | 2016-08-10 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 测量时间交织模数转换***的失配误差的方法和*** |
CN106209103A (zh) * | 2016-07-29 | 2016-12-07 | 电子科技大学 | 基于频谱分析的tiadc增益和时间误差的校正方法 |
CN106341132A (zh) * | 2016-08-08 | 2017-01-18 | 中国工程物理研究院电子工程研究所 | 时间交织采样adc的误差盲校正方法 |
CN106374920A (zh) * | 2016-09-05 | 2017-02-01 | 中山大学 | 一种基于多项式模型的tiadc***的估计与补偿实现方法 |
CN106375252A (zh) * | 2016-09-05 | 2017-02-01 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种iq与tiadc频率联合失真的修正方法及*** |
CN107046425A (zh) * | 2017-04-18 | 2017-08-15 | 中国电子科技集团公司第二十四研究所 | 基于采样间隔差值统计的采样时间误差检测*** |
CN112636754A (zh) * | 2020-11-06 | 2021-04-09 | 中山大学 | 宽带交叉采样***通道失配估计方法、***、装置及介质 |
CN113595553A (zh) * | 2021-07-22 | 2021-11-02 | 华中科技大学 | 一种基于时间交织并行采样的数据采集***及方法 |
CN114267407A (zh) * | 2022-03-03 | 2022-04-01 | 合肥悦芯半导体科技有限公司 | 一种精度校正方法、装置、***以及精度校正设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110063149A1 (en) * | 2008-12-29 | 2011-03-17 | Intersil Americas Inc. | Calibration of offset, gain and phase errors in m-channel time-interleaved analog-to-digital converters |
CN102075464A (zh) * | 2011-01-18 | 2011-05-25 | 电子科技大学 | Tiadc***通道误差的联合估计及实时校正方法 |
CN102761335A (zh) * | 2011-04-29 | 2012-10-31 | 英特希尔美国公司 | 在双信道时间交错模拟至数字转换器中的取样时间和增益不匹配误差估计的无乘法器算法 |
CN103812506A (zh) * | 2012-11-13 | 2014-05-21 | 邹月娴 | 一种基于信号频域稀疏性的tiadc时间失配参数盲测量方法 |
US20140266823A1 (en) * | 2013-03-15 | 2014-09-18 | Integrated Device Technology, Inc. | Systems and methods for estimation of offset and gain errors in a time-interleaved analog-to-digital converter |
US20150070198A1 (en) * | 2013-09-06 | 2015-03-12 | Broadcom Corporation | Flexible ADC Calibration Technique Using ADC Capture Memory |
US9088293B1 (en) * | 2014-09-04 | 2015-07-21 | Semtech Corporation | Calibration of a time-interleaved analog-to-digital converter (ADC) circuit |
US9143147B1 (en) * | 2014-07-03 | 2015-09-22 | Keysight Technologies, Inc. | Calibration of inter-slice gain and offset errors in time-interleaved analog-to- digital converter |
US20150303934A1 (en) * | 2014-04-17 | 2015-10-22 | The Board Of Regents, The University Of Texas System | System and method for dynamic path-mismatch equalization in time-interleaved adc |
CN105024696A (zh) * | 2015-07-02 | 2015-11-04 | 大唐微电子技术有限公司 | 多通道并行模数转换***采样时间误差的校准装置及方法 |
-
2016
- 2016-01-22 CN CN201610044047.7A patent/CN105720983A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110063149A1 (en) * | 2008-12-29 | 2011-03-17 | Intersil Americas Inc. | Calibration of offset, gain and phase errors in m-channel time-interleaved analog-to-digital converters |
CN102075464A (zh) * | 2011-01-18 | 2011-05-25 | 电子科技大学 | Tiadc***通道误差的联合估计及实时校正方法 |
CN102761335A (zh) * | 2011-04-29 | 2012-10-31 | 英特希尔美国公司 | 在双信道时间交错模拟至数字转换器中的取样时间和增益不匹配误差估计的无乘法器算法 |
CN103812506A (zh) * | 2012-11-13 | 2014-05-21 | 邹月娴 | 一种基于信号频域稀疏性的tiadc时间失配参数盲测量方法 |
US20140266823A1 (en) * | 2013-03-15 | 2014-09-18 | Integrated Device Technology, Inc. | Systems and methods for estimation of offset and gain errors in a time-interleaved analog-to-digital converter |
US20150070198A1 (en) * | 2013-09-06 | 2015-03-12 | Broadcom Corporation | Flexible ADC Calibration Technique Using ADC Capture Memory |
US20150303934A1 (en) * | 2014-04-17 | 2015-10-22 | The Board Of Regents, The University Of Texas System | System and method for dynamic path-mismatch equalization in time-interleaved adc |
US9143147B1 (en) * | 2014-07-03 | 2015-09-22 | Keysight Technologies, Inc. | Calibration of inter-slice gain and offset errors in time-interleaved analog-to- digital converter |
US9088293B1 (en) * | 2014-09-04 | 2015-07-21 | Semtech Corporation | Calibration of a time-interleaved analog-to-digital converter (ADC) circuit |
CN105024696A (zh) * | 2015-07-02 | 2015-11-04 | 大唐微电子技术有限公司 | 多通道并行模数转换***采样时间误差的校准装置及方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105846822A (zh) * | 2016-03-24 | 2016-08-10 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 测量时间交织模数转换***的失配误差的方法和*** |
CN106209103A (zh) * | 2016-07-29 | 2016-12-07 | 电子科技大学 | 基于频谱分析的tiadc增益和时间误差的校正方法 |
CN106209103B (zh) * | 2016-07-29 | 2019-09-24 | 电子科技大学 | 基于频谱分析的tiadc增益和时间误差的校正方法 |
CN106341132A (zh) * | 2016-08-08 | 2017-01-18 | 中国工程物理研究院电子工程研究所 | 时间交织采样adc的误差盲校正方法 |
CN106341132B (zh) * | 2016-08-08 | 2019-05-24 | 中国工程物理研究院电子工程研究所 | 时间交织采样adc的误差盲校正方法 |
CN106375252B (zh) * | 2016-09-05 | 2019-05-03 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种iq与tiadc频率联合失真的修正方法及*** |
CN106375252A (zh) * | 2016-09-05 | 2017-02-01 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种iq与tiadc频率联合失真的修正方法及*** |
CN106374920A (zh) * | 2016-09-05 | 2017-02-01 | 中山大学 | 一种基于多项式模型的tiadc***的估计与补偿实现方法 |
CN107046425A (zh) * | 2017-04-18 | 2017-08-15 | 中国电子科技集团公司第二十四研究所 | 基于采样间隔差值统计的采样时间误差检测*** |
CN107046425B (zh) * | 2017-04-18 | 2020-09-01 | 中国电子科技集团公司第二十四研究所 | 基于采样间隔差值统计的采样时间误差检测*** |
CN112636754A (zh) * | 2020-11-06 | 2021-04-09 | 中山大学 | 宽带交叉采样***通道失配估计方法、***、装置及介质 |
CN113595553A (zh) * | 2021-07-22 | 2021-11-02 | 华中科技大学 | 一种基于时间交织并行采样的数据采集***及方法 |
CN113595553B (zh) * | 2021-07-22 | 2023-10-20 | 华中科技大学 | 一种基于时间交织并行采样的数据采集***及方法 |
CN114267407A (zh) * | 2022-03-03 | 2022-04-01 | 合肥悦芯半导体科技有限公司 | 一种精度校正方法、装置、***以及精度校正设备 |
CN114267407B (zh) * | 2022-03-03 | 2022-06-10 | 合肥悦芯半导体科技有限公司 | 一种精度校正方法、装置、***以及精度校正设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105720983A (zh) | 用于时间交织模数转换***的误差估计方法和装置 | |
Kang et al. | Phase difference correction method for phase and frequency in spectral analysis | |
CN104007316B (zh) | 一种欠采样速率下的高精度频率测量方法及其测量仪 | |
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN106483374A (zh) | 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法 | |
CN107102255A (zh) | 单一adc采集通道动态特性测试方法 | |
CN101674087A (zh) | 一种时间交替adc***通道失配误差的获取方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
CN103647554A (zh) | 一种测量时间交替adc采集***偏置与增益误差的方法 | |
CN105915221A (zh) | 一种tiadc***及方法 | |
CN105510706B (zh) | 一种高精度欠采样测频方法 | |
Li et al. | A background correlation-based timing skew estimation method for time-interleaved ADCs | |
Liu et al. | Adaptive semiblind background calibration of timing mismatches in a two-channel time-interleaved analog-to-digital converter | |
CN106341132A (zh) | 时间交织采样adc的误差盲校正方法 | |
CN109274372A (zh) | 一种tiadc***通道间采样时刻失配误差提取方法 | |
CN105044769B (zh) | 提高地震信号的分辨率的方法 | |
CN101308175A (zh) | 相位谱分析仪 | |
Ramos et al. | Comparison of impedance measurements in a DSP using ellipse-fit and seven-parameter sine-fit algorithms | |
CN103812506A (zh) | 一种基于信号频域稀疏性的tiadc时间失配参数盲测量方法 | |
CN115001494B (zh) | 一种交织采样的后台自适应自校准方法 | |
CN110285881A (zh) | 一种基于全相位滤波的密集谱频率估计法 | |
US10498348B1 (en) | Time error and gain offset estimation in interleaved analog-to-digital converters | |
CN105846822A (zh) | 测量时间交织模数转换***的失配误差的方法和*** | |
Asami et al. | Technique to improve the performance of time-interleaved AD converters with mismatches of non-linearity | |
Song et al. | Digital estimation and compensation of analog errors in frequency-interleaved ADCs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination |