本发明属于信号处理
技术领域：
，具体涉及一种双通道多指数衰减正弦信号欠采样方法。
背景技术：
：如何高效地测量多指数衰减正弦(multipleexponentiallydampedsinusoids，meds)信号的参数是很多实际应用中多次出现的问题，如语音分析，核磁共振成像等其他应用。为了能够方便快捷地测量meds信号的参数，许多研究者提出了各种各样的方法。umesh等学者提出了快速最大似然(fml)估计来进行参数估计。fotinea等学者针对谱估计提出了一种状态空间方法。但是，以上方法均建立在著名的奈奎斯特采样理论的基础之上，该理论要求采样速率是信号带宽的两倍。面对宽带信号，由于需要大量的采样值以及繁杂的参数测量，这些方法就需要非常高的采样速率。考虑如下形式的连续时间meds信号z(t)：其中，ck,k＝1,2,…k(分量个数k是已知的)是复振幅，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率。假设以速率fs对信号z(t)进行采样，则采样值表示为：其中，m∈z 。利用零化滤波器方法，我们可以从2k个z[m]的样本中估算出参数ck和uk。设为估计值，由于三角函数的周期性，欠采样会导致频谱混叠现象。近来，由于许多应用中meds信号的参数测量问题多次出现，很多研究者对该问题给予了不少关注，同时也提出了很多测量方法。众所周知，如果采用欠奈奎斯特对待测信号的采样值进行参数测量，我们需要解决的关键问题就是频谱混叠。r.g.baraniuk等学者提出了一种基于压缩感知(compressedsensing，cs)理论的测量方法，但是cs要求对待测信号进行随机采样，该要求对于硬件实现难度较大。s.qin等人提出了一个具有互质速率的双通道欠奈奎斯特采样和参数测量方案，但是，该方案被s.huang等人证明了在某些情况下不能获得唯一解。s.huang等人提出了一种具有成对互质速率的三通道欠奈奎斯特采样和参数测量方案，虽然该方案能够解决频谱混叠问题，但是需要很多的采样样本才能获得唯一解。可以看出这些系统并非为meds信号参数测量而专门设计的，而且这些系统需要知道频率个数。然而，现有方法大多数是基于奈奎斯特采样理论的，需要大量的采样值，因此，如何在通过较少的采样值测量meds信号的参数的同时解决欠奈奎斯特采样带来的频谱混叠问题仍然是一个关键难题。技术实现要素：为了解决meds信号的欠采样引发的频谱混叠和参数测量问题，为了克服已有技术的不足，本发明提供了一种双通道多指数衰减正弦信号欠采样方法，可以从少至3k个样本中测量出已知分量个数为k的meds信号。本发明解决其技术问题所采用的技术方案是：一种双通道多指数衰减正弦信号欠采样方法，包括以下步骤：步骤一，生成待测meds信号：用上位机labview软件控制nipxie-6363任意波形发生器生成待测信号，待测meds信号形式表示为如下形式：其中，k为分量个数，并且k为已知量，ck是复振幅，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率；步骤二，模拟预处理过程：采用模拟预处理模块将待测信号分流至两个通道中，模拟预处理模块包括功分器、噪声滤波器和信号放大器；步骤三，低速adc采样过程：采用上位机labview软件控制nipxi-6251模块板卡，使用两个ai通道分别对模拟信号进行采样，其中一个通道的采样时钟与另一个通道的采样时钟跳变频率相同，但是上升沿相差一个较小的时间差te，且满足te≤1/fmax，fmax表示信号的最大频率，所以该部分称为时钟交错采样系统；步骤四，运行零化滤波器谱估计算法：在上位机中，首先对采集到的样本进行离散傅立叶变换，从而获取fourier系数；然后运行谱估计算法，从获取到的fourier系数中估计出meds信号的参数步骤五，运行复频率联合估计算法，通过参数采样速率fs，主、次通道采样时钟上升沿相差的时间间隔te以及n′≥k个连续采样值ze[n′]联合估计出阻尼因子rk以及极点频率fk的估计值算法的输入为参数采样速率fs、交错时间te以及n′≥k个连续采样值ze[n′]，输出为阻尼因子rk和极点频率fk的估计值进一步，所述步骤三中，信号采样过程的数学模型描述如下：3.1，在主采样通道中，信号z(t)根据时钟clk1被均匀采样，并且采样可表示为如下形式：其中，k为分量个数，并且k为已知量，ck是复振幅，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率，m∈z ；3.2，在次采样通道中，信号z(t)根据时钟clk2被均匀采样，并且采样由下式给出：其中，是复振幅ck的估计值，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率，ts是采样周期，te是主、次通道采样时钟上升沿相差的时间间隔，n′∈z 。再进一步，所述步骤五的算法过程如下：5.1，估计阻尼因子rk。阻尼因子rk可以很容易地从中估算出来，估算过程如下：由与sk＝rk j2πfk得则得阻尼因子rk的估计值如下式：同时得到极点频率fk的估计值如下式：其中，∠(·)表示(·)的幅角主值，且0≤∠(·)＜2π；5.2，由公式(3)和公式(4)得其中，是已知的，是未知的；5.3，求解bk的估计值：利用n′≥k个连续采样值ze[n′](1≤n′≤k)构建一个n′×1的向量，该向量形式如下：ze＝[xe[0],xe[1],…,xe[n′-1]]t(8)利用的幅角主值创建一个k×n′的范德蒙矩阵，该矩阵形式如下：利用和的幅角主值创建一个n′×k的对角矩阵通过求解得到系数矩阵5.4，求解极点频率fk的估计值，由于三角函数的周期性，估计的归一化频率与的参数相差2πmk，即：其中，∠(·)是(·)的幅角主值，0≤∠(·)＜2π，mk∈z,k∈{1,2,…,k}，因为fmax表示信号的最大频率，得到如下不等式：进一步化简公式(12)得如下不等式：记公式(13)可转换为：因为q∈z ，则有如果设置主、次通道采样时钟上升沿相差的时间间隔则有：与公式(14)相结合，则有：0≤mkfste＜1(16)所以对于mk有唯一解，其解如下式：其中，极点频率fk的估计值由唯一确定。本发明设计了一种双通道时钟交错欠采样系统对生成的meds信号进行采样，然后通过零化滤波器法对该信号进行谱估计，估算出参数其次通过复频率联合估计算法对该信号进行频率估计，从而完成参数测量进行重构，最后在上位机中显示重构效果。本发明的有益效果主要表现在：采用本发明方法，可以从少至3k个样本中测量出已知分量个数为k的meds信号；同时，本发明还提供了硬件原型来实现所提系统，通过仿真和实验结果验证了方案的有效性。附图说明图1是双通道meds信号欠采样硬件系统的框图。具体实施方式下面结合附图对本发明作进一步描述。参照图1，一种双通道多指数衰减正弦信号欠采样方法，包括以下步骤：步骤一，生成待测meds信号：用上位机labview软件控制nipxie-6363任意波形发生器生成待测信号，待测meds信号形式表示为如下形式：其中，k为分量个数，并且k为已知量，ck是复振幅，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率；步骤二，模拟预处理过程：采用模拟预处理模块将待测信号分流至两个通道中，模拟预处理模块包括功分器、噪声滤波器和信号放大器；步骤三，低速adc采样过程：采用上位机labview软件控制nipxi-6251模块板卡，使用两个ai通道分别对模拟信号进行采样，其中一个通道的采样时钟与另一个通道的采样时钟跳变频率相同，但是上升沿相差一个较小的时间差te，且满足te≤1/fmax，fmax表示信号的最大频率，所以该部分称为时钟交错采样系统，信号采样过程的数学模型描述如下：3.1，在主采样通道中，信号z(t)根据时钟clk1被均匀采样，并且采样可表示为如下形式：其中，k为分量个数，并且k为已知量，ck是复振幅，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率，m∈z ；3.2，在次采样通道中，信号z(t)根据时钟clk2被均匀采样，并且采样由下式给出：其中，是复振幅ck的估计值，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率，ts是采样周期，te是主、次通道采样时钟上升沿相差的时间间隔，n′∈z ；步骤四，运行零化滤波器谱估计算法：在上位机中，首先对采集到的样本进行离散傅立叶变换，从而获取fourier系数；然后运行谱估计算法，从获取到的fourier系数中估计出meds信号的参数步骤五，运行复频率联合估计算法。通过参数采样速率fs，主、次通道采样时钟上升沿相差的时间间隔te以及n′≥k个连续采样值ze[n′]联合估计出阻尼因子rk以及极点频率fk的估计值本过程是关键的一步，算法的输入为参数采样速率fs、交错时间te以及n′≥k个连续采样值ze[n′]，输出为阻尼因子rk和极点频率fk的估计值算法过程如下：5.1，估计阻尼因子rk。阻尼因子rk可以很容易地从中估算出来，估算过程如下：由与sk＝rk j2πfk得则得阻尼因子rk的估计值如下式：同时得到极点频率fk的估计值如下式：其中，∠(·)表示(·)的幅角主值，且0≤∠(·)＜2π；5.2，由公式(3)和公式(4)得其中，是已知的，是未知的；5.3，求解bk的估计值：利用n′≥k个连续采样值ze[n′](1≤n′≤k)构建一个n′×1的向量，该向量形式如下：ze＝[xe[0],xe[1],…,xe[n′-1]]t(8)利用的幅角主值创建一个k×n′的范德蒙矩阵，该矩阵形式如下：利用和的幅角主值创建一个n′×k的对角矩阵通过求解得到系数矩阵5.4，求解极点频率fk的估计值，由于三角函数的周期性，估计的归一化频率与的参数相差2πmk，即：其中，∠(·)是(·)的幅角主值，0≤∠(·)＜2π，mk∈z,k∈{1,2,…,k}，因为fmax表示信号的最大频率，得到如下不等式：进一步化简公式(12)得如下不等式：记公式(13)可转换为：0≤mk≤(q-1)(14)因为q∈z ，则有如果设置主、次通道采样时钟上升沿相差的时间间隔则有：与公式(14)相结合，则有：0≤mkfste＜1(16)所以对于mk有唯一解，其解如下式：其中，极点频率fk的估计值由唯一确定。本实施例的双通道meds信号欠采样硬件系统。该系统由两个具有相同采样率的欠奈奎斯特采样通道组成，其中一个采样通道的采样时钟与另一个采样通道的采样时钟具有相同的频率，但是上升沿错开一个很小的时段。实验对比：为了验证本发明方法的优势，我们将本发明方法与s.qin提出的延时采样系统和s.huang提出的三通道互质采样系统进行了比较，比较结果如表1所示。方法频谱混叠采样数qin所提系统适用4khuang所提系统适用需要大量采样数本发明方法适用3k表1可以看出，本发明所提出的系统仅用3k个样本就可以解决频谱混叠，所需样本个数比其他采样系统所需要的采样数少得多。本发明还进行了其他的相关实验，从而检测了系统的有效性。对于meds信号本发明仍然能够有效地完成参数测量。假设meds信号x(t)含有k个频率分量且k已知，具有最大频率fmax＝10khz，复振幅ck在(0,1 i)内随机取值，阻尼因子rk在(0,1)上随机分布，并且极点频率fk在(0,fmax)上分布。为了方便比较，我们尽可能的设置相似的系统参数：qin的系统和本发明方法的信道采样率设置为fs＝fmax/4＝2.5khz；huang的系统的三条信道采样率分别设置为fs＝2,2.5,3.3khz。首先，在无噪声的情况下，我们用少量样本来验证本发明所用采样系统的有效性。设置待测meds信号x(t)的分量个数k＝5。假设k是已知的，并且从qin的系统中选取4k＝20个采样值，同时，从huang的系统中选取188个采样值恢复信号。在本发明方法中，只要3k＝15个采样值(k是已知的)就能恢复信号。无噪声情况下的恢复结果对比如表2所示。表2从表2可以看出，三个都系统能够估计所有的极点频率，但是qin和huang的系统不能估计复振幅ck和阻尼因子rk，而本发明能够估计这两个参数，更重要的是，本发明方法的采样值个数远少于qin和huang的系统所需要的采样值。当前第1页1 2 3
技术特征：1.一种双通道多指数衰减正弦信号欠采样方法，其特征在于，所述方法包括以下步骤：
步骤一，生成待测meds信号：用上位机labview软件控制nipxie-6363任意波形发生器生成待测信号，待测meds信号形式表示为如下形式：
其中，k为分量个数，并且k为已知量，ck是复振幅，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率；
步骤二，模拟预处理过程：采用模拟预处理模块将待测信号分流至两个通道中，模拟预处理模块包括功分器、噪声滤波器和信号放大器；
步骤三，低速adc采样过程：采用上位机labview软件控制nipxi-6251模块板卡，使用两个ai通道分别对模拟信号进行采样，其中一个通道的采样时钟与另一个通道的采样时钟跳变频率相同，但是上升沿相差一个较小的时间差te，且满足te≤1/fmax，fmax表示信号的最大频率，所以该部分称为时钟交错采样系统；
步骤四，运行零化滤波器谱估计算法：在上位机中，首先对采集到的样本进行离散傅立叶变换，从而获取fourier系数；然后运行谱估计算法，从获取到的fourier系数中估计出meds信号的参数
步骤五，运行复频率联合估计算法，通过参数采样速率fs，主、次通道采样时钟上升沿相差的时间间隔te以及n′≥k个连续采样值ze[n′]联合估计出阻尼因子rk以及极点频率fk的估计值输入为参数采样速率fs、交错时间te以及n′≥k个连续采样值ze[n′]，输出为阻尼因子rk和极点频率fk的估计值
2.如权利要求1所述的一种双通道多指数衰减正弦信号欠采样方法，其特征在于，所述步骤三中，信号采样过程的数学模型描述如下：
3.1，在主采样通道中，信号z(t)根据时钟clk1被均匀采样，并且采样可表示为如下形式：
其中，k为分量个数，并且k为已知量，ck是复振幅，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率，m∈z ；
3.2，在次采样通道中，信号z(t)根据时钟clk2被均匀采样，并且采样由下式给出：
其中，是复振幅ck的估计值，sk＝rk j2πfk是复频率，rk是阻尼因子，fk是极点频率，ts是采样周期，te是主、次通道采样时钟上升沿相差的时间间隔，n′∈z 。
3.如权利要求1或2所述的一种双通道多指数衰减正弦信号欠采样方法，其特征在于，所述步骤五的算法过程如下：
5.1，估计阻尼因子rk，阻尼因子rk可以很容易地从中估算出来，估算过程如下：由与sk＝rk j2πfk得
则得阻尼因子rk的估计值如下式：
同时得到极点频率fk的估计值如下式：
其中，∠(·)表示(·)的幅角主值，且0≤∠(·)＜2π；
5.2，由公式(3)和公式(4)得
其中，是已知的，是未知的；
5.3，求解bk的估计值：利用n′≥k个连续采样值ze[n′]构建一个n′×1的向量，1≤n′≤k，该向量形式如下：
利用的幅角主值创建一个k×n′的范德蒙矩阵，该矩阵形式如下：
利用和的幅角主值创建一个n′×k的对角矩阵
通过求解得到系数矩阵
5.4，求解极点频率fk的估计值，由于三角函数的周期性，估计的归一化频率与的参数相差2πmk，即：
其中，∠(·)是(·)的幅角主值，0≤∠(·)＜2π，mk∈z,k∈{1,2,…,k}，因为fmax表示信号的最大频率，得到如下不等式：
进一步化简公式(12)得如下不等式：
记公式(13)可转换为：
0≤mk≤(q-1)(14)
因为q∈z ，则有如果设置主、次通道采样时钟上升沿相差的时间间隔则有：
与公式(14)相结合，则有：
0≤mkfste＜1(16)
所以对于mk有唯一解，其解如下式：
其中，极点频率fk的估计值由唯一确定。
技术总结一种双通道多指数衰减正弦信号欠采样方法，对生成的MEDS信号进行采样，然后通过零化滤波器法对该信号进行谱估计，估算出参数其次通过复频率联合估计算法对该信号进行频率估计，从而完成参数测量进行重构，最后在上位机中显示重构效果。采用本发明方法，可以从少至3K个样本中测量出已知分量个数为K的MEDS信号；同时，本发明还提供了硬件原型来实现所提系统，通过仿真和实验结果验证了方案的有效性。
技术研发人员：黄国兴;倪安;卢为党;彭宏
受保护的技术使用者：浙江工业大学
技术研发日：2020.11.09
技术公布日：2021.03.12
