CN103197278A - 一种基于多普勒频移变化率的舰船线谱噪声源定位方法 - Google Patents
一种基于多普勒频移变化率的舰船线谱噪声源定位方法 Download PDFInfo
- Publication number
- CN103197278A CN103197278A CN2013100685517A CN201310068551A CN103197278A CN 103197278 A CN103197278 A CN 103197278A CN 2013100685517 A CN2013100685517 A CN 2013100685517A CN 201310068551 A CN201310068551 A CN 201310068551A CN 103197278 A CN103197278 A CN 103197278A
- Authority
- CN
- China
- Prior art keywords
- signal
- noise source
- doppler
- omega
- 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.)
- Granted
Links
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种基于多普勒频移变化率的舰船线谱噪声源定位方法,在现有的单个水听器测量舰船辐射噪声的平台上,利用测量时舰船与测量水听器之间相对运动产生的多普勒效应,结合首次提出的L-相位分布,提取多普勒信号的多普勒频移变化率,实现了舰船低频线谱噪声源的定位,该方法实施简单,测量精度比通过特性法要高。
Description
技术领域
本发明是一种舰船线谱噪声源定位方法,涉及信号处理、水声测量等领域。
背景技术
舰船辐射噪声的总声级取决于所有辐射噪声源,主要的噪声源来源于舰船的主机推进系统、螺旋桨、辅机机械系统和航行时产生的水动力效应。其中,由辅机机械的往复运动引起的线谱噪声具有频率低、功率高和稳定性强等特点,其携带着舰船重要的特征信息,可被远距离探测,是水声设备探测舰船的关键。降低舰船噪声不仅可以有效地减小舰船航行对海洋生物的噪声影响,同时也有利于提高装备在舰船上的各种水声设备的性能。因此,舰船降噪成为舰船声学设计的迫切任务,而精确识别舰船上主要噪声源的位置为舰船的声学设计提供了参考,同时其也可为探寻研究舰船内部机械的振动机理提供重要依据。
以往对舰船辐射噪声的测量主要利用其通过特性来实现。在进行辐射噪声测量时,被测舰船做匀速直线航行,由远及近接近测量水听器并由近及远离开测量水听器,在这个过程中测量水听器记录被测舰船辐射噪声的声压级(或功率谱级)与测量水听器相对位置的对应关系,一般称之为舰船辐射噪声的通过特性。利用水下测距信息与辐射噪声同步记录的方法可得到舰船部位与辐射噪声信号对应的通过特性曲线。舰船辐射噪声通过特性包括总声级通过特性,1/3oct频带级通过特性和线谱通过特性。总声级通过特性曲线的最大值对应于舰船总声级最大的位置,而1/3oct频带声压级和线谱通过特性则反映不同频率的噪声级和被测舰船不同位置的关系。把测量水听器与船体上某一个参考位置的关系确定,那么在空间(或时间)上,测量数据与船体部位就形成了一一对应关系,它反映了舰船通过测量水听器时,声压沿艇体的分布。通过对纵向通过特性分析,可得到不同频带噪声级与艇体部位的对应关系,为确定舰船的噪声源位置提供了一种方法。利用舰船通过特性来进行舰船线谱噪声源定位方法,测试简单但是噪声源定位精度不理想。
徐灵基(Xu Ling-ji,Yang Yi-xin,and Tian Feng.Doppler method for identification ofnoise sources on underwater moving target[C].In:Asia-Pacific Signal and InformationProcessing Association Annual Summit and Conference2011,Xi’an,2011,83-88.)提出利用多普勒频移方法来识别水下运动目标低频线谱声源的位置,取得了一定的效果。时洁(时洁,杨德森,时胜国.基于最差性能优化的运动声源稳健聚焦定位识别方法研究[J].物理学报,2011,60(6):1-11.)、王志伟(王志伟,徐灵基,杨益新等.直线阵识别水下运动目标噪声源方法研究[J].振动与冲击,2012;31(19):118-122.)利用声聚焦阵列的高分辨技术来测量水下运动目标的噪声源,定位精度较高。然而对于舰船上的低频线谱声源,上述方法需要的阵列孔径太大,且对阵列布放的形状要求严格,很难应用于实际的舰船辐射噪声测试中。
发明内容
为了克服现有技术的不足,本发明采用的舰船通过特性法对舰船线谱噪声源定位精度低的不足,本发明利用被测舰船与测量水听器相对运动引起的多普勒频移变化率特性,在仅使用单个水听器的测量条件下,实现对舰船线谱噪声源的定位。
本发明解决其技术问题所采用的技术方案包括以下步骤:
1)在舰船上任意固定位置Lc处布放一个已知的水声同步定位仪作为信标源,舰船以速度ν做匀速直线航行,由远及近接近测量水听器并由近及远离开测量水听器,测量水听器与舰船运动轨迹的正横距离为R0,在这个过程中测量水听器接收舰船辐射的声信号,并将其转换成电压信号s(t),0<t≤T,t表示时间,T表示数据记录的时间长度;
2)把采集得到的水听器接收信号s(t)做快速傅立叶变换分析频谱,找到信号中存在的低频线谱,其频率选择为m=1,2,…,M,M表示低频线谱的个数,通过带通滤波器分离出存在线谱的频带宽度为各个小区域信号,得到各个低频线谱噪声源多普勒信号sm(t),将分离出的线谱噪声源多普勒信号sm(t)复数化得到多普勒解析信号zm(t),其中zm(t)=sm(t)+jH[sm(t)],H[sm(t)]表示信号sm(t)的希尔伯特变换;
3)多普勒信号zm(t)瞬时频率的变化率称作多普勒频移变化率,表达式为 fm0表示舰船上各辐射线谱噪声源信号的频率,tm0为各线谱噪声源通过正横位置离测量水听器最近时的正横时刻;利用L-相位分布提取多普勒信号zm(t)的多普勒频移变化率,多普勒信号zm(t)的LPD为: 其中Ω为频率变化率,τ为时间延迟,L为常数,取为大于等于8的整数;通过一维峰值搜索得到信号zm(t)的多普勒频移变化率
4)获得已知位置Lc信标源通过测量水听器时的正横时刻tc0、舰船运动轨迹的正横距离R0,根据步骤(3)中信号zm(t)的多普勒频移变化率的表达式,利用非线性最小二乘估计法得到舰船上各线谱噪声源通过测量水听器时的正横时刻和线谱声源的频率估计值
本发明的有益效果是:在现有的单个水听器测量舰船辐射噪声的平台上,利用测量时舰船与测量水听器之间相对运动产生的多普勒效应,结合首次提出的L-相位分布,提取多普勒信号的多普勒频移变化率,实现了舰船低频线谱噪声源的定位,该方法实施简单,测量精度比通过特性法要高。
附图说明
图1是舰船辐射噪声测量模型示意图;
图2是舰船线谱噪声源定位方法的总体流程框图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
在对舰船辐射噪声进行测量时,测量海区可人为确定,一般选择海面较平静、水深较深的海区,这样可以降低水声信道对测量的影响,水中声速c可测,为已知量。
本发明包括以下步骤:
1)在舰船上任意固定位置Lc处布放一个已知的水声同步定位仪作为信标源,舰船以速度ν做匀速直线航行,由远及近接近测量水听器并由近及远离开测量水听器,测量水听器与舰船运动轨迹的正横距离为R0,在这个过程中测量水听器接收舰船辐射的声信号,并将其转换成电压信号s(t),0<t≤T,通过前放后用数据采集仪记录下来,T表示数据记录的时间长度。
2)由于被测舰船与测量水听器之间存在相对运动,会形成多普勒效应,测量水听器接收到的线谱噪声源信号称为多普勒信号。接收水听器记录的信号需要预处理,目的为:一是确定存在的线谱噪声源信号;二是滤波去除多普勒信号的带外噪声;三是复数化得到解析形式的多普勒信号。先把采集得到的水听器接收信号s(t)做快速傅立叶变换(FFT)分析频谱,找到信号中存在的低频线谱,其频率一般选择为m=1,2,…,M,M表示低频线谱的个数,即待分析的低频线谱噪声源个数。通过带通滤波器分离出存在线谱的频带宽度为各个小区域信号,其中v可通过水声同步定位仪获得,这便得到各个低频线谱噪声源多普勒信号sm(t),m=1,2,…,M,t表示时间。将分离出的线谱噪声源多普勒信号sm(t)复数化得到多普勒解析信号zm(t),其中zm(t)=sm(t)+jH[sm(t)],H[sm(t)]表示信号sm(t)的希尔伯特变换。
3)多普勒信号zm(t)的瞬时频率是变化的,其瞬时频率的变化率称作多普勒频移变化率,表达式为 fm0表示舰船上各辐射线谱噪声源信号的频率,tm0为各线谱噪声源通过正横位置离测量水听器最近时的正横时刻。根据Stone-Weierstrass理论,多普勒信号zm(t)可以表示为高阶多项式相位信号。由于多普勒信号近似的高阶相位信号的阶数一般都大于3,可以利用L-相位分布(L-phase distribution,LPD)提取多普勒信号zm(t)的多普勒频移变化率。多普勒信号zm(t)的LPD为: 其中t为时间,Ω为频率变化率,τ为时间延迟,L为常数,一般取为大于等于8的整数。在得到多普勒信号的LPD后,通过一维峰值搜索可以估计得到信号zm(t)的多普勒频移变化率
4)水声同步定位仪可以获得已知位置Lc信标源通过测量水听器时的正横时刻tc0、舰船运动轨迹的正横距离R0,根据步骤(3)中信号zm(t)的多普勒频移变化率的表达式,利用非线性最小二乘估计法得到舰船上各线谱噪声源通过测量水听器时的正横时刻和线谱声源的频率估计值
图1为舰船辐射噪声的测量模型,测量水听器位于S处,被测舰船匀速运动,速度为ν,轨迹为MN,正横位置O为其运行轨迹上离测量水听器距离最近点,在舰船上固定位置Lc处布放一个水声同步定位仪作为信标,舰船通过测量水听器的过程中,测量水听器接收舰船辐射的声信号,并将其转换成电压信号s(t),t为时间,通过前置放大后用数据采集仪以采样率FS记录时间长度为T的数据,得到信号s(t)的离散时间序列s(n),n=1,…,N,N=TFS,n与t的对应关系为t=nFS。
图2为舰船线谱噪声源定位方法的总体流程,具体实施如下:
1)对水听器接收信号预处理,先将采集到的信号s(n)做快速傅立叶变换(FFT)分析频谱,确定信号中存在的低频线谱,其频率一般选择为m=1,2,…,M,M表示低频线谱的个数,即待分析的低频线谱噪声源个数。通过FIR带通滤波器分离出存在线谱且频带宽度为的各个小区域信号,其中v可通过水声同步定位仪获得,这便得到各个线谱噪声源多普勒信号sm(n),m=1,2,…,M。将分离出的线谱噪声源多普勒信号sm(n)复数化得到zm(n),其中zm(n)=sm(n)+jH[sm(n)],H[sm(n)]表示信号sm(n)的希尔伯特变换。
2)利用LPD的离散形式,获得多普勒解析信号zm(n)的时间-频率变化率分布: θ为时间延迟点数,Ω为频率变化率,L一般取为大于等于8的整数,时间序列点n′=0.25FS:0.25FS:0.25KFS,K是T0.25向下取整的整数值,即分析的时间间隔为0.25秒。再通过一维峰值搜索得到多普勒解析信号zm(n)的瞬时频率变化率曲线: 其中t′=n′FS,共K个值。
3)水声同步定位仪可以获得已知位置Lc信标源通过测量水听器时的正横时刻tc0、舰船运动轨迹的正横距离R0。多普勒频移变化率的模型为 通过非线性最小二乘估计法可以得到各噪声源通过测量水听器时正横时刻的估计值 k=1,2,…,K为时间t′的标记序号。令向量变量x=[tm0;fm0],假设最小化目标函数 下面给出和的具体计算步骤如下:(a)对提取出的多普勒频移变化率曲线进行一维峰值搜索得到正横时刻的初值(b)对于步骤(a)给定的确定线谱噪声源频率的初值 (c)容易获得函数Q(x)对向量变量x的一阶偏导数为再利用Levenberg-Marquardt迭代算法求解得到参数估计值和其迭代关系式为 x(m)表示第m次迭代变量x的值,μ为收敛因子,I为单位矩阵,迭代初始值为步骤(a)、(b)所得。
Claims (1)
1.一种基于多普勒频移变化率的舰船线谱噪声源定位方法,其特征在于包括下述步骤:
1)在舰船上任意固定位置Lc处布放一个已知的水声同步定位仪作为信标源,舰船以速度ν做匀速直线航行,由远及近接近测量水听器并由近及远离开测量水听器,测量水听器与舰船运动轨迹的正横距离为R0,在这个过程中测量水听器接收舰船辐射的声信号,并将其转换成电压信号s(t),0<t≤T,t表示时间,T表示数据记录的时间长度;
2)把采集得到的水听器接收信号s(t)做快速傅立叶变换分析频谱,找到信号中存在的低频线谱,其频率选择为m=1,2,…,M,M表示低频线谱的个数,通过带通滤波器分离出存在线谱的频带宽度为各个小区域信号,得到各个低频线谱噪声源多普勒信号sm(t),将分离出的线谱噪声源多普勒信号sm(t)复数化得到多普勒解析信号zm(t),其中zm(t)=sm(t)+jH[sm(t)],H[sm(t)]表示信号sm(t)的希尔伯特变换;
3)多普勒信号zm(t)瞬时频率的变化率称作多普勒频移变化率,表达式为 fm0表示舰船上各辐射线谱噪声源信号的频率,tm0为各线谱噪声源通过正横位置离测量水听器最近时的正横时刻;利用L-相位分布提取多普勒信号zm(t)的多普勒频移变化率,多普勒信号zm(t)的LPD为: 其中Ω为频率变化率,τ为时间延迟,L为常数,取为大于等于8的整数;通过一维峰值搜索得到信号zm(t)的多普勒频移变化率
4)获得已知位置Lc信标源通过测量水听器时的正横时刻tc0、舰船运动轨迹的正横距离R0,根据步骤(3)中信号zm(t)的多普勒频移变化率的表达式,利用非线性最小二乘估计法得到舰船上各线谱噪声源通过测量水听器时的正横时刻和线谱声源的频率估计值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310068551.7A CN103197278B (zh) | 2013-03-04 | 2013-03-04 | 一种基于多普勒频移变化率的舰船线谱噪声源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310068551.7A CN103197278B (zh) | 2013-03-04 | 2013-03-04 | 一种基于多普勒频移变化率的舰船线谱噪声源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103197278A true CN103197278A (zh) | 2013-07-10 |
CN103197278B CN103197278B (zh) | 2014-12-31 |
Family
ID=48719955
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310068551.7A Active CN103197278B (zh) | 2013-03-04 | 2013-03-04 | 一种基于多普勒频移变化率的舰船线谱噪声源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103197278B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777453A (zh) * | 2015-04-23 | 2015-07-15 | 西北工业大学 | 舰船线谱噪声源定位的波束域时频分析方法 |
CN105589066A (zh) * | 2015-12-14 | 2016-05-18 | 西北工业大学 | 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法 |
CN105807272A (zh) * | 2016-04-29 | 2016-07-27 | 中国科学院声学研究所 | 一种基于多普勒效应的电动水下航行器航速测量方法及装置 |
KR20170054175A (ko) * | 2015-11-09 | 2017-05-17 | 소나테크 주식회사 | 미식별 수중신호 위치의 추정 방법 및 시스템 |
KR101783822B1 (ko) * | 2016-03-09 | 2017-10-10 | 서울대학교산학협력단 | 수상 이동체 내의 소음원 위치 추정 방법 |
CN110068331A (zh) * | 2018-01-24 | 2019-07-30 | 北京致感致联科技有限公司 | 水下导航定位设备及系统 |
CN111679296A (zh) * | 2020-05-28 | 2020-09-18 | 哈尔滨工业大学 | 一种基于gp迭代外推的单星定位方法 |
CN113514145A (zh) * | 2020-04-09 | 2021-10-19 | 中国船舶重工集团公司第七六0研究所 | 水面船舶推进系统与辅机系统水下噪声贡献分离测试方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0253277A2 (de) * | 1986-07-12 | 1988-01-20 | STN ATLAS Elektronik GmbH | Passives Verfahren zur Schätzung von Zieldaten eines im Wasser sich bewegenden, zeitkontinuierliche Wasserschallsignale abstrahlenden Ziels |
DE19857760C1 (de) * | 1998-12-15 | 2000-03-23 | Stn Atlas Elektronik Gmbh | Verfahren zur passiven akustischen Peilung eines Schall ins Wasser abstrahlenden Ziels |
CN102243302A (zh) * | 2011-04-15 | 2011-11-16 | 东南大学 | 一种水声目标辐射噪声的线谱时间累积特征提取方法 |
-
2013
- 2013-03-04 CN CN201310068551.7A patent/CN103197278B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0253277A2 (de) * | 1986-07-12 | 1988-01-20 | STN ATLAS Elektronik GmbH | Passives Verfahren zur Schätzung von Zieldaten eines im Wasser sich bewegenden, zeitkontinuierliche Wasserschallsignale abstrahlenden Ziels |
DE19857760C1 (de) * | 1998-12-15 | 2000-03-23 | Stn Atlas Elektronik Gmbh | Verfahren zur passiven akustischen Peilung eines Schall ins Wasser abstrahlenden Ziels |
CN102243302A (zh) * | 2011-04-15 | 2011-11-16 | 东南大学 | 一种水声目标辐射噪声的线谱时间累积特征提取方法 |
Non-Patent Citations (4)
Title |
---|
LINGJI XU等: "Doppler Method for Identification of Noise Sources on Underwater Moving Target", 《APSIPA ASC 2011 XI’AN》 * |
Y. T. CHAN等: "Passive Localization from Doppler-Shifted Frequency Measurements", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
徐灵基等: "舰船辐射噪声线谱多普勒分析的ESPRIT方法", 《声学技术》 * |
王志伟等: "直线阵识别水下运动目标噪声源方法研究", 《振动与冲击》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777453A (zh) * | 2015-04-23 | 2015-07-15 | 西北工业大学 | 舰船线谱噪声源定位的波束域时频分析方法 |
KR20170054175A (ko) * | 2015-11-09 | 2017-05-17 | 소나테크 주식회사 | 미식별 수중신호 위치의 추정 방법 및 시스템 |
KR102374304B1 (ko) | 2015-11-09 | 2022-03-16 | 대우조선해양 주식회사 | 미식별 수중신호 위치의 추정 방법 및 시스템 |
CN105589066A (zh) * | 2015-12-14 | 2016-05-18 | 西北工业大学 | 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法 |
CN105589066B (zh) * | 2015-12-14 | 2018-01-09 | 西北工业大学 | 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法 |
KR101783822B1 (ko) * | 2016-03-09 | 2017-10-10 | 서울대학교산학협력단 | 수상 이동체 내의 소음원 위치 추정 방법 |
CN105807272A (zh) * | 2016-04-29 | 2016-07-27 | 中国科学院声学研究所 | 一种基于多普勒效应的电动水下航行器航速测量方法及装置 |
CN105807272B (zh) * | 2016-04-29 | 2018-04-17 | 中国科学院声学研究所 | 一种基于多普勒效应的电动水下航行器航速测量方法及装置 |
CN110068331A (zh) * | 2018-01-24 | 2019-07-30 | 北京致感致联科技有限公司 | 水下导航定位设备及系统 |
CN113514145A (zh) * | 2020-04-09 | 2021-10-19 | 中国船舶重工集团公司第七六0研究所 | 水面船舶推进系统与辅机系统水下噪声贡献分离测试方法 |
CN111679296A (zh) * | 2020-05-28 | 2020-09-18 | 哈尔滨工业大学 | 一种基于gp迭代外推的单星定位方法 |
CN111679296B (zh) * | 2020-05-28 | 2023-04-25 | 哈尔滨工业大学 | 一种基于gp迭代外推的单星定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103197278B (zh) | 2014-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103197278B (zh) | 一种基于多普勒频移变化率的舰船线谱噪声源定位方法 | |
CN103176163B (zh) | 基于相位模型的舰船线谱噪声源位置识别方法 | |
CN103293515B (zh) | 一种舰船线谱噪声源纵向分布特性测量方法 | |
CN105589066B (zh) | 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法 | |
CN103076594B (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
CN103344961B (zh) | 舰船速度和距离联合测量的被动声多普勒相位方法 | |
CN103675819A (zh) | 一种可用于被动合成孔径阵列偏航的目标检测方法及系统 | |
CN103529441A (zh) | 一种被动合成孔径目标信号检测和分辨方法及系统 | |
CN104407340A (zh) | 拖曳线列阵阵形标定装置及方法 | |
CN105204026A (zh) | 一种基于声场干涉条纹的单水平阵被动测速和测距装置及其方法 | |
CN103076590A (zh) | 一种基于频率预估的水声脉冲信号的定位方法 | |
CN113687308B (zh) | 一种基于弯曲波的冰上震源定位方法 | |
CN115166817A (zh) | 一种基于冰层模态群慢度差特征的冰声定位方法 | |
CN108318876A (zh) | 一种利用单个水听器估计水下目标深度和距离的方法 | |
RU103616U1 (ru) | Устройство определения местоположения течи продукта из подводных трубопроводов | |
CN116908855A (zh) | 一种极地水下脉冲简正波声学测距方法 | |
CN114757240B (zh) | 一种水声目标运动参数联合估计方法 | |
CN108646248A (zh) | 一种针对低速运动声源的被动声学测速测距方法 | |
KR101446439B1 (ko) | Fmcw 레이더의 고정밀 주파수 추정 방법 및 고정밀 주파수 추정 장치 | |
JP5708018B2 (ja) | アクティブソーナー装置 | |
CN115236592B (zh) | 一种基于单阵元时频曲线匹配的冰声定位方法 | |
CN117169816B (zh) | 一种深海声影区宽带声源被动定位方法、介质及系统 | |
JP6610224B2 (ja) | バイスタティックアクティブソーナー装置およびその受信器 | |
Yang et al. | Study on line-spectrum stability of ship underwater radiated noise |
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 |