一种稳健的基于线谱多普勒频移估计运动声源速度的方法技术

本发明专利技术实施例公开了一种稳健的基于线谱多普勒频移估计运动声源速度的方法，方法包括：初始化计算参数，所述计算参数包括窗函数长度、傅里叶变换长度以及重叠长度；根据计算参数执行短时傅里叶变换，以得到接收信号的STFT时频谱图；根据STFT时频谱图计算出未迭代时的多普勒有关参数，所述多普勒有关参数包括线谱频率、速度、正横时刻和正横距离；根据STFT时频谱图计算时频能量聚集度；根据DopplerCT时频分析算法对多普勒有关参数进行迭代处理；当达到迭代停止要求后，输出最终的多普勒有关参数。本发明专利技术显著提升了低信噪比环境下运动声源速度、正横距离等多普勒有关参数估计过程的稳健性。稳健性。稳健性。

【技术实现步骤摘要】
一种稳健的基于线谱多普勒频移估计运动声源速度的方法


[0001]本专利技术涉及信号处理
，更具体地说是一种稳健的基于线谱多普勒频移估计运动声源速度的方法。

技术介绍

[0002]由于机械结构转动的原因，汽车、飞机、舰船等类型的目标在运动过程中不可避免的向外辐射噪声。一般地，目标辐射噪声信号可分为线谱噪声和宽带噪声两种类型，其中线谱噪声的声源级显著高于宽带噪声的声源级，能够携带目标声源的特征信息传播较远距离，因此线谱噪声信息常用于目标识别和方位跟踪。当线谱声源与接收点之间存在相对运动时，接收线谱信号的瞬时频率将会产生多普勒频移，根据多普勒频移可以估计运动声源速度、正横距离等多普勒有关参数。
[0003]目前较为普遍的运动声源速度估计方式是提取线谱瞬时频率曲线后，依据多普勒频移的理论公式，基于最小二乘曲线拟合技术进行估计。早期，学者们通过相位插值方法估计接收信号中的多普勒频移变化曲线。随着时频分析技术的兴起，1994年，Ferguson和Quinn提出利用短时傅里叶变换(STFT，short
‑
time Fourier transform，也称为spectrogram)和Wigner
‑
Ville分布(WVD，Wigner
–
Ville distribution)从接收信号中提取多普勒频移变化曲线。随后，为了进一步提高多普勒频移变化曲线的估计精度，学者们又尝试了众多高分辨时频分析算法，包括多项式Wigner
‑
Ville分布(polynomial Wigner
–
Ville distribution)、多项式线性调频小波变换(PCT，polynomial chirplet transform)、匹配Wigner变换(MWT，matched Wigner transform)、扩展样条调频小波变换(ESCT，enhanced spline
‑
kernelled chirplet transform)、多普勒调频小波变换(DopplerCT,Doppler chirplet transform)、Wigner
‑
Ville分布交叉项抵消方法(WVD
‑
CTR，Wigner
–
Ville distribution cross
‑
term rejection)等。相比于经典的STFT方法，高分辨时频分析算法给出的时频图中线谱信号的时频能量聚集度较高，一定程度上提高了多普勒频移变化曲线的提取精度，理论上有助于在低信噪比条件下实现稳健的运动声源速度估计。但是在实际应用中发现，仅仅提高时频图中线谱信号的时频能量聚集度对于在低信噪比条件下参数估计稳健性的提升有限。
[0004]高分辨时频分析算法通常基于迭代的方式不断提高信号的时频能量聚集度以及参数估计的精度，即不断重复“计算时频图、提取线谱瞬时频率曲线、基于曲线拟合方式估计运动声源速度”的过程，直至迭代收敛，这个流程设计中最容易出错的地方是初次迭代。初次迭代时由于时频图中线谱信号的时频能量聚集度还比较低，当信噪比较低时，受噪声影响，很难准确提取线谱瞬时频率曲线，使得曲线拟合后给出错误的运动声源速度估计结果，基于错误的参数估计结果继续下一次迭代，只会导致参数估计结果难以收敛。根据以上分析可知，这类算法实际上对噪声较为敏感，一般很难适应低信噪比环境。但在实际应用中，由于接收点距离声源较远、背景噪声较高等原因，低信噪比环境难以避免。
[0005]因此，有必要设计一种适用于低信噪比环境的运动声源速度的估计方法。

技术实现思路

[0006]本专利技术的目的在于克服现有技术的不足，提供一种稳健的基于线谱多普勒频移估计运动声源速度的方法，旨在提升低信噪比环境下运动声源速度、正横距离等多普勒有关参数估计过程的稳健性。
[0007]为实现上述目的，本专利技术采用以下技术方案：
[0008]第一方面，一种稳健的基于线谱多普勒频移估计运动声源速度的方法，包括：
[0009]初始化计算参数，所述计算参数包括窗函数长度、傅里叶变换长度以及重叠长度；
[0010]根据计算参数执行短时傅里叶变换，以得到接收信号的STFT时频谱图；
[0011]根据STFT时频谱图计算出未迭代时的多普勒有关参数，所述多普勒有关参数包括线谱频率、速度、正横时刻和正横距离；
[0012]根据STFT时频谱图计算时频能量聚集度；
[0013]根据DopplerCT时频分析算法对多普勒有关参数进行迭代处理；
[0014]当达到迭代停止要求后，输出最终的多普勒有关参数。
[0015]其进一步技术方案为：所述初始化计算参数，所述计算参数包括窗函数长度、傅里叶变换长度以及重叠长度，包括：
[0016]设定计算STFT时频图所需要的计算参数；
[0017]设定计算DopplerCT时频图所需要的计算参数。
[0018]其进一步技术方案为：所述DopplerCT时频图所需的窗函数长度和傅里叶变换长度均大于STFT时频图所需的窗函数和傅里叶变换长度。
[0019]其进一步技术方案为：所述根据STFT时频谱图计算出未迭代时的多普勒有关参数，所述多普勒有关参数包括线谱频率、速度、正横时刻和正横距离中，计算未迭代时的多普勒有关参数的公式为：
[0020][0021]式中，P(t,f)表示时频图，t表示时间，f表示频率，σ表示瞬时频率曲线的粗细程度，f
u
为频率区间的上界，f
d
为频率区间的下界。
[0022]其进一步技术方案为：所述根据STFT时频谱图计算时频能量聚集度中，计算时频能量聚集度的公式为：
[0023][0024]式中，H
α
(P)表示P(t,f)的Renyi熵，Renyi熵用于度量P(t,f)的时频能量聚集度，α表示Renyi熵阶数，∫∫P(u,v)dudv表示通过积分计算P(t,f)的总能量。
[0025]其进一步技术方案为：所述根据DopplerCT时频分析算法对多普勒有关参数进行迭代处理，包括：
[0026]根据多普勒有关参数计算得到DopplerCT时频图；
[0027]根据DopplerCT时频图计算出迭代更新后的多普勒有关参数；
[0028]根据DopplerCT时频图计算时频能量聚集度；
[0029]根据时频能量聚集度判断迭代是否终止。
[0030]其进一步技术方案为：所述根据多普勒有关参数计算得到DopplerCT时频图，计算得到DopplerCT时频图的公式如下：
[0031][0032][0033]式中，DopplerCT(t,f,p)表示参数p条件下计算得到的DopplerCT时频图，中τ表示t被积分量，y表示接收线谱信号，w表示窗函数，c表示声速，参数p＝[f0,c,v,τ
c
,d
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种稳健的基于线谱多普勒频移估计运动声源速度的方法，其特征在于，包括：初始化计算参数，所述计算参数包括窗函数长度、傅里叶变换长度以及重叠长度；根据计算参数执行短时傅里叶变换，以得到接收信号的STFT时频谱图；根据STFT时频谱图计算出未迭代时的多普勒有关参数，所述多普勒有关参数包括线谱频率、速度、正横时刻和正横距离；根据STFT时频谱图计算时频能量聚集度；根据DopplerCT时频分析算法对多普勒有关参数进行迭代处理；当达到迭代停止要求后，输出最终的多普勒有关参数。2.根据权利要求1所述的一种稳健的基于线谱多普勒频移估计运动声源速度的方法，其特征在于，所述初始化计算参数，所述计算参数包括窗函数长度、傅里叶变换长度以及重叠长度，包括：设定计算STFT时频图所需要的计算参数；设定计算DopplerCT时频图所需要的计算参数。3.根据权利要求2所述的一种稳健的基于线谱多普勒频移估计运动声源速度的方法，其特征在于，所述DopplerCT时频图所需的窗函数长度和傅里叶变换长度均大于STFT时频图所需的窗函数和傅里叶变换长度。4.根据权利要求1所述的一种稳健的基于线谱多普勒频移估计运动声源速度的方法，其特征在于，所述根据STFT时频谱图计算出未迭代时的多普勒有关参数，所述多普勒有关参数包括线谱频率、速度、正横时刻和正横距离中，计算未迭代时的多普勒有关参数的公式为：s.t.f
d
＜f＜f
u
；式中，P(t,f)表示时频图，t表示时间，f表示频率，σ表示瞬时频率曲线的粗细程度，f
u
为频率区间的上界，f
d
为频率区间的下界。5.根据权利要求1所述的一种稳健的基于线谱多普勒频移估计运动声源速度的方法，其特征在于，所述根据STFT时频谱图计算时频能量聚集度中，计算时频能量聚集度的公式为：式中，H
α
(P)表示P(t,f)的Renyi熵，Renyi熵用于度量P(t,f)的时频能量聚集度，α表示Renyi熵阶数，∫∫P(u,v)dudv表示通过积分计算P(t,f)的总能量。6.根据权利要求1所述的一种稳健的基于线谱多普勒频移估计运动声源速度的方法，其特征在于，所述根据DopplerCT时频分析算法对多普勒有关参数进行迭代处理，包括：根据多普勒有关参数计算...

【专利技术属性】
技术研发人员：杨益新，梁宁宁，周建波，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：