CN108809308B - 一种tiadc采集系统的时间误差估计及校正方法 - Google Patents

一种tiadc采集系统的时间误差估计及校正方法 Download PDF

Info

Publication number
CN108809308B
CN108809308B CN201810600898.4A CN201810600898A CN108809308B CN 108809308 B CN108809308 B CN 108809308B CN 201810600898 A CN201810600898 A CN 201810600898A CN 108809308 B CN108809308 B CN 108809308B
Authority
CN
China
Prior art keywords
adc
value
time error
tiadc
phase
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.)
Active
Application number
CN201810600898.4A
Other languages
English (en)
Other versions
CN108809308A (zh
Inventor
杨扩军
叶芃
潘志翔
张沁川
黄武煌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810600898.4A priority Critical patent/CN108809308B/zh
Publication of CN108809308A publication Critical patent/CN108809308A/zh
Application granted granted Critical
Publication of CN108809308B publication Critical patent/CN108809308B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/10Calibration or testing
    • H03M1/1009Calibration
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/12Analogue/digital converters
    • H03M1/124Sampling or signal conditioning arrangements specially adapted for A/D converters
    • H03M1/1245Details of sampling arrangements or methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种TIADC采集系统的时间误差估计及校正方法,通过给TIADC系统一个高频标准正弦激励获得各ADC采样数据,再对各ADC采样数据做无频谱泄漏的FFT获得其频谱,并利用频域折叠理论精确获得各ADC相位谱值,然后利用各ADC相位谱值获得时间误差估计值,并对估计值进行修正以获得时间误差校正值,最后根据时间误差校正值和ADC内部相位调节单元对TIADC时间误差进行校正。

Description

一种TIADC采集系统的时间误差估计及校正方法
技术领域
本发明属于时域测试技术领域,更为具体地讲,涉及一种TIADC采集系统的时间误差估计及校正方法。
背景技术
随着科技的快速发展,以GHz频段为代表,超高频信号得到了越来越多的应用。尤其在通信,航空航天,空间探测等领域,超高频信号的应用迅速增长。为了能够准确识别和测量这些信号,对数据采集系统的采样率提出了更高的要求。现有的单片ADC由于材料、工艺的限制已经无法满足对这些信号的采样需求,因此并行采集技术作为一种对任意信号的实时采集技术得到了广泛的研究。
多片ADC并行交替采样系统结构如图1所示。N个通道的ADC以并行方式交替采样,第i通道的ADC的编号为ADCi,其采样时钟为SCLKi,i=0,1,...,M-1。设系统的总采样率为fs,根据TIADC的理论,单片ADC的采样率为fsingle=fs/N,且采样时钟的相位关系为:
Figure BDA0001693187020000011
在采样时钟的驱动下,每片ADC对信号调理通道处理后的信号x(t)进行采样,并将采样数据xi(n)发送给数据拼合模块进行数据重组,获得原始信号的重构y(n)。
这种并行结构虽然提升了最大采样率,但在实际应用中会由于不同通道间器件参数、布局的不一致性产生偏置、增益和时间三种通道失配误差,其中增益和时间误差与输入信号的频率f0有关。失配误差的引入大大降低了系统的整体性能,使高频信号的识别和测量受到了很大限制。为了使系统能够尽可能真实地还原出原始高频信号,必须要对通道失配误差进行估计和校正。然而现有的误差估计方法无法在单通道为欠采样的模式下估计时间失配误差,因此设计一种新的时间误差估计和校正方法就成为高速并行采集系统中亟需解决的问题。
发明内容
本发明的目的在于克服现有技术的不足,提供一种TIADC采集系统的时间误差估计及校正方法,模拟延时单元延迟值进行时间误差的估计和校正。
为实现上述发明目的,本发明一种TIADC采集系统的时间误差估计及校正方法,其特征在于,包括以下步骤:
(1)、设TIADC采集系统的ADC个数为M,单通道ADC的采样率为fsingle
(2)、设输入频率为f0的标准正弦信号为x(t),将x(t)输入至TIADC采集系统,并对x(t)进行采样,得到各ADC的采样数据yk[n],k=0,1,…,M-1;
(3)、求各ADC采样数据yk[n]的均值ak,再从各ADC采样数据yk[n]中减去均值ak,从而去除直流偏置,得到
Figure BDA0001693187020000021
Figure BDA0001693187020000022
(4)、对
Figure BDA0001693187020000023
做无频谱泄露的FFT变换,然后取频点f0'=f0-μfsingle处的相位谱值∠Xk(jω),其中,μ为常数;
(5)、依次计算各∠Xk(jω)-∠X0(jω)的值,k=1,2,…,M-1,若其值小于0,则在∠Xk(jω)上加2π,反之则保持∠Xk(jω)不变,从而得到调整后的相位谱值ΔΦk
(6)、计算各ADC的时间误差估计值δk
Figure BDA0001693187020000024
其中,
Figure BDA0001693187020000025
Figure BDA0001693187020000026
表示向下取整;Ts为TIADC系统的采样间隔,Ts=1/Mfsingle
(7)、对各ADC的时间误差估计值δk进行修正;
Figure BDA0001693187020000027
则将δk修正为:
Figure BDA0001693187020000028
Figure BDA0001693187020000029
则将δk修正为:
Figure BDA00016931870200000210
Figure BDA00016931870200000211
则保持δk不变,δk'=δk
(8)、计算各ADC的时间误差校正值tk:tk=-δk';
(9)、根据各ADC的时间误差校正值tk对各ADC进行校正;
(9.1)、设各ADC内部相位调节寄存器的调节步进延时为S,则第k个ADC的寄存器调节值为Rk=[tk/S],其中,[·]为四舍五入操作;
(9.2)、设第k个ADCk的相位调节寄存器的初值为R0k,第0个ADC相位调节寄存器的初值为R00,计算调整后第k个ADCk的相位调节寄存器值R'k,R'k=R0k+Rk
(9.3)、判断R'k是否超过ADCk内部相位寄存器的调节范围:
若R'k大于可调范围的最大值,则减小第0通道相位调节寄存器初值R00,再返回步骤(2);
若R'k小于可调范围的最小值,则增大第0通道相位调节寄存器初值R00,再返回步骤(2);若R'k在可调范围之内,则调节R0k至R'k,完成校正。
本发明的发明目的是这样实现的:
本发明一种TIADC采集系统的时间误差估计及校正方法,通过给TIADC系统一个高频标准正弦激励获得各ADC采样数据,再对各ADC采样数据做无频谱泄漏的FFT获得其频谱,并利用频域折叠理论精确获得各ADC相位谱值,然后利用各ADC相位谱值获得时间误差估计值,并对估计值进行修正以获得时间误差校正值,最后根据时间误差校正值和ADC内部相位调节单元对TIADC时间误差进行校正。
同时,本发明一种TIADC采集系统的时间误差估计及校正方法还具有以下有益效果:
(1)、通过使用频域折叠理论,解除了原有TIADC时间误差估计方法只能适用于低频信号的限制,将时间误差估计应用于单ADC为欠采样下的高频场合。
(2)、给出了时间误差估计值的具体修正过程,解决了由于随机环境误差导致时间误差计算错误的情况,增加了校正算法的稳定性和正确性。
附图说明
图1是本发明一种TIADC采集系统的时间误差估计及校正方法原理图;
图2是图1所示TIADC采集系统的一种等效结构模型示意图;
图3是输入校正信号频谱与TIADC系统采样率关系示意图;
图4是第一种情况下输入校正信号频谱折叠后与单ADC采样率关系示意图;
图5是第二种情况下输入校正信号频谱折叠后与单ADC采样率关系示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明一种TIADC采集系统的时间误差估计及校正方法原理图。
在本实施例中,如图1所示,本发明一种TIADC采集系统的时间误差估计及校正方法,包括以下步骤:
S1、设TIADC采集系统的ADC个数为M=16,单通道ADC的采样率为fsingle;如图2所示,每个通道等效为一个代表通道频率响应的滤波器和ADC组成,而时间误差由通道的频响滤波器Hi(jω)表示。
S2、设输入频率为f0的标准正弦信号为x(t),将x(t)输入至TIADC采集系统,并对x(t)进行采样,得到各ADC的采样数据yk[n],k=0,1,…,15;
S3、求各ADC采样数据yk[n]的均值ak,再从各ADC采样数据yk[n]中减去均值ak,从而去除直流偏置,得到
Figure BDA0001693187020000041
Figure BDA0001693187020000042
S4、对
Figure BDA0001693187020000043
做无频谱泄露的FFT变换,然后取频点f0'=f0-μfsingle处的相位谱值∠Xk(jω),其中,μ为常数,且满足:μfsingle<f0<(μ+1)fsingle
在本实施例中,设输入的校正信号x(t)为高频正弦信号,即:
x(t)=sin(2πf0t+θ)
经过频响滤波器后,输出为:
Figure BDA0001693187020000046
其中,f0为x(t)的频率,
Figure BDA0001693187020000047
为时间误差系数,其与时间误差δk的关系为
Figure BDA0001693187020000048
θ为x(t)的相位;
将其改写为:
Figure BDA0001693187020000044
Figure BDA0001693187020000045
则xk(n)的连续傅里叶变换为
Figure BDA0001693187020000051
其中δ表示dirac函数。由dirac函数的性质,有:
Figure BDA0001693187020000052
因此有:
ω=2πf0时,Xk(jω)的相位为
Figure BDA0001693187020000058
ω=-2πf0时,Xk(jω)的相位为
Figure BDA0001693187020000059
即f0和-f0处频谱的相位相反。
由于采样是对原始信号频谱的周期延拓,故分析时只考虑频谱范围为0~fsingle内即可。
原始信号频谱如图3所示,其中正频率f0处谱线用实线表示,负频率-f0处谱线用虚线表示。根据f0和fsingle的关系,分两种情况讨论:
Figure BDA0001693187020000053
如图4所示,此时有:
Figure BDA0001693187020000054
即采样后信号频点f0'=f0-μfsingle处谱线与原信号f0处谱线对应。
Figure BDA0001693187020000055
如图5所示,此时有:
Figure BDA0001693187020000056
即采样后信号频点f0'=f0-μfsingle处谱线与原信号f0处谱线对应。实际应用中FFT常取频率范围
Figure BDA0001693187020000057
此时若使用(μ+1)fsingle-f0频点处值,则根据前述讨论,其相位与原始相位相反,为了简化分析,使用f0'=f0-μfsingle处谱线数据;
S5、以第0通道为参考通道,其他通道以第0通道为基准进行校正,第0通道相对自身无通道失配误差,有δ0=0,依次计算各∠Xk(jω)-∠X0(jω)的值,k=1,2,…,15,若其值小于0,则在∠Xk(jω)上加2π,反之则保持∠Xk(jω)不变,从而得到调整后的相位谱值ΔΦk
S6、计算各ADC的时间误差估计值δk
Figure BDA0001693187020000061
其中,
Figure BDA0001693187020000062
Figure BDA0001693187020000063
表示向下取整;Ts为TIADC系统的采样间隔,Ts=1/Mfsingle
S7、对各ADC的时间误差估计值δk进行修正;
Figure BDA0001693187020000064
则将δk修正为:
Figure BDA0001693187020000065
Figure BDA0001693187020000066
则将δk修正为:
Figure BDA0001693187020000067
Figure BDA0001693187020000068
则保持δk不变,δk'=δk
S8、计算各ADC的时间误差校正值tk:tk=-δk';
S9、根据各ADC的时间误差校正值tk对各ADC进行校正;
S9.1、设各ADC内部相位调节寄存器的调节步进延时为S,则第k个ADC的寄存器调节值为Rk=[tk/S],其中,[·]为四舍五入操作;
S9.2、设第k个ADCk的相位调节寄存器的初值为R0k,第0个ADC相位调节寄存器的初值为R00,计算调整后第k个ADCk的相位调节寄存器值R'k,R'k=R0k+Rk
S9.3、判断R'k是否超过ADCk内部相位寄存器的调节范围:
若R'k大于可调范围的最大值,则减小第0通道相位调节寄存器初值R00,此时S5中的∠Xk(jω)-∠X0(jω)值发生改变,再返回步骤S2;
若R'k小于可调范围的最小值,则增大第0通道相位调节寄存器初值R00,此时S5中的∠Xk(jω)-∠X0(jω)值发生改变,再返回步骤S2;
若R'k在可调范围之内,则调节R0k至R'k,完成校正。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。

Claims (2)

1.一种TIADC采集系统的时间误差估计及校正方法,其特征在于,包括以下步骤:
(1)、设TIADC采集系统的ADC个数为M,单通道ADC的采样率为fsingle
(2)、设输入频率为f0的标准正弦信号为x(t),将x(t)输入至TIADC采集系统,并对x(t)进行采样,得到各ADC的采样数据yk[n],k=0,1,…,M-1;
(3)、求各ADC采样数据yk[n]的均值ak,再从各ADC采样数据yk[n]中减去均值ak,从而去除直流偏置,得到
Figure FDA0003172998870000011
Figure FDA0003172998870000012
(4)、对
Figure FDA0003172998870000013
做无频谱泄露的FFT变换,然后取频点f0'=f0-μfsingle处的相位谱值∠Xk(jω),其中,μ为常数;
(5)、依次计算各∠Xk(jω)-∠X0(jω)的值,k=1,2,…,M-1,若其值小于0,则在∠Xk(jω)上加2π,反之则保持∠Xk(jω)不变,从而得到调整后的相位谱值ΔΦk
(6)、计算各ADC的时间误差估计值δk
Figure FDA0003172998870000014
其中,
Figure FDA0003172998870000015
Figure FDA0003172998870000016
表示向下取整;Ts为TIADC系统的采样间隔,Ts=1/Mfsingle
(7)、对各ADC的时间误差估计值δk进行修正;
Figure FDA0003172998870000017
则将δk修正为:
Figure FDA0003172998870000018
Figure FDA0003172998870000019
则将δk修正为:
Figure FDA00031729988700000110
Figure FDA00031729988700000111
则保持δk不变,δk'=δk
(8)、计算各ADC的时间误差校正值tk:tk=-δk';
(9)、根据各ADC的时间误差校正值tk对各ADC进行校正;
(9.1)、设各ADC内部相位调节寄存器的调节步进延时为S,则第k个ADC的寄存器调节值为Rk=[tk/S],其中,[·]为四舍五入操作;
(9.2)、设第k个ADCk的相位调节寄存器的初值为R0k,第0个ADC相位调节寄存器的初值为R00,计算调整后第k个ADCk的相位调节寄存器值R'k,R'k=R0k+Rk
(9.3)、判断R'k是否超过ADCk内部相位寄存器的调节范围:
若R'k大于可调范围的最大值,则减小第0通道相位调节寄存器初值R00,再返回步骤(2);
若R'k小于可调范围的最小值,则增大第0通道相位调节寄存器初值R00,再返回步骤(2);若R'k在可调范围之内,则调节R0k至R'k,完成校正。
2.根据权利要求1所述 的一种高速并行采集系统的时间误差估计及校正方法,其特征在于,所述常数μ满足:μfsingle<f0<(μ+1)fsingle
CN201810600898.4A 2018-06-12 2018-06-12 一种tiadc采集系统的时间误差估计及校正方法 Active CN108809308B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810600898.4A CN108809308B (zh) 2018-06-12 2018-06-12 一种tiadc采集系统的时间误差估计及校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810600898.4A CN108809308B (zh) 2018-06-12 2018-06-12 一种tiadc采集系统的时间误差估计及校正方法

Publications (2)

Publication Number Publication Date
CN108809308A CN108809308A (zh) 2018-11-13
CN108809308B true CN108809308B (zh) 2021-10-19

Family

ID=64085259

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810600898.4A Active CN108809308B (zh) 2018-06-12 2018-06-12 一种tiadc采集系统的时间误差估计及校正方法

Country Status (1)

Country Link
CN (1) CN108809308B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109379080A (zh) * 2018-09-21 2019-02-22 电子科技大学 用于时间交替采样的时间误差自适应消除方法
EP3872994A4 (en) * 2018-11-16 2021-11-03 Huawei Technologies Co., Ltd. ERROR CORRECTION PROCESS AND ANALOG / DIGITAL CONVERTER INTERLACED IN TIME
CN110557122B (zh) * 2019-09-25 2022-04-19 电子科技大学 一种tiadc系统频响非一致性误差的校正方法
CN113114241B (zh) * 2021-03-10 2022-04-19 电子科技大学 一种时间交替架构采集系统中频响失配误差的校正方法
CN114355284B (zh) * 2022-01-04 2023-05-05 电子科技大学 一种利用频谱主分量的时差估计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103647554A (zh) * 2013-12-16 2014-03-19 电子科技大学 一种测量时间交替adc采集系统偏置与增益误差的方法
CN104993827A (zh) * 2015-07-08 2015-10-21 中国电子科技集团公司第二十四研究所 模数转换器误差估计校正的装置及其方法
CN105656485A (zh) * 2015-12-30 2016-06-08 北京坤驰科技有限公司 一种多通道时间交错adc测量校准方法和装置
CN105933005A (zh) * 2016-04-19 2016-09-07 中国电子科技集团公司第四十研究所 一种基于等效采样的时域交替模数转换器失配校准方法
CN106209103A (zh) * 2016-07-29 2016-12-07 电子科技大学 基于频谱分析的tiadc增益和时间误差的校正方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI282216B (en) * 2005-04-13 2007-06-01 Realtek Semiconductor Corp Correlation circuit for time-interleaved ADC and method thereof
US8890739B2 (en) * 2012-12-05 2014-11-18 Crest Semiconductors, Inc. Time interleaving analog-to-digital converter
US9401726B2 (en) * 2014-11-26 2016-07-26 Silicon Laboratories Inc. Background calibration of time-interleaved analog-to-digital converters

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103647554A (zh) * 2013-12-16 2014-03-19 电子科技大学 一种测量时间交替adc采集系统偏置与增益误差的方法
CN104993827A (zh) * 2015-07-08 2015-10-21 中国电子科技集团公司第二十四研究所 模数转换器误差估计校正的装置及其方法
CN105656485A (zh) * 2015-12-30 2016-06-08 北京坤驰科技有限公司 一种多通道时间交错adc测量校准方法和装置
CN105933005A (zh) * 2016-04-19 2016-09-07 中国电子科技集团公司第四十研究所 一种基于等效采样的时域交替模数转换器失配校准方法
CN106209103A (zh) * 2016-07-29 2016-12-07 电子科技大学 基于频谱分析的tiadc增益和时间误差的校正方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A subband FFT-based method for static errors compensation in Time-Interleaved ADCs;Maher Jridi;《 2011 IEEE 54th International Midwest Symposium on Circuits and Systems (MWSCAS)》;20110923;第1-4页 *
Blind timing skew estimation based on spectra sparsity and all phase FFT for time-interleaved ADCs;S.J. Liu;《2015 IEEE International Conference on Digital Signal Processing (DSP)》;20150910;第926-930页 *
数理统计和频谱分析的TIADC误差校正方法;杨扩军;《电子科技大学学报》;20180131;第43-50页 *

Also Published As

Publication number Publication date
CN108809308A (zh) 2018-11-13

Similar Documents

Publication Publication Date Title
CN108809308B (zh) 一种tiadc采集系统的时间误差估计及校正方法
CN108494402B (zh) 一种基于正弦拟合的tiadc系统误差估计和补偿方法
CA2929382C (en) Method of removing decaying dc component from power system fault signal
CN103647554A (zh) 一种测量时间交替adc采集系统偏置与增益误差的方法
CN107645770B (zh) 一种相位校准方法及装置
CN103957009A (zh) 一种对压缩采样系统低通滤波器进行补偿的方法
CN113965217B (zh) 一种双通道、单通道s/c/x三波段宽带单比特数字测频接收机
CN112162152B (zh) 基于相位直线拟合的正弦波相参脉冲串信号频率估计方法
CN105259410A (zh) 一种强噪声干扰下的欠采样波形的频率估计方法及其装置
CN103823177A (zh) 基于窗函数设计的滤波器性能检测方法和系统
US9748967B1 (en) Periodic signal averaging with a time interleaving analog to digital converter
CN107425853B (zh) 基于fft的双通道tiadc系统失配误差盲校正方法
CN106777505A (zh) 基于频偏识别的欠采样信号的鲁棒的频率估计方法及装置
CN111611686A (zh) 一种通信信号时频域的检测方法
Bai et al. Blind calibration method for two-channel time-interleaved analog-to-digital converters based on FFT
CN107132511B (zh) 一种精确的雷达线性调频源预失真方法
CN115597641A (zh) 一种光纤法布里-珀罗传感器的动态腔长解调方法及系统
CN111555995B (zh) 一种宽带采集系统的相频响应测量及补偿方法
CN110297199B (zh) 一种基于全相位fft的铯光泵磁力仪频率测量方法及系统
CN114563769A (zh) 数字相控阵接收通道相位非线性的测量方法及装置
CN107124183B (zh) 双通道tiadc系统失配误差盲校正方法
CN107576842B (zh) 一种宽频带同步采样方法
WO2016062033A1 (zh) 频率补偿处理方法及装置
US10298249B1 (en) Time error and gain offset estimation in interleaved analog-to-digital converters
JPH04328471A (ja) ディジタル信号測定装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant