CN106595835B - 一种基于错位叠加算法的声信标信号提取方法及系统 - Google Patents
一种基于错位叠加算法的声信标信号提取方法及系统 Download PDFInfo
- Publication number
- CN106595835B CN106595835B CN201611104475.0A CN201611104475A CN106595835B CN 106595835 B CN106595835 B CN 106595835B CN 201611104475 A CN201611104475 A CN 201611104475A CN 106595835 B CN106595835 B CN 106595835B
- Authority
- CN
- China
- Prior art keywords
- signal
- superposition
- acoustic marker
- algorithm
- staggered
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于错位叠加算法的声信标信号提取方法,包括以下步骤:搭建声信标信号采集系统;对采集的原始信号随机选取叠加起始点;根据按相邻周期(或相同时间间隔)叠加的方式确定各截取信号段的长度;将各截取的信号段进行线性叠加,以实现对声信标信号的提取;将该算法编写程序导入单片机中,完成声信标信号提取系统的搭建。本发明采用的错位叠加算法提取的声信标信号无需进行时频转换,方法简便易行,提取的信号波形准确性高。
Description
技术领域
本发明涉及一种信号处理领域的信号提取方法及系统,特别是一种基于错位叠加算法的声信标信号提取方法和相应的系统。
背景技术
当飞机在海上失事时,为了能快速找到飞行数据记录器(即俗称的黑匣子),飞行数据记录器上都装有声信标。当飞行数据记录器落入水中后,声信标便自动激活,开始向水中发出特定信号形式的声信号,供搜寻设备对其进行声学定位。声信标在海底工作,其发射出来的声波在海洋中传播时,随着距离的增加会被减弱(即传播损失);海洋中的噪声不受人们控制而依赖于海洋环境本身(即海洋环境噪声),其产生的原因有海水中分子的热运动、海洋生物、船只航行以及如地震、风暴、下雨的自然现象等。由于传播损失和海洋环境噪声的干扰,必须在信号处理手段上进行优化。
针对确定的被提取信号以及检测频带近似白噪声的海洋背景干扰特性,常用的信号提取方法主要有:人耳听闻、窄带滤波、匹配滤波、能量检测、FFT滤波及自适应滤波等。在水声信号处理领域中,短时傅里叶变换、小波变换、Hilbert-Huang变换是最常用的三种方法。1946年Gabor提出了短时傅里叶变换(STFT),通过特定的平移窗函数来分解信号的频谱,提取出它的局部信息,提高时间分辨能力。但短时傅里叶变换是以Fourier变换为基础的,故不能从根本上克服Fourier变换的局限性。
1981年,Morlet首先提出了小波分析(Wavelet Analysis)这一概念,通过一种伸缩和平移小波对信号作变换达到了时频局部化分析的目的,用以测量声音信号的频率定位。但是,小波变换本质上是一种窗口可调的Fourier变换,其小波窗内的信号必须是平稳的,因而也没有摆脱Fourier变换的局限。1998年,Norden E.Huang等提出了Hilbert-Huang变换(Hilbert-Huang Transform,HHT)。它是一种分析非线性、非平稳信号的新方法,该方法不同于Fourier变换分析信号的思想,从根本上克服了Fourier变换的局限性。
但是该方法过程较为复杂,且对待处理信号的特征有一定限制。这些方法在信号提取方面已经有了成熟的应用,为了更加快速而精确的提取信号,本发明介绍了一种基于随机减量技术提出的新的信号提取方法——错位叠加法,通过信号的叠加使待提取信号增强到一定程度,从而提取出所需要的信号。
发明内容
为解决现有技术存在的上述问题,本发明的目的在于提供一种方法简便易行、处理结果准确性高的基于错位叠加算法的声信标信号提取方法。
—获取待提取信号的声信标模拟脉冲信号;
—采用错位叠加算法对采集的声信标模拟信号进行相邻周期叠加和/或相同时间间隔叠加;
—将由相邻周期叠加和相同时间间隔叠加算法计算得出的结果与37.5kHz正弦脉冲信号做相关性分析,计算相关性系数,选取相关性系数较大的处理结果,即更为接近正弦波的结果;提取该结果作为输出错位叠加后的信号,完成声信标信号的提取。
作为优选的实施方式,所述的相邻周期叠加过程如下:
设采集信号的正弦脉冲部分的信号周期为T,将整体信号a向左平移距离T得到信号a1,再与原始信号a叠加得到信号b1;
将信号a1向左平移距离T得到信号a2,再与原始信号a叠加得到信号b2。通过对上述过程的重复叠加,使得正弦部分的能量得到增强;
叠加计算公式为:
其中:S1[n]为原始信号;M为叠加次数,T为叠加步长,k=0,1,2,…。
作为优选的实施方式,所述的相同时间间隔叠加过程如下:
任取一个采样点作为叠加起始点,在离该点1个脉冲距离处往后,选取一段叠加长度l进行叠加,计算公式为:
其中,C[t]为原始信号;M为叠加次数;k=0,1,2,…。
一种基于错位叠加算法的声信标信号提取系统,其特征在于包括:
接收声信标发出的模拟脉冲信号的信号采集单元;
采用错位叠加算法对采集的声信标模拟信号进行相邻周期叠加和/或相同时间间隔叠加处理的计算单元;
将由相邻周期叠加和相同时间间隔叠加算法计算得出的结果与37.5kHz正弦脉冲信号做相关性分析,计算相关性系数,选取相关性系数较大的处理结果,即更为接近正弦波的结果;提取该结果作为输出错位叠加后的信号,完成声信标信号的提取。作为优选的实施方式,所述的相邻周期叠加过程如下:
设采集信号的正弦脉冲部分的信号周期为T,将整体信号a向左平移距离T得到信号a1,再与原始信号a叠加得到信号b1;
将信号a1向左平移距离T得到信号a2,再与原始信号a叠加得到信号b2。通过对上述过程的重复叠加,使得正弦部分的能量得到增强;
叠加计算公式为:
其中:S1[n]为原始信号;M为叠加次数,T为叠加步长,k=0,1,2,…。
作为优选的实施方式,所述的相同时间间隔叠加过程如下:
任取一个采样点作为叠加起始点,在离该点1个脉冲距离处往后,选取一段叠加长度l进行叠加,计算公式为:
其中,C[t]为原始信号;M为叠加次数;k=0,1,2,…。
与现有技术相比,本发明公开的一种基于错位叠加算法的声信标信号提取方法及系统具有以下有益效果:1、简单易行。通过软件编程将算法写入硬件平台中,即可测得在不同位置采集的声信标信号经过处理后的结果。2、结果准确性高。经过多次叠加,合成信号与被关注的正弦信号之间的相关性达到80%以上,提取效果显著。
附图说明
本发明共有附图7张,其中:
图1是三个正弦信号及其合成信号;
图2是信号D1在不同叠加次数下的波形;
图3是按相同时间间隔叠加的计算过程示意图;
图4是声信标信号提取系统框图;
图5为本发明实施例中的原始信号;
图6为本发明实施例中按相邻周期叠加20次的效果示意图;
图7位本发明实施例中按相同时间间隔叠加20次的效果示意图。
具体实施方式
下面结合附图对本发明进行进一步地说明。如图1-4所示,一种基于错位叠加算法的声信标信号提取方法包括以下步骤:
A、搭建声信标信号采集系统,包括声信标、水听器、采集卡等;
B、利用错位叠加算法对采集的信号进行处理;
C、将算法编写程序导入单片机,完成声信标信号提取系统的搭建。
本发明步骤A所述的搭建声信标信号采集系统过程可描述为:
A1、将声信标放入水中,将水听器置于离声信标约1000米处接收声信标发射的37.5kHz的脉冲信号;
A2、水听器接采集卡并连在电脑上,采集卡将水听器接收的信号转换成数字信号并储存在电脑上。
本发明步骤B所述的利用错位叠加算法对采集的信号进行处理过程可描述为:
利用错位叠加算法对信号进行处理,这里提出两种叠加方式:第一种是按照信号自身周期进行叠加,第二种是按照相同时间间隔进行叠加。两种方式的目的都是为了增强所关注的信号,削弱其他成分,从而提取有效信号。
1、相邻周期的叠加
设采集信号的正弦脉冲部分的信号周期为T,将整体信号a向左平移距离T得到信号a1,再与原始信号a叠加得到信号b1;将信号a1向左平移距离T得到信号a2,再与原始信号a叠加得到信号b2。通过对上述过程的重复叠加,可使正弦部分的能量得到增强。
叠加计算公式为:
其中:S1[n]为原始信号;M为叠加次数;T为叠加步长。
2、相同时间间隔的叠加
与第一种方式相似,取一个采样点作为叠加起始点,在离该点1个脉冲距离处往后选取一段叠加长度l进行叠加试验。计算公式为:
其中:C[t]为原始信号;M为叠加次数;k=0,1,2,…。
本发明步骤C所述的将算法编写程序导入单片机过程可描述为:
该算法有两种叠加方式:
比较分别由相邻周期叠加和相同时间间隔叠加算法计算得出的结果;选取更为接近正弦波的结果;提取该结果作为输出错位叠加后的信号,完成声信标信号的提取。
将其编写程序分别导入单片机,然后接入信号采集系统里。整个声信标信号提取系统框图如图4所示。
实施例
在本实施例中采用的航空用声信标的信号形式满足水下声信标最低性能标准,其发射信号主要参数如表1所示。
表1声信标信号参数
将采集的一段声信标信号作为试验信号,如下图5-7所示。按错位叠加算法的两种叠加方式分别对试验信号进行20次错位叠加,可得到如下图所示的信号。右边的信号为局部放大图。两种叠加结果均有接近正弦的趋势,且第一种叠加方式的效果比第二种好。因此最终选择按相邻周期叠加的方式,这种方式处理的效果更好。
将处理后的信号分别与37.5kHz正弦脉冲信号做相关性分析,此分析可借助MATLAB中的corrcoef函数实现。通过计算,第一种方式的相关性系数为0.82,而第二种方式的相关性系数为0.76,可见按相邻周期叠加的效果更好,更接近声信标的原始信号。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (4)
1.一种基于错位叠加算法的声信标信号提取方法,其特征在于包括以下步骤:
—获取待提取信号的声信标模拟脉冲信号;
—采用错位叠加算法对采集的声信标模拟信号进行相邻周期叠加和相同时间间隔叠加;
—将由相邻周期叠加和相同时间间隔叠加算法计算得出的结果与37.5kHz正弦脉冲信号做相关性分析,计算相关性系数,选取相关性系数较大的处理结果,即更为接近正弦波的结果;提取该结果作为输出错位叠加后的信号,完成声信标信号的提取;
所述的相同时间间隔叠加过程如下:
任取一个采样点作为叠加起始点,在离该点1个脉冲距离处往后,选取一段叠加长度l进行叠加,计算公式为:
其中,C[t]为原始信号;M为叠加次数;k=0,1,2,…,M。
2.根据权利要求1所述的基于错位叠加算法的声信标信号提取方法,其特征还在于所述的相邻周期叠加过程如下:
设采集信号的正弦脉冲部分的信号周期为T,将整体信号a向左平移距离T得到信号a1,再与原始信号a叠加得到信号b1;
将信号a1向左平移距离T得到信号a2,再与原始信号a叠加得到信号b2;通过对上述过程的重复叠加,使得正弦部分的能量得到增强;
叠加计算公式为:
其中:S1[n]为原始信号;M为叠加次数,T为叠加步长,k=0,1,2,…,M。
3.一种基于错位叠加算法的声信标信号提取系统,其特征在于包括:
接收声信标发出的模拟脉冲信号的信号采集单元;
采用错位叠加算法对采集的声信标模拟信号进行相邻周期叠加和相同时间间隔叠加处理的计算单元;
将由相邻周期叠加和相同时间间隔叠加算法计算得出的结果与正弦脉冲信号做相关性分析,计算相关性系数,选取相关性系数较大的处理结果,即更为接近正弦波的结果;提取该结果作为输出错位叠加后的信号,完成声信标信号的提取;
所述的相同时间间隔叠加过程如下:
任取一个采样点作为叠加起始点,在离该点1个脉冲距离处往后,选取一段叠加长度l进行叠加,计算公式为:
其中,C[t]为原始信号;M为叠加次数;k=0,1,2,…,M。
4.根据权利要求3所述的基于错位叠加算法的声信标信号提取系统,其特征还在于:特征还在于所述的相邻周期叠加过程如下:
设采集信号的正弦脉冲部分的信号周期为T,将整体信号a向左平移距离T得到信号a1,再与原始信号a叠加得到信号b1;
将信号a1向左平移距离T得到信号a2,再与原始信号a叠加得到信号b2;通过对上述过程的重复叠加,使得正弦部分的能量得到增强;
叠加计算公式为:
其中:S1[n]为原始信号;M为叠加次数,T为叠加步长,k=0,1,2,…,M。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611104475.0A CN106595835B (zh) | 2016-12-05 | 2016-12-05 | 一种基于错位叠加算法的声信标信号提取方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611104475.0A CN106595835B (zh) | 2016-12-05 | 2016-12-05 | 一种基于错位叠加算法的声信标信号提取方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106595835A CN106595835A (zh) | 2017-04-26 |
CN106595835B true CN106595835B (zh) | 2019-03-12 |
Family
ID=58595294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611104475.0A Expired - Fee Related CN106595835B (zh) | 2016-12-05 | 2016-12-05 | 一种基于错位叠加算法的声信标信号提取方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106595835B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110793623A (zh) * | 2019-11-15 | 2020-02-14 | 大连海事大学 | 一种基于自适应性错位叠加算法的发动机故障成分提取方法 |
CN112435692B (zh) * | 2020-11-12 | 2023-09-05 | 南京世海声学科技有限公司 | 一种基于多脉冲相干累加的信标信号高增益检测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013088275A1 (en) * | 2011-12-16 | 2013-06-20 | Eads Singapore Pte. Ltd. | A new design of underwater locator beacon with integrated pressure sensor |
CN204086566U (zh) * | 2014-11-06 | 2015-01-07 | 安徽理工大学 | 一种基于阻抗分析的飞机黑匣子探测实验装置 |
CN104535169A (zh) * | 2014-12-03 | 2015-04-22 | 北京神州普惠科技股份有限公司 | 一种基于光纤水听器阵列的噪声测量装置及测量方法 |
CN204731415U (zh) * | 2015-07-07 | 2015-10-28 | 大连测控技术研究所 | 一种多频率水声信标远距离快速探测系统 |
CN105656564A (zh) * | 2016-01-17 | 2016-06-08 | 中国船舶重工集团公司第七六○研究所 | 一种不依赖水介质的水声信标信号提取装置 |
CN105675121A (zh) * | 2016-01-17 | 2016-06-15 | 中国船舶重工集团公司第七六○研究所 | 一种便携式多功能水声信标检测装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160192185A1 (en) * | 2014-09-03 | 2016-06-30 | Peeramyd, Llc | Beacon Activated Software |
-
2016
- 2016-12-05 CN CN201611104475.0A patent/CN106595835B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013088275A1 (en) * | 2011-12-16 | 2013-06-20 | Eads Singapore Pte. Ltd. | A new design of underwater locator beacon with integrated pressure sensor |
CN204086566U (zh) * | 2014-11-06 | 2015-01-07 | 安徽理工大学 | 一种基于阻抗分析的飞机黑匣子探测实验装置 |
CN104535169A (zh) * | 2014-12-03 | 2015-04-22 | 北京神州普惠科技股份有限公司 | 一种基于光纤水听器阵列的噪声测量装置及测量方法 |
CN204731415U (zh) * | 2015-07-07 | 2015-10-28 | 大连测控技术研究所 | 一种多频率水声信标远距离快速探测系统 |
CN105656564A (zh) * | 2016-01-17 | 2016-06-08 | 中国船舶重工集团公司第七六○研究所 | 一种不依赖水介质的水声信标信号提取装置 |
CN105675121A (zh) * | 2016-01-17 | 2016-06-15 | 中国船舶重工集团公司第七六○研究所 | 一种便携式多功能水声信标检测装置 |
Non-Patent Citations (2)
Title |
---|
《基于异响信号的发动机故障成分提取》;宁大勇等;《光学 精密工程》;20151030;第23卷(第10期);第314页第2-4节及图1 |
《随机减量技术的方法与理论》;陈德成等;《振动与冲击》;19840430(第12期);第31-40页 |
Also Published As
Publication number | Publication date |
---|---|
CN106595835A (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109815942B (zh) | 基于海洋环境噪声信号的简正波特征提取方法 | |
CN110852201B (zh) | 一种基于多脉冲包络谱匹配的脉冲信号检测方法 | |
CN107461611B (zh) | 基于小波和emd重构相结合的泄漏检测方法及泄漏检测装置 | |
CN106503336B (zh) | 一种海豚嘀嗒声信号建模与合成的方法 | |
CN103076594B (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN107340055B (zh) | 一种基于多测度融合的随机共振微弱信号检测方法 | |
CN101893698B (zh) | 噪声源测试分析方法及其装置 | |
CN106595835B (zh) | 一种基于错位叠加算法的声信标信号提取方法及系统 | |
CN104168232A (zh) | 一种水声信道中多径时延与多普勒频移的测定方法 | |
CN107886078A (zh) | 一种基于分层自适应阈值函数的小波阈值降噪方法 | |
CN110376575A (zh) | 一种基于阻尼参数匹配随机共振的低频线谱检测方法 | |
CN113049080A (zh) | 一种舰船辐射噪声的gdwc听觉特征提取方法 | |
CN114167423A (zh) | 基于深度回归网络的雷达海浪参数测量方法 | |
CN105158734A (zh) | 一种基于阵不变量的单矢量水听器被动定位方法 | |
CN103557925B (zh) | 水下目标gammatone离散小波系数听觉特征提取方法 | |
CN103559893B (zh) | 一种水下目标gammachirp倒谱系数听觉特征提取方法 | |
CN104280776A (zh) | 一种自适应小波阈值求取方法 | |
CN107332642B (zh) | 基于海洋哺乳类动物信号合成的低截获概率波形设计方法 | |
Zeng et al. | Underwater sound classification based on Gammatone filter bank and Hilbert-Huang transform | |
CN107635181A (zh) | 一种基于信道学习的多址感知源反馈优化方法 | |
CN104407341A (zh) | 利用鲸声及仿其发声特点的伪装隐蔽水下声源装置及方法 | |
Liu et al. | Speech enhancement based on Hilbert-Huang transform | |
CN107907279A (zh) | 基于小波系数幅值的多相流管道泄漏声波信号分析方法 | |
Li et al. | Feature extraction method for ship-radiated noise based on extreme-point symmetric mode decomposition and dispersion entropy | |
Veeraiyan et al. | Frequency domain based approach for denoising of underwater acoustic signal using EMD |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190312 Termination date: 20191205 |