CN114488100A - 一种齿鲸回声定位单脉冲信号提取方法 - Google Patents

一种齿鲸回声定位单脉冲信号提取方法 Download PDF

Info

Publication number
CN114488100A
CN114488100A CN202210066053.8A CN202210066053A CN114488100A CN 114488100 A CN114488100 A CN 114488100A CN 202210066053 A CN202210066053 A CN 202210066053A CN 114488100 A CN114488100 A CN 114488100A
Authority
CN
China
Prior art keywords
signal
pulse
echo
correlation peak
correlation
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.)
Pending
Application number
CN202210066053.8A
Other languages
English (en)
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.)
Harbin Engineering University
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 CN202210066053.8A priority Critical patent/CN114488100A/zh
Publication of CN114488100A publication Critical patent/CN114488100A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种齿鲸回声定位单脉冲信号提取方法,包括对齿鲸回声定位脉冲串预处理;使用相关法检测实测信号,选定回声定位单脉冲信号作为参考信号,得到脉冲串信号与参考信号的互相关函数。计算互相关函数局部相关峰;多次迭代得到局部相关峰幅值方差满足阈值条件的局部相关峰;通过设置阈值对局部相关峰进行消除重复点、滤除野点,并滤除无信号段产生的判定为相关幅度低的点;利用局部相关峰位置在时域中的对应,对信号进行时域截取,同一脉冲串中单脉冲信号持续时间取为d,选择长为d的时间窗将回声定位单脉冲信号从脉冲串中进行提取。本发明具有运算量小、检测结果准确的优势,能够精确提取时域中回声定位脉冲串中的单脉冲信号位置。

Description

一种齿鲸回声定位单脉冲信号提取方法
技术领域
本发明属于信号处理领域,涉及一种齿鲸回声定位单脉冲信号提取方法。
背景技术
众所周知,齿鲸经历千万年的进化具备非常灵活的生物声呐系统,可以在不同的环境、针对不同的目标实现探测、跟踪、并捕获猎物,对不同的目标与水下环境具有良好的适应性。探究齿鲸生物声呐系统的机理对仿生技术的研究与发展有着重要的指导意义。
齿鲸以脉冲串的形式发射回声定位脉冲信号,通过对回波的处理从而实现各种生理活动。每个脉冲串所包含的脉冲个数并不相等,少的仅仅几个,多的几百个。在目标探测与识别的过程中,脉冲串可以从信号检测理论的角度深入理解。
中国专利说明书CN112712047A中公开了一种基于图像处理的海洋哺乳动物回声定位信号检测方法,根据动物回声定位信号在时频图中近似直线的特征采用图像滤波和直线检测进行处理,并进一步采用随机森林检测模型提升检测准确率,在低信噪比条件下,不仅保持了较高准确率而且有更高的召回率,但该方法是从时频域的角度进行检测,时频分辨力受限,同时计算量较大。中国专利说明书CN106503336A中公开了一种海豚嘀嗒声信号建模与合成的方法,采用瑞利脉冲作为合成海豚嘀嗒声信号的基本单元,设计双脉冲结构模型来建模海豚嘀嗒声信号,有效地实现了不同类型海豚嘀嗒声信号的合成,但并未给出嘀嗒声脉冲串的检测方法。
发明内容
针对上述现有技术,本发明要解决的技术问题是提供一种齿鲸回声定位单脉冲信号提取方法,采用多次迭代的局部相关峰搜索,并且根据脉冲串内参数连续性的生物特征对提取后的漏检与错检情况进行修正。
为解决上述技术问题,本发明的一种齿鲸回声定位单脉冲信号提取方法,包括以下步骤:
步骤一:对齿鲸回声定位脉冲串进行预处理;
步骤二:使用相关法检测实测信号,选定回声定位单脉冲信号作为参考信号f(n),得到采集到的脉冲串信号与参考信号的互相关函数。
步骤三:计算步骤二互相关函数的局部相关峰;
步骤四:判断局部相关峰幅值的方差是否小于设定的阈值T,若满足,则提取最终的局部相关峰,并执行步骤五,否则在局部相关峰中再次计算得到局部相关峰并再次执行步骤四;
步骤五:通过设置阈值对步骤四提取的局部相关峰进行消除重复点、滤除野点的操作,并滤除无信号段产生的判定为相关幅度低的点;
步骤六:利用步骤五获得的局部相关峰位置在时域中的对应,对信号进行时域截取,同一脉冲串中单脉冲信号持续时间取为d,选择长为d的时间窗将回声定位单脉冲信号从脉冲串中进行提取。
进一步的,在步骤六后执行:根据步骤六中得到的局部相关峰位置,计算预处理后的脉冲串内的单脉冲间隔,当脉冲间隔为正常间隔的设定倍数,则判定存在漏检的单脉冲信号,将脉冲间隔修正为若干正常间隔,然后选择长为d的时间窗在脉冲间隔之内提取回声定位单脉冲信号并用于修正步骤六提取结果。
进一步的,预处理包括:将接收到的回声定位信号经过一个f0~f1的带通滤波处理,经过滤波后,消除直流分量与相对于回声定位信号的低频脉冲干扰。
进一步的,互相关函数为:
Figure BDA0003480091520000021
g(n)为待检测的回声定位脉冲串。
本发明的有益效果:海洋哺乳动物的声呐系统能力是目前人工声呐无法比拟的,研究齿鲸的回声定位声呐系统对人工声呐信号设计、仿生学等学科领域有着极大的意义。齿鲸中白鲸的回声定位信号常以脉冲串形式出现,检测并提取脉冲串中的单脉冲信号,是白鲸声呐分析的关键步骤。与前人的单脉冲信号检测方法不同,本发明提出了一种多次迭代搜索相关峰局部最大值的方法,并通过计算相邻脉冲间隔来反馈修正漏检和错检的脉冲信号位置,并取得了很好的效果。本发明针对齿鲸回声定位脉冲串内单脉冲信号提取的问题,提出了多次迭代的局部相关峰搜索方法,并以脉冲串内参数连续性的生物特征为依据,对提取后的漏检与错检情况进行修正。相较于传统的信号检测方法,本发明提出的通过多次迭代筛选局部相关峰方法的数学运算更加简单,思路更加明确,计算量较低;同时,利用单脉冲信号脉冲间隔连续的生物学特性,通过对n次迭代筛选后的局部相关峰位置计算脉冲间隔,并对脉冲间隔不连续变化的前后脉冲信号间隔之内再次进行脉冲检测与提取,实现了对漏检脉冲位置的反馈修正,大大提高了信号检测概率,能够精确提取时域中回声定位脉冲串中的单脉冲信号位置,提升了系统性能。
附图说明
图1为白鲸回声定位脉冲串中单脉冲信号提取方法的过程示意图;
图2为待检测提取单脉冲信号位置的白鲸回声定位脉冲串时域波形图;
图3为白鲸回声定位脉冲串的脉冲间隔计算结果图;
图4为经过脉冲间隔反馈修正后的单脉冲信号提取结果图。
具体实施方式
下面结合说明书附图和具体实施方式对本发明做进一步说明。
结合图1,白鲸回声定位脉冲串内单脉冲信号的提取方法,具体包括:
步骤一:白鲸回声定位脉冲串的预处理;
将接收到的白鲸回声定位信号经过一个f0~f1的带通滤波处理,一方面减少实测场地内的低频环境噪声;另一方面将研究的频率范围限制在f0~f1范围内,对不同种类的回声定位信号进行统一研究频率范围的处理。
经过滤波后,消除直流分量与相对于回声定位信号的低频脉冲干扰,其结果如图2所示。
步骤二:选定单脉冲参考信号,并进行互相关计算;
互相关描述的是两个不同时间序列之间取值的相关程度。互相关函数的定义为
Figure BDA0003480091520000031
式中,x(n)和y(n)均为能量有限的信号,Φx,y(k)为x(n)和y(n)的互相关计算结果。
海豚回声定位信号具有宽带特性,即相关函数具有锐利的相关峰,便于检测信号的初始时刻。本发明使用相关法检测实测信号,利用任意白鲸回声定位单脉冲信号作为计算互相关函数的参考信号,设其中参考信号为f(n),待检测的回声定位脉冲串为g(n),其互相关函数为:
Figure BDA0003480091520000032
根据上式可以获得采集到的脉冲串信号与参考信号的互相关函数,即得到脉冲序列中与参考信号的相关程度。
步骤三:多次迭代筛选局部相关峰;
本发明选择计算步骤二中得到的互相关结果的局部最大值(局部相关峰),并将局部相关峰位置作为判断回声定位单脉冲位置的依据。从时间轴的角度观察步骤二中得到的相关结果(f*g)(τ),在具有高相关计算值的时域位置意味着待测信号该局部与标准信号的相关性强、相似度高,大概率可以将该位置判定为单脉冲存在的位置。
本方法利用回声定位信号幅度连续变化的生物发声规律作为基础,在第一次迭代得到的局部相关峰计算结果中,存在很多的干扰脉冲(如池壁反射、池底反射、其他白鲸回声定位信号等),同时干扰脉冲的相关峰值明显低于相邻的回声定位脉冲产生的相关峰。为消除干扰脉冲的影响,通过对相关峰再次求解局部最大值以过滤掉不满足幅值连续变化的相关峰,反复筛选相关峰获得回声定位脉冲的真实单脉冲起始位置。
在经过n次的迭代筛选后,相邻局部相关峰幅值的方差小于设定的阈值T,说明非感兴趣信号干扰已经被基本抑制,同时帧内单脉冲信号的起始位置被初步提取出。
步骤四:设定阈值过滤,并提取单脉冲信号位置;
经过步骤三的n次迭代筛选局部相关峰,剩下的相关峰中仍旧存在极少数的低相关性噪点。通过设置合适的阈值,对第n次筛选结果做消除重复点、滤除野点的操作,并滤除无信号段产生的个别相关幅度较低的点,消除剩余低相关性脉冲噪声的干扰。
利用获得的局部相关峰位置在时域中的对应,对信号进行时域截取。同一脉冲串中单脉冲信号持续时间近似相等,取为d。选择长为d的时间窗将回声定位脉冲从脉冲串中进行提取,足以覆盖单脉冲信号的完整时间脉宽。
步骤五:计算回声定位脉冲信号参数,反馈修正单脉冲信号提取过程;
步骤三中的多次迭代滤波,导致提取出的脉冲串相关峰存在一定的漏检情况,进而导致脉冲串中的单脉冲间隔出现整数倍的误差。利用回声定位脉冲串脉冲间隔连续变化的生物特性,在单脉冲检测器的后端加上一个基于脉冲间隔反馈修正的部分。在时域上,对帧内脉冲间隔不连续变化的前后单脉冲信号间隔之内再次进行脉冲检测与提取,以修正漏检的问题。
下面结合具体参数给出实施例:
结合图1,白鲸回声定位脉冲串内单脉冲信号的提取方法,具体包括:
步骤一:白鲸回声定位脉冲串的预处理;
将接收到的白鲸回声定位信号经过一个10-150kHz的带通滤波处理,一方面减少实测场地内的低频环境噪声;另一方面将研究的频率范围限制在10-150kHz范围内,对不同种类的回声定位信号进行统一研究频率范围的处理。
经过滤波后,消除直流分量与相对于回声定位信号的低频脉冲干扰。
步骤二:选定单脉冲参考信号,并进行互相关计算;
互相关描述的是两个不同时间序列之间取值的相关程度。互相关函数的定义为
Figure BDA0003480091520000041
式中,x(n)和y(n)均为能量有限的信号,Φx,y(k)为x(n)和y(n)的互相关计算结果。
海豚回声定位信号具有宽带特性,即相关函数具有锐利的相关峰,便于检测信号的初始时刻。本发明使用相关法检测实测信号,利用任意白鲸回声定位单脉冲信号作为计算互相关函数的参考信号,设其中参考信号为f(n),待检测的回声定位脉冲串为g(n),其互相关函数为:
Figure BDA0003480091520000051
根据上式可以获得采集到的脉冲串信号与参考信号的互相关函数,即得到脉冲序列中与参考信号的相关程度。
步骤三:多次迭代筛选局部相关峰;
本发明选择计算步骤二中得到的互相关结果的局部最大值(局部相关峰),并将局部相关峰位置作为判断回声定位单脉冲位置的依据。从时间轴的角度观察步骤二中得到的相关结果(f*g)(τ),在具有高相关计算值的时域位置意味着待测信号该局部与标准信号的相关性强、相似度高,大概率可以将该位置判定为单脉冲存在的位置。
本方法利用回声定位信号幅度连续变化的生物发声规律作为基础,在第一次迭代得到的局部相关峰计算结果中,存在很多的干扰脉冲(如池壁反射、池底反射、其他白鲸回声定位信号等),同时干扰脉冲的相关峰值明显低于相邻的回声定位脉冲产生的相关峰。为消除干扰脉冲的影响,通过对相关峰再次求解局部最大值以过滤掉不满足幅值连续变化的相关峰,反复筛选相关峰获得回声定位脉冲的真实单脉冲起始位置。
在经过6次的迭代筛选后,相邻局部相关峰幅值的方差小于设定的阈值T,说明非感兴趣信号干扰已经被基本抑制,同时帧内单脉冲信号的起始位置被初步提取出。
步骤四:设定阈值过滤,并提取单脉冲信号位置;
经过步骤三的6次迭代筛选局部相关峰,剩下的相关峰中仍旧存在极少数的低相关性噪点。通过设置合适的阈值,对第6次筛选结果做消除重复点、滤除野点的操作,并滤除无信号段产生的个别相关幅度较低的点,消除剩余低相关性脉冲噪声的干扰。
利用获得的局部相关峰位置在时域中的对应,对信号进行时域截取。同一脉冲串中单脉冲信号持续时间近似相等,取为100微秒。选择长为100微秒的时间窗将回声定位脉冲从脉冲串中进行提取,足以覆盖单脉冲信号的完整时间脉宽。
步骤五:计算回声定位脉冲信号参数,反馈修正单脉冲信号提取过程;
步骤三中的多次迭代滤波,导致提取出的脉冲串相关峰存在一定的漏检情况,进而导致脉冲串中的单脉冲间隔出现整数倍的误差。利用回声定位脉冲串脉冲间隔连续变化的生物特性,在单脉冲检测器的后端加上一个基于脉冲间隔反馈修正的部分。具体而言,根据步骤四中得到的相关峰位置,计算脉冲串内的单脉冲间隔参数。其结果中包含少量二倍或三倍于正常间隔范围的野点,这显然不满足脉冲间隔连续变化的生物特性,属于漏检的单脉冲信号。将这些野点对应的脉冲间隔除以二或三,即可得到连续变化的脉冲间隔,如图3所示。
根据修正后的脉冲间隔,在步骤四中相邻单脉冲位置之间再次提取漏检的信号,提高检测的准确性,单脉冲信号的最终提取结果如图4所示。

Claims (4)

1.一种齿鲸回声定位单脉冲信号提取方法,其特征在于,包括以下步骤:
步骤一:对齿鲸回声定位脉冲串进行预处理;
步骤二:使用相关法检测实测信号,选定回声定位单脉冲信号作为参考信号f(n),得到采集到的脉冲串信号与参考信号的互相关函数。
步骤三:计算步骤二互相关函数的局部相关峰;
步骤四:判断局部相关峰幅值的方差是否小于设定的阈值T,若满足,则提取最终的局部相关峰,并执行步骤五,否则在局部相关峰中再次计算得到局部相关峰并再次执行步骤四;
步骤五:通过设置阈值对步骤四提取的局部相关峰进行消除重复点、滤除野点的操作,并滤除无信号段产生的判定为相关幅度低的点;
步骤六:利用步骤五获得的局部相关峰位置在时域中的对应,对信号进行时域截取,同一脉冲串中单脉冲信号持续时间取为d,选择长为d的时间窗将回声定位单脉冲信号从脉冲串中进行提取。
2.根据权利要求1所述的一种齿鲸回声定位单脉冲信号提取方法,其特征在于:在步骤六后执行:根据步骤六中得到的局部相关峰位置,计算预处理后的脉冲串内的单脉冲间隔,当脉冲间隔为正常间隔的设定倍数,则判定存在漏检的单脉冲信号,将所述脉冲间隔修正为若干正常间隔,然后选择长为d的时间窗在所述脉冲间隔之内提取回声定位单脉冲信号并用于修正步骤六提取结果。
3.根据权利要求1所述的一种齿鲸回声定位单脉冲信号提取方法,其特征在于:所述预处理包括:将接收到的回声定位信号经过一个f0~f1的带通滤波处理,经过滤波后,消除直流分量与相对于回声定位信号的低频脉冲干扰。
4.根据权利要求1所述的一种齿鲸回声定位单脉冲信号提取方法,其特征在于:所述互相关函数为:
Figure FDA0003480091510000011
g(n)为待检测的回声定位脉冲串。
CN202210066053.8A 2022-01-20 2022-01-20 一种齿鲸回声定位单脉冲信号提取方法 Pending CN114488100A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210066053.8A CN114488100A (zh) 2022-01-20 2022-01-20 一种齿鲸回声定位单脉冲信号提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210066053.8A CN114488100A (zh) 2022-01-20 2022-01-20 一种齿鲸回声定位单脉冲信号提取方法

Publications (1)

Publication Number Publication Date
CN114488100A true CN114488100A (zh) 2022-05-13

Family

ID=81472202

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210066053.8A Pending CN114488100A (zh) 2022-01-20 2022-01-20 一种齿鲸回声定位单脉冲信号提取方法

Country Status (1)

Country Link
CN (1) CN114488100A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115078845A (zh) * 2022-07-26 2022-09-20 深圳华创电科技术有限公司 脉冲重复周期的采集方法、设备及计算机存储介质
CN115150009A (zh) * 2022-06-24 2022-10-04 自然资源部第三海洋研究所 一种应用于pam系统的海豚回声定位信号检测方法
CN116660876A (zh) * 2023-05-26 2023-08-29 吉林大学 一种中华白海豚回声定位信号自动检测和定位方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102253361A (zh) * 2011-06-10 2011-11-23 哈尔滨工程大学 一种超短基线垂直运动多目标信号检测方法
CN105676169A (zh) * 2016-01-19 2016-06-15 中国电子科技集团公司第三研究所 一种脉冲声目标的定向方法及装置
CN111175729A (zh) * 2020-01-18 2020-05-19 中国科学院水生生物研究所 一种基于齿鲸类高频声呐信号的实时在线监测预警系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102253361A (zh) * 2011-06-10 2011-11-23 哈尔滨工程大学 一种超短基线垂直运动多目标信号检测方法
CN105676169A (zh) * 2016-01-19 2016-06-15 中国电子科技集团公司第三研究所 一种脉冲声目标的定向方法及装置
CN111175729A (zh) * 2020-01-18 2020-05-19 中国科学院水生生物研究所 一种基于齿鲸类高频声呐信号的实时在线监测预警系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘凇佐;乔钢;尹艳玲;: "一种利用海豚叫声的仿生水声通信方法", 物理学报, no. 14, 23 July 2013 (2013-07-23) *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115150009A (zh) * 2022-06-24 2022-10-04 自然资源部第三海洋研究所 一种应用于pam系统的海豚回声定位信号检测方法
CN115078845A (zh) * 2022-07-26 2022-09-20 深圳华创电科技术有限公司 脉冲重复周期的采集方法、设备及计算机存储介质
CN115078845B (zh) * 2022-07-26 2022-11-29 深圳华创电科技术有限公司 脉冲重复周期的采集方法、设备及计算机存储介质
CN116660876A (zh) * 2023-05-26 2023-08-29 吉林大学 一种中华白海豚回声定位信号自动检测和定位方法

Similar Documents

Publication Publication Date Title
CN114488100A (zh) 一种齿鲸回声定位单脉冲信号提取方法
CN110865357B (zh) 一种基于参数优化vmd的激光雷达回波信号降噪方法
CN104569948B (zh) 海杂波背景下子带自适应glrt‑ltd检测方法
CN110852201B (zh) 一种基于多脉冲包络谱匹配的脉冲信号检测方法
CN106154332B (zh) 一种井中微地震纵横波事件初至识别方法
CN113221631B (zh) 基于卷积神经网络的序列脉冲抗干扰目标检测方法
CN110988834B (zh) 一种基于自适应阈值双门限的脉冲到达时间测量方法
CN110879389B (zh) 基于多基地ir-uwb生物雷达信号的多人体目标识别定位方法
CN110879388B (zh) 基于ir-uwb生物雷达信号的人与动物非接触探测区分方法
US11525889B2 (en) Non-contact method for detectiing and distinguishing human and animal based on IR-UWB bio-radar signal
CN111024207B (zh) 一种矢量水听器线谱自动检测与判决方法
CN106772299B (zh) 一种基于距离匹配的pd雷达微弱目标动态规划检测方法
CN104483668A (zh) 一种高精度雷达信号检测和跟踪系统及其方法
Rose A review of problems and new directions in the application of fisheries acoustics on the Canadian East Coast
CN112578435B (zh) 岩石超声测试初至拾取方法及系统
CN106443792A (zh) 一种水下多缆声学网络定位方法
CN112462321A (zh) 基于脉冲匹配累加的声信标信号截获处理方法
CN117153168A (zh) 鲸类回声定位信号的检测与提取方法、系统、设备及介质
McCauley et al. Migratory patterns and estimated population size of pygmy blue whales (Balaenoptera musculus brevicauda) traversing the Western Australian coast based on passive acoustics.
CN111413702A (zh) 用于宽带探鱼仪的高效目标分割方法
CN105548994B (zh) 一种基于压缩感知的水下多亮点目标的检测方法
CN111812608B (zh) 基于mtd脉冲积累及模态分解的雷达目标方位角估计方法
Merkens et al. Comparison of high-frequency echolocation clicks (likely Kogia) in two simultaneously collected passive acoustic data sets sampled at 200 kHz and 320 kHz
CN116502152A (zh) 雷达天线扫描方式的识别方法、装置、设备和存储介质
Hamilton et al. Separating overlapping echolocation: An updated method for estimating the number of echolocating animals in high background noise levels

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