远距离支援干扰下目标跟踪的PF-MHT融合策略:提高20%探测范围

需积分: 9 0 下载量 65 浏览量 更新于2024-08-11 收藏 910KB PDF 举报
本文档主要探讨了"远距离支援干扰下的目标跟踪技术",发表于2011年11月的《北京航空航天大学学报》第37卷第11期。在现代战争环境中,远距离支援干扰(standoff jamming, SOJ)对雷达的性能构成挑战,尤其是在干扰机掩护区域内,雷达难以有效检测目标,而其他区域的目标航迹也变得不稳定。针对这一问题,研究者徐海全、王国宏和关成斌提出了两种关键的算法:噪声干扰下的粒子滤波(Particle Filter, PF)检测前跟踪(Track Before Detect, TBD)算法以及门限自适应的多假设跟踪(Multiple Hypothesis Tracking, MHT)算法。 PF-TBD算法通过利用粒子滤波原理,在目标可能存在的概率空间中搜索,即使在干扰环境下也能提前预测和跟踪目标,从而减少丢失目标的风险。MHT算法则是一种多模式跟踪策略,它能够处理多个潜在目标的跟踪,增强了算法的鲁棒性和准确性。 作者将这两种算法进行了融合,设计了一种智能策略,使得算法可以根据雷达探测区域所受到的干扰程度自动选择合适的跟踪方法。这种策略优化了目标航迹管理,显著提高了雷达探测区域的有效覆盖范围,提升了约20%以上,同时保持了算法的高效执行,确保了实时性和响应性。 此外,论文强调了这种方法对干扰情况的高度适应性,不仅适用于远距离支援干扰,也适用于其他类型的压制干扰,对于复杂战场环境中的目标跟踪具有重要的理论参考价值和实用意义。论文的关键词包括远距离支援干扰、粒子滤波检测前跟踪、多假设跟踪和区域划分,这些概念和技术在现代军事和电子战领域中是前沿研究的重点。 这篇文章提供了对抗远程干扰环境下目标跟踪的有效解决方案,展示了在技术复杂度与实际效能之间找到平衡的能力,对于提高雷达系统的生存能力和作战效率具有显著的贡献。

%% 远距离支援干扰仿真 dread=pi/180; %角度转弧度 radeg=180/pi; %弧度转角度 Pt=1000e3; %%雷达发射功率w Gt=40; %%雷达天线主瓣方向上的增益db Gt_coef=10^(Gt/10); Pj=2000; %%干扰机发射功率w Gj=20; %%干扰机天线主瓣方向上的增益db Gj_coef=10^(Gj/10); gamaj=0.5; %%干扰信号对雷达天线的极化损失 sigma=5; %%被掩护目标的雷达有效反射面积 H=8e3; %%干扰机的天线与雷达的高度差m h=8e3; %%被保护目标相对干扰机高度m Dj=300e3; %%干扰机的天线与雷达的水平距离m Kj=10; dfjddfr=2; %%△fj和△fr分别干扰信号和雷达接收机的中频带宽 theta_mid=5; k=0.1; A=sqrt((Pt*Gt_coef*sigma*Kj*dfjddfr)/(4*pi*Pj*Gj_coef*gamaj)); Ht=H-h; Dt0=zeros(1,360); for i=1:length(Dj) Rj=sqrt(Dj(i)^2+H^2); for theta=-179:1:180 if(0<=abs(theta)&&abs(theta)<=theta_mid/2) Dt0(:,theta+180)=sqrt(A*Rj-Ht^2); elseif(theta_mid/2<abs(theta)&& abs(theta)<=90) Dt0(:,theta+180)=sqrt((A*Rj/(sqrt(k)*theta_mid))*abs(theta)-Ht^2); else Dt0(:,theta+180)=sqrt(90*(A*Rj/(sqrt(k)*theta_mid))-Ht^2); end end figure theta=(-179:1:180)*dread; polarplot(theta,Dt0); rticks(50e3:50e3:200e3); rticklabels({'r = 50km','r = 100km','r = 150km'}) end sigma1=[5 10 20 30 40]; figure for i=1:length(sigma1) A=sqrt((Pt*Gt_coef*sigma1(i)*Kj*dfjddfr)/(4*pi*Pj*Gj_coef*gamaj)); Djj=(50:10:600)*1e3; Rmax=zeros(1,length(Djj)); for j=1:length(Djj) Rj=sqrt(Djj(j)^2+H^2); Rmax(:,j)=sqrt(90*(A*Rj/(sqrt(k)*theta_mid))-Ht^2); end scatter(Djj/1e3,Rmax/1e3); hold on end title('干扰机雷达水平距离与最大探测距离的关系'); xlabel('干扰机距雷达水平距离/km');ylabel('雷达最大探测距离/km') xlim([50,600]) legend('5','10','20','30','40') %%% fig 5 hh=0:0.2:8; Rmax1=zeros(1,length(hh)); for l=1:length(hh) Ht=H-hh(l)*1e3; Rmax1(:,l)=sqrt(90*(A*Rj/(sqrt(k)*theta_mid))-Ht^2); end figure scatter(hh,Rmax1/1e3)转换成C代码

2023-07-11 上传