本发明属于时间误差估算技术领域，更为具体地讲，涉及一种估计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。
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。