一种时域信号的数字时频测量方法及相应的目标识别方法技术

本发明专利技术提供一种时域信号的数字时频测量方法，包括：1)根据待测信号的采样率v设定时间窗口长度ΔT，确定时频测量的频谱范围，该频谱范围在1/ΔT到v/2之间；2)用时间窗口截取待测信号；3)在所确定的频谱范围内设定离散频点序列，对于每个频点，用频率值等于当前频点，且相位差恒定为90度的两个正弦信号作为参考信号，分别对当前待处理信号分片进行相关计算，将两个相关计算结果分别作为实部和虚部计算模与幅角，得到当前时间点、当前频点的幅度值和相位值；4)重复执行步骤2)至3)，直到得到被测信号每个时间点和频点组合所对应的幅度值和相位值。本发明专利技术能够在准确测出幅度时频谱的同时准确测出相位时频谱，且抗噪能力强，响应速度快，时间分辨力高。

【技术实现步骤摘要】
一种时域信号的数字时频测量方法及相应的目标识别方法
本专利技术涉及电子与时频测量
，具体地说，本专利技术涉及一种时域信号的数字时频测量方法。
技术介绍
频谱仪是一种典型的频谱测量装置，它广泛应用于电工电子、物理化学、生物医学与国防安全等各领域。随着应用水平和要求的逐步提高，人们对频谱仪功能的要求也越来越高。例如超声波诊断、心电图和脉冲雷达方面的应用，频谱仪不仅需要具有优良的频谱分辨能力，还应当具备优良的时间分辨能力。目前国内外市场流行的频谱仪是基于傅里叶变换的频谱仪，对于幅度谱，它具备良好的频谱分辨能力和时间分辨能力，能够在一定程度上处理复杂的瞬态响应的信号及其衍生信号的演变。然而，傅里叶变换存在时域的窗口截断效应，它需要利用有限的时间窗口内的信号来近似表示无限时长的信号，这造成相位测量与幅度测量存在明显偏差。即现有的基于傅里叶变换的频谱仪在准确测出信号的幅度谱的同时，其同时所测出的相位谱存在较大的偏差甚至错误。而相位谱能够反映被测信号的大量特征，如果不能准确地同步测出幅度谱和相位谱，就难以完整地分析信号(尤其是复杂的瞬态响应的信号及其衍生信号)中所携带的信息。为克服上述缺陷，中国专利申请CN101308175A提出了一种在傅里叶变换基础上的改进方案，该方案在傅里叶变换过程中引入了一些参数对其相位谱进行修正，从而在一定程度上缩小了相位偏差。然而这种方案不能从根本上改变傅里叶变化中的时域窗口截断效应，并且这种方案有可能因人为引入参数带来额外的偏差，因此其相位谱的准确度仍然存在不足，有待提高。因此，当前迫切需要一种能够准确地同步测出时域信号的幅度谱和相位谱的解决方案。
技术实现思路
本专利技术的任务是提供一种能够准确地同步测出时域信号的幅度谱和相位谱的解决方案。为实现上述专利技术目的，本专利技术提供了一种数字时频测量方法，包括下列步骤：1)设定时间窗口长度ΔT，根据所设定的时间窗口长度ΔT和待测信号的采样率v确定能够测量的频谱范围，该频谱范围是：1/ΔT到v/2；2)用时间窗口截取待测信号，得到当前时间点对应的待处理信号分片；3)在步骤1)确定的范围内设定离散频点序列，对于频点序列中的每个频点，用频率值等于当前频点，且相位差恒定为90度的两个正弦参考信号，分别对当前待处理信号分片进行相关计算；将所述两个参考信号对应的相关计算结果分别作为实部和虚部组成复数，计算所组成的复数模与幅角，将所述模与幅角分别作为当前时间点、当前频点的幅度值和相位值；4)将下一个时间点设定为当前时间点，重复执行步骤2)至3)，直到得到被测信号每个时间点和频点组合所对应的幅度值和相位值。其中，所述步骤1)中，所确定的所述频谱范围是：2/ΔT到v/5。其中，所述步骤3)还包括：对于每个当前频点f，计算其所对应的当前周期k/f，从当前的所述待处理信号分片中舍去待处理信号分片末尾的一段数据，以保证参与相关计算的待处理信号的时间长度是当前频点对应周期长度1/f的整数倍；并且，每个所述参考信号的时间长度均与所述参与相关计算的待处理信号的时间长度一致。其中，所述步骤3)还包括：通过在所述步骤1)确定的范围内线性取点或非线性取点设定所述的离散频点序列。其中，所述步骤3)中，所述非线性取点包括：对数均匀取点，多项式函数均匀取点或者倒数均匀取点。即在步骤1)确定的范围内，沿着对数，多项式或者倒数函数在横坐标轴上进行均匀取点，所取点的纵坐标就是所取的频点值，这样所得到的离散频点序列就会沿着对数，多项式或者倒数函数排列，从而更加灵活地获得所需的时频谱。其中，所述步骤3)中，当设定所述的离散频点序列采用线性取点时，对于离散频点序列中的任意一个频点f，使时间窗口长度ΔT是该频点f对应的周期长度1/f的整数倍。与现有技术相比，本专利技术具有下列技术效果：1、能够在准确测出幅度时频谱的同时准确测出相位时频谱。2、抗噪能力强。3、响应速度快。4、时间分辨力高。附图说明图1示出了本专利技术技术原理的示意图；图2示出了本专利技术一个实施例的流程图；图3示出了本专利技术一个实施例中的一个待分析信号的时域波形图；图4示出了对于图3的被分析信号，根据本专利技术一个实施例所测的幅度与相位频谱(灰线表示)与传统的傅里叶变换方案所测的幅度与相位频谱(黑色方框连线表示)的对照示意图；图5示出了本专利技术一个实施例中的另一个待分析信号的时域波形图；图6示出了对于图5的被分析信号，根据本专利技术一个实施例所测的幅度频谱；图7示出了对于图5的被分析信号，根据本专利技术一个实施例所测的时频联合分析幅度谱图；图8示出了对于图5的被分析信号，根据本专利技术一个实施例所测的时频联合分析相位谱图；图9示出了本专利技术一个实施例中的又一个待分析信号的时域波形图；图10示出了对于图9的被分析信号，根据本专利技术一个实施例所测的幅度谱图；图11示出了对于图9的被分析信号，在较强噪声情形下根据本专利技术一个实施例所测的116微秒到124微秒之间的信号分片的信号时域分布图和频谱分布图。其中a)b)分别为不含目标信号时背景噪声的时域和频域谱；c)d)分别为单次测量的带噪声目标信号的时域和频域谱；e)f)分别为10次带噪声目标信号测量的平均时域和频域谱；g)h)分别为100次带噪声目标信号测量的平均时域和频域谱；i)j)分别为1000次带噪声目标信号测量的平均时域和频域谱。具体实施方式下面，结合附图和实施例对本专利技术做进一步的描述。为便于理解，首先介绍本专利技术的测量原理。图1示出了本专利技术技术原理的示意图，本专利技术提出锁相环与相关器组合的方法来实现高密度频谱的幅度和相位的快速准确测量，此时相位的测量直接通过锁相环的两路信号分别与被测信号进行相关而获得，而无需反复扫描相位，大大节省处理时间。具体来讲，本专利技术使用的锁相时频测量的基本原理是：首先设定待测频率ω＝2πf，由锁相环负责产生相位差恒定为90度的两路正弦波信号，如cos(ωt)与sin(ωt)，即复数信号eiωt；然后让相位被锁定的两路信号分别与被测信号进行如公式(1a)或(1b)的相关运算，获得信号在该设定频率下的复数(实部和虚部，或幅度和相位)。原理公式表达如下：(连续情形)(1a)(离散情形)(1b)其中复数函数S(ω)是复数信号S(t)的频谱，Sk为复数离散信号序列，i为虚数单位即-1的平方根，n为每周期数字样本个数，N为n的自然数倍，信号长度T为1/f的自然数倍。对于通常的信号，复数函数S(t)的虚部为零，此时复数函数S(t)等于实数函数S(t)，而复数S(ω)的实部和虚部分别表示，即其中，根据公式(1a)，将复数S(ω)按实部和虚部分别表示，得到公式(2a)和(3a)，根据公式公式(1b)，将复数S(ω)按实部和虚部分别表示，得到公式(2b)和(3b)(连续情形)(2a)(离散情形)(2b)(连续情形)(3a)(离散情形)(3b)需要指出的是，本专利技术中使用的基本原理，即公式(1a)和(1b)，与傅里叶变换的积分或求和在形式上非常相似，这是因为本专利技术的公式具备时域到频域或实空间向密度空间变换所普遍具有的特征。实际上，本专利技术使用的时间分辨锁相原理在本质上不同于傅里叶变换。具体来说，傅里叶变换假设原始信号在正负无穷大区间内，只有满足该理想假设时任意两个频率之间的正交性才能保证，进而使得频谱准确性可以保证。但该理想本文档来自技高网...

【技术保护点】
一种时域信号的数字时频测量方法，包括下列步骤：1)接收待测数字信号，根据待测信号的采样率v设定时间窗口长度ΔT，确定时频测量的频谱范围，所述频谱范围在1/ΔT到v/2之间，其中1/ΔT<v/2；2)用时间窗口截取待测信号，得到当前时间点对应的待处理信号分片；3)在步骤1)确定的范围内设定离散频点序列，对于离散频点序列中的每个频点，用频率值等于当前频点，且相位差恒定为90度的两个正弦信号作为参考信号，分别对当前待处理信号分片进行相关计算；将两个所述参考信号对应的相关计算结果分别作为复数的实部和虚部，然后计算所述复数的模与幅角，并将所述模与幅角分别作为当前时间点、当前频点的幅度值和相位值；4)将下一个时间点设定为当前时间点，重复执行步骤2)至3)，直到得到被测信号每个时间点和频点组合所对应的幅度值和相位值。

【技术特征摘要】
1.一种时域信号的数字时频测量方法，包括下列步骤：1)接收待测数字信号，根据待测信号的采样率v设定时间窗口长度ΔT，确定时频测量的频谱范围，所述频谱范围在1/ΔT到v/2之间，其中1/ΔT<v/2；2)用时间窗口截取待测信号，得到当前时间点对应的待处理信号分片；3)在步骤1)确定的范围内设定离散频点序列，对于离散频点序列中的每个频点，用频率值等于当前频点，且相位差恒定为90度的两个正弦信号作为参考信号，分别对当前待处理信号分片进行相关计算；将两个所述参考信号对应的相关计算结果分别作为复数的实部和虚部，然后计算所述复数的模与幅角，并将所述模与幅角分别作为当前时间点、当前频点的幅度值和相位值；4)将下一个时间点设定为当前时间点，重复执行步骤2)至3)，直到得到被测信号每个时间点和频点组合所对应的幅度值和相位值。2.根据权利要求1所述的时域信号的数字时频测量方法，其特征在于，所述步骤1)中，所确定的所述频谱范围是：2/ΔT到v/5。3.根据权利要求2所述的时域信号的数字时频测量方法，其特征在于，所述步骤3)还包括：对于每个当前频点f，计算其所对应的当前周期k/f，从当前的所述待处理信号分片中舍去待处理信号分片末尾的一段数据，以保证参与相关计算的待处理信号的时间长度是当前频点对应周期长度1/f的整...

【专利技术属性】
技术研发人员：陆俊，沈保根，邵晓萍，
申请(专利权)人：中国科学院物理研究所，
类型：发明
国别省市：北京;11