CN112269962B - 一种辐射噪声线谱稀疏反演方法 - Google Patents
一种辐射噪声线谱稀疏反演方法 Download PDFInfo
- Publication number
- CN112269962B CN112269962B CN202011064518.3A CN202011064518A CN112269962B CN 112269962 B CN112269962 B CN 112269962B CN 202011064518 A CN202011064518 A CN 202011064518A CN 112269962 B CN112269962 B CN 112269962B
- Authority
- CN
- China
- Prior art keywords
- noise
- line spectrum
- radiation noise
- signal
- representing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 164
- 230000005855 radiation Effects 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 72
- 230000008569 process Effects 0.000 claims description 35
- 230000003595 spectral effect Effects 0.000 claims description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000000737 periodic effect Effects 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 230000021615 conjugation Effects 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005309 stochastic process Methods 0.000 claims 1
- 230000009467 reduction Effects 0.000 abstract description 8
- 238000005259 measurement Methods 0.000 abstract description 6
- 239000007787 solid Substances 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 8
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开一种辐射噪声线谱稀疏反演方法,该方法针对现有辐射噪声声源分辨及声源强度估计精度不高的问题,通过利用目标辐射噪声谱的频域稀疏特性,建立压缩感知理论框架下的目标辐射噪声谱的稀疏表示模型,反演目标辐射噪声线谱在L1范数意义下的频域噪声系数,旨在提取高分辨的目标辐射噪声谱精细特征,获取噪声源强度的高精度估计结果,从而提升目标辐射噪声测量水平,为减振降噪的工程实现提供坚实的理论基础与技术支撑。
Description
技术领域
本发明涉及一种方法,具体的是涉及一种辐射噪声线谱稀疏反演方法,主要是用于目标分离及目标强度估计,属于噪声信号处理技术领域。
背景技术
通常来说,在减振降噪工程中,如何准确识别出起主要作用的噪声源,并从频率和空间域上获取各主要噪声源的空间分布和频率特征,定量识别各主要噪声源的贡献大小,从而有针对性地采取有效减振降噪措施是一项关键技术。低频段辐射噪声经常被用于目标检测、跟踪与识别、目标参数估计之中。线谱主要位于辐射噪声的低频段,所以,高精度地估计低频辐射噪声的参数是提升噪声测量水平的重要环节。
目前辐射噪声强度获取及方位的高精度高分辨测量主要采用基于声压阵和声矢量阵的阵列信号处理方法,存在无法识别声源的绝对强度大小、无法分辨毗邻入射信号、分辨精度低、信号处理成本高等问题。近年来,稀疏信号分解与信号重构作为一种新的理论逐渐被引入到信号参数估计领域,逐步成为信号处理、图像恢复和无线通信等领域研究的热点问题。稀疏信号分解与信号重构理论指出当信号满足稀疏或在某个变换域稀疏的条件时,能够通过少量的信号投影测量值准确或近似重构原始信号。稀疏反演方法通过最优化问题可以在稀疏模型下精确反演多峰辐射噪声谱,峰值和峰的位置分别对应噪声的声源强度与噪声源方位,在信息提取和参数分析等方面具有复杂信号稀疏化、特征差异明显和噪声影响小等优势。
综上,需提供一种辐射噪声声源分辨及声源强度估计精度高的辐射噪声线谱稀疏反演方法是解决上述技术问题的关键所在。
发明内容
针对上述背景技术中存在的诸多缺陷与不足,本发明对此进行了改进和创新,目的在于提供一种辐射噪声线谱稀疏反演方法,该方法能精确反演多峰稀疏线谱,是确定目标辐射噪声线谱强度的一种切实可行的方法。在信噪比满足一定要求时,声源识别的峰值平均误差小于0.5dB,保证了良好的反演精度。
本发明的方法针对现有辐射噪声声源分辨及声源强度估计精度不高的问题,通过利用目标辐射噪声谱的频域稀疏特性,建立压缩感知理论框架下的目标辐射噪声谱的稀疏表示模型,反演目标辐射噪声线谱在L1范数意义下的频域噪声系数,旨在提取高分辨的目标辐射噪声谱精细特征,获取噪声源强度的高精度估计结果,从而提升目标辐射噪声测量水平,为减振降噪的工程实现提供坚实的理论基础与技术支撑。
为解决上述问题并达到上述的发明目的,本发明一种辐射噪声线谱稀疏反演方法是通过采用下列的设计结构以及采用下列的技术方案来实现的:
一种辐射噪声线谱稀疏反演方法,利用目标辐射噪声谱的频域稀疏特性,建立压缩感知理论框架下的目标辐射噪声谱的稀疏表示模型,反演目标辐射噪声线谱在L1范数意义下的频域噪声系数,提取高分辨的目标辐射噪声谱精细特征,获取噪声源强度的高精度估计结果。
优选的,步骤1,根据线谱特征,建立水下大型结构辐射噪声在噪声窄带线谱上的稀疏表示模型;
步骤2,假设用一个宽平稳随机过程拟合,目标辐射噪声用随机过程表示,具体目标辐射噪声为:
式中,F(t)表示目标辐射噪声随机过程,H(t)表示宽带平稳随机过程,Sl(t) 表示具有随机初始相位的周期信号,t表示时间变量,表示求和,L是周期数,辐射噪声功率谱表示为
式中,F(f)表示辐射噪声功率谱,f表示频率变量,表示求极限,N是每段信号的傅里叶变换点数,E{·}表示数学期望值,Fk,N(f)表示采样获取的辐射噪声功率,k是信号段的数目;
步骤3,根据步骤2得出,若噪声的有效频带为[flow,fhigh],flow表示频带范围下限,fhigh表示频带范围上限,则可以用频率采样间隔Δf从频带[flow,fhigh]内的中取L个不同的频率点来离散频谱方程,且包含一些带有辐射噪声信息的离散频率分量,则(1)可表示为:
式中,F(t)表示目标辐射噪声随机过程,H(t)表示宽带平稳随机过程,表示求和,al为频率fl信号对应的分量振幅,值一定;/>表示不同频率信号间的相位误差,fl对应第l个线谱频率,j表示虚数单位,π表示圆周率;式(3) 中的傅里叶变换具体为:
式中,Fk,N(f)表示目标辐射噪声的傅立叶变换,F(t)表示目标辐射噪声随机过程,为傅里叶变换因子,Hk,N(f)表示宽带平稳随机过程的傅里叶变换,/>表示求和,al为频率fl信号对应的分量振幅,/>表示sinc函数,N 是每段信号的傅里叶变换点数,将结果代入式(4),假定连续谱分量和不同频率的线谱分量相互独立,令:
式中,G(f)表示等式右边求极限的结果,{}*表示求共轭,表示求极限, Hk,N(f)表示宽带平稳随机过程的傅里叶变换,N是每段信号的傅里叶变换点数, E{·}表示数学期望值,1/N{sin[Nπ(fl-f)]/π/(fl-f)}2,当N→∞,若f≠fl,结果趋向于0;若f=fl,结果趋向于∞,且拥有δ函数的性质可以得到:
式中,F(f)表示辐射噪声功率谱,G(f)为求极限的结果,δ(f-fl)是δ函数, E{·}表示求数学期望,al为频率fl信号对应的分量振幅;
步骤4,利用辐射噪声线谱集中稳定的能量来提高检测性能,线谱携带的频率信息可用于目标参数估计和类型识别;其中,
根据信号处理基本理论中时域与频域的对应关系,若频域是线谱,则时域信号为正弦信号,δ函数与三角函数为傅立叶变换对,则线谱表示为:
式中,s(t)表示线谱信号,L为线谱个数;al、fl和分别表示第l个线谱的振幅、频率和初始相位;
根据欧拉公式,将上述表达式中的三角函数转化为指数函数,统一公式中的常数项,使初始相位等于0度,则式(7)重新表示为:
式中,s(t)表示线谱信号,表示求和,L表示谱线的个数,al和fl分别表示第l个谱线谱强度和对应频率,/>为欧拉公式表示因子,fl对应第l个线谱频率,j表示虚数单位,π表示圆周率;令t=2nτ(n=1,2,3,…,Nn),τ为半回波间隔时间,Nn为辐射噪声回波信号个数且Nn>L;构造适当的稀疏字典矩阵/>其中/>根据稀疏理论,目标辐射噪声线谱的稀疏表示模型为:
s=Da (9)
式中,为辐射噪声的理想回波信号,a=[a1,a2,…,aL]T为与窄带线谱振幅相对应的矩阵振幅系数,向量D表示字典矩阵,通过求解式(9) 得到目标辐射噪声的线谱。
优选的,所述步骤4中,基于模型实现辐射噪声谱反演获取目标参数,并基于L1范数最小化的稀疏反演算法将目标辐射噪声幅值的L1范数约束为一个凸优化问题,在辐射噪声谱稀疏模型下求解得到最优解;其中,
范数最小约束问题表示为:
式中,||·||0表示L0范数,为/>的L0范数,即向量/>中非零元素的个数;是字典D的稀疏系数组成的向量;s.t.表示“受限于…”;/>表示观测到的目标辐射噪声信号,D为字典:L0范数最小化约束保证了只有少数元素不为零,非零元素由目标辐射噪声谱线的幅度组成,具体为噪声的绝对强度。
优选的,所述L1范数是L0范数的最优凸近似,并考虑到实际观测信号中必然存在噪声,根据式(10),建立了基于L1范数最小化约束的噪声稀疏模型,该噪声稀疏模型为:
式中,是/>的L1范数;||·||1表示L1范数,||·||2表示L2范数;s.t.表示“受限于…”;D为字典,/>表示观测到的目标辐射噪声信号;ε>0表示允许的边界误差,实际取经验值。
本发明与现有技术相比所产生的有益效果是:
1、本发明针对现有辐射噪声声源分辨及声源强度估计精度不高的问题,通过利用目标辐射噪声谱的频域稀疏特性,建立压缩感知理论框架下的目标辐射噪声谱的稀疏表示模型,反演目标辐射噪声线谱在L1范数意义下的频域噪声系数,旨在提取高分辨的目标辐射噪声谱精细特征,获取噪声源强度的高精度估计结果,从而提升目标辐射噪声测量水平,为减振降噪的工程实现提供坚实的理论基础与技术支撑;
2、本发明由于线谱只占整个辐射噪声谱的一部分,因而在频域上具有稀疏性,因此辐射噪声的线谱反演问题可以转化为基于逆Fourier的稀疏频率系数估计问题,首先在压缩感知理论框架下建立目标辐射噪声线谱的稀疏模型,然后反演L1范数意义下的线谱频域系数,该系数与目标辐射噪声强度相对应;线谱反演结果表明,该方法能精确反演多峰稀疏线谱,是确定目标辐射噪声线谱强度的一种切实可行的方法。在信噪比满足一定要求时,声源识别的峰值平均误差小于0.5dB,保证了良好的反演精度;
3、本发明不存在表达冗余,充分利用辐射噪声线谱的频域稀疏性,将辐射噪声的线谱反演问题可以转化为基于逆Fourier的稀疏频率系数估计问题;
4、与现有线谱检测方法相比,本发明基于稀疏信号分解理论,估计线谱更精准,并且能获取线谱对应的准确的声源强度信息,对实际中利用辐射噪声线谱集中稳定的能量提高检测性能及利用线谱携带的频率信息进行目标参数估计和类型识别具有重要意义;
5、本发明通过仿真计算,证明本发明能精确反演多峰稀疏线谱,是确定目标辐射噪声线谱强度的一种切实可行的方法。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明,其中:
图1为本发明目标辐射噪声线谱图;
图2为本发明分辨2个线谱信号的多峰反演图(其中,下图系上图的强度方向的局部放大);
图3为本发明分辨3个线谱信号的多峰反演图(其中,下图系上图的强度方向的局部放大);
图4为本发明在信噪比等于30dB条件下反演两个线谱信号的反演结果;
图5为本发明在信噪比等于20dB条件下反演两个线谱信号的反演结果;
图6为本发明在信噪比等于10dB条件下反演两个线谱信号的反演结果;
图7为本发明在信噪比等于5dB条件下反演两个线谱信号的反演结果;
图8为本发明反演线谱强度峰值相对误差随信号强度变化的性能仿真图;
图9为本发明反演线谱强度峰值精度指数随信号强度变化的性能仿真图。
具体实施方式
为了使本发明实现的技术手段、创造特征、达成目的与功效易于明白了解,下面结合附图以及具体实施方式对本发明的技术方案作更进一步详细的说明,需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
为了使本发明实现的技术手段、创造特征、达成目的与功效易于明白了解,下面结合附图以及具体实施方式对本发明的技术方案作更进一步详细的说明,需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
如图2~9所示的一种辐射噪声线谱稀疏反演方法,
进一步的,步骤1,根据线谱特征,建立水下大型结构辐射噪声在噪声窄带线谱上的稀疏表示模型;
在上述步骤1中,目标辐射噪声主要由机械噪声、螺旋桨噪声和水动力噪声组成,在100~6000Hz的频率范围内包含了所有的能量,目标辐射噪声的功率谱具有两个分量,一个是由离散的窄带单频点组成的噪声线谱,另一个是噪声幅值随频率连续变化的宽带连续谱;其中,
线谱反映了噪声信号的周期分量,主要由机械往复运动、螺旋桨桨叶和水动力产生,在辐射噪声的低频范围内存在离散的线谱,螺旋桨振动形成的线谱分量是一个稳定的离散线谱,与目标的速度无关,通常是100~1000Hz的低频强线谱,低频强线谱由各种机械振动和螺旋桨桨叶切断的水流形成;桨叶速度谱由桨叶基频和谐频的低频线谱分量组成,与船舶的航速和螺旋桨的桨叶数有关,满足关系式fl=lNpv1,其中fl是第l次谐波的频率,Np是螺旋桨叶片的数量, v1是螺旋桨的速度;
随着下潜深度的增加,线谱越来越明显,尤其是在目标低速航行时。
步骤2,假设用一个宽平稳随机过程拟合,目标辐射噪声用随机过程表示,具体目标辐射噪声为:
式中,F(t)表示目标辐射噪声随机过程,H(t)表示宽带平稳随机过程,Sl(t) 表示具有随机初始相位的周期信号,L是周期数,辐射噪声功率谱表示为
式中,F(f)表示辐射噪声功率谱,表示求极限,N是每段信号的傅里叶变换点数,E{·}表示数学期望值,k是信号段的数目;
上述式(1)与式(2)方程是数学模型下噪声谱的严格定义,即信号段的数目和每个信号段的时间长度趋于无穷大,在实际应用中,我们只能求有限长和有限个数的集合平均。
步骤3,根据步骤2得出,若噪声的有效频带为[flow,fhigh],则可以用频率采样间隔Δf从频带[flow,fhigh]内的中取L个不同的频率点来离散频谱方程,且包含一些带有辐射噪声信息的离散频率分量,则(1)可表示为:
式中,F(t)表示目标辐射噪声随机过程,H(t)表示宽带平稳随机过程,al为频率fl信号对应的分量振幅,值一定,式(3)中的傅里叶变换具体为:
式中,F(t)表示目标辐射噪声随机过程,Fk,N(f)表示目标辐射噪声的傅立叶变换,Hk,N(f)表示宽带平稳随机过程的傅里叶变换,al为频率fl信号对应的分量振幅,N是每段信号的傅里叶变换点数,将结果代入式(4),假定连续谱分量和不同频率的线谱分量相互独立,令:
式中,{}*表示共轭,表示求极限,N是每段信号的傅里叶变换点数,E{·} 表示数学期望值,1/N{sin[Nπ(fl-f)]/π/(fl-f)}2,当N→∞,若f≠fl,结果趋向于0;若f=fl,结果趋向于∞,且拥有δ函数的性质可以得到:
上述式(3)、式(4)、式(5)和式(6)方程表明,时域信号的离散频率分量在功率谱中产生线谱。
步骤4,利用辐射噪声线谱集中稳定的能量来提高检测性能,线谱携带的频率信息可用于目标参数估计和类型识别;其中,
根据信号处理基本理论中时域与频域的对应关系,若频域是线谱,则时域信号为正弦信号,δ函数与三角函数为傅立叶变换对,则线谱表示为:
式中,s(t)表示线谱信号,L为线谱个数,在具体实施中,假设待反演获取的线谱个数分别为L=2个、L=3个和L=4个;al、fl和分别表示第l个线谱的振幅、频率和初始相位,在实例中假设不同信号振幅不同,令al取3.5dB、 5dB,频率fl取100Hz~6000Hz之间的任意频率,初始相位/>
因此,利用某些特定振幅、频率和初相周期正弦函数的线性组合作为模拟线谱的模型是合适的。
根据欧拉公式,将上述表达式中的三角函数转化为指数函数,统一公式中的常数项,使初始相位等于0度,则式(7)重新表示为:
式中,s(t)表示线谱信号,L表示谱线的个数,al和fl分别表示第l个谱线谱强度和对应频率。令t=2nτ(n=1,2,3,…,Nn),τ为半回波间隔时间,Nn为辐射噪声回波信号个数且Nn>L;构造适当的稀疏字典矩阵其中/>根据稀疏理论,目标辐射噪声线谱的稀疏表示模型为:
s=Da (9)
式中,为辐射噪声的理想回波信号,a=[a1,a2,…,aL]T为与窄带线谱振幅相对应的矩阵振幅系数,D是字典矩阵,通过求解式(9)得到目标辐射噪声的线谱。
优选的,所述步骤2中,基于模型实现辐射噪声谱反演获取目标参数,并基于L1范数最小化的稀疏反演算法将目标辐射噪声幅值的L1范数约束为一个凸优化问题,在辐射噪声谱稀疏模型下求解得到最优解;其中,
范数最小约束问题表示为:
式中,为/>的L0范数,即向量/>中非零元素的个数;/>是字典D的稀疏系数组成的向量;/>表示观测到的目标辐射噪声信号,L0范数最小化约束保证了只有少数元素不为零,非零元素由目标辐射噪声谱线的幅度组成,具体为噪声的绝对强度。
优选的,所述L1范数是L0范数的最优凸近似,并考虑到实际观测信号中必然存在噪声,根据式(10),建立了基于L1范数最小化约束的噪声稀疏模型,该噪声稀疏模型为:
式中,是/>的L1范数;||·||2表示L2范数;ε>0表示允许的边界误差,一般取经验值。
在本发明中,由于L0范数极小化问题是一个NP难问题,L1范数是L0范数的最优凸近似,更容易获得最优解。考虑到实际观测信号中必然存在噪声,从而根据式(10)建立了上述的基于L1范数最小化约束的噪声稀疏模型。
本发明的效果可以通过以下3组具体仿真实例进一步说明:
设计仿真实验充分验证所提出的谱反演方法的性能。为了便于观察所提算法在不同频段的稀疏反演效果,将整个频段分为[1Hz,10Hz]、[10Hz,100Hz]、 [100Hz,1000Hz]和[1000Hz,10000Hz]四个子带。考虑到目标辐射噪声主要包含 10Hz~6000Hz的大部分能量,而线谱主要集中在1000Hz以下的频率范围内,因此仿真时合成信号主要分布在6000Hz以下。在所有的实验中,信噪比定义为
式中,s是根据公式(1)获得的模拟理想回波信号,是噪声回波信号。假设在模拟的辐射噪声线谱模型中,非零线谱个数假设为M=2,3,4,分别对应应不同峰值为h,利用所提出算法对信号/>进行反演得到的非零线谱个数为M′,对应峰值为h′。当M=M′,则峰值相对误差比Δhr和峰值平均误差Δh获取如下
其中|·|表示绝对值。如果M≠M′,则存在峰值数误差ΔM=M′-M。
1、为了验证所提出的方法对不同谱线峰值数的反演效果,仿真了辐射噪声峰值数为两个和三个的稀疏反演结果,图中蓝色线为理想峰值个数,红色为所提方法还原峰值谱结果。
由图2和图3可以看出,所给出的辐射噪声谱稀疏反演方法对不同频率的多个源具有良好的反演性能。用所提谱反演算法处理不同峰值数的稀疏模型时,得到的辐射噪声谱与仿真模型接近。其对应值与仿真模型基本一致,谱估计峰值与仿真模型峰值之间存在微小误差。在图2中,使用所提算法获得的辐射噪声线谱的两个峰值的峰值相对误差比仅为0.74%。在图3中,使用所提稀疏算法得到的辐射噪声线谱的三个峰值的峰值相对误差比为1.41%,也就是说,多个峰值反演得到的最大峰值平均误差也小于1dB。因此利用该方法估计辐射噪声源强度具有准确性。
2、为了验证所提稀疏算法在不同信噪比下的反演效果,在模拟的稀疏模型中,将不同信噪比为30dB、20dB、10dB和5dB的高斯白噪声加入到辐射噪声的回波曲线中,然后用本发明所提算法反演不同信噪比下辐射噪声的线谱。反演结果如图4~7所示。
通过对比图4~7中各信噪比下的线谱反演结果和辐射噪声仿真模型,可以看出当信噪比较高时,所提算法得到的稀疏线谱更接近仿真模型。随着信噪比的降低,噪声对反演结果的影响逐渐增大。当信噪比小于10dB时,所提算法的反演结果存在伪线谱目标。当信噪比小于5dB时,所提算法得到的噪声线谱趋于平坦,稀疏线谱无法准确恢复。
3、为了评价辐射噪声源强度对谱反演和声源大小估计精度的影响,引入精度指数(ACI)来评价线谱分量的提取是否正确。ACI表示为ACI=p′/p,其中是p实际线谱数,p′是反向谱线数。如果峰值相对误差比Δhr≤5%,则认为对应该频率的谱线被成功地反演。根据ACI的定义,ACI越大,反演精度越高,即线谱特征提取越好,噪声源强度估计越准确。
图8和图9显示了峰值相对误差和ACI与噪声源特性的关系曲线。噪声特性由噪声功率定义。当噪声强度从0.1dB增加到1dB时,算法估计的噪声强度越来越接近实际值,峰值相对误差比逐渐减小,如图8所示。当噪声特性大于 0.3dB时,ACI指数等于1,如图9所示,这意味着所有稀疏线谱都被精确地反演。
当然,本发明还可有其他多种实施例,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (3)
1.一种辐射噪声线谱稀疏反演方法,其特征在于,利用目标辐射噪声线谱的频域稀疏特性,建立压缩感知理论框架下的目标辐射噪声线谱的稀疏表示模型,反演目标辐射噪声线谱在L1范数意义下的频域噪声系数,提取高分辨的目标辐射噪声线谱精细特征,获取噪声源强度的高精度估计结果;其中,具体包括以下步骤:
步骤1,根据线谱特征,建立水下大型结构辐射噪声线谱在噪声窄带线谱上的稀疏表示模型;
步骤2,假设用一个宽平稳随机过程拟合,目标辐射噪声线谱用随机过程表示,具体地用于表示目标辐射噪声线谱的随机过程的表达式为:
式中,F(t)表示目标辐射噪声随机过程,用于表示目标辐射噪声线谱;H(t)表示宽带平稳随机过程,表示具有随机初始相位的周期信号,t表示时间变量,L0表示总周期数,l0表示范围区间[1,L0]内对应周期信号的第l0个周期,若l0=1,则代表周期信号的第1个周期,辐射噪声功率谱表示为
式中,F(f)表示辐射噪声功率谱,f表示频率变量,表示求极限,L1表示单个周期长度,N是每段信号的傅里叶变换点数,E{·}表示数学期望值,Fk,N(f)表示采样获取的辐射噪声功率,k是信号段的数目;
步骤3,根据步骤2得出,若噪声的有效频带为[flow,fhigh],flow表示频带范围下限,fhigh表示频带范围上限,则可以用频率采样间隔Δf从频带[flow,fhigh]内的中取Lq个不同的频率点来离散频谱方程,且包含一些带有辐射噪声信息的离散频率分量,则式(1)可表示为:
式中,F(t)表示目标辐射噪声随机过程,H(t)表示宽带平稳随机过程,t表示时间变量,L0表示总周期数,l0表示范围区间[1,L0]内对应周期信号的第l0个周期,为频率fl信号对应的分量的振幅,值一定;/>表示不同频率信号间的相位误差,fl对应第l个线谱频率,j表示虚数单位,π表示圆周率;式(3)中的傅里叶变换具体为:
式中,Fk,N(f)表示F(t)的傅里叶变换,F(t)表示目标辐射噪声随机过程,为傅里叶变换因子,相较于/>多了一个负号,因此表示含义不同;Hk,N(f)表示宽带平稳随机过程的傅里叶变换,al表示第l个线谱的振幅,/>表示sinc函数,N是每段信号的傅里叶变换点数,将结果代入式(4),假定连续谱分量和不同频率的线谱分量相互独立,令:
式中,G(f)表示等式右边求极限的结果,表示求极限,Hk,N(f)表示宽带平稳随机过程的傅里叶变换,N是每段信号的傅里叶变换点数,E{·}表示数学期望值;为了获取线谱能量特征,需要对式(4)进行求共轭运算,式(4)由Hk,N(f)和/>两部分组成,对Hk,N(f)求共轭得到G(f),对/>求共轭得到/>由数学原理可知,对于式子1/N{sin[Nπ(fl-f)]/π/(fl-f)}2,当N→∞时,若f≠fl,则结果趋向于0;若f=fl,则结果趋向于∞;再根据δ函数的性质即在负无穷到正无穷时间区间内积分为1且在t≠0它的值为0,那么式(4)取共轭后,就可以得到:
式中,F(f)表示辐射噪声功率谱,f表示频率变量,δ(f-fl)是δ函数,E{·}表示求数学期望,al为频率fl信号对应的分量振幅,Al表示功率;
步骤4,利用辐射噪声线谱集中稳定的能量来提高检测性能,线谱携带的频率信息可用于目标参数估计和类型识别;其中,
根据信号处理基本理论中时域与频域的对应关系,若频域是线谱,则时域信号为正弦信号,δ函数与三角函数为傅立叶变换对,则线谱信号表示为:
式中,s(t)表示线谱信号,L为线谱个数;al、fl和分别表示第l个线谱的振幅、频率和初始相位;
根据欧拉公式,将式(7)中的三角函数转化为指数函数,统一公式中的常数项,使初始相位等于0度,则式(7)重新表示为:
式中,s(t)表示线谱信号,L表示谱线的个数,为不同频率信号间的相位误差,fl对应第l个线谱频率,j表示虚数单位,π表示圆周率;令t=2nτ(n=1,2,3,…,Nn),τ为半回波间隔时间,Nn为辐射噪声回波信号个数且Nn>L;构造适当的稀疏字典矩阵/>其中/>根据稀疏理论,目标辐射噪声线谱的稀疏表示模型为:
s=Da (9)
式中,为辐射噪声的理想回波信号,a=[a1,a2,…,aL]T为与窄带线谱振幅相对应的矩阵振幅系数,向量D表示字典矩阵,通过求解式(9)得到目标辐射噪声的线谱。
2.根据权利要求1所述一种辐射噪声线谱稀疏反演方法,其特征在于,所述步骤4中,基于模型实现辐射噪声谱反演获取目标参数,并基于L1范数最小化的稀疏反演算法将目标辐射噪声幅值的L1范数约束为一个凸优化问题,在辐射噪声谱稀疏模型下求解得到最优解;其中,
范数最小约束问题表示为:
式中,||·||0表示L0范数,为/>的L0范数,即向量/>中非零元素的个数;是字典D的稀疏系数组成的向量;/>表示观测到的目标辐射噪声信号,D为字典;L0范数最小化约束保证了只有少数元素不为零,非零元素由目标辐射噪声谱线的幅度组成,具体为噪声的绝对强度。
3.根据权利要求2所述一种辐射噪声线谱稀疏反演方法,其特征在于,所述L1范数是L0范数的最优凸近似,并考虑到实际观测信号中必然存在噪声,根据式(10),建立了基于L1范数最小化约束的噪声稀疏模型,该噪声稀疏模型为:
式中,是/>的L1范数;||·||1表示L1范数,||·||2表示L2范数;D为字典,表示观测到的目标辐射噪声信号;ε表示允许的边界误差,其中,ε>0;Nn表示辐射噪声回波信号个数,n表示第n个回波信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011064518.3A CN112269962B (zh) | 2020-09-30 | 2020-09-30 | 一种辐射噪声线谱稀疏反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011064518.3A CN112269962B (zh) | 2020-09-30 | 2020-09-30 | 一种辐射噪声线谱稀疏反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112269962A CN112269962A (zh) | 2021-01-26 |
CN112269962B true CN112269962B (zh) | 2024-05-31 |
Family
ID=74338769
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011064518.3A Active CN112269962B (zh) | 2020-09-30 | 2020-09-30 | 一种辐射噪声线谱稀疏反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112269962B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539597A (zh) * | 2009-04-29 | 2009-09-23 | 哈尔滨工程大学 | 一种分辨具有相同频带辐射噪声的多目标方法 |
JP2014174149A (ja) * | 2013-03-13 | 2014-09-22 | Techno Science Japan Co Ltd | 放射ノイズ推定装置 |
CN108919240A (zh) * | 2018-04-23 | 2018-11-30 | 东南大学 | 一种基于组稀疏结构的水声目标辐射噪声调制谱重构方法 |
CN108918964A (zh) * | 2018-07-11 | 2018-11-30 | 国网江苏省电力有限公司宿迁供电分公司 | 一种稀疏性增强的谐波分析方法 |
CN109239712A (zh) * | 2018-07-24 | 2019-01-18 | 哈尔滨工程大学 | 基于水下声场和声能流的噪声探测方法 |
CN110118962A (zh) * | 2019-04-30 | 2019-08-13 | 东南大学 | 一种水声目标机动状态的辐射噪声仿真方法 |
-
2020
- 2020-09-30 CN CN202011064518.3A patent/CN112269962B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539597A (zh) * | 2009-04-29 | 2009-09-23 | 哈尔滨工程大学 | 一种分辨具有相同频带辐射噪声的多目标方法 |
JP2014174149A (ja) * | 2013-03-13 | 2014-09-22 | Techno Science Japan Co Ltd | 放射ノイズ推定装置 |
CN108919240A (zh) * | 2018-04-23 | 2018-11-30 | 东南大学 | 一种基于组稀疏结构的水声目标辐射噪声调制谱重构方法 |
CN108918964A (zh) * | 2018-07-11 | 2018-11-30 | 国网江苏省电力有限公司宿迁供电分公司 | 一种稀疏性增强的谐波分析方法 |
CN109239712A (zh) * | 2018-07-24 | 2019-01-18 | 哈尔滨工程大学 | 基于水下声场和声能流的噪声探测方法 |
CN110118962A (zh) * | 2019-04-30 | 2019-08-13 | 东南大学 | 一种水声目标机动状态的辐射噪声仿真方法 |
Non-Patent Citations (1)
Title |
---|
李耀波 ; 任志良 ; 刘恩洋 ; .一种鱼雷辐射噪声的线谱估计方法.指挥控制与仿真.2006,第28卷(第01期),63-70. * |
Also Published As
Publication number | Publication date |
---|---|
CN112269962A (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108919240B (zh) | 一种基于组稀疏结构的水声目标辐射噪声调制谱重构方法 | |
CN109840495B (zh) | 一种低信噪比条件下的轴频线谱增强方法 | |
CN107085140B (zh) | 基于改进的SmartDFT算法的非平衡系统频率估计方法 | |
Li et al. | Period-assisted adaptive parameterized wavelet dictionary and its sparse representation for periodic transient features of rolling bearing faults | |
Qin et al. | The order tracking of rotating machinery based on instantaneous frequency estimation | |
Hinich | A statistical theory of signal coherence | |
CN111665489A (zh) | 一种基于目标特性的线谱提取方法 | |
CN102538768A (zh) | 一种基于双频高频地波雷达的浅海水深测量方法 | |
CN113011093B (zh) | 基于LCWaveGAN的船舶航行噪声仿真生成方法 | |
Gao et al. | Automatic detection of underwater propeller signals using cyclostationarity analysis | |
CN101709997A (zh) | 一种振动信号处理的谐波窗函数 | |
CN106483520B (zh) | 一种船舶辐射噪声调制系数估计方法 | |
CN112269962B (zh) | 一种辐射噪声线谱稀疏反演方法 | |
CN1588124A (zh) | 利用多普勒变换的被动测速测距法及装置 | |
CN104954298B (zh) | 冲击噪声下且带有数据丢失的信号的频率估计方法 | |
CN106533394A (zh) | 一种基于自适应滤波器幅频响应的高精度频率估计方法 | |
Lan et al. | Wavelet and improved Hilbert–Huang transform method are used to study the spectrum distribution and energy of turbine pressure pulsation | |
CN115902791A (zh) | 基于s波段测波雷达时间多普勒谱的海浪反演方法及系统 | |
Liu et al. | Acoustic spectrum and signature analysis on underwater radiated noise of a passenger ship target based on the measured data | |
CN107450058A (zh) | 基于FrFT和HT的雷达信号时频参数估计方法 | |
CN110320018B (zh) | 一种基于二阶循环平稳特性的旋转机械复合故障诊断方法 | |
Li et al. | Sparse Bayesian learning for horizontal wavenumber retrieval in underwater acoustical signal processing | |
CN110275157B (zh) | 基于遗传算法自适应随机共振的矢量声定向方法 | |
CN112949163A (zh) | 一种基于解析函数理论的海浪谱和波高反演方法 | |
Liu et al. | Target depth extraction based on the character of sound field vertical wavenumber spectrum |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |