一种水声信号谱峰检测方法技术

本发明专利技术公开了一种水声信号谱峰检测方法，对包含谱峰的信号波形消除极性得到B(n)，求增益系数c(n)，c(n)与B(n)相乘得到输出信号C(n)，对C(n)中值滤波平滑得到D(n)，利用高斯窗对D(n)拟合得到E(n)，对E(n)提取极大值点以及对应的极大值。本发明专利技术通过绝对值去极性、双滑动窗谱峰增强、滑动平均滤波、高斯拟合等步骤能对信号波形中存在的谱峰进行检测，有效减少水声信道对原始谱峰的影响，实现谱峰特征的增强。强。强。

【技术实现步骤摘要】
一种水声信号谱峰检测方法


[0001]本专利技术属于水声信号的特征提取
，涉及一种水声信号谱峰检测方法。

技术介绍

[0002]在进行水声信号识别或参数估计时经常需要对谱峰进行检测，比如对水声频移键控(FSK)识别时会对频谱的谱峰进行检测，对水声正交频分复用(OFDM)信号识别时会对其自相关谱峰进行检测，对FSK、相移键控(PSK)以及直接序列扩频(DSSS)等通信信号进行载频估计时需要对频域的谱峰位置进行检测。
[0003]在水声中信号的谱峰检测方法大多都来源于无线电中。无线电中有许多经典的谱峰提取方法，比如直接谱峰检测算法、三点谱峰检测算法、高斯拟合算法、多项式拟合算法等。直接谱峰检测虽然计算复杂度低，但是抗噪性能差。三点谱峰检测、高斯拟合、多项式拟合算法对谱峰的形状要求较高。水声信道相比于常见的无线信道更为复杂，具有复杂的时
‑
空
‑
频变特性，在频域表现为频率选择性衰落，并且当信噪比较低或信道时变明显时，利用相关得到的谱峰也会变得不稳定。所以在水声信道下利用信号频域或相关波形进行谱峰检测具有较高的挑战性。

技术实现思路

[0004]针对上述现有技术，本专利技术要解决的技术问题是提供一种针对水声信道下谱峰特征不明显的水声信号谱峰检测方法，能有效地实现对谱峰特征的增强，并且还能对干扰进行抑制。
[0005]为解决上述技术问题，本专利技术的一种水声信号谱峰检测方法，包括以下步骤：
[0006]S1、用A(n)表示包含谱峰的信号波形，n为信号波形第n个采样点的位置，n＝1,2,
···
,N，N为采样点个数，当谱峰波形存在负极性时，令B(n)＝|A(n)|，当谱峰波形不存在负极性时，令B(n)＝A(n)；
[0007]S2、以B(n)的第n个采样点为中心，建立两个滑动窗，短滑动窗长度为2s
win
+1，长滑动窗长度为2l
win
+1，并且l
win
＝2s
win
；通过所述两个滑动窗的能量差得到增益系数c(n)，；通过所述两个滑动窗的能量差得到增益系数c(n)，是第n个采样点的短滑动窗能量，是第n个采样点的长滑动窗能量，其中输出信号为C(n)＝B(n)
·
c(n)，其中n＝1,2,
···
,N，当c(n)值为负数、n
‑
s
win
≤0、n
‑
l
win
≤0、n+s
win
＞N或n+l
win
＞N时，令c(n)＝0；
[0008]S3、对S2中得到的C(n)波形进行滑动平均滤波，经过滑动平均滤波后的输出信号为：其中2P+1是滑动平均滤波长度，n＝1,2,
···
,N，当n
‑
P≤0或n+P＞N时，令D(n)＝C(n)；
[0009]S4、对S3中得到的D(n)做高斯拟合，高斯窗函数为：其中L是高斯窗函数长度，
‑
(L
‑
1)/2≤n≤(L
‑
1)/2，α是高斯窗函数的宽度系数，α与高斯函数的标准差σ＝(L
‑
1)/2α成反比；经过高斯拟合后的输出信号为：其中L＝2Q+1是高斯窗函数的长度，n＝1,2,
···
,N，然后对E(n)做归一化得到优化后的谱峰波形，当n
‑
Q≤0或n+Q＞N时，令E(n)＝D(n)；
[0010]S5、提取E(n)中所有极大值以及对应的极大值点，并按照极大值从大到小排序，如果是单个谱峰，则提取E(n)中最大的极大值以及对应的极大值点作为谱峰，如果有多个谱峰，则提取E(n)中前m个极大值以及对应的极大值点作为谱峰，m等于谱峰数量。
[0011]进一步的，S3中滑动平均滤波长度等于S2中短滑动窗长度。
[0012]进一步的，S4中高斯窗函数长度等于S2中短滑动窗长度。
[0013]进一步的，S2中短滑动窗长度、S3中滑动平均滤波长度和S4高斯窗函数长度相等。
[0014]本专利技术的有益效果：本专利技术通过绝对值去极性、双滑动窗谱峰增强、滑动平均滤波、高斯拟合等步骤能对信号波形中存在的谱峰进行检测，有效减少水声信道对原始谱峰的影响，实现谱峰特征的增强。本专利技术主要有两个方面的应用：
[0015]第一个方面，我们可以利用本专利技术进行非合作水声通信信号进行谱峰特征的提取，通过这些特征实现非合作水声通信信号的识别。比如水声OFDM一般存在循环前缀的结构，我们对信号的自相关波形进行谱峰检测，确定是否存在两个相关峰的情况。又比如水声FSK信号在频域有多个谱峰，我们对信号频域波形进行谱峰检测，确定存在多少个谱峰，根据检测到的谱峰数量确定FSK的调制阶数。
[0016]第二个方面，我们可以利用本专利技术进行水声FSK信号各载波频率的估计。利用本专利技术进行谱峰检测最终会得到极大值点与对应的极大值，在频域上极大值点对应的就是信号的频率。
附图说明
[0017]图1为本专利技术所提供的水声信号谱峰提取方法的流程图；
[0018]图2为本专利技术所提供的进行测试的信道冲激响应图；
[0019]图3为仿真的8FSK进行傅里叶变换得到的原始频域波形；
[0020]图4为在图3原始频域波形的基础上得到增益系数；
[0021]图5为图4增益系数与图3原始频域波形乘积的结果；
[0022]图6为对图5中波形进行滑动平均滤波的结果；
[0023]图7为对图6中波形进行高斯拟合的结果；
[0024]图8为将图7中横坐标限定在10kHz至14kHz内的结果。
具体实施方式
[0025]下面结合说明书附图和具体实施例对本专利技术做进一步说明。
[0026]结合图1，本专利技术包括以下步骤：
[0027]S1、用A(n)表示包含谱峰的信号波形，n为信号波形第n个采样点的位置。为了消除
双极性信号的不确定性，我们首先对A(n)取绝对值得到B(n)，即B(n)＝|A(n)|，当谱峰波形不存在负极性时，令B(n)＝A(n)；
[0028]S2、针对全部采样点逐个进行以下操作：以B(n)的采样点n为中心，建立一长一短两个滑动窗，短滑动窗的长度为2s
win
+1，长滑动窗的长度为2l
win
+1，并且l
win
＝2s
win
。通过这两个滑动窗的能量差得到增益系数c(n)，c(n)可以表示为其中当c(n)值为负数，n
‑
s
win
≤0，n
‑
l
win
≤0，n+s
win
＞N或n+l
win
＞N时，令c(n)＝0，输出信号为C(n)＝B(n)
·
c(n)。通过乘以增益系数c(n)，当信号波形本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种水声信号谱峰检测方法，其特征在于，包括以下步骤：S1、用A(n)表示包含谱峰的信号波形，n为信号波形第n个采样点的位置，n＝1,2,
···
,N，N为采样点个数，当谱峰波形存在负极性时，令B(n)＝|A(n)|，当谱峰波形不存在负极性时，令B(n)＝A(n)；S2、以B(n)的第n个采样点为中心，建立两个滑动窗，短滑动窗长度为2s
win
+1，长滑动窗长度为2l
win
+1，并且l
win
＝2s
win
；通过所述两个滑动窗的能量差得到增益系数c(n)，；通过所述两个滑动窗的能量差得到增益系数c(n)，是第n个采样点的短滑动窗能量，是第n个采样点的长滑动窗能量，其中输出信号为C(n)＝B(n)
·
c(n)，其中n＝1,2,
···
,N，当c(n)值为负数、n
‑
s
win
≤0、n
‑
l
win
≤0、n+s
win
＞N或n+l
win
＞N时，令c(n)＝0；S3、对S2中得到的C(n)波形进行滑动平均滤波，经过滑动平均滤波后的输出信号为：其中2P+1是滑动平均滤波长度，n＝1,2,
···<...

【专利技术属性】
技术研发人员：刘凇佐，方涛，王虔，青昕，乔钢，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：