一种避碰声呐时延估计及防跨周期的方法

(整期优先)网络出版时间:2021-09-16
/ 2

一种避碰声呐时延估计及防跨周期的方法

黄 凡

(海装广州局驻昆明地区第三军事代表室 , 云南 昆明 )


摘 要: 分别利用匹配滤波算法与短时傅里叶算法对测距脉冲及其回波进行时延检算,对比二者在窄脉冲(低功耗)条件下的仿真结果,针对测距脉冲“跨周期”现象,采取不同中心频率下正调频脉冲—负调频脉冲交替发射的方式,仿真结果验证了该探测方式的有效性。

关键词: 避碰声呐; 时延估计; 匹配滤波

1 前言

避碰声呐一般安装于水下自主航行器的头部,为其提供前方有无障碍物及障碍物距离信息。避碰声呐通过主动发射脉冲、检测回波信号时延对障碍目标完成距离估计,在产品低功耗要求下,探测脉冲的脉宽受到限制,短时傅里叶算法的时延估计性能随脉宽的减小迅速下降。本文基于匹配滤波算法,通过包络提取和正负调频信号交替发射,提高了避碰声呐在窄脉宽条件下的测距精度,解决了测距脉冲的跨周期问题。

2 两种时延估计算法

2.1 短时傅里叶

短时傅里叶算法的基本思想是用窗函数来截取信号,它是一种线性变换,满足叠加原理,概念明确,是目前研究非平稳信号应用最为广泛的一种方法[1]。设信号为6142e5223b55b_html_1d16150c677c5351.gif ,分析窗为6142e5223b55b_html_3ffebbb28b2f77c9.gif ,可定义6142e5223b55b_html_1d16150c677c5351.gif 的短时傅里叶变换为:

6142e5223b55b_html_f396ef61ff888c5.gif (1)

时间6142e5223b55b_html_ea65dd98fd3aac4e.gif 处的短时傅里叶变换的计算过程如下:

将窗函数6142e5223b55b_html_3ffebbb28b2f77c9.gif 由时间零点处平移至时间6142e5223b55b_html_ea65dd98fd3aac4e.gif 处,得到6142e5223b55b_html_1a386b025f860c25.gif

信号被平移后的分析窗加窗隔断后,得到短时信号:

6142e5223b55b_html_ea3c751d2469c0c2.gif (2)

利用傅里叶变换得到短时信号6142e5223b55b_html_1f3ebe65b12203c6.gif 的频谱。

设分析窗的持续时间为6142e5223b55b_html_94bf5042c8c56c4c.gif ,则加窗截断相当于取6142e5223b55b_html_b8140527989f5f00.gif 范围内的非平稳信号进行傅里叶变换。若期望时间分辨率高,则分析窗的持续时间必须尽量短。

2.2 匹配滤波

匹配滤波器是白噪声背景下有规信号检测的最佳接收机[2]。设接收机的输入是确知的有限能量信号6142e5223b55b_html_19315f8fb1abc63d.gif (其复数形式为6142e5223b55b_html_84b8ae54ecabd278.gif ),接收系统响应函数是6142e5223b55b_html_d14d66ef6c344654.gif ,传递函数是6142e5223b55b_html_75994c3130f70bf8.gif ,在6142e5223b55b_html_34e96eacc65c58ca.gif 时的输出是:

6142e5223b55b_html_13a0cb35df4162c9.gif (3)

6142e5223b55b_html_974bc63c851eb0fc.gif (4)

使6142e5223b55b_html_855056019e641f9c.gif 达到最大值的系统定义为与信号6142e5223b55b_html_84b8ae54ecabd278.gif 匹配的匹配滤波器系统简称匹配滤波器。

由Schwartz不等式

6142e5223b55b_html_d5d4c1380b9a882.gif (5)

可以确定式(4)最大的条件为:

6142e5223b55b_html_cdb301c9e84c8da.gif (6)

6142e5223b55b_html_46d1621d387a9918.gif 是任一常数,通常取1,相应的响应函数是

6142e5223b55b_html_bfc18009802b93ad.gif (7)

匹配滤波器的基本特性如下[3]

(1)响应函数是被匹配信号的延迟共轭镜像波形;

(2)匹配滤波器相当于一个能计算互相关函数的互相关器;

(3)匹配滤波器对延迟了的匹配信号也是匹配的,只是最大值出现时间向延迟方向平移了一个同样的量。

3 两种时延估计方法仿真对比

相对于短时傅里叶测距算法,匹配滤波算法的优势主要体现在窄脉宽(低功耗)条件下测距结果的抗干扰能力,下面分别在不同探测脉宽、不同信噪比下对二者的测距结果进行仿真分析,归一化仿真条件见表1。

表1 时延估计方法仿真参数

参数

短时傅里叶测距算法

匹配滤波测距算法

信号形式

200kHz单频信号

中心频率200kHz、带宽20kHz chirp信号

采样频率(kHz)

500

500

目标距离(m)

75

75

脉宽(ms)

5

2

5

2

信噪比(dB)

0

-10

-16

0

-10

-16

6142e5223b55b_html_e9f48950dabb5bed.gif

图1 脉宽=5ms,SNR=0dB时两种算法测距结果对比


6142e5223b55b_html_64a3bc051a6cf398.gif

图2 脉宽=5ms,SNR=-10dB时两种算法测距结果对比

6142e5223b55b_html_2117ecdbc05fb8fd.gif

图3 脉宽=5ms,SNR=-16dB时两种算法测距结果对比

当探测脉冲为5ms、信噪比为0dB时,两种算法的测距结果均较为理想,在目标距离75m处形成较为明显的峰值,峰值相对背景有约20dB的增益空间,门限值设置易于把握;对比图1、图2及图3可以看出,随着信噪比降低,两种算法的测距性能呈现出不同程度的下降,但对于距离检算效果影响不大。


6142e5223b55b_html_1263f8bc25f3eb9.gif

图4 脉宽=2ms,SNR=0dB时两种算法测距结果对比

6142e5223b55b_html_f3c43d114bba11eb.gif

图5 脉宽=2ms,SNR=-10dB时两种算法测距结果对比


6142e5223b55b_html_cc0d778fb96a9e58.gif

图6 脉宽=2ms,SNR=-16dB时两种算法测距结果对比

图4、图5、图6为探测脉宽2ms时各信噪比下的测距结果仿真。对比图4和图1可知,脉宽变窄对两种算法的抗干扰能力均有一定影响,但在0dB信噪比下基于以上两种算法仍可得到较为准确的目标距离;对比图2和图5,当信噪比下降至-10dB时,短时傅里叶算法对2ms探测脉宽的计算结果出现伪峰,对于门限设定要求十分苛刻,而匹配滤波算法在不同探测脉宽下的测距结果相差不大,其抗干扰能力明显优于短时傅里叶算法;当信噪比下降至-16dB时,短时傅里叶算法在2ms探测脉宽下已无法实现对于目标距离的准确估计,峰值出现在190m左右,而匹配滤波算法依然有6dB左右的门限冗余空间,只要探测门限设置恰当依然可以准确得出障碍目标距离。

4 “跨周期”现象抑制效果仿真

“跨周期”是主动声呐在探测目标过程中较为常见的误判现象[4,5,6]。当探测距离大于脉冲发射周期与介质声速乘积的一半时,回波脉冲与相对应的发射脉冲不在一个探测周期内,由此产生的距离计算失准即为“跨周期”现象。避碰声呐正常探测过程中的发射脉冲和回波脉冲序列如图7所示。

6142e5223b55b_html_e2a0dcd766d13fbc.gif 图7 避碰声呐正常探测时发射脉冲及回波脉冲示意图

图中6142e5223b55b_html_31f2476c5a8a5c57.gif 为探测周期,6142e5223b55b_html_ea65dd98fd3aac4e.gif 为回波延时,因此可得目标距离6142e5223b55b_html_adcb83f369b4dae0.gif 为:

6142e5223b55b_html_92944c5c7d353405.gif (8)

式(7)中6142e5223b55b_html_499e5adf91c5a607.gif 为水中声速。当避碰声呐出现“跨周期”现象时,其发射脉冲和回波脉冲序列如图8所示。

6142e5223b55b_html_d36c222c6edde053.gif 图8 避碰声呐“跨周期”探测时发射脉冲及回波脉冲示意图

当探测距离大于脉冲发射周期与介质声速乘积的一半时,回波脉冲将会落入下一个探测周期内,如图8中的1号脉冲回波在2号脉冲的探测周期内,2号脉冲在3号脉冲的探测周期内并以此类推。上图中的真实距离6142e5223b55b_html_adcb83f369b4dae0.gif 与式(8)定义相同,而探测距离6142e5223b55b_html_c63b76c446a791c9.gif 为:

6142e5223b55b_html_e7711614a69fc81e.gif (9)

由此导致探测距离比实际距离要大,探测误差为一个探测周期的半声程6142e5223b55b_html_74d4f3f0f722cc92.gif 。以匹配滤波算法为例,若相邻两脉冲的频率、相位等信号形式相同,那么必然会出现回波脉冲“跨周期”现象;相反,邻近探测脉冲若采用不同中心频率探测,且将调频方式分为正、负调频,“跨周期”现象将得到有效改善。

6142e5223b55b_html_2d81d071a6084a20.gif

图9 避碰声呐探测脉冲信号形式及回波分布

图9为模拟避碰声呐的探测脉冲,其中1-3号脉冲是中心频率分别为210kHz、170kHz、190kHz的正调频信号,4-6号脉冲是中心频率为210kHz、170kHz、190kHz的负调频信号,分别以1号、4号探测脉冲的信号形式作为基准信号对各回波信号进行匹配滤波,所得结果见图10.

6142e5223b55b_html_83ee7252db613a6.gif

图10 中心频率210kHz的正负调频信号与各脉冲的匹配滤波结果

从图10中可以看出,匹配滤波算法对不同中心频率回波信号的抑制效果较好,相应脉冲与相邻干扰脉冲的匹配结果相差在10倍以上,中心频率相同、调频形式不同的脉冲信号虽然匹配结果相对大一些,然而与相应脉冲的匹配结果相比仍有5倍左右的差距,因此只要判别门限设置得当,“跨周期”现象即可通过该方法得到有效解决。

5 结束语

通过仿真结果可知,相比于短时傅里叶算法,匹配滤波算法的时延估计精度对回波信噪比要求不高,尤其是在低功耗-窄脉冲的探测要求下,匹配滤波算法优势体现的更加明显。针对避碰声呐探测过程中的“跨周期”现象,采取不同中心频率信号结合正负调频的手段,将探测信号分为6种并交替发射,通过匹配滤波算法仿真验证了该探测方式的有效性。



参考文献

[1] 石海泉.波导环境下矢量声场特性分析及特征提取技术[D],哈尔滨:哈尔滨工程大学出版社,2013.

[2] 鞠德航,林可祥.信号检测理论导论[M],北京:科学出版社,1975.

[3] 朱埜.主动声呐检测信息原理[M],北京:科学出版社,2014.

[4] 侯自强,李贵斌.声呐信号处理-原理和设备[M],北京:海洋出版社,1986.

[5] 李启虎.进入21世纪的声呐技术[D],信号处理,2012.

[6] 孙超,刘雄厚.MIMO声纳:概念与技术特点探讨[D],声学技术,2012.


作者简介:黄凡(1987-),男,工学博士,工程师,主要从事目标特性及辨识技术研究工作。