CN105680858A - 一种估计tiadc并行采集系统时间偏移误差的方法 - Google Patents
一种估计tiadc并行采集系统时间偏移误差的方法 Download PDFInfo
- Publication number
- CN105680858A CN105680858A CN201610023333.5A CN201610023333A CN105680858A CN 105680858 A CN105680858 A CN 105680858A CN 201610023333 A CN201610023333 A CN 201610023333A CN 105680858 A CN105680858 A CN 105680858A
- Authority
- CN
- China
- Prior art keywords
- phase
- passage
- tiadc
- time offset
- error
- 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.)
- Granted
Links
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
Abstract
本发明公开了一种估计TIADC并行采集系统时间偏移误差的方法,首先通过各通道ADC采样得到各通道采样值,然后分别对各通道采样值进行希尔伯特变换,接着进行反正切、相位补偿以及相位解旋操作得到各通道采样输出离散信号的相位。假设通道0为参考通道,即不存在时间偏移误差,则可根据前文离散信号相位与时间偏移误差的关系求出各个通道存在的时间偏移误差。
Description
技术领域
本发明属于时间误差估算技术领域,更为具体地讲,涉及一种估计TIADC并行采集系统时间偏移误差的方法。
背景技术
当前,典型TIADC并行采集系统用M片采样率为fs/M的ADC进行前端采样,通过精确控制相邻两个通道ADC之间的采样时钟相位差为Ts=1/fs来实现M片ADC的并行逐次采样,再通过后端合成拼接实现以相同分辨率,而等效采样率提高为原来单片ADC采样率M倍,即fs对模拟输入信号进行采样。
如图1所示,正常情况下,理想的ADC并行采样不存在通道失配误差,然而,在实际工程应用中,单片ADC芯片的制作工艺以及印制电路板的布局布线不可能做到完全相同,这些差异将会引入通道失配误差。
其中,TIADC并行采集系统通道失配误差主要有三种:增益、偏置及时间偏移误差。三种误差在各个通道中如图2所示,其中,gi,oi(i=0,1,...,M-1)分别表示增益,偏置误差,Ti表示各通道ADC采样时刻,即
Ti=kMTs+iTs+δiTs
其中,δi表示通道i中的时间偏移相对误差,k表示每个通道抽样值序号。
理想情况下,gi、oi和δi三种通道失配误差值都应该为0,否则由失配误差引起的谐波分量将会出现在系统采样输出信号的频谱中,影响整个采样系统信噪比以及无杂散动态范围等指标,降低系统采样输出信号质量。
发明内容
本发明的目的在于克服现有技术的不足,提供一种估计TIADC并行采集系统时间偏移误差的方法,在假设整个TIADC并行采集系统的偏置与增益误差gi=oi=0的条件下,通过估算时间偏移误差来提高系统采样输出信号质量。
为实现上述发明目的,本发明一种估计TIADC并行采集系统时间偏移误差的方法,其特征在于,包括以下步骤:
(1)、设频率及初相已知的正弦信号x(t)为:
x(t)=Acos(2πf0t+θ)
其中,A为正弦信号x(t)的幅度,f0为正弦信号x(t)的频率,θ为正弦信号x(t)的初相;
(2)、将正弦信号x(t)经过TIADC并行采集系统采样后,得到各通道的输出信号,可表示为:
xi(n)=Acos(2πf0nMTs+2πf0(i+δi)Ts+θ)
其中,M为TIADC并行采集系统中ADC的总个数,i=0,1,2,…,M-1,δi为时间偏移误差,n为离散信号抽样值序号,Ts为TIADC并行采集系统等效采样周期;
设0<f0<fs/2且n≥0,fs为TIADC并行采集系统等效采样频率;用表示xi(n)的相位,即
(3)、将各个通道的输出信号进行希尔伯特变换
(3.1)、将输入的正弦信号x(t)的频率范围f0∈(0,fs/2)作如下的子区间划分
(3.2)、将通道i的输出信号xi(n)进行希尔伯特变换,得到信号为xih(n);
当m为偶数时:
则有:xih(n)=Asin(2πf0nMTs+2πf0(i+δi)Ts+θ);
当m为奇数时:
则有:xih(n)=-Asin(2πf0nMTs+2πf0(i+δi)Ts+θ);
(4)、计算的估计值
(5)、根据点(xi(n),xih(n))具体所在象限对进行补偿,使[-π/2,π/2]区间的相位角映射到[0,2π]区间范围内,即:
(6)、利用相位解旋算法去除的相位模糊现象
设修正项c(n)满足:
且c(0)=0;
将修正项c(n)与相加,得到相位估计值即
设rtk(k=0,1,..,M-2)表示通道k+1与通道k之间的时间偏置误差差值,即
rtk=δk+1-δk
结合相位估计值如果相位估计值的样本长度为N,则rtk可以表示为
设通道0为参考通道,不存在时间偏移误差,即δ0=0,则每个通道的时间偏移误差估计值可以表示为
其中,k=0,1,..,M-2。
本发明的发明目的是这样实现的:
本发明一种估计TIADC并行采集系统时间偏移误差的方法,首先通过各通道ADC采样得到各通道采样值,然后分别对各通道采样值进行希尔伯特变换,接着进行反正切、相位补偿以及相位解旋操作得到各通道采样输出离散信号的相位。假设通道0为参考通道,即不存在时间偏移误差,则可根据前文离散信号相位与时间偏移误差的关系求出各个通道存在的时间偏移误差。
同时,本发明一种估计TIADC并行采集系统时间偏移误差的方法还具有以下有益效果:
本发明是基于数字信号处理中运用普遍的希尔伯特变换,原理简单。其次,借助一个已知频率及初相的正弦信号,通过估计各相邻通道采样输出信号之间相位差计算得到各通道之间存在的采样时间偏移误差,具有较高的估计精度,同时可有效降低硬件成本。
附图说明
图1是TIADC并行采集系统原理框图;
图2是TIADC并行采集系统通道失配误差示意图;
图3是本发明一种估计TIADC并行采集系统时间偏移误差的方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图3是本发明一种估计TIADC并行采集系统时间偏移误差的方法流程图。
在本实施例中,如图3所示,本发明一种估计TIADC并行采集系统时间偏移误差的方法,包括以下步骤:
S1、设频率及初相已知的正弦信号x(t)为:
x(t)=cos(2π×1000×t+π/2)
即,输入正弦信号频率为f0=1000Hz,初相为π/2;
S2、将正弦信号x(t)经过M=4通道TIADC并行采集系统(设系统等效采样周期Ts=10-4s)采样后,得到各通道的输出信号,可表示为:
xi(n)=cos(2π×1000×4n×10-4+2π×1000×(i+δi)×10-4+π/2)
其中,i=0,1,2,3,δi为时间偏移误差,n为离散信号抽样值序号;
设n≥0,用表示xi(n)的相位,即
从上式可以看出,只要求出就可以计算得到时间偏移误差δi。
S3、将各个通道的输出信号进行希尔伯特变换
首先,我们对希尔伯特变换进行简单介绍:理想希尔伯特变换的频率响应如下
式所示
其周期为2π。从上式可以看出,希尔伯特变换器是幅频特性为1的全通滤波器,信号通过希尔伯特变换器之后,其负频率部分作π/2相移,而正频率成分作-π/2相移。可以求得希尔伯特变换的冲激响应为:
其中,n表示离散信号抽样值序号,则hhil(n)与Hhil(ω)之间有如下的关系:
图1中,各通道ADC的采样率为fs/M,对单独某一个通道而言,输入信号频率f0有可能会大于通道采样率fs/M。
因此,需要先将输入的正弦信号x(t)的频率范围f0∈(0,fs/2)作如下的子区间划分
再将通道i的输出信号xi(n)进行希尔伯特变换,得到信号为xih(n);
当m为偶数时:
则有:xih(n)=Asin(2πf0nMTs+2πf0(i+δi)Ts+θ);
当m为奇数时:
则有:xih(n)=-Asin(2πf0nMTs+2πf0(i+δi)Ts+θ);
由前文,f0=1000Hz,fs=1/Ts=10000Hz,可求得m=0,为偶数,因此可得到xih(n)=sin(2π×1000×4n×10-4+2π×1000×(i+δi)×10-4+π/2);
S4、计算的估计值
上式中,由于反正切函数的值域为[-π/2,π/2],因此需要根据点(xi(n),xih(n))具体位于哪个象限对上式进行补偿;
S5、根据点(xi(n),xih(n))具体所在象限对进行补偿,使[-π/2,π/2]区间的相位角映射到[0,2π]区间范围内,即:
经过上式所示的相位补偿后,其相位角区间为[0,2π),然而是关于n的单调函数,而上式中经过修正后的取值范围为[0,2π),因此非单调,这一现象称为相位模糊。
S6、利用相位解旋算法去除的相位模糊现象
设修正项c(n)满足:
且c(0)=0;
将修正项c(n)与相加,得到相位估计值即
设rtk(k=0,1,..,M-2)表示通道k+1与通道k之间的时间偏置误差差值,即
rtk=δk+1-δk
结合相位估计值如果相位估计值的样本长度为N,则rtk可以表示为
设通道0为参考通道,不存在时间偏移误差,即δ0=0,则每个通道的时间偏移误差估计值可以表示为
其中,k=0,1,..,M-2。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (2)
1.一种估计TIADC并行采集系统时间偏移误差的方法,其特征在于,包括以下步骤:
(1)、设频率及初相已知的正弦信号x(t)为:
x(t)=Acos(2πf0t+θ)
其中,A为正弦信号x(t)的幅度,f为正弦信号x(t)的频率,θ为正弦信号x(t)的初相;
(2)、将正弦信号x(t)经过TIADC并行采集系统采样后,得到各通道的输出信号,可表示为:
xi(n)=Acos(2πf0nMTs+2πf0(i+δi)Ts+θ)
其中,M为TIADC并行采集系统中ADC的总个数,i=0,1,2,…,M-1,δi为时间偏移误差,n为离散信号抽样值序号,Ts为TIADC并行采集系统等效采样周期;
设0<f0<fs/2且n≥0,fs为TIADC并行采集系统等效采样频率;用表示xi(n)的相位,即
(3)、将各个通道的输出信号进行希尔伯特变换
(3.1)、将输入的正弦信号x(t)的频率范围f0∈(0,fs/2)作如下的子区间划分
(3.2)、将通道i的输出信号xi(n)进行希尔伯特变换,得到信号为xih(n);
当m为偶数时:
则有:xih(n)=Asin(2πf0nMTs+2πf0(i+δi)Ts+θ);
当m为奇数时:
则有:xih(n)=-Asin(2πf0nMTs+2πf0(i+δi)Ts+θ);
(4)、计算的估计值
(5)、根据点(xi(n),xih(n))具体所在象限对进行补偿,使[-π/2,π/2]区间的相位角映射到[0,2π]区间范围内,即:
(6)、利用相位解旋算法去除的相位模糊现象
设修正项c(n)满足:
且c(0)=0;
将修正项c(n)与相加,得到相位估计值即
设rtk(k=0,1,..,M-2)表示通道k+1与通道k之间的时间偏置偏置误差差值,即
rtk=δk+1-δk
结合相位估计值如果相位估计值的样本长度为N,则rtk可以表示为
设通道0为参考通道,不存在时间偏移误差,即δ0=0,则每个通道的时间偏移误差估计值可以表示为
其中,k=0,1,..,M-2。
2.根据权利要求1所述的估计TIADC并行采集系统时间偏移误差的方法,其特征在于,所述的希尔伯特变换是利用幅频特性为1的全通滤波器对各通道的输出信号进行变换处理,使其负频率部分作π/2相移,而正频率部分作-π/2相移,其希尔伯特变换换的冲激响应可表示为:
其中,n表示离散信号抽样值序号,hhil(n)与Hhil(ω)之间有如下的关系:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610023333.5A CN105680858B (zh) | 2016-01-14 | 2016-01-14 | 一种估计tiadc并行采集系统时间偏移误差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610023333.5A CN105680858B (zh) | 2016-01-14 | 2016-01-14 | 一种估计tiadc并行采集系统时间偏移误差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105680858A true CN105680858A (zh) | 2016-06-15 |
CN105680858B CN105680858B (zh) | 2018-08-07 |
Family
ID=56300682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610023333.5A Expired - Fee Related CN105680858B (zh) | 2016-01-14 | 2016-01-14 | 一种估计tiadc并行采集系统时间偏移误差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105680858B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107453755A (zh) * | 2017-07-11 | 2017-12-08 | 电子科技大学 | 一种基于混合架构的高速高精度多通道并行采集系统 |
CN107994902A (zh) * | 2017-12-22 | 2018-05-04 | 北京工业大学 | 一种基于机器学习的双通道模数转换器子转换器间时间误差校正方法 |
CN115913231A (zh) * | 2023-01-06 | 2023-04-04 | 上海芯炽科技集团有限公司 | 一种tiadc的采样时间误差数字估计方法 |
CN116840549A (zh) * | 2023-07-03 | 2023-10-03 | 广州致远仪器有限公司 | 时间交织adc的误差确定方法、装置、设备以及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080291066A1 (en) * | 2006-09-29 | 2008-11-27 | Optichron, Inc. | Composite analog to digital receiver with adaptive self-linearization |
CN101783683A (zh) * | 2008-12-29 | 2010-07-21 | 英特赛尔美国股份有限公司 | 双通道时间交错型模数转换器中的误差估计与校正 |
CN102457277A (zh) * | 2010-10-27 | 2012-05-16 | 英特希尔美国公司 | 用于时间交错式模拟至数字转换器的稳健增益和相位校准方法 |
CN103647554A (zh) * | 2013-12-16 | 2014-03-19 | 电子科技大学 | 一种测量时间交替adc采集系统偏置与增益误差的方法 |
-
2016
- 2016-01-14 CN CN201610023333.5A patent/CN105680858B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080291066A1 (en) * | 2006-09-29 | 2008-11-27 | Optichron, Inc. | Composite analog to digital receiver with adaptive self-linearization |
CN101783683A (zh) * | 2008-12-29 | 2010-07-21 | 英特赛尔美国股份有限公司 | 双通道时间交错型模数转换器中的误差估计与校正 |
CN102457277A (zh) * | 2010-10-27 | 2012-05-16 | 英特希尔美国公司 | 用于时间交错式模拟至数字转换器的稳健增益和相位校准方法 |
CN103647554A (zh) * | 2013-12-16 | 2014-03-19 | 电子科技大学 | 一种测量时间交替adc采集系统偏置与增益误差的方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107453755A (zh) * | 2017-07-11 | 2017-12-08 | 电子科技大学 | 一种基于混合架构的高速高精度多通道并行采集系统 |
CN107453755B (zh) * | 2017-07-11 | 2019-12-27 | 电子科技大学 | 一种基于混合架构的高速高精度多通道并行采集系统 |
CN107994902A (zh) * | 2017-12-22 | 2018-05-04 | 北京工业大学 | 一种基于机器学习的双通道模数转换器子转换器间时间误差校正方法 |
CN115913231A (zh) * | 2023-01-06 | 2023-04-04 | 上海芯炽科技集团有限公司 | 一种tiadc的采样时间误差数字估计方法 |
CN116840549A (zh) * | 2023-07-03 | 2023-10-03 | 广州致远仪器有限公司 | 时间交织adc的误差确定方法、装置、设备以及存储介质 |
CN116840549B (zh) * | 2023-07-03 | 2024-03-22 | 广州致远仪器有限公司 | 时间交织adc的误差确定方法、装置、设备以及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105680858B (zh) | 2018-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105680858A (zh) | 一种估计tiadc并行采集系统时间偏移误差的方法 | |
US10291245B2 (en) | Device and method for correcting error estimation of analog-to-digital converter | |
TWI413780B (zh) | The phase correction circuit of the encoder signal | |
CN101674087B (zh) | 一种时间交替adc系统通道失配误差的获取方法 | |
CN102751991B (zh) | 基于δς调制原理的旋转变压器数字转换装置及方法 | |
CN102075464A (zh) | Tiadc系统通道误差的联合估计及实时校正方法 | |
Kang et al. | Phase difference correction method for phase and frequency in spectral analysis | |
US20120262318A1 (en) | System and Method for Background Calibration of Time Interleaved Analog to Digital Converters | |
US20060232460A1 (en) | Estimation circuit for time-interleaved adc and method thereof | |
CN103067006A (zh) | 一种针对时间交替模数转换系统时间误差的实时校正方法 | |
US9337853B2 (en) | Interleaved A/D converter | |
CN102761335A (zh) | 在双信道时间交错模拟至数字转换器中的取样时间和增益不匹配误差估计的无乘法器算法 | |
CN106209103B (zh) | 基于频谱分析的tiadc增益和时间误差的校正方法 | |
CN103647554A (zh) | 一种测量时间交替adc采集系统偏置与增益误差的方法 | |
WO2021114939A1 (zh) | 时间交织逐次逼近型模数转换器及其校准方法 | |
CN103941089A (zh) | 基于dft的正弦信号频率估计方法 | |
CN106027051A (zh) | 一种适用于流水线模数转换器的后台校准电路及校准方法 | |
CN105720983A (zh) | 用于时间交织模数转换系统的误差估计方法和装置 | |
US9748967B1 (en) | Periodic signal averaging with a time interleaving analog to digital converter | |
CN103957009A (zh) | 一种对压缩采样系统低通滤波器进行补偿的方法 | |
CN108809308A (zh) | 一种tiadc采集系统的时间误差估计及校正方法 | |
CN108445735B (zh) | 一种采用延迟链结构的层次式tdc的校正方法 | |
CN105141312A (zh) | 一种n通道时域交织模数转换器时钟偏差的数字后台校准算法 | |
CN111865331B (zh) | 一种mwc扩展系统传递矩阵的相位校准方法 | |
CN104122439B (zh) | 提高相位校正精度的电能表 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180807 Termination date: 20210114 |