CN101539597B - 一种分辨具有相同频带辐射噪声的多目标方法 - Google Patents

一种分辨具有相同频带辐射噪声的多目标方法 Download PDF

Info

Publication number
CN101539597B
CN101539597B CN2009100719165A CN200910071916A CN101539597B CN 101539597 B CN101539597 B CN 101539597B CN 2009100719165 A CN2009100719165 A CN 2009100719165A CN 200910071916 A CN200910071916 A CN 200910071916A CN 101539597 B CN101539597 B CN 101539597B
Authority
CN
China
Prior art keywords
target
axle
frequency
demon
line spectrum
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.)
Expired - Fee Related
Application number
CN2009100719165A
Other languages
English (en)
Other versions
CN101539597A (zh
Inventor
王逸林
殷敬伟
惠俊英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Rugao Productivity Promotion Center
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN2009100719165A priority Critical patent/CN101539597B/zh
Publication of CN101539597A publication Critical patent/CN101539597A/zh
Application granted granted Critical
Publication of CN101539597B publication Critical patent/CN101539597B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明提供的是一种基于矢量DEMON信号处理的分辨具有相同频带辐射噪声的多目标方法。(1)将矢量水听器声压信号与振速信号进行互谱处理,完成包络检测,获得DEMON线谱;(2)对DEMON线谱进行融合与净化;(3)采用基于最大公约数算法的方位加权的串行迭代法对于DEMON线谱提取轴频;(4)通过轴频信息分辨多目标数量,并给出各目标的方位信息。本发明的优点是可以区分多目标,即使多个目标的辐射噪声具有相同频带,但由于各目标的轴频不同,仍然可以分辨。

Description

一种分辨具有相同频带辐射噪声的多目标方法
(一)技术领域
本发明涉及的是水声领域,更确切地说,涉及一种实现单矢量传感器分辨多目标的方法。
(二)背景技术
矢量传感器利用所接收到目标的声压和振速信息可以测得目标方位,平均声强器可以通过目标连续谱声辐射来测向。当多个目标具有频带重叠的连续谱特性时,单矢量传感器将无法分辨多目标,而只能测得声场中合成声强流的方向。但由于DEMON线谱可以反映舰船目标的螺旋桨叶片数和主轴转速等物理参数,是舰船辐射噪声的不变特征,因此,利用矢量传感器的对舰船辐射噪声DEMON线谱提取并帮助区分多目标轴频,再利用目标的DEMON线谱的轴频来对多目标测向与识别。可以解决单矢量传感器多目标分辨和测向问题。
文献“用单个矢量水听器基于DEMON线谱的双目标分辨”(《应用声学》,2005,3)中也是利用DEMON来分辨多个目标,但其是利用人工分簇的方法来分辨两个目标的,并未利用方位信息。
文献“基于DEMON线谱的轴频提取方法研究”(《应用声学》,2005,6)提出了基于标量水听器的最大公约数算法提取DEMON线谱轴频,但只应用于单目标情况。
(三)发明内容
本发明的目的在于提供一种可以大大提高轴频提取的准确率,同时可以对每个目标测向的一种分辨具有相同频带辐射噪声的多目标方法。
本发明的目的是这样实现的:
(1)将矢量水听器测得的声压信号与振速信号进行互谱处理,完成包络检测,获得DEMON线谱;
(2)对DEMON线谱进行融合与净化;
(3)采用基于最大公约数算法的方位加权的串行迭代法对于DEMON线谱提取轴频;
(4)通过轴频信息分辨多目标数量,并给出各目标的方位信息。
为实现本发明的目的,单矢量传感器分辨具有相同频带辐射噪声的多目标的方法,本发明利用目标DEMON线谱是由舰船螺旋桨运动产生的,是舰船的特征谱线,提出通过单个矢量传感器提取多目标DEMON线谱的轴频,来分辨多目标并测向的方法。
本发明的优点是可以在具有相同频带辐射噪声的多目标的情况下,通过对DEMON线谱轴频的提取,实现多目标分辨。与常规的DEMON谱分析方法相比,矢量传感器提供的方位信息可以大大提高轴频提取的准确率,同时可以对每个目标测向。
(四)附图说明
图1是矢量DEMON谱分析框图;
图2是常规DEMON第一次提取时差频的品质因数;
图3是常规DEMON第二次提取时差频的品质因数;
图4是矢量DEMON线谱方位分布图;
图5是矢量DEMON谱多目标轴频提取流程图;
图6是矢量DEMON第一次提取时差频的品质因数;
图7是矢量DEMON第二次提取时差频的品质因数。
(五)具体实施方式
下面结合附图举例对本发明做更详细地描述:
首先是对DEMON谱的提取:本发明这个步骤与传统矢量DEMON分析法相同,见文献“用单个矢量水听器基于DEMON线谱的双目标分辨”。
第二步是DEMON线谱的融合与净化:本发明这个步骤与传统DEMON分析法相同,见文献“基于DEMON线谱的轴频提取方法研究”。
第三步对融合和净化后的DEMON线谱进行测向:
设有舰船辐的DEMON线谱:
Figure G2009100719165D00021
其中,f1为基频,
Figure G2009100719165D00022
为初始相位。
平面波条件下,矢量水听器接收到的声压和振速分别为:
p=cos(2πf0t)·s1(t)          (2)
vx=cos(2πf0t)·s1(t)·cos(θ)(3)
vy=cos(2πf0t)·s1(t)·sin(θ)(4)
以上3式中,f0是被调制的高频信号的频率,θ是目标方位。
为获得DEMON线谱,需要对被调制的信号进行解调,为此将(2)分别与(3)和(4)式相乘,得到声强流:
Ix=p·vx=cos2(2πf0t)·s2 1(t)·cos(θ)(5)
Iy=p·vy=cos2(2πf0t)·s1 2(t)·sin(θ)(6)
上两式的声强流,通过一个低通滤波器后就可以检出包络,包络中含有DEMON谱和目标方位信息。
然后用如下公式测向:
θ = tan - 1 Iy ( f ) Ix ( f ) - - - ( 7 )
第四步是本发明的创新之处,也是关键所在。最大公约数算法提取轴频法只能用于单目标轴频提取,本发明以其为基础,加入方位加权的串行迭代法提取多个目标轴频,流程图如图5。
(1)求出所有线谱之间的差频,应用最大公约数算法计算所有差频的品质因数,设定一个品质因数门限,品质因数最大且过门限的差频是第一个目标的轴频;
(2)当确定第一个目标的轴频后,根据串行迭代法,应将该目标的轴频和其所有的倍频剔除,否则无法使用最大公约数算法提取下一个目标的轴频。但同时应将多个目标的公共倍频保留,否则会降低其他目标的品质因数。因此在去除第一个目标的轴频及其倍频时参照方位信息,如果倍频方位与轴频方位相同或在一个误差范围内,说明该倍频只有该目标的能量,不是其它目标轴频的倍频,可以剔除;如果倍频方位与轴频方位不同说明该谱线合成了其它目标的能量,这根谱线也是其它目标的倍频,因此不能剔除需要保留。
(3)对剩下的线谱重复以上算法提取第二个目标的轴频。如此反复直到所有差频的品质因数都不过门限,表示所有目标的轴频都被提取。
第五步通过轴频信息分辨多目标数量。
是舰船的特征谱线,即使是同一型号的舰船轴频也会有所不同,轴频中只包含该舰船的能量。因此用第三步中的测向方法测得的方位为该目标的真实方位,而不是声场中合成声强流的方向。对各个目标的轴频用第三步中的测向方法测向,并给出各目标的方位信息。
下面是对本方法的仿真,以双目标为例,两个目标的轴频分别为10Hz和13Hz,信噪比为6dB。
图2为第一次提取时差频的品质因数,图3为第二次提取时差频的品质因数,由图可见由于第一次将不该剔除的线谱剔除,导致第二目标的品质因数下降无法提取轴频。
图4是矢量DEMON线谱方位分布图。根据图4的方位对数据作处理,结果如图6和图7,由图中可见,可以完整提取两个目标的轴频。
基于矢量的DEMON线谱分析与提取技术可以很好的分辨具有相同频带辐射噪声的多目标,即使多个目标的轴频具有公共倍频。

Claims (1)

1.一种分辨具有相同频带辐射噪声的多目标方法,其特征是:
(1)将矢量水听器测得的声压信号与振速信号进行互谱处理,完成包络检测,获得DEMON线谱;
(2)对DEMON线谱进行融合与净化;
(3)对融合和净化后的DEMON线谱进行测向;
设有舰船辐射的DEMON线谱:
Figure FSB00000269842500011
其中,f1为基频,
Figure FSB00000269842500012
为初始相位;
平面波条件下,矢量水听器接收到的声压和振速分别为:
p=cos(2πf0t)·s1(t)
vx=cos(2πf0t)·s1(t)·cos(θ)
vy=cos(2πf0t)·s1(t)·sin(θ)
其中,f0是被调制的高频信号的频率,θ是目标方位;
将p分别与vx和vy相乘,得到声强流:
Ix=p·vx=cos2(2πf0t)·s2 1(t)·cos(θ)
Iy=p·vy=cos2(2πf0t)·s1 2(t)·sin(θ)声强流,通过一个低通滤波器后就检出包络,包络中含有DEMON线谱和目标方位信息;
然后用如下公式测向:
θ = tan - 1 Iy ( f ) Ix ( f ) ;
(4)采用基于最大公约数算法的方位加权的串行迭代法对于DEMON线谱提取轴频;其中方位加权的串行迭代法提取多个目标轴频的方法为:
(a)求出所有线谱之间的差频,应用最大公约数算法计算所有差频的品质因数,设定一个品质因数门限,品质因数最大且过门限的差频是第一个目标的轴频;
(b)当确定第一个目标的轴频后,去除第一个目标的轴频及其倍频时参照方位信息,如果倍频方位与轴频方位相同或在一个误差范围内,可以剔除;如果倍频方位与轴频方位不同说明该线谱合成了其它目标的能量,不能剔除需要保留;
(c)对剩下的线谱重复步骤(a)、(b)提取第二个目标的轴频,如此反复直到所有差频的品质因数都不过门限,表示所有目标的轴频都被提取;
(5)通过轴频信息分辨多目标数量,并给出各目标的方位信息。
CN2009100719165A 2009-04-29 2009-04-29 一种分辨具有相同频带辐射噪声的多目标方法 Expired - Fee Related CN101539597B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100719165A CN101539597B (zh) 2009-04-29 2009-04-29 一种分辨具有相同频带辐射噪声的多目标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100719165A CN101539597B (zh) 2009-04-29 2009-04-29 一种分辨具有相同频带辐射噪声的多目标方法

Publications (2)

Publication Number Publication Date
CN101539597A CN101539597A (zh) 2009-09-23
CN101539597B true CN101539597B (zh) 2011-01-26

Family

ID=41122894

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100719165A Expired - Fee Related CN101539597B (zh) 2009-04-29 2009-04-29 一种分辨具有相同频带辐射噪声的多目标方法

Country Status (1)

Country Link
CN (1) CN101539597B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102255671B (zh) * 2011-07-06 2014-01-29 哈尔滨工程大学 一种单矢量传感器水声多址通信方法
CN104091048B (zh) * 2014-06-19 2017-02-22 哈尔滨工程大学 被动垂直双矢量水听器的水下目标深度测定方法
CN105182345B (zh) * 2015-09-26 2017-10-03 哈尔滨工程大学 基于单矢量潜标的连续谱信号目标自动检测方法
CN105424167A (zh) * 2015-11-09 2016-03-23 哈尔滨工程大学 一种水下自容式声强仪及声强测量方法
CN107422211A (zh) * 2017-08-21 2017-12-01 江苏益邦电力科技有限公司 宽带载波通讯设备的辐射噪声诊断方法
CN108921014B (zh) * 2018-05-21 2021-05-14 西北工业大学 一种基于改进噪声包络信号识别的螺旋桨轴频搜索方法
CN109489799B (zh) * 2018-12-25 2020-07-14 西北工业大学 一种基于双矢量水听器的海底声速分步反演方法
CN109443516B (zh) * 2018-12-25 2020-07-14 西北工业大学 一种基于噪声场垂直振速信号的海底声速被动获取方法
CN110221307B (zh) * 2019-05-28 2022-12-13 哈尔滨工程大学 一种多被动声纳非合作多目标线谱信息融合方法
CN110146846B (zh) * 2019-06-06 2021-04-13 青岛理工大学 一种声源位置估计方法、可读存储介质及计算机设备
CN110596458B (zh) * 2019-07-16 2021-02-02 西北工业大学 Demon谱谐波线谱和基频自动估计方法
CN110716172B (zh) * 2019-10-22 2021-07-09 哈尔滨工程大学 基于频率挑选的矢量水听器包络谱估计方法
CN111735525B (zh) * 2020-05-28 2023-03-31 哈尔滨工程大学 一种适用于无人声纳的demon谱特征提取方法
CN112269962B (zh) * 2020-09-30 2024-05-31 中国船舶重工集团公司七五0试验场 一种辐射噪声线谱稀疏反演方法
CN113887450B (zh) * 2021-10-09 2022-05-20 中国人民解放军91977部队 基于调制谱特征的水中目标在线筛选方法
CN116400337B (zh) * 2023-06-08 2023-08-18 中国人民解放军国防科技大学 基于线段检测的舰船噪声调制线谱提取与轴频估计方法

Also Published As

Publication number Publication date
CN101539597A (zh) 2009-09-23

Similar Documents

Publication Publication Date Title
CN101539597B (zh) 一种分辨具有相同频带辐射噪声的多目标方法
CN104268883B (zh) 一种基于边缘检测的时频谱曲线提取方法
CN106023996B (zh) 基于十字形声阵列宽带波束形成的声识别方法
EP1973101B1 (en) Pitch extraction with inhibition of harmonics and sub-harmonics of the fundamental frequency
CN111665489B (zh) 一种基于目标特性的线谱提取方法
CN110852201A (zh) 一种基于多脉冲包络谱匹配的脉冲信号检测方法
CN106052849A (zh) 一种汽车车内非平稳异常噪声源识别方法
CN110389325B (zh) 一种旋翼无人机的雷达微多普勒信号提取方法
CN110855374B (zh) 一种水声目标辐射噪声调制特征提取方法
RU2718144C1 (ru) Способ классификации, определения координат и параметров движения шумящего в море объекта в инфразвуковом диапазоне частот
CN103344947A (zh) 基于微多普勒效应的微动目标特征提取方法
CN110133580B (zh) 一种随机共振增强的声矢量信号定向方法
CN103176163A (zh) 基于相位模型的舰船线谱噪声源位置识别方法
CN105785346B (zh) 一种基于相位方差加权的未知目标线谱检测方法及系统
CN106100769B (zh) 一种多个不同体制卫星下微弱回波信号联合检测方法
CN103637821A (zh) 基于超声RF时间序列的肝纤维化程度Fisher识别方法
CN110346752A (zh) 基于互质稀疏阵的无模糊测向方法
Huang et al. A practical fundamental frequency extraction algorithm for motion parameters estimation of moving targets
CN106483520B (zh) 一种船舶辐射噪声调制系数估计方法
EP1686561B1 (en) Determination of a common fundamental frequency of harmonic signals
CN113887450B (zh) 基于调制谱特征的水中目标在线筛选方法
CN103267568B (zh) 一种汽车电子控制单元的声音在线检测方法
CN110716172B (zh) 基于频率挑选的矢量水听器包络谱估计方法
Liu et al. Replay attacks detection using phase and magnitude features with various frequency resolutions
Wang et al. Doppler-shift invariant feature extraction for underwater acoustic target classification

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: RUGAO PRODUCTIVITY PROMOTION CENTER

Free format text: FORMER OWNER: HARBIN ENGINEERING UNIV.

Effective date: 20130923

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 150001 HARBIN, HEILONGJIANG PROVINCE TO: 226503 NANTONG, JIANGSU PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20130923

Address after: 226503 Deng yuan community, Rugao Economic Development Zone, Jiangsu 15

Patentee after: Rugao Productivity Promotion Center

Address before: 150001 Heilongjiang, Nangang District, Nantong street, building No. 145, Harbin Engineering University, Intellectual Property Office

Patentee before: Harbin Engineering Univ.

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110126

Termination date: 20180429

CF01 Termination of patent right due to non-payment of annual fee