微小卫星飞行姿态动环境规划算法研究

周梦晓，蔡倩，周建江，蒋晖，胡容

（南京航空航天大学，江苏南京210016）

微小卫星飞行姿态动环境规划算法研究

周梦晓，蔡倩，周建江，蒋晖，胡容

（南京航空航天大学，江苏南京210016）

为了提高微小卫星的使用效率和生存能力，关键是提高微小卫星的低可观测特性，设计了一种以混沌算法和粒子群算法为基础的微小卫星低可观测飞行姿态动环境规划算法。通过对微小卫星飞行中的俯仰角和方位角等姿态角进行实时调整，可以有效降低微小卫星在威胁雷达方向上的RCS值，提高微小卫星的低可观测特性。粒子群优化（PSO）算法可以降低计算复杂度以提高规划的实时特性，加入混沌（chaos）运动可以提高算法的精确程度。通过对工作频率在VHF波段的威胁雷达对微小卫星的威胁性进行仿真，结果显示规划后微小卫星的低可观测性能明显改善，满足飞行姿态规划的需求。

微小卫星；姿态规划；低可观测性；混沌算法；粒子群优化算法

0 引言

随着微小卫星在侦察、导航、数据中继及早期预警等领域得到越来越广泛的应用，微小卫星的安全性能也越来越得到重视。空间探测设备，反卫星武器的快速发展越来越多地限制了微小卫星的使用[1⁃3]。为提高在轨微小卫星的使用效能及生存能力，对其进行低可观测性设计是一个有效的方案。目前，对中、低轨卫星的探测主要由雷达完成，因此低可观测性设计中低雷达探测特性是关键。除了对微小卫星的外形进行低可观测外形设计外，对其进行适当的飞行姿态规划以提高其低可观测特性也很有必要[4⁃5]。

对微小卫星的飞行姿态规划可从静环境和动环境两个方面进行规划。静环境主要是指环境中的威胁分布和各威胁设备的威胁性均已知，可以精确地进行规划，有效地提高卫星的低可观测特性；动环境规划则是规划前并不知道威胁分布，也不知道有哪些威胁设备及威胁设备的威胁性，只有当卫星处于当前位置时才获得该位置所对应的威胁设备的分布情况，所以需要对卫星在飞行过程中进行实时调整来降低卫星的低可观测性能。本文主要对微小卫星动环境规划进行研究。

微小卫星姿态规划与无人机等的航迹规划的基本思想都是在一定的约束条件下寻找一种飞行方式，使无人机等所受到威胁最低，所以对微小卫星的姿态规划可以参考航迹规划方面的算法。目前国内外对路径规划算法的研究很多，常见的有Dijkstra算法、A*算法及相关的改进A*算法[6]、概率地图方法（Probabilistic Roadmap Method，PRM）[7]、Voronoi图法[8]、快速扩展随机树算法（Rapidly Random Exploring Trees，RRT）[9]等以及智能优化算法如粒子群算法[10]，遗传算法[11]，混沌算法[12]等，前几种方法都属于图搜索算法，算法的实现需要知道威胁分布情况，对于动环境规划并不适用，所以本文采用智能优化算法对微小卫星低可观测飞行姿态进行动环境规划。

1 混沌优化算法

混沌优化算法（Chaos Optimization Algorithm）的基本思想是把混沌变量从混沌空间映射到解空间，然后利用混沌变量具有遍历性、随机性和规律性的特点进行搜索。混沌优化算法具有对初值敏感、易跳出局部极小、搜索速度快、计算精度高、全局渐近收敛的特点。

混沌算法首先应使优化变量在空间中处于混沌状态，一般选用式（1）所示的Logisic映射来产生混沌变量，其中μ是控制参量。设0≤xn≤1,n=0,1,2,…，μ=4时，Logisic映射完全处于混沌状态[12]。混沌优化算法需要在混沌区间任意设定i个初值（不能为方程（1）的不动点0.25，0.5，0.75）作为i个初始混沌变量。

设一类连续对象的优化问题为求函数最小值，如式（2）所示：

采用混沌优化算法的基本步骤如下：

（1）算法初始化：给式（1）中的xn赋予i个大小在[0，1]之间的初值，可以得到i个变量xi.n+1，将这i个变量作为初始混沌变量；

（2）首先初始化相应的性能指标，随机设定一个最优解xopt，fopt作为初始最优解，因为优化的是最小值，所以先初始化一个较大的fopt，保证fopt在接下来的迭代搜索中能取到当前优化问题的最优解；

（3）用混沌变量进行迭代搜索；

Step1：因为优化问题的变量区间为[a,b]，所以需要将i个混沌变量xi.n+1，按式（3）分别转化为i个混沌变量x′i.n+1，x′i.n+1在[a,b]之间：

Step2：用混沌变量x′i.n+1进行迭代搜索，k是迭代次数，fi(k)表示变量为x′i.n+1时的函数值。若fi(k)≤fopt，则fopt=fi(k)，xopt=x′i.n+1，否则不执行任何操作；

Step3：k=k+1；

Step4：将x′i.n+1通过式（3）的逆运算转换回混沌变量区间[0，1]，即将x′i.n+1转换回xi.n+1；

Step5：将xi.n+1代入式（1）中的xn，继续进行混沌映射，得到新的混沌变量；

Step6：重复步骤（3），当fopt保持不变时，迭代搜索完成，输出xopt，fopt。

图1为初始值分别为0.501和0.502两个点的混沌演化轨道，初始距离仅为0.001，采用Logistic映射迭代后，两个点逐渐分离。可以看出，前6次迭代两个点的距离还比较近，在图中表现为两个点基本重合，当迭代次数超过7次后，两个点迅速分离，分别按照各自的混沌轨道运行。由此例可以看出混沌运动的初值敏感性和遍历性。

图1 0.501和0.502两点的混沌演化轨道

2 粒子群优化算法

粒子群优化算法（Particle Swarm Optimization，PSO）是基于群体智能的一种进化计算方法。PSO算法中每个粒子就是解空间中的一个解，它根据自己的经历和整个粒子群的经历来调整自己。每个粒子在飞行过程的最好位置就是该粒子本身所找到的最优解。整个粒子群在飞行过程中经历的最好位置就是整个种群目前所找到的最优解。粒子本身的最优解叫做个体最优解pbest，整个粒子群体找到的最优解叫做全局最优解gbest。每个粒子根据上述两个最优解结合更新公式不断更新自己，从而产生新一代粒子群体，粒子的“好坏”程度由适应度函数来评价。与一般的进化算法相比，PSO概念简单、容易实现并且需要调整的参数少，目前广泛应用于各种优化领域。

在PSO算法中，每个粒子可以看作是解空间中的一个点，假设粒子种群规模为N，则第i(i=1,2,…,N)个粒子的位置可表示为Xi。粒子的位置即是适应度函数的变量，根据粒子位置及适应度函数可以计算出粒子的适应度，根据适应度判断粒子的“好坏”程度。粒子所经历过的“最好”位置记为pbest(i)，该粒子的更新速度用Vi来表示，粒子群体所经历的“最好”的位置的用gbest表示，第i个粒子的速度和位置更新公式为：

式中：c1,c2为常数，称为学习因子；r1,r2是[0，1]上的随机数；ω称为惯性权重，同时粒子在更新自己的速度和位置的时候还受最大更新速度Vmax和最小更新速度Vmin的约束[13⁃15]，即Vi∈[Vmax,Vmin]。

粒子群算法的步骤如下：

（1）首先在搜索空间里初始化粒子种群。假设粒子种群规模为M，种群中的粒子记为Xi，粒子维数为n，表示为(ai1,ai2,…,ain)，即粒子在搜索空间中位置的坐标。同时初始化每个粒子的飞行速度V，也是一个n维向量(Vi1,Vi2,…,Vin)，其中i=1,2,…,M；

（2）根据每个粒子的位置和适应度函数，计算出每个粒子的适应度；

（3）比较适应度的大小，将每个粒子的当前适应度与该粒子的个体最优解相比较，若当前适应度优于个体最优解，则用当前位置取代个体最优解的位置；否则，个体最优解保持不变；

（4）首先比较得出最优粒子，然后将该粒子的适应度与种群最优解比较，若当前适应度优于种群最优解，则用当前位置取代个体最优解的位置；否则，种群最优解保持不变；

（5）根据更新公式（4），（5）来更新粒子的速度和位置；

（6）重复步骤（2）～（5），直到适应度函数达到最优或者满足迭代条件。

3 动环境规划算法设计

3.1 雷达分布及威胁水平计算模型

由于雷达是低可观测微小卫星最主要的威胁，所以在威胁水平评估建模中，主要考虑雷达作为威胁设备。当卫星处于某一位置时，并不是地球上所有的雷达都对卫星具有威胁性，而只是在某一角度内的雷达才对卫星具有威胁性，如图2所示。为了计算出卫星在某一位置时的威胁性水平，首先应分辨出哪些雷达对卫星具有威胁性[16]。

图2 威胁区、威胁设备及威胁方向定义

建立如图2所示的直角坐标系。取微小卫星位置S为坐标系原点，x轴指向为星体飞行方向，z轴指向为地球质心方向。图2中，Oe为地球质心，▲表示威胁雷达或激光设备，△表示无威胁雷达或激光设备。

自卫星位置S向威胁分布球面引切线，由切点可确定一个平面AOcB垂直于z轴，平面AOcB球面划分曲面ACB、曲面ADB两个部分。由雷达的工作原理可知仅分布在曲面ACB上的雷达能够探测到微小卫星，因此，可将该曲面定义为威胁区，威胁区内的设备定义为威胁设备，威胁设备相对于微小卫星的方向称为威胁方向。

威胁区俯仰角θz为切线SA与z轴的夹角，(θT,φT)为威胁方向的俯仰角及方位角，(θs,φs,γs)为微小卫星低可观测外形的姿态，θs，φs，γs分别表示微小卫星相应的俯仰角、方位角及横滚角[4]。

威胁区对应俯仰角θz可根据卫星飞行高度h及地球半径Re按公式（6）求得：

3.2 雷达威胁性定义

根据雷达方程及二元假设检验理论[4⁃5]，可得到雷达检测概率p随σ/R4单调递增。因此雷达的威胁性可由雷达监测概率来表示，本文将雷达威胁性T定义为

式中：R为雷达作用距离；σ为微小卫星的雷达散射截面积（Radar Cross⁃Section，RCS）；k为雷达威胁等级加权系数。当一个区域内有多部雷达时，只要被一部雷达探测到，则微小卫星就被发现，所以只要保证威胁性最大的雷达探测不到微小卫星即可。将威胁性最大的雷达对微小卫星的威胁性定义为卫星在该区域面临的威胁水平T，即：

式中：n为威胁雷达数；Ti（i=1，2，…,n）为第i部雷达的威胁性。

3.3 规划空间压缩

（1）对卫星产生威胁的雷达是威胁区内的雷达，所以微小卫星俯仰角θ的规划只需要在威胁区俯仰角θz内即可，即：

（2）根据姿控最大速度，前一时刻卫星的飞行姿态为(θ0,φ0,γ0)，则t时刻卫星外形姿态(θ,φ,ν)应满足：

式中：ωθ,ωφ,ωγ分别为卫星俯仰角θ，方位角φ，横滚角γ单位时间内的最大调整角度。

（3）目前，卫星低可观测外形多为锥形轴对称设计，所以对横滚角的调整对规划性能的影响很小，因此，可将4维规划空间{θ,φ,γ,t}简化为3维{θ,φ,t}，既可降低规划空间规模和规划计算复杂度，又可以降低规划算法设计难度[16]。

3.4 规划代价

卫星规划性能主要考虑两个方面：威胁水平和姿控能耗。提高微小卫星的低可观测性能首先要保证卫星拥有较低的威胁水平，其次，由于微小卫星的工作特点，对微小卫星的姿态的控制需要在有限的姿控能耗下完成。

假设威胁雷达的威胁等级相同，则可取威胁等级系数k=1。根据式（8）则可以得到，微小卫星在t时刻面临的威胁水平T(t)为：

在t时刻，定义卫星的姿控能耗代价C(t)为：

式中：aθ,aφ,aγ分别为θ,φ,γ的能耗加权系数；Δθ,Δφ, Δγ分别为前一时刻到t时刻θ,φ,γ的变化角度，且满足：

卫星的规划代价则为：

3.5 混沌粒子群（chaos PSO）规划算法

上述分析可知，t时刻雷达对微小卫星的威胁性可由T及C决定。对微小卫星飞行姿态规划就是寻找一组最优飞行姿态(θt,φt,γt)，使式（18）所确定的f(t)的值最小，因此，可以采取智能优化算法进行规划。在粒子群算法的基础上引入混沌运动，增加了粒子在规划空间的遍历性，从而提高了规划算法的精确性。

根据粒子群优化算法特征[13⁃15]，惯性权值ω可定义为：

式中：si为当前迭代次数；sI为最大迭代次数；W1，W2分别为惯性权值的初始值和终止值。粒子的适应度函数即为卫星的规划代价，粒子的位置即为卫星飞行姿态。

结合混沌优化算法和粒子群优化算法的基于局部近似最优解的混沌粒子群优化规划算法的步骤如下[17⁃18]：

（1）以微小卫星之前的飞行姿态为基础，根据式（20）～（25）初始化规模为N的初始粒子种群：

式中：θinitial，φinitial，γinitial分别为动环境规划中的初始规划姿态（俯仰角、方位角、横滚角），即微小卫星前一时刻的飞行姿态，rθ，rφ，rγ为[-1，1]内的随机数；tc为规划时长，pθ，pφ，pγ为初始化的粒子飞行姿态，i表示第i个粒子，vθ,φ,γ(i)表示第i个粒子的初始更新速度；fopt为个体最优解；gopt为种群最优解。

（2）计算出当前位置内的威胁雷达数目，根据威胁雷达分布计算出卫星规划代价。

（3）将当前的规划代价与个体最优解相比较，若小于个体最优解，则用当前姿态代替粒子最优姿态，保存粒子个体最优姿态，同样的，将当前种群中的最优解与先前保存的种群最优解相比较，若小于先前种群最优解，则用当前种群最优个体的信息取代先前的种群最优个体的信息，并保存种群最优个体。

（4）更新粒子信息速度更新公式：

姿态更新公式：

式中：c1，c2是常数，称为学习因子；r1，r2是[0，1]上的随机数，ω也为一个常数，称为惯性权值；粒子移动速度v∈(vmin,vmax)，vmin，vmax是常数，用来设定粒子移动速度。粒子移动速度与惯性权值都是用来维护全局最优解与局部最优解的平衡，本文引入动态权值ω，即减小了粒子群算法对粒子移动速度的依赖，所以本文中粒子移动速度的设定对粒子群算法性能影响较小。

若更新后的姿态不满足式（10）～（12）的范围，若θ＜θ0-ωθt，则θ=θ0-ωθt；若θ＞θ0+ωθt，则θ=θ0+ωθt。同理判断出φ，γ的取值。

（5）如果粒子有重叠现象，即粒子聚集在某个极值附近的情况，则保留其中一个粒子不变，其他粒子赋予混沌运动，首先将粒子逆运算到混沌区间即[0，1]，若粒子区间为[a,b]，根据逆运算公式x′=(x-a)/(b-a)转换，然后再采用Logistic映射进行迭代，再根据式（3）变换到优化变量空间中，其中需要判断当更新后的姿态不满足式（10）～（12）的范围时，若θ＜θ0-ωθt，则θ=θ0-ωθt；若θ＞θ0+ωθt，则θ=θ0+ωθt。同理判断出φ，γ的取值。

（6）重复算法步骤（2）～（5）直到满足迭代数或者最优解稳定，输出最优粒子的信息。

4 算法仿真结果及分析

本文的仿真条件如下：地球半径取6 371 km，卫星的飞行高度为500 km，俯仰角θ范围为0°～69°（69°为仿真条件下的威胁区俯仰角θz），方位角φ是0°～360°。威胁设备为地面雷达。根据微小卫星面临雷达威胁特征，规划VHF波段进行仿真，雷达工作频率取300 MHz。因为是动环境规划，所以微小卫星实时计算出其所在位置所面临的雷达威胁水平，粒子的群体规模为30，W1=0.9，W2=0.1，因为提高卫星低可观测性能更主要的是降低雷达对卫星的威胁性，能耗代价对卫星规划代价影响相对较小，所以姿控能耗系数取较小的值，将姿控能耗系数俯仰角θ取值为0.05，方位角φ取值为0.03，卫星在飞行过程的姿态变化一般较缓慢，所以取θ，φ的最大姿控速度均为0.2(°)/s，c1=0.2，c2=0.8，vmin=-10，vmax=10，vmin，vmax用来限制卫星在一个规划时长tc内卫星姿态变化角度范围在[-10°,10°]。

首先对卫星在处于某一位置时，分别采用粒子群算法和混沌粒子群算法进行规划仿真。当两种算法的迭代次数均为1～50次时，规划得到的最小威胁代价水平如图3所示。

图3 迭代次数为1～50次时chaos PSO和PSO的规划结果

当算法迭代次数均为50次时，将两种算法均重复运行50次，规划得到的最小威胁水平如图4所示。

从仿真结果可以看出，随着迭代次数的增加，混沌粒子群算法规划处的结果趋于稳定，而粒子群算法则有一定的波动性，因为粒子群算法本身容易陷入局部最优解，混沌运动的遍历性和随机性可以使其跳出局部最优解，增强了粒子群算法的搜索性能；同时迭代次数均为50次时，混沌粒子群算法也比粒子群算法更稳定，所以可以看出混沌粒子群算法性能更优。

图4 迭代次数均为50次时chaos PSO和PSO重复运行50次的规划结果

接下来对卫星绕地球运行一周进行规划，每个位置算法规划迭代次数均为50次，将两种算法进行比较的同时与无规划（θ，φ均为0°）情况进行比较。卫星绕球心每飞行1°进行飞行姿态规划，即tc=1°。

图5所示为雷达在300 MHz工作频率时对卫星进行飞行姿态规划后卫星在飞行过程中的俯仰角和方位角的变化，实线和点线分别表示采用chaos PSO规划后的卫星俯仰角θ，方位角φ的变化，点划线和虚线分别表示采用PSO规划后的卫星俯仰角θ，方位角φ的变化。由于设定的俯仰角的能耗系数比方位角的大，所以在图中俯仰角的变化比较小，主要通过调整方位角来减小卫星威胁代价。

图5 300 MHz频率采用chaos PSO和PSO规划时卫星飞行俯仰角θ，方位角φ变化

图6所示是进行飞行姿态规划后的卫星在飞行过程中的威胁代价图。实线（chaos PSO）采用chaos PSO算法规划后的威胁代价水平分布。点线（PSO）表示的是采用PSO算法进行规划后的威胁代价水平分布，虚线（noplan）表示的则是无规划时的卫星威胁代价水平分布。

表1300 MHz频率卫星无规划时，PSO规划时及chaos PSO规划时的仿真数据

由表1和图5，图6可以看出不管是PSO算法还是chaos PSO算法，都可以在很大程度上降低卫星的威胁代价，图中在某些地方规划后的威胁代价大于未规划时的威胁代价，是因为图中的威胁水平曲线是微小卫星在当前位置面临的威胁水平，而本文的规划代价为初始位置到当前位置的总体威胁水平，所以会出现在某些位置的代价值偏大的情况，但是总体上的代价值会远远小于未规划时的代价。由于规划是实时规划，因为最大调整角度的约束，所以卫星某一时刻某个位置的规划结果很依赖于前一时刻的卫星飞行姿态。如图中所示，有些地方PSO的规划结果会优于chaos PSO就是因为前一时刻的飞行姿态不同。但是从整体规划来看，chaos PSO的规划结果还是优于PSO，因为chaos PSO对最优解的搜索能力比PSO更强更稳定。

图6 300 MHz频率采用chaos PSO和PSO规划时微小卫星的威胁代价水平分布

5 结论

本文以粒子群算法和混沌算法为基础，对微小卫星进行了动环境的飞行姿态规划的算法设计。采用粒子群算法和混沌粒子群算法分别对微小卫星在某一位置和飞行一个周期两种情况进行规划，规划结果显示混沌粒子群算法比粒子群算法的规划性能更好，同时通过与无规划时微小卫星所面临的威胁水平相比较，可以看出两种算法均能有效提高微小卫星的低可观测性能。

微小卫星飞行姿态规划是个复杂的多学科课题，对卫星进行姿控规划的同时也对能源系统设计提出新的要求，同时进行姿态规划可能会对卫星功能及寿命产生一定影响。这些问题需要通过多学科的交流合作进行进一步研究。

[1]黄汉文.卫星隐身概念研究[J].航天电子对抗，2010，26（6）：22⁃24.

[2]FULGHUM D A，WALL R.Russia′s top designers claim anti⁃stealth skills[J].Avitio Week＆Space Technology，2001，8：82⁃85.

[3]张建民，孙健.国外隐身技术的应用与发展分析[J].舰船电子工程，2012，32（4）：18⁃21.

[4]苏抗，周建江.有限姿控能力的低可观测微小卫星姿态实时规划[J].航空学报，2010，31（9）：1841⁃1846.

[5]苏抗.微小卫星低可观测关键技术研究[D].南京：南京航空航天大学，2011.

[6]史辉，曹闻，朱述龙，等.A*算法的改进及其在路径规划中的应用[J].测绘与空间地理信息，2009，32（6）：208⁃211.

[7]刘贤敏，王琪，许书诚.改进型概率地图航迹规划方法[J].火力与指挥控制，2012，37（4）：121⁃124.

[8]吴海彬，林宜.基于改进Voronoi图的移动机器人在线路径规划[J].中国工程机械学报，2007，5（1）：117⁃121.

[9]崔挺，李俨，张明庄.基于改进RRT算法的无人机航迹规划[J].电子设计工程，2013，21（12）：50⁃53.

[10]王新增，慈林林，李俊山，等.基于改进粒子群优化算法的无人机实时航迹规划[J].微电子学与计算机，2011，28（4）：87⁃90.

[11]鲁艺，吕跃，罗燕，等.基于改进遗传算法的UAV航迹规划[J].电光与控制，2012，19（1）：29⁃33.

[12]李兵，蒋慰孙.混沌优化方法及其应用[J].控制理论与应用，1997，14（4）：613⁃615.

[13]MASEHIAN E，SEDIGHIZADEH D.A multi⁃objective PSO⁃based algorithm for robot path planning[C]//2010 IEEE Inter⁃national Conference on Industrial Technology（ICIT）.[S.l.]：IEEE，2010：465⁃470.

[14]吴晓军，李峰，马悦，等.均匀搜索粒子群算法的收敛性分析[J].电子学报，2012，40（6）：1115⁃1120.

[15]NICKABADI A，EBADZADEH M M，SAFABAKHSH R.A novel particle swarm optimization algorithm with adaptive iner⁃tia weight[J].Applied Soft Computing，2011，11（4）：3658⁃3670.

[16]苏抗，周建江.微小卫星低可观测外形飞行姿态规划[J].航空学报，2011，32（4）：720⁃728.

[17]杨俊杰，周建中，喻菁，等.基于混沌搜索的粒子群优化算法[J].计算机工程与应用，2005，41（16）：69⁃71.

[18]陈如清，俞金寿.混沌粒子群混合优化算法的研究与应用[J].系统仿真学报，2008，20（3）：685⁃688.

Research of dynamic⁃environment planning algorithm about flying attitude of micro⁃satellite

ZHOU Meng⁃xiao，CAI Qian，ZHOU Jian⁃jiang，JIANG Hui，HU Rong
（Nanjing University of Aeronautics and Astronautics，Nanjing 210016，China）

In order to improve the efficiency and viability of micro⁃satellite，especially the low observability of the micro⁃satellite，a dynamic⁃environment attitude planning algorithm based on the PSO algorithm and chaos algorithm was designed.The RCS value of the radar which threaten the micro⁃satellite can be reduced，and the low observability of the micro⁃satellite can be im⁃proved by real⁃time adjustment of the pitching angle and the azimuth angle of micro⁃satellite.The PSO algorithm can reduce the cal⁃culation complexity and improve the efficiency of real⁃time planning，and when the chaos movement is added into it，the accuracy of the algorithm can be improved.The simulation result of work frequency threatened in VHF wave band from radars shows that the low observability of the micro⁃satellite is obviously improved and the algorithm can meet the demands of flying attitude planning.

micro⁃satellite；attitude planning；low observability；chaos algorithm；PSO algorithm

TN927⁃34

A

1004⁃373X（2015）07⁃0001⁃06

周梦晓（1989—），男，江苏南通人，硕士研究生。主要研究方向为计算机测控、微小卫星低可观测飞行姿态规划等。

2014⁃10⁃14

国家“863”计划（4015⁃187001）

蔡倩，女，江苏常州人，副教授，硕士生导师。研究方向为计算机测控技术、测试信号处理、嵌入式系统及应用。

周建江，男，江苏南通人，教授，博士生导师。研究方向为雷达目标特性分析、特征控制与目标识别、射频隐身、机载电子信息系统。