一种用于时间交织ADC采样时间失配的提取方法与流程

一种用于时间交织adc采样时间失配的提取方法
技术领域
1.本发明属于集成电路技术领域，具体涉及一种用于时间交织adc采样时间失配的提取方法。


背景技术：

2.时间交织模数转换器的原理如图1所示，将多个子adc(time
‑
interleaved，模数转换器)组合到一个adc系统中交错输出，从而成倍提高adc系统的采样率。然而，子adc通道间采样间隔的不匹配即时间失配会极大程度降低时间交织adc的动态性能。因此，对于采样时间失配(以下简称时间失配)的提取和补偿是必不可少的。
3.对于有n个通道的时间交织adc，定义t
s
为整个adc系统的采样周期，k表征第k个子adc通道(k＝1,...,n)，
m
为该通道adc的第
m
次输出(m＝0,1,..., ∞)。该情况下采样时间可以表示为t
m,k
＝(n
·
m k)
·
t
s
。该adc的输出可以表示为
4.y[n]＝y(t
m,k
)＝y((n
·
m k)
·
t
s
)
ꢀꢀꢀꢀ
(1)
[0005]
如果k通道存在时间失配
△
t
k
，实际采样时间为t
m,k
＝(n
·
m k)
·
t
s

△
t
k
，选第一通道为参考，有
△
t1＝0，此时adc输出表示为
[0006]
y[n]＝y(t
m,k
)＝y((n
·
m k)
·
t
s

△
t
k
)
ꢀꢀꢀꢀ
(2)
[0007]
存在量化噪声e
q
和热噪声e
th
，adc输出表示为
[0008]
y[n]＝y(t
m,k
) e
q
e
th
＝y((n
·
m k)
·
t
s

△
t
k
) e
q
e
th
ꢀꢀꢀꢀ
(3)
[0009]
对于时间交织adc的时间失配，常用的传统提取算法采用相关函数求得两个交织通道的输出乘积的均值：
[0010]
r(τ)＝e(y[n]
·
y[n 1])＝e(y(t
m
,k)
·
y(t
m,k
τ))
ꢀꢀꢀꢀ
(4)
[0011]
并用其表征adc的采样间隔，定义
τ
为相关函数的采样间隔。(4)中常用相邻通道进行计算，即理想情况下τ＝t
s
。如图2所示，第k
‑
1通道与第k通道的实际间隔值为τ＝t
s

△
t
k
‑△
t
k
‑1，此时相关函数r(τ)表示为r(t
s

△
t
k
‑△
t
k
‑1)。将r(t
s

△
t
k
‑△
t
k
‑1)简化为r
k
‑
1,k
表示，通过一阶泰勒展开：
[0012]
r
k
‑
1,k
＝r(t
s
) r'(t
s
)
·
△
t
k
‑
r'(t
s
)
·
△
t
k
‑1ꢀꢀꢀꢀ
(5)
[0013]
其中r'(t
s
)为τ＝t
s
时相关函数的斜率。此外，对于一个n通道时间交织adc，通常第一个通道作为参考通道，因此有
△
t1＝0。然后计算输出乘积均值之差，用于表征时间失配：
[0014][0015]
其中d
cor
‑
k
为自相关函数之差，用于表征时间失配的大小，从(6)可知理想情况下d
cor
‑
k
正比于
△
t
k
，完成d
cor
‑
k
即完成了
△
t
k
的提取。定义nc为计算(4)中均值的样本量，根据统计特性可知nc越大(6)中正比例关系受其他因素干扰越小。d
cor
‑
k
与
△
t
k
的正比例关系如
图3所示(图3模型中n＝4，f
s
＝1ghz，nc＝2
15
，rms jitter为200fs，且热噪声使输入信号的snr保持在74.5db左右)。该提取方法具有同时提取每个通道时间失配的优点，计算出d
cor
‑
k
可通过估计r'(t
s
)从而计算出
△
t
k
的准确值从而进行后续补偿，也可以通过调节模拟端可变延时线(variable delay line,vdl)将d
cor
‑
k
调整为0从而对时间失配
△
t
k
进行直接补偿。
[0016]
然而，上述提取方法具有以下两个主要缺陷：
[0017]
1.在样本量nc有限且输入频率接近n
f
·
f
s
2(n
f
＝0,1,2,...)时，(3)中量化噪声e
q
和热噪声e
th
的影响会等效为时间失配的提取误差，且nc越小该误差越大。如图4所示，噪声影响等效为时间失配的提取误差，从而影响精度。图4(a)中图4(b)中可以看到量化噪声e
q
和热噪声e
th
的影响较图3更明显，且更小的nc对应更大的误差。图4模型中n＝4，f
s
＝1ghz，rms jitter为200fs，且热噪声使输入信号的snr保持在74.5db左右。
[0018]
2.在输入频率接近n
f
·
f
s
2(n
f
＝0,1,2,...)时，在较大
△
t
k
范围d
cor
‑
k
与
△
t
k
的正比例关系中存在较大的非线性，如图5所示，该非线性效应可能导致时间失配的错误提取，从而使时间失配的补偿过程会收敛到错误的值。图5模型中n＝4，f
s
＝1ghz，nc＝2
14
，rms jitter为200fs，且热噪声使输入信号的snr保持在74.5db左右。


技术实现要素：

[0019]
针对上述传统校准算法受噪声以及非线性影响的问题，本发明提出了一种基于通道间隔拉长的自相关函数的时间交织adc时间失配提取方法。该方法的流程框图如图6所示，包括以下步骤：
[0020]
s1、时间交织adc将输入信号量化为数字信号；
[0021]
s2、定义采样间隔τ＝c
ci
·
t
s
，其中c
ci
为通道间隔系数，计算自相关函数：
[0022][0023]
s3、计算每个通道对应的自相关函数之差d
cor
‑
k,cci
：
[0024][0025]
其中，r'(c
ci
·
t
s
)为自相关函数的斜率；
[0026]
s4、基于以下正比例关系，提取或直接补偿各通道对应的采样时间失配
△
t
k
：
[0027]
[0028]
其中，e
k
为量化噪声e
q
和热噪声e
th
的影响在计算d
cor
‑
k,cci
时的等效误差，nl(
△
t
k
)为等效的非线性误差。
[0029]
计算相关函数时采样间隔采用τ＝c
ci
·
t
s
，取c
ci
＝n
c
·
n 1(n
c
＝1,2,3...)，此时(5)中r(τ)可表示为r(c
ci
·
τ)，r(c
ci
·
t
s

△
t
k
‑△
t
k
‑1)简化为r
k
‑
1,k,cci
。此时有
[0030]
r
k
‑
1,k,cci
＝r(c
ci
·
t
s
) r'(c
ci
·
t
s
)
·
△
t
k
‑
r'(c
ci
·
t
s
)
·
△
t
k
‑1ꢀꢀ
(10)
[0031]
根据(6)计算自相关函数之差d
cor
‑
k
，此时(6)中d
cor
‑
k
由相隔c
ci
个周期后的相邻通道之间的相关函数相加减得到，将通道间隔系数引入d
cor
‑
k
中得到d
cor
‑
k,cci
，从而得到式(8)，基于d
cor
‑
k,cci
与
△
t
k
的正比例关系可实现时间失配的提取，通过每个通道对应的d
cor,cci
的值调整至0后时间失配即被补偿完成，可通过数字滤波器、可变延时线等方法进行补偿。
[0032]
考虑误差与非线性效应，定义e
k
为量化噪声e
q
和热噪声e
th
的影响在计算d
cor
‑
k,cci
时的等效误差，nl(
△
t
k
)为等效的非线性误差。由于d
cor
‑
k,cci
总是相邻通道的自相关函数做加减计算求得，e
k
与nl(
△
t
k
)的值与通道间隔相关性极小，因此得到式(9)。
[0033]
未拉长通道之前c
ci
＝1，等效的时间失配提取误差e
△
t
为
[0034][0035]
在输入频率接近n
f
·
f
s
/2(n
f
＝0,1,2,...)时，r'(c
ci
·
t
s
)接近于0，等效误差不可忽略。对于单音信号、多音信号以及宽带信号等通信系统中的信号，在r'(c
ci
·
t
s
)接近于0时都具有以下特性
[0036]
r'(c
ci
·
t
s
)≈c
ci
·
r'(t
s
)
ꢀꢀ
(12)
[0037]
因此拉长间隔后有
[0038][0039]
此时等效的提取误差被抑制，从而实现高精度的时间失配提取。
[0040]
如图7(a)所示，该方法可有效减小输入频率接近n
f
·
f
s
/2(n
f
＝0,1,2,...)时量化噪声e
q
和热噪声e
th
的影响等效的提取误差(c
ci
＝1为背景技术提取效果，c
ci
＝5为拉长间隔为5个周期后的提取效果，可以发现拉长间隔后量化噪声e
q
和热噪声e
th
的影响等效的提取误差变小)。同时，如图7(b)所示，该方法可以明显抑制较大
△
t
k
范围内d
cor
‑
k
与
△
t
k
的正比例关系中存在的非线性，此时在该
△
t
k
范围d
cor
‑
k
不再与x轴相交第二次
[0041]
本发明的有益效果为：
[0042]
本发明提出的方法可用于在保证与背景技术相同的样本量情况下，通过抑制非线性和噪声影响从而提高提取精度；或者保证与背景技术相同提取精度的同时减小计算所需的样本量，从而对时间失配进行快速提取，实现时间失配校准的快速收敛。
附图说明
[0043]
图1为时间交织adc原理框图
[0044]
图2为采样时间失配在n通道时间交织adc中的影响
[0045]
图3为理想情况下d
cor
‑
k
与
△
t
k
的正比例关系图
[0046]
图4为输入频率分别为(a)以及(b)时在较小
△
t
k
范围内
d
cor
‑
k
与
△
t
k
的正比例关系在nc＝2
12
和nc＝2
16
时的对比图
[0047]
图5输入频率为时在较大
△
t
k
范围内d
cor
‑
k
与
△
t
k
存在非线性关系的说明图
[0048]
图6为本发明提出的基于通道间隔拉长的相关函数的时间交织adc时间失配提取方法的流程框图
[0049]
图7为采用本发明的提取方法与传统方法进行提取时精度提高的效果说明图，(a)背景技术方法，(b)为本发明方法。
具体实施方式
[0050]
在发明内容部分已经对本发明的技术方案进行了详细描述，下面结合附图和仿真示例说明本发明的实用性。
[0051]
以四通道(n＝4)时间交织adc为例，且只考虑采样时间失配，假设其他失配已被校准，详细步骤如下：
[0052]
1.对于有4个通道的时间交织adc，输入正弦信号定义t
s
为整个adc系统的采样周期，k表征第k个子adc通道(k＝1,...,4)，
m
为该通道adc的第
m
次输出(m＝0,1,..., ∞)。该情况下采样时间可以表示为t
m,k
＝(4m k)
·
t
s
。选第一通道为参考，有
△
t1＝0，此时adc输出表示为
[0053]
y[n]＝y(t
m,k
)＝y((4m k)
·
t
s

△
t
k
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0054]
2.求出nc＝2
12
且基于c
ci
＝5拉长通道间隔的通道间自相关函数值：
[0055][0056]
3.求出基于c
ci
＝5拉长通道间隔的相关函数之差：
[0057][0058]
此时，d
cor
‑
k，cci＝5
和采样时间失配
△
t
k
在(16)中呈正比例关系，根据(16)，可通过估计r'(5
·
t
s
)从而计算出
△
t
k
的准确值从而进行补偿，也可以直接通过模拟端vdl的补偿将d
cor
‑
k，cci＝5
调整为0从而补偿
△
t
k
。