舰载预警机对海侦察预警阵位配置研究
2024-10-10董苇青张毅戴梓琴李雪菲
摘 要:为提高舰载预警机对敌水面舰艇的搜索能力,保证舰载预警机能够安全持续地跟踪监视预警,提出了舰载预警机对海面目标侦察预警阵位配置方法。基于阵位配置的基本要求和作战背景的符合条件,分析了影响舰载预警机提供航母编队预警距离的因素,建立了舰载预警机巡逻线配置模型,充分考虑了舰载预警机在巡逻线端点的情况,并计算了最小无源探测距离。根据来袭水面舰艇的位置和运动状态及其作战能力对舰载预警机的阵位配置进行定量研究,利用Matlab仿真实验,验证了模型的有效性,为舰载预警机阵位配置提出建议,满足作战使用的需要。
关键词:舰载预警机;海面侦察;战术使用;阵位配置
引用格式:
董苇青,张毅,戴梓琴,等.
舰载预警机对海侦察预警阵位配置研究
.指挥控制与仿真,2024,46(4):77-81.
DONG W Q,ZHANG Y,DAI Z Q,et al.
Research on sea reconnaissance and early warning position configuration of carrier-based early warning aircraft
.Command Control & Simulation,2024,46(4):77-81.
中图分类号:E926.37 文献标志码:A DOI:10.3969/j.issn.1673-3819.2024.04.010
Research on sea reconnaissance and early warning position
configuration of carrier-based early warning aircraft
DONG Weiqing, ZHANG Yi, DAI Ziqin, LI Xuefei
(Naval Aviation University, Yantai 264001, China)
Abstract:In order to improve the search ability of carrier-based early warning aircraft to enemy surface ships and ensure the safe and continuous tracking, monitoring and early warning of carrier-based early warning aircraft, a method for reconnaissance and early warning position configuration of carrier-based early warning aircraft for sea surface targets is proposed. Based on the basic requirements of its position configuration and the conditions that meet the combat background, the factors affecting the early warning distance of the carrier formation provided by the carrier-based early warning aircraft are analyzed, and the patrol line configuration model of the carrier-based early warning aircraft is established. The situation of the carrier-based early warning aircraft at the end of the patrol line is fully considered, and the minimum passive detection distance is calculated. According to the position and motion state of the enemy attacking the surface ship and its combat capability, the quantitative research on the position configuration of our carrier-based early warning aircraft is carried out. The effectiveness of the model is verified by Matlab simulation experiments, and suggestions are made for the position configuration of the carrier-based early warning aircraft to meet the needs of combat use.
Key words:carrier-based early warning aircraft; sea surface reconnaissance; tactical operation; position configuration
收稿日期: 2023-06-14修回日期: 2023-07-17
作者简介:
董苇青(1997—),男,本科,研究方向为预警侦察。
1 小波阈值算法
1.1 小波阈值去噪
假设一维噪声信号模型可以表示如下:
x(t)=h(t)+ε(t)(1)
其中,x(t)是混合信号,h(t)是原始信号,ε(t)是噪声。在最简单的情况下,ε(t)可以假设为高斯白噪声。实际工程中,h(t)通常被认为是低频信号或相对稳定的信号,而噪声信号通常是高频信号。小波变换的目的就是抑制高频的无用信号并恢复低频的有用部分。小波阈值的去噪过程包括小波变换分解、阈值处理和信号重构。其原理如图1所示。
对于确定的分解水平,可以使用离散小波变换将信号分解为近似系数序列和细节系数序列。近似系数代表低频分解,通常相当于有用信号,细节系数代表中高频分解,通常被视作噪声。图2显示了一维信号的典型3级小波分解示意图。其中,Y(t)为原始信号,CA为低频信息,CD为高频小波系数。
经过小波分解处理后,有用信号的小波系数幅度远大于噪声的幅度。当小波分解水平增加时,有用信号的小波系数的幅度很少改变,但噪声的小波系数幅度迅速衰减。小波分解的这个特性保证了其可以构造阈值和阈值函数来处理分解系数,从而达到消除噪声的目的。最后,通过应用逆小波变换重建去噪信号。
小波阈值去噪的过程包括小波分解、阈值收缩过程和小波重构。小波分解涉及小波基和分解水平的选择,阈值收缩处理需要考虑阈值计算和阈值函数。
1.2 改进阈值函数
小波阈值是将小波系数与给定阈值进行比较。如果小波系数小于阈值,则小波系数置零。如果小波系数大于该阈值,则应保留或修改小波系数。常用的阈值函数有硬阈值函数和软阈值函数。
软阈值函数为
J,K=Sgn(WJ,K)(WJ,K-λ), |WJ,K|≥λ0,|WJ,K|<λ(2)
硬阈值函数为
J,K=WJ,K, |WJ,K|≥λ0,|WJ,K|<λ(3)
在式(2)和(3)中,WJ,K是处理前的小波系数;J,K是处理后的小波系数;λ是设置的阈值;Sgn(·)是符号函数;J是最大分解尺度,j=1,2,…,J;K是小波系数的长度,k=1,2,…,nj;nj=N/2J-j+1;N是信号的长度。
软阈值函数通常会对小波系数进行平滑处理来减小噪声,然而,平滑性过强会导致信号的细节信息丢失,从而导致过度去噪。并且,实际的小波系数和处理后的小波系数始终存在恒定的偏差,这会降低重建信号的质量。硬阈值函数会将小于阈值的小波系数置零,从而导致信号中的小幅度信号被完全消除,造成有效信息的丢失。此外,由于硬阈值函数直接将小波系数截断,在信号重建过程中会出现振荡现象,进一步产生伪吉布斯现象。
为了克服软阈值函数和硬阈值函数的上述缺点,本文提出了一种改进的阈值函数,这种改进的阈值函数可以通过改变参数因子,进行灵活调整,使函数值位于硬阈值和软阈值函数值之间。改进的阈值函数如下:
J,K=WJ,K-Sgn(WJ,K)×βλeα(|WJ,K|-λ), |WJ,K|≥λ0,|WJ,K|<λ(4)
式中,α和β是调节因子,用于调节J,K的变化速度,其中,α用于调节函数向渐近线的逼近程度,取值范围为[0,∞],β用于调节阈值λ阶段点的y值,取值范围为[0,1]。图3展示了当β为1时,函数曲线随α的变化情况,随着α的增大,改进的阈值函数更加逼近软阈值函数。图4展示了当α为1时,函数曲线随β的变化情况,随着β的增大,改进的阈值函数在截止点的阈值更接近硬阈值函数。
本文提出的改进阈值函数通过调节α和β两个调节因子,就可以更便捷地控制阈值函数,在消除阈值函数不连续性的同时,使函数更迅速靠近期望的位置。改进方法保证了阈值函数的连续性,不会导致突变的产生,从而避免了伪吉布斯现象的发生,使得降噪后信号能够保持原始信号的平滑性。同时,改进方法可以减小由于小波分解重构对信号压缩带来的偏差,从而提高重构信号与真实信号的逼近程度。
1.3 基于改进模拟退火算法优化的阈值计算
阈值常用的选择方法主要有固定阈值、极值阈值、无偏似然估计阈值、极大极小阈值和启发式阈值。但上述方法需要根据信号和噪声的性质及应用需求进行调整,对于小波系数而言,全局最优的阈值并不一定局部最优。为了克服这一缺点,本文采用了一种基于改进模拟退火算法的阈值计算方法来优化阈值计算。
模拟退火算法是一种启发式优化算法,用于解决各种组合优化和搜索问题。但是传统的模拟退火算法会存在丢失最优解、容易陷入局部最优解等问题。
为了解决上述传统模拟退火算法存在的缺点,进一步提高算法的搜索能力,笔者结合目标小波阈值的特点,对算法做出以下改进:
1) 记录全局最优解,解决原始算法由于随机扰动造成的最优解丢失的问题。
2) 采用正态分布作为扰动的邻域分布,选择当前解和最优解的差值作为方差。在远离最优解时,大方差会带来大扰动,当靠近最优解时,小方差避免跳出最优解。
3) 采用方差阈值的方式,连续多次方差小于阈值则增大方差跳出当前状态,防止陷入局部最优解。增大方差时回火升温,避免因为温度过低导致搜索能力变差。
改进算法的流程图如图5所示。
该研究通过改进模拟退火算法,提高了算法的寻优能力,减少搜索陷入局部最优解的情况;通过将阈值的计算转化为基于信噪比的优化,解决了阈值选择问题;通过对优化算法的改进,增强了其搜索能力。
2 仿真实验及分析
2.1 噪声分析
本次实验将真实声波数据添加噪声,针对改进算法进行仿真验证,对实验结果进行分析,为真实信号的处理提供了参考。为了更好地近似光纤声波传感的真实信号,在仿真实验中以微软可定制语音数据库MS-SNSD中的无噪声波信号作为模拟原始信号。向模拟信号添加噪声,噪声选择数据集NoiseX-92中的白噪声以及分布式光纤传感声波信号稳态值模拟的真实噪声。
白噪声常常用来模拟自然环境中的背景噪声,它包含了所有频率的噪声成分。真实噪声包含系统使用过程中存在的噪声。添加不同类型的噪声有助于验证去噪算法的鲁棒性等性能。不同类型的噪声反映了实际环境中存在的干扰,使算法能够应对各种复杂的去噪场景。通过在仿真实验中添加这些噪声,可以评估去噪算法在不同噪声条件下的效果。
2.2 评价指标
实验中,使用信噪比(SNR)、均方误差(MSE)、波形互相关系数(CC)定量地比较降噪算法的性能。SNR可以反映降噪算法的去噪效果,MSE和CC可以反映去噪后信号和真实信号之间的接近程度。信噪比、均方误差和波形相似度的计算公式如下:
SNRdB=10log10PsignalPnoise=10log10∑ni=1Y2i∑ni=1(Yi-i)2(5)
MSE=1n∑ni=1(Yi-i)2(6)
CC=Rxy(m)=1N∑Nk=1x(k)y(k-m)(7)
其中,Yi是原始信号幅值,i是去噪信号的幅值,x(k)为原始信号序列,y(k-m)为去噪信号序列,m为两个序列的移位长度,N为信号序列的长度。
2.3 实验参数设置
1)小波基选择
不同的小波基具有不同的时频特性。为选择最适合声波信号去噪的小波基,本文对仿真信号选取了不同的小波基函数,计算了降噪后的信噪比,基于上述降噪评价规则对比了最优小波基函数。通过添加高斯白噪声(SNR为5 dB)对信号进行仿真,以测试小波基haar、Sym、db、coif和bior函数对声波信号降噪的影响。各小波基去噪后的信噪比如表1所示。表1中展示了各小波基在相同条件下对仿真信号的去噪结果,各小波系仅记录去噪效果最好的小波。其中,sym小波系中的sym8小波对仿真信号去噪效果最好,因此后续的仿真实验采用sym8小波作为小波阈值算法的小波基。
2)分解层数选择
为选择最适合声波信号去噪的分解层数,本文对仿真信号选取了不同的分解层数,计算了降噪后的信噪比,通过添加高斯白噪声(SNR为5 dB)对信号进行仿真,以测试不同分解层数对声波信号降噪的影响。表2中展示了各分解层数在sym8小波基下对仿真信号的去噪效果。从表2可以看出,随着分解层数的提高,信噪比先增大后减小,这是因为分解层数过少时,信号去噪不彻底,去噪效果不理想。更高的分解层数能够带来更好的去噪效果,但是分解层数过高时,信号的有效信息丢失严重,反而会导致信噪比的降低,还会使得运算量增大,算法的处理速度变慢。对本文仿真信号而言,选择3层分解层数可以取得最高的信噪比,因此后续的仿真实验采用3层分解层数。
2.4 实验结果分析
为了验证本文所提方法对小波阈值去噪算法的提升效果,笔者利用数据集中的声波信号进行了噪声消除的仿真实验,对实验结果进行了分析。为了更好地接近光纤声波传感信号的实际信号,在仿真实验中以数据集声波信号作为模拟信号,向模拟信号添加噪声,将具有高斯分布的白噪声和真实光纤声波传感系统噪声作为噪声添加到模拟信号中。信噪比设置为 5 dB。原始信号及包含噪声的信号波形图如图6所示。原始信号及包含噪声的信号频谱图如图7所示。
与硬阈值函数小波变换和软阈值函数小波变换去噪方法相比,采用改进的阈值函数及改进模拟退火算法优化阈值的小波变换方法进行去噪的对比效果如图8、图9所示。
表3和表4中分别展示了加入5 dB信噪比的不同类型噪声后,本文改进方法和硬、软小波阈值去噪算法的信噪比、均方误差和波形互相关系数的值。从表中可以看出,本文所提方法在白噪声和真实噪声的仿真实验中三项数据的表现均优于原始算法。从信噪比的值可以看出,本文所提方法对小波阈值算法的去噪效果提升显著。
从表4中可以看出,在加入真实噪声后,本文所提算法依然能取得更好的效果,说明算法具有更强的鲁棒性,能够应对更加真实的噪声情况。
图10展示了不同信噪下两种噪声去噪后的趋势变化。
从图10可以看出,随着信噪比的提高,算法的去噪能力逐渐变弱,这是由于信噪比较低时,信号中含有的噪声分量较多,去噪带来的信噪比提升较大;而信噪比高时,信号中含有的噪声分量较少,去噪时容易过度去噪,造成有效信息的损失,从而使信噪比提升幅度变小,其中的小波硬阈值去噪方法在20 dB信噪比进行去噪处理时反而会导致信噪比的下降。
3 结束语
本文提出了一种新的阈值函数和阈值计算方法用于改进小波阈值去噪算法,旨在对光纤传感声波数据取得更好的去噪效果。新的阈值函数解决了软、硬阈值函数存在的缺点,保证了阈值函数的连续性。该函数结合经典的硬阈值函数和软阈值函数,通过形状调整因子调整整体速度,即用
α和β调整J,K的增长率来使其逼近软阈值函数和硬阈值函数,并且可以指定阈值处的截断点,具有更好的灵活性和可控性。本文提出的新阈值计算方法自适应地求解最优阈值,减小了小波阈值去噪方法对于阈值选择的依赖性,解决了小波阈值去噪算法过度依赖参数选择的情况。综合数据测试结果表明,所提改进方法能够抑制高频和低频噪声,有利于声波信号的准确识别,提高了声波信号的可用性。此外,与传统的小波阈值去噪方法相比,本文所提方法能够显著提高两种仿真声波信号的信噪比。在真实光纤声波传感信号噪声的仿真中进一步证实了所提出的改进小波阈值去噪方法可以有效对真实信号进行处理,证明了所提算法的鲁棒性。
参考文献:
[1] GONZALEZ-HERRAEZ M. Distributed acoustic sensing for seismic monitoring[C]//2021 Optical Fiber Communications Conference and Exhibition (OFC), San Francisco, CA, USA, 2021.
[2] SUN W H, ZHU S, LI W, et al. Noise suppression of distributed acoustic sensing based on f-x deconvolution and wavelet transform[J]. IEEE Photonics Journal, 2020, 12(1): 1-8.
[3] HICKE K, HUSSELS M-T, EISERMANN R, et al. Condition monitoring of industrial infrastructures using distributed fiber optic acoustic sensors[C]//2017 25th Optical Fiber Sensors Conference (OFS), Jeju, Korea (South), 2017.
[4] 李柯,龚岩栋,张卓.φ-OTDR降噪处理的研究进展[J].激光与光电子学进展,2022,59(23):25-33.
LI K, GONG Y D, ZHANG Z. φ- Research progress of φ-OTDR noise reduction processing[J]. Laser & Optoelectronics Progress, 2022,59(23): 25-33.
[5] 张芳. 扩频光时域反射技术理论与空间定位应用研究[D].秦皇岛:燕山大学,2021.
ZHANG F. Research on the theory and spatial positioning application of spread spectrum optical time domain reflection technology[D].Qinhuangdao: Yanshan University, 2021.
[6] 汪豫. 基于Ф-OTDR的光纤振动传感关键技术研究[D].成都:电子科技大学,2018.
WANG Y. Study on the optical fiber vibration sensing technology based on Ф-OTDR[D]. Chengdu: University of Electronic Science and Technology of China, 2018.
[7] 蒋浩. 基于Φ-OTDR的管道第三方破坏监测技术研究[D].杭州:浙江大学,2023.
JIANG H. Base on Φ- Research on third-party pipeline damage monitoring technology for OTDR[D].Hangzhou: Zhejiang University, 2023.
[8] 吕龙宝,庞拂飞,刘奂奂,等.基于奇异值分解的Φ-OTDR相位提取[J].光通信技术,2018,42(9):51-54.
LYU L B, PANG F F, LIU H H, et al. Phase-extraction Ф-OTDR based on SVD[J]. Optical Communication Technology, 2018,42(9): 51-54.
[9] 齐晓辉,宋宛鸿,李文召.一种改进的小波去噪方法在OTDR中的应用[J].光通信技术,2021,45(4):54-58.
QI X H, SONG W H, LI W Z. Application of an improved wavelet denoising method in OTDR[J]. Optical Communication Technology, 2021,45(4): 54-58.
[10]WANG P F, LV Y J, WANG Y,et al., Adaptability and anti-noise capacity enhancement for Φ-OTDR with deep learning[J]. Journal of Lightwave Technology, 2020,38,(23): 6 699-6 706.
[11]陈霄,徐慨,王潋,等.基于小波去噪和三阶全息图的时延估计算法[J].指挥控制与仿真,2015,37(4):69-73.
CHEN X, XU K, WANG L, et al. Time delay estimation algorithm based on wavelet denoising and third order hologram[J]. Command Control & Simulation, 2015,37(4): 69-73.
[12]WANG C, DENG C X, HU Z B. An improved wavelet threshold function and its application in image edge detection[C]//2019 International Conference on Wavelet Analysis and Pattern Recognition (ICWAPR), Kobe, 2020: 1-7.
[13]LI C, CAO D G YUAN Y Y. Research on improved wavelet denoising method for sEMG signal[C]//2019 Chinese Automation Congress (CAC), Hangzhou,2020: 5 221-5 225.
(责任编辑:张培培)
引用格式:师雪玮,徐大林,刘志成,等.基于改进小波阈值的声波信号去噪算法与仿真[J].指挥控制与仿真,2024,46(5):69-76.
SHI X W,XU D L,LIU Z C,et al.Research on algorithm and simulation of acoustic signaldenoising based on improved wavelet threshold[J].Command Control & Simulation,2024,46(5):69-76.
作者简介:
师雪玮(2000—),男,硕士研究生,研究方向为光纤传感和信号去噪。
徐大林(1964—),男,研究员。