基于先验知识的机载雷达自适应训练样本选取的方法.pdf
《基于先验知识的机载雷达自适应训练样本选取的方法.pdf》由会员分享,可在线阅读,更多相关《基于先验知识的机载雷达自适应训练样本选取的方法.pdf(6页珍藏版)》请在咨信网上搜索。
1、第37卷 第2期2023年4月空 天 预 警 研 究 学 报Journal of Air&Space Early Warning ResearchVol.37 No.2Apr.2023收稿日期:2023-04-30基金项目:国防科技卓越青年科学基金项目(2019-JCJQ-ZQ-006)作者简介:李虎(1990-),男,硕士生,主要从事空时自适应处理研究基于先验知识的机载雷达自适应训练样本选取的方法李虎,谢文冲,熊元燚,侯铭(空军预警学院,武汉 430019)摘要:针对实际复杂环境下机载雷达面临的均匀样本数量严重不足的问题,提出了一种基于先验知识的机载雷达自适应训练样本选取方法;针对选取得到的
2、训练样本距离间隔较远导致明显功率差异的问题,本文方法对训练样本进行了进一步修正该方法从距离-多普勒单元角度出发,在考虑地形遮蔽影响的同时,基于地貌类型和地貌坡度构建了各训练单元的加权欧氏距离测度,实现了对均匀训练样本的有效选取仿真结果表明,相对于现有KA-STAP方法,基于本文所提样本选取方法的STAP处理后的杂波抑制性能得到显著改善关键词:先验知识;机载雷达;训练样本选取;空时自适应处理;非均匀杂波抑制中图分类号:TN958文献标识码:A文章编号:2097-180X(2023)02-0089-05空时自适应处理(STAP)技术对于机载雷达在下视工作时面临的强杂波信号具有良好的抑制效果,能够有
3、效提高预警系统的检测性能STAP 的核心在于构造出待检测单元的杂波协方差矩阵在实际中,杂波协方差矩阵是未知的,通常的做法是利用待检测单元两侧的训练样本做最大似然估计,用估计出的样本协方差矩阵代替未知的杂波协方差矩阵然而,外部环境的变化,比如地形变化与起伏、天气状况、季节变化、人为开发、离散强散射点等造成杂波分布非均匀,用来估计的训练样本与待检测单元的杂波分布不一致,使得估计出的样本协方差矩阵存在偏差,导致 STAP 杂波抑制性能严重下降为了解决非均匀杂波抑制问题,国内外学者从两方面展开了研究一方面,为了筛选出均匀训练样本,研究者提出了多种非均匀检测器来检测和剔除奇异样本,其中广义内积(GIP)
4、非均匀检测器1是典型代表但是 GIP 用待检测单元的估计协方差矩阵来表征杂波特性,并不能直接表征待检测单元自身的杂波特性,所选出的训练样本仅是与待检测单元杂波协方差矩阵结构上相近,有可能与待检测单元自身的杂波特性并不相符另一方面,知识辅助 STAP(KA-STAP)方法2-9利用地理环境的先验知识改善了非均匀环境下的杂波抑制性能,是当前 STAP 领域的研究热点在 KA-STAP 方法中,文献7提出的根据先验地形数据信息选取样本的方法是间接利用先验知识类 STAP 方法的代表本文在文献7方法基础上,提出一种基于先验知识的自适应训练样本选取方法该方法首先利用机载雷达平台参数、地形地貌参数、地理高
5、程参数构建先验知识信息数据库;其次确定各距离-多普勒(RD)单元对应的信息点;再次构建基于地貌类型和地貌坡度的加权欧式距离测度,完成均匀样本的选取;最后对选取的样本进行功率差异性修正1经典的根据地貌类型属性选取训练样本方法经典的根据地貌类型属性选取训练样本方法7的思路是,首先将单个距离单元按照多普勒划分为若干个 RD 单元;然后将 RD 单元的地理信息量化为地形矢量,地形矢量中归类每种地形所占的比例,基于训练 RD 单元与待检测 RD 单元的地形矢量求二者之间的欧氏距离;最后利用距离测度挑选样本图 1 为文献7的 RD 单元划分区域经度abcd纬度O图 1文献7的 RD 单元划分区域文献7方法
6、的具体步骤如下:DOI:10.3969/j.issn.2097-180X.2023.02.002空 天 预 警 研 究 学 报2023年90Step 1 找出 RD 单元中 4 个顶点的最大最小值经纬度,而后划定矩形范围作为该 RD 单元的范围,矩形方框中的圆点为该 RD 单元区域的信息点(本文的信息点是指以经纬栅格形式存储的地表覆盖信息与高程信息的点)Step 2 根据中国科学院空天信息创新研究院提供 的 全 球 30m 精细 地 表 覆 盖 数 据(GLC_FCS30-2020)10-11统计 RD 单元的地貌类型属性由于该数据库包含有 29 种地貌类型,则该 RD 单元的地貌类型属性可以
7、用一个包含 29 个元素的向量来表示(每一个元素对应一种地貌类型),将该地貌类型的信息点数除以 RD 单元内总的信息点数,得到该地貌类型在 RD 单元中的占比,从而形成含有 29 个元素的归一化向量第 l 个待检测 RD 单元的地貌类型归一化向量可以表示为tl=tl1tl2 tl29TStep 3 计算第l个训练单元与第 l 个待检测RD 单元的地貌类型归一化向量的欧氏距离为dl=i=129(tli-tli)2;选取距离最小的前 L 个样本作为训练样本,L 需满足 RMB 准则12Step 4 利用选取的训练样本估计杂波噪声协方差矩阵经典的根据地貌类型属性选取训练样本的方法存在以下不足:RD
8、单元的划分区域粗糙,将矩形虚线框内的所有信息点默认为 RD 单元的信息点,而且这种划分方式只适合地形平坦的区域,当地形出现较大起伏时,矩形框中的点由于具有不同的海拔高度导致到载机雷达的视距发生明显变化,可能不在同一距离单元内;选取的地貌类型属性分类太细,如阔叶林、针叶林、混叶林的经验后向散射系数差别不大,将其单独分类时,选取的训练样本会放大相近地貌类型间的细微差异,不能得到足够数量的有效样本;在求取欧氏距离时,各地貌类型赋予相同的权重,忽视了不同地貌类型间后向散射系数区别大导致回波功率区别大的特点,从而选取的训练样本不够准确;该方法是从地貌类型角度出发,考虑待检测 RD 单元与训练单元的杂波特
9、性,忽略了地貌坡度变化对后向散射系数的影响;未考虑距离、天线方向图调制对 RD 单元回波功率造成的影响2基于先验知识的自适应训练样本选取方法针对以上经典的根据地貌类型属性选取训练样本方法不足,本文提出一种基于先验知识的自适应训练样本选取方法2.1信息点与杂波块的映射关系将地理信息数据库中的每一个信息点精准对应到杂波信号模型划分的杂波块,是基于地理信息进行训练样本筛选的重要前提首先完成数据库中各信息点和载机的大地坐标系数值到平面直角坐标系的转换13,然后统一在平面直角坐标系下将各信息点映射到不同距离单元考虑到地球曲率的影响,为了使数据库与雷达杂波模型相匹配,计算信息点对应的俯仰角时应考虑地球曲率
10、的影响假设地球等效半径为Re,载机坐标为(xayaHa),信息点坐标为(xcycHc),有=arcsin(Re+Ha)2+R2-(Re+Hc)2)/(2R(Re+Ha)Rcos=(xa-xc)2+(ya-yc)2)1/2(1)通过上述方程组即可求得信息点相对于载机的距离 R 和俯仰角信息点对应的方位角为=arccos(xc/(xa-xc)2+(ya-yc)2)1/2)(2)根据得到的各信息点的距离和方位角参数,即可将信息点映射到对应的距离单元内2.2基于先验知识的自适应训练样本选取方法本文所提的知识辅助训练样本选取方法基本思路为:首先确定各 RD 单元对应的信息点,并将受地形遮蔽影响的信息点进
11、行剔除;然后基于地貌类型和地貌坡度构建各 RD 单元的加权归一化向量,通过欧式距离选取出距离较近的若干个 RD 单元作为训练样本;最后对选取得到的训练样本进行俯仰天线增益和距离导致的功率差异性修正具体方法步骤如下Step 1 RD 单元划分本文提出的基于地理信息的样本选取方法是针对当前实际装备中采用的多普勒后处理STAP 方法该类方法是针对每一个待检测多普勒单元,在训练样本中选取独立同分布的多普勒单元作为训练样本因此将相同距离环内的信息点按照多普勒频率重新进行划分,划分后的单元称为 RD 单元假设第 l 个距离单元对应的距离为Rl,俯仰角为l,则l=arcsin(Ha/Rl+(R2l-H2a)
12、/(Rl(Re+Ha)(3)该距离环上杂波块的最大归一化多普勒频率为fl=2Vrcosl/(fr)(4)式中,Vr为载机速度,fr为雷达脉冲重复频率,第2期李虎,等:基于先验知识的机载雷达自适应训练样本选取的方法91为雷达工作波长因此,第 l 个距离单元的归一化多普勒频率取值区间为-fl fl假设划分的多普勒单元数为 J,则第 j 个多普勒频点的频率公式为flj=-fl+2fl(j-1)/Jj=1J+1(5)对应的方位角为lj=arccos(fljfr/(2Vrcosl)(6)则第 l 个距离单元的方位角集合为l=l1l2lJ+1以J+1个方位角组成边界点,形成 J个多普勒单元根据信息点对应的
13、方位角信息,即可将信息点映射到对应的多普勒单元Step 2 统计地貌类型分类占比,形成地貌类型归一化向量经过合并后各 RD 单元的地貌类型属性可以用一个包含 8 个元素的向量表示14每一个元素代表一个地貌类型,剔除受遮蔽的信息点15后,将该地貌类型的信息点数除以 RD 单元内总的信息点数,得到该地貌类型在 RD 单元中的占比,从而形成含有 8 个元素的归一化向量,则第l 个 RD 单元的地貌类型归一化向量表示为Il=i1i8TStep 3 统计地貌坡度分类占比,形成地貌坡度归一化向量坡度系数是表示地表单元陡缓的程度,反映了地面在该点的倾斜度信息点坡度模型如图 2 所示载机位置dkBB1B2X3
14、Y3Z3O3图 2信息点坡度模型设B1、B2表示与B相邻的信息点,k为由3个信息点组成的面元的法线向量,d为平面O3X3Y3的法线向量,为向量k和d的夹角则信息点B的坡度为=arccos(kd)(7)经过大量实验对数据统计分析得到的经验值发现16,对于同一地貌类型,不同的坡度范围对应有不同的经验散射系数通常将坡度值分为 3 个区间,即1、12因此,RD 单元的地表坡度属性可以用一个包含 3 个元素的向量表示,每一个元素代表一个坡度范围剔除受遮蔽的信息点后,将该坡度范围的信息点数除以 RD 单元内总的信息点数,得到该坡度范围在 RD 单元中的占比,从而形成含有 3 个元素的归一化向量则第 l 个
15、 RD 单元的地貌坡度归一化向量表示为Ql=q1 q2 q3TStep 4 合并 RD 单元的地貌类型与地貌坡度归一化向量形成新归一化向量,计算对应欧氏距离,选取均匀训练样本同时考虑地貌类型和坡度后的第 l 个待检测 RD 单元对应的归一化向量为Wl=Il;Ql,式中合并后的归一化向量的维数为 111第l个训练单元与第 l 个待检测距离多普勒单元之间的欧氏距离为El=i=111(Wli-Wli)2(8)式中,Wl=WlAl,表示 Hardmard 积;Al=a1 a2a81/31/31/3T,ai为第 i 种地貌类型相对于所有 8 种地貌类型的归一化散射强度通过权值Al,增大了散射系数较大的地
16、貌类型的权重由式(8)得到各训练单元的欧氏距离,将其从小到大排序,选出距离最小的前 L 个样本作为训练样本,其中 L 需满足 RMB 准则Step 5 修正影响因素,估计杂波噪声协方差矩阵当根据地形地貌等先验信息进行样本选取时,选取的训练样本距离待检测单元有可能较远,导致训练样本与待检测单元之间存在较大的功率差异因此需要对选取的训练样本数据进行修正来补偿功率非均匀的影响从雷达距离方程可知,影响不同 RD 单元杂波功率的因素包括距离、所包含信息点对应的擦地角、俯仰天线方向图等假设经过空时滤波后的待检测 RD 单元的回波功率为Pj,则对待检测距离单元的 J 个多普勒频点的功率值求平均,即P=1Jj
17、=1JPj(9)同 样 可 得 到 第 l 个 训 练 单 元 的 平 均 功 率Pl通过待检测距离单元回波功率对第 l 个训练单元进行归一化处理,即cl=Pl/P(10)经过修正后的第 l 个距离单元第 k 个多普勒单元的回波信号为xlk=xlk/cl(11)式中,xlk为第 l 个距离单元、第 k 个多普勒单元的原始回波信号Step 6 利用修正后的训练样本数据估计杂空 天 预 警 研 究 学 报2023年92波噪声协方差矩阵,形成空时自适应权值,完成空时自适应杂波抑制处理3仿真实验假设载机在某地向正北飞行,波束指向为正东方向,载机飞行中的偏航角、俯仰角、横滚角均为 0,雷达采用正侧视均匀
18、线阵,具体仿真参数:载机高度为 8km,载机速度为 140m/s,波束指向0,目标俯视角为 1,峰值发射功率为 230kW,工作波长为 0.23 m,脉冲重复频率为 2 434.8 Hz,接收机带宽为 5MHz,相参脉冲数为 16,接收子阵数为 16机载雷达探测区域位于岛屿与海域交界处为了便于比较本文所提样本选取方法的有效性,本文基于 1DT 方法进行验证,由 RMB 准则可知,仅需要 32 个训练样本即可设定待检测单元距离门编号为 101,仅在待检测单元左右两侧相邻的各 100 个距离门内选取训练样本1)地貌坡度情况图 3 给出了机载雷达探测区域内地貌坡度情况(其中蓝紫色区域坡度值较小,地形
19、平坦;黄绿色区域坡度值较大,地形陡峭;海域区域没有坡度起伏)从图 3 可以看出:该区域地表坡度分布范围为 060,坡度分布不均匀;待检测距离单元的第 13 号多普勒通道处为海岛、水域交错分布区域,该待检测 RD 单元的地貌坡度取值最不均匀;本文方法中利用地貌坡度分类占比进行样本选取是合理的待检测距离单元1 34 08 01 2 01 6 02 0 00距离单元00.5 1.0 1.5-0.5-1.0-1.5方位角/()1 02 03 004 05 06 0图 3机载雷达探测区域内地貌坡度情况2)影响因素修正效果图 4 给出了第 10 号多普勒通道经过功率影响因素修正前后的效果对比由图 4 可知
20、,修正后的数据较修正前功率更均匀,体现在 2 个方面:待检测多普勒通道的功率为 117.375dB,修正前多普勒通道的均值为 120.154dB,修正后多普勒通道的均值为 118.308dB,修正后的均值更接近待检测通道功率;修正前数据方差为 81.2,修正后数据方差为 52.8,修正后的数据方差显著减小04 08 01 2 01 6 02 0 0距离单元8 01 0 01 2 01 4 01 6 0多普勒单元功率/d B(a)修正前04 08 01 2 01 6 02 0 0距离单元8 01 0 01 2 01 4 01 6 0多普勒单元功率/d B(b)修正后图 4功率影响因素修正前后效果
21、对比3)信杂噪比(SCNR)损失图 5 从 SCNR 损失角度给出了分别选用 32个训练样本和 60 个训练样本情况下的修正效果(fd/fr为归一化多普勒频率)从图 5 可以看出,选取 60 个样本后的修正效果更为明显,其原因是经本文方法选取的 60 个样本中包含的离待检测单元距离远的样本相对更多,功率差异较大,因此修正后的补偿效果更好-0.3-0.5-0.1 0.10.30.5fd/fr-1 2 0-1 0 0-8 0-6 0-4 0-2 00S C N R 损失/d B修正前修正后(a)32 个样本-0.3-0.5-0.1 0.10.30.5fd/fr-1 2 0-1 0 0-8 0-6
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 先验 知识 机载 雷达 自适应 训练 样本 选取 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。