第2期 2013年4月 雷达科学与技术 Radar SCience and TechnoIogy Vo1.n No.2 April 2O13 基于极化不变量特征的雷达目标识别技术 王福友 ,罗钉 ,刘宏伟 (1.中航工业雷达与电子设备研究院,江苏无锡214063; 2.西安电子科技大学雷达信号处理国家重点实验室,陕西西安710071) 摘 要:极化是雷达目标具有的特性之一。以电磁散射计算仿真的圆锥形弹头模型、球形和圆柱形诱 饵模型为研究对象,在极化不变量理论基础上对这些简单目标的极化特性进行了试验分析研究,提出了一 种新的组合极化不变量特征(功率矩阵迹与行列式的比值)用于雷达目标识别,并给出了其对应实际的物理 意义。文中以SVM为分类器,提出基于功率矩阵迹、去极化系数和功率矩阵迹与行列式的比值特征进行分 类识别,结果表明,该方法可以有效地将弹头和诱饵进行分类识别。 关键词:雷达目标识别;极化;功率矩阵迹;去极化系数;功率矩阵迹与行列式的比值 中图分类号:TN958 文献标识码:A 文章编号:1672—2337(2013)02—0165—08 Radar Target Classification Based on Some Invariant Properties of the Polarization WANG Fu—you -一.LUO Ding .LIU Hong—wei。 (1.AVIC Radar and Avionic Institute,Wuxi 214063,China;2.National Key Lab of Radar Signal Processing,Xidian University,Xi’an 710071,China) Abstract:Polarization is one of the main target features.This paper uses electromagnetic scattering computation to model cone warhead,bait(sphere and column).Experiment analysis of the invariant proper— ties of the polarization scattering matrix(IPPSMs)for the simple targets is given.A new IPPSM(the ratio for trace of power scattering matrix to determinant)is first presented and its physical characteristic is also given. SVM is used as classifier,and trace of power scattering matrix,depolarization and TPSMD are presented for target classification.The recognition results validate the effectiveness of theoretical analysis. Key words:radar target recognition;polarization;trace of power scattering matrix;depolarization; trace of power scattering matrix to determinant 重要的军事意义L6。]。 1 引言 作为战术打击和战略威慑武器的中、远程各 型弹道导弹的研发一直受到各主要军事强国和新 兴国家的重视l_1]。而从弹道导弹研制成功不久后 雷达回波的极化特征是继幅度、相位和多普 勒频移之后的另一个重要特征。由于极化信息处 理使用了雷达回波的矢量信息,使得许多在时域、 频域、空域用常规手段难以解决的问题有了解决 办法。利用极化信息能提高目标检测信噪比、抗 干扰能力和目标识别能力 等。因此,作为时域、 的1950年开始,弹道导弹防御技术就成为美苏等 国的研究热点。目前,美国等国已经建立起实用 的弹道导弹防御系统 ]。 频域和空域信号处理技术的补充,极化信息处理 技术对目标识别具有潜在应用价值。 极化是电磁波的固有属性,被认为是与目标 几何外形直接相关的一种特征,而全极化雷达往 往被认为可以接收雷达回波的所有能量,得到了 极化散射矩阵,更与目标特性紧密联系。极化用 于目标识别的一个最简单的例子就是极化可以有 弹道导弹防御面临的环境和技术日益复杂, 如何实现在导弹自由飞行段和再入段的预警和识 别,并引导我方导弹进行拦截,是弹道导弹防御的 关键,雷达目标识别技术的研究对于导弹防御系 统和战略、战术弹道导弹突防技术的发展均具有 收稿日期:2012—09—20;修回日期:2012—11-20 166 雷达科学与技术 第11卷第2期 效识别球形诱饵,其RCS可以和弹头迎头方向的 RCS近似相等,运用传统的单一极化雷达很难区 分弹头和球形诱饵,而全极化雷达可以有效区分 弹头和球形诱饵,因为球在不同方位的极化散射 矩阵近似相同,而弹头在不同方位的极化散射矩 阵是和诱饵球不同的。 从公开的文献上看,当前国内外研究得比较多 的是极化SAR,主要是用于遥感等,取得了较成功 的研究 ]。国外鲜有报道极化用于弹道导弹目标识 别,但据预测美国反导雷达具有极化识别的能力,目 前国内具有极化能力的雷达较少,主要是成本较大, 但随着科技的发展以及工艺的改进,相信具有全极 化测量能力的极化雷达会得到进一步的应用。 Bickel较早地提出极化不变量的理论l】 ,但 没有给出试验验证和分析,本文在电磁散射计算 仿真数据基础上,对锥形弹头模型、球形和圆柱形 诱饵模型的极化不变量特征进行了分析,提出了 新的极化不变量特征,得出了一些实用的结论。 2极化不变量目标识别原理 极化与目标的几何外形直接相关,因此本文 主要是依据极化不变量特征进行目标分类识别。 分时极化或同时极化情况下,可以采用测量极 化散射矩阵的方法进行目标识别。在H—V基下, 信号的发射电场和散射电场存在以下线性关系: E =[s]E. (1) 式中,Ei和E 分别为人射波和散射波的极化状 态,Es]为2×2的复矩阵(Sinclair矩阵),得到极 化散射矩阵: s一[Ss、H, .喜、H,V ] (2) 3对于互易性目标,可以近似地认为 SHv≈SvH (3) 通过极化散射矩阵进行相应变换可得到极化 不变量: (1)行列式模值l△l l△l—l desES]l—l s S 、,一S 、,I (4) (2)功率矩阵迹P P 一l SHH l +l S、,、,I。+2 l SHv I。 (5) (3)去极化系数D D一1一 (6) (4)本征极化方向角 。 。一 _ar1 ctarctan an ㈩ ) 式中, fS1一SHH+Svv S 2===SHH—Svv (8) lS12一SHv (5)本征极化椭圆率r。 1ro 一 一 arctan n ■ (9)) 式中, S'l2一S12 cos(2q ̄。)一÷s 2 sin(2q ̄。) (1o) 由于功率矩阵迹P 和行列式△均为不变量, 本文提出一个组合的不变量系数 一 P I S f +1SHHSvv—S v S、,、,1 +2 1 S 、,1 (11) 表1给出了以上极化特征对应的物理意义,具 有明确的物理意义,可能更能有效地利用极化特 征进行目标识别。 表1 极化散射矩阵提取的特征对应的物理意义 极化特征 物理意义 行列式模l△l 行列式值j△粗略反映了目标的粗细或“胖瘦” 功率矩阵迹P1 P 代表了一对正交极化天线所接收到的总功 率,它大致反映了目标的大小 D大致反映了目标散射中心的个数: 去极化系数D ①0<D<0.5时对应的是一孤立的散射中心 目标;②0.5<D≤1时对应的是多散射中心 的组合体目标 本征极化方向角伽 90反映了目标极化方向 本征极化椭圆率ro ro反映了目标对称性差异 组合的不变量系 表示目标的纵横比 仿真试验结果 由于极化涉及复杂的电磁散射机理和目标特 性,本文采用电磁散射计算仿真数据进行分析,对 于纯金属电磁计算仿真的结果和微波暗室测量的 结果相当。根据文献[63和文献[12]所述,导弹弹 头的外形可能是圆锥体,因此本文选取圆锥体作 为弹头模型,而由文献[11]和文献E12]可知,诱饵 可能是充气球体,此外考虑长度等因素,本文还选 择了圆柱体,圆柱体与圆锥体等高,底面直径也相 等。图1给出了本文电磁散射计算所用的三种目 标模型示意图,分别为圆锥、球、圆柱。极化方式 采用线极化HH、VV、HV、VH极化。 2013年第2期 王福友:基于极化不变量特征的雷达目标识别技术 167 (a) 圆锥体(仿真弹头) (b)球 (c) 圆柱 图1简单目标模型示意图 3.1 电磁计算仿真试验模型及参数 3.1.1 圆锥体 圆锥体高为1.5 m(包含锥尾5 cm高圆柱部 分,由于该部分相对较小,本文将该目标统称为圆 锥体),锥顶半角为9。,底面半径为0.24 m,频率为 8.5~10.5 GHz,频率步进为2 MHz,方位角为 0。~180。,角度步进为0.2。,O。对应于锥尖,180。对 应于锥尾。 3.1.2球体 在各种极化模式下,金属球的RCS与视角无 关,其仿真随频率变化的RCS数据。球的RCS 结果是一致的,因此本文测量的球目标尺寸是 以波长归一化的目标特征尺寸大小的参数 ka_l ,ka∈(0,500),忌为波数,愚一27【 一 2 ̄f/c,a为球的半径,尼n步进为0.1。 3.1.3 圆柱体 圆柱和圆锥目标高和底面直径相同,圆柱体 高为1.5 m,底面半径为0.24 m。频率为8.5~ 10.5 GHz,频率步进为2 MHz,方位角为O。~90。, 角度步进为0.1。,0。对应于底面,9O。为垂直对 圆柱。 3.2 电磁计算仿真RCS结果与分析 图2(a)~(c)给出的是图1(a)~(c)所对应的 简单目标的不同极化下的RCS。由于圆锥体是标 准体,满足互易性,所以图2(a)给出的是HH、 VV、HV极化下频率为10 GHz的RCS; 40 … HH极化 20 I ……vv极化 鲁 奄 0 八 .HV极.. 化 8 一20 —柏 l l 1唧 『 .60 0 45 90 135 180 方位角/(。) (a) 圆锥体的RCS 。 -2。 g 墓-6o HH极化 8。 一loo g∞∞ 、8《 -120 0 l00 ∞ 200 0 300 ∞ 400 加 500 /ca (b)球体的RCS 一 由于球体为各向同性,其同极化HH和VV 极化近似相等,HV和VH极化相等,所以图2(b) 只给出球的HH极化和HV极化模式下频率为 10 GHz的RCS回波,由图2(b)可知,交叉极化 RCS接近于零(接近一80 dB)。需要说明的是当 忌口<0.5时,目标处于瑞利区,当0.5≤ka≤20 时,目标处于谐振区,由图2(b)可知,在谐振区内, 目标的RCS随频率变化震荡,当忌n>20时,此时 目标的RCS主要取决于其形状和表面的粗糙程 度,由于是球体,在光学区类似于一条直线,同时 也验证了文献[12]对瑞利区、谐振区和光学区的 划分理论。 168 雷达科学与技术 第l1卷第2期 图2(c)给出的是圆柱全极化模式下的RCS, 可以看出HV极化和VH极化的RCS近似相等; 由图2(a)和(c)可知,在全方位变化下,任意一种 目标的同极化的RCS起伏趋势一致,RCS值近似 相等,交叉极化也具有相同的结论,同极化与交叉 极化的差值在2O dB左右。 由图2可知,除了个别方位外,圆锥体、球体、 圆柱体的同极化的RCS接近,这对于传统的单一 极化雷达是很难通过RCS对这三类目标进行分类 识别的,此外由于圆柱的长度和圆锥体的长度相 同,传统的一维像通过测量目标的长度也将很难 区分这两类目标,这也说明用球体和圆柱体作诱 饵是可能的。 4特征提取与识别 4.1特征提取与分析 依据电磁计算仿真数据提取极化散射矩阵, 由极化散射矩阵提取极化不变量特征,图3~图5 分别给出了圆锥体、球、圆柱极化不变量特征:行 列式模值l△l、功率矩阵迹尸。、去极化系数D、本 征极化方向角 。、本征极化椭圆率r。和本文提出 的新的组合极化不变量特征功率矩阵迹与行列式 的比值玑 由图3~图5(a)、(b)可知,I△l和P 特征 除了量值的不同,其量值的变化趋势是一致 的。从图3(c)的D特征除了个别方位,在全方 位0。~180。范围内,0.6<D≤1,说明圆锥体 是多散射中心目标,散射中心个数≥2,这也验证 了文献[12]的观点:一般认为平底锥弹头有3个主 要散射中心。由图4(c)可知,D≈0,依据文献 Elo]可知目标是单散射中心,说明球体是各向同 性,与雷达视角无关。由图5(c)可知,圆柱体的D 特征在0~1的震荡变化,说明圆柱的散射中心随 着方位角是变化的,有的方位角圆柱的散射中心 是1个,有的方位角的散射中心的个数≥2,并没 有规律性。 由图3~图5(d)、(e)可知,三类目标的 。和 r。特征也具有较明显的差异,但在实际应用当中, 大型极化雷达很难保证极化相位的测量精度,而 这两个特征与相位关联较明显,因此本文对于这 两个特征不作深入分析,也不作为本文分类识别 所用特征。 由图3~图5(f)可知,三个目标的'7特征也存 在差异,球体目标近似是一条直线,说明其纵横比 不变,而圆柱和圆锥体有一定的差异,尤其是在方 位角(45。~9O。)区间。 40 20 0 ∞ -口 司-20 f 6 —— -40 《”f 『 _6o -80 0 45 90 l35 180 方位角,(。) (a)行列式模值l△l 60 40 20 I 葛0 \ —20 ^ _40 1} 嗍 .60 —80 0 45 90 135 180 方位角,(。) (b) 功率矩阵迹P … Q 45 90 135 l80 方位角,(。) (c)去极化系数D 方位角,(。) (d)本征极化方向角 o 17O 20 雷达科学与技术 _ 第11卷第2期 0 1 f ∞ 、一 \ —20 n 曲 ∞ 0 珈 枷 _口 =一40 厂\ 8 司 —— 一60 I f —80 懒 1 —100 30 60 90 方位角,(。/ 行列式模值l△l { 1 ^/ I 。, Ev 30 一 一 6O 0 方位角,(。)枷 功率矩阵迹P 方位角,(。) (c)去极化系数D I — r 如 p l ’I 。l 叫4 11 30 60 方位角,(。) 本征极化方向角po 岬 . 1 I30 60 9O 方位角/(。) 本征极化椭圆率r。 30 24 ∞ l8 12 6 j f ..¨ 上l L 0 30 60 90 方位角,(。) (f)行列式与功率矩阵迹比值日 图5 圆柱体的极化不变量特征 4.2识别结果及分析 图6给出了极化目标分类识别的流程图。根 据式(7)~(11)可知, 。特征和r。特征容易受到 相位的影响,而大型极化雷达很难保证极化相位 测量精度,因此本文提出采用与幅度相关性较大 的P 、D以及叼特征进行分类识别,而目标的几何 外形特征主要取决于极化的幅度信息。 结果 图6极化目标分类识别流程图 4.2.1 SVM分类器 特征提取后接下来就需要进行分类识别,而 目标分类识别离不开稳定的分类器,SVM(Sup— port Vector Machine)相比其他分类器泛化能力较 2013年第2期 王福友:基于极化不变量特征的雷达目标识别技术 171 强,收敛速度较快,且对小样本训练效果也较好, 故本文采用SVM分类器,其核函数采用RBF。 4.2.2基于双特征分类识别 由图3~图5(a)、(b)可知,l△}特征和P 特 征相关性较近,因此这两个特征选取P 特征,再 加一个D特征,形成双特征。在本文分类识别过 程中,SVM训练样本为在特征数据中随机抽取 150个,再用SVM训练后的分类线对全部数据进 行分类。图7给出基于SVM分类器的P 和D双 特征的分类识别结果,分类线可以较好地将两类 目标进行分类。表2给出了基于P 和D双特征的 目标分类识别结果。 分类正确率(圆柱:93.17%,球:95.24%) l 0 = 0.8 Q 0 6 雾 I .04 麓 0.2 / 。~\ \ 胜 0 120 —9O—60 —30 0 30 P1/dB (a) 球体和圆柱体的分类识别结果 分类正确率(圆锥体:98.71%,球:95.49%) ’。 秘 { +圆锥体 Q 1 球 f\ Pl/dB (b) 圆锥体和球体的分类识别结果 分类正确率(圆锥体:84.92o/o,圆柱:76.05%) ? Q 籀 J+圆锥柱体【 .{ + {毒 + + Pl/dB (c) 圆锥体和圆柱体的分类识别结果 图7基于P 和D双特征SVM分类识别结果 表2几种简单目标的分类识别结果 4.2.3基于三特征分类识别 由图7(c)可知,圆锥体和圆柱的分类效果相 对于前两组效果差一些,也为了验证本文提出的 新的极化不变量叼特征的有效性,图8给出了基于 P 、D和田三特征的分类识别效果图。SVM训练 样本为在特征数据随机抽取150个,由于是三维 图,无法画出分类界限。表3给出了基于这三个特 征的分类识别结果。 分类正确率(圆锥体:89.7%.圆柱:83.37%) 50 笔0 一50 —100 1.0 图8基于P 、D和口三特征的SVM分类识别结果 表3几种简单目标的分类识别结果 4.2.4分析和讨论 由图4和表2可知,球体和圆柱体、圆锥体和 球体的分类识别率较好,平均识别率分别为 94.21%和97.1O ,分析这是因为球体和圆柱体、 球体和圆锥体存在较大的几何差异,而极化与目标 的几何外形有关。而相对前两组,圆椎体和圆柱 的分类识别效果就差了些,平均识别率为80.49, 这主要是因为圆锥体与圆柱长度和底面直径都相 等,几何外形相近程度比球体的相近的更大。 由表3可知,加入本文提出的组合极化不变量 172 雷达科学与技术 第11卷第2期 钾特征,目标识别的结果是,圆锥体为89.7 9/6,圆柱 体为83.37 ,都有不同程度的提高,平均识别率为 86.54 ,相比基于双特征识别平均识别率80.49 , 提高了6个百分点,验证本文提出该特征的有效性, 此外考虑稳健性,建议采用三特征分类识别。 5 结束语 本文在电磁计算仿真模型圆锥形仿真弹头、 球形和圆柱形诱饵数据基础上对极化不变量特征 进行了提取,并提出了一种新的极化不变量,并给 出了其对应的物理意义,文章利用SVM分类器, 提出基于双特征和三特征对这几种目标进行分类 识别,分类试验结果验证了本文方法的有效性。 得出的主要结论如下: (1)在全方位下,圆锥体和圆柱体的任何一种 目标的同极化的RCS起伏趋势一致,RCS值近似 相等,交叉极化也具有相同的结论,同极化和交叉 极化的差值在20 dB以上; (2)通过RCS很难将弹头目标和球形诱饵及 圆柱诱饵进行分类,尤其是圆柱诱饵更具有迷惑性; (3)在全方位模式下,锥形弹头的散射中心个 数≥2,球体诱饵的个数散射中心个数为1个,圆柱 形诱饵的个数随着方位角是变化的,有时是1个有 时是多个(≥2); (4)本文提出的双特征法(功率矩阵迹、去极 化系数)和三特征法(功率矩阵迹、去极化系数、功 率矩阵迹与行列式的比值)可以有效地将弹头与 球形诱饵、圆柱形诱饵进行分类识别,平均识别率 分别为90 和80 以上,考虑到稳健性及识别率 建议采用三特征联合分类识别效果会更好;本文 方法可能用于弹头和诱饵的分类识别中。 参考文献: E1]Pridotkas J J.ICBMS—Their Relevance to US Security and the United States Air Force in the 2 1 st Century, ADA424619JR].Alabama:Air University Press Max— well AFB,2003:1-40. -12]Isby D C.No Go-Ahead Yet for Russian Heavy ICBM [J].Jane’S Missiles and Rockets,2010,14(12):14—17. [32 Aydin A T.Strategic Missile Defense:A Threat to Future Nuclear Arms Reductions[R].Monterey, USA:Arms Control Association,2O11:1-11. [4]Ingwersen P A,Lemnios W z.Radars for Ballistic Missile Defense Research[J].Lincoln Laboratory Journal,2002,12(2):245—266. E52 Crawford J F,Reed E W,Hines J J,et a1.Ground Based Radar-Prototype Antenna[C]//lEE National Conference on Antennas and Propagation,York,UK.[S.n.], 1999:249—252. E6]王璐,刘宏伟.基于时频图的微动目标运动参数提取 和特征识别的方法[J].电子与信息学报,2010,32 (8):1812-1817. WANG Lu。LIU Hong—wei.Method for Micro—Mo— tion Target Recognition and Motion Parameter Extrac— tion Based on Time-Frequency Analysisl-J].Journal of Electronics 8L Information Technology,2010,32(8): 1812—1817.(in Chinese) [7]陶勇,胡卫东,徐劲.基于轨道匹配和改进的空间目标 识别方法口].雷达科学与技术,2006,4(3):134—139. TAO Yong。HU Wei—dong,XU Jin.A Method of Space Object Identification Based on Orbit Matching and Improving[J].Radar Science and Technology, 2006,4(3):134—139.(in Chinese) -183戴幻尧,李永祯,王雪松,等.有源多假目标干扰的极化 识别新方法I-J].航天电子对抗,2010,26(5):19—22. Dai Huanyao,LI Yong—zhen,Wang Xue—song,et a1. Polarization Discrimination of Multi——False-Target Jam—— ming Based on Pattern Difference Structure of Sum and Difference Beam FJ].Aerospace Electronic war- fare,2010,26(5):19—22.(in Chinese) [9]Arii M,Zyl J J,Yunjin K.Adaptive Model—Based De— composition of Polarimetric SAR Covariance Matrices I-J].IEEE Trans on Geoscience and Remote Sensing, 2011,49(3):1104—1113. ri0]Bickel S H.Some Invariant Properties of the Polari— zation Scattering Matrix[J].Proceedings of the IEEE,1965,53(8):1070—1072. [11]Lisbeth G,George N L,David C W.The Continuing Debate on National Missile Defenses[J].Physics To— day,2000,53(12):36—42. -112]黄培康,殷红成,许小剑.雷达目标特性l-M].北京: 电子工业出版社,2005:16—18. 作者简介 - 王河体理会I《雷达学报》的审稿专家。S市T福技议P友术发,博检表、,士雷在索男论后达国,文核1目内主9心3标8核O要期识心余年研刊别期篇生究《和刊,方于信 多雷和黑向号篇达国龙为处被信江雷际理E省达号学I》讷总处术和 E—mail:hellofew@163.corn