0引言在现代战争中，大量的空中机动目标出现在战场，并且以其高机动性、高灵活性、全天候作战和有效的攻击火力，承担着越来越重要的作用，使战争由平面发展到立体空间，对战略战术和军队组成等产生了重大影响。
但是，随着地面雷达和低空导弹的日益完善以及反空火力的增强，空中机动目标在作战中的危险性不断增加。
由于雷达对低空目标的快速反应能力很差、波束宽等原因，基于雷达体制收稿日期：2017-05-01修回日期：2017-06-23基金项目：国家自然科学基金资助项目（61201391）作者简介：鲁鹏威（1992-），男，江苏泰州人，在读硕士。
研究方向：角度跟踪和轨迹预测。
*摘要：在现代战争中，空中目标的低空作战方式严重威胁了地面的作战单位。
针对轻型防空武器系统，提出了一种基于角度的跟踪算法，对低空机动目标的飞行特性进行分析，采用“当前”统计模型，利用扩展卡尔曼滤波理论，同时通过最小均方算法对机动频率进行自适应调整，实现了对目标的角跟踪和预测。
仿真结果表明，与固定机动频率的角跟踪算法相比，该角跟踪算法具有更好的跟踪性能，误差明显减小，并且该算法具有计算速度快、稳定性好等特点。
关键词：角跟踪，“当前”统计模型，自适应滤波，最小均方算法中图分类号：TN953+.5；TJ0文献标识码：ADOI ：10.3969/j.issn.1002-0640.2018.06.006引用格式：鲁鹏威，贾方秀，王晓鸣，等.自适应低空机动目标角跟踪算法研究［J ］.火力与指挥控制，2018，43（6）：32-36.自适应低空机动目标角跟踪算法*鲁鹏威1，贾方秀1，王晓鸣1，刘铭2（1.南京理工大学智能弹药技术国防重点学科实验室，南京210094；2.西安交通大学人居环境与建筑工程学院，西安710049）文章编号：1002-0640（2018）06-0032-05Vol.43，No.6Jun ，2018火力与指挥控制Fire Control &Command Control 第43卷第6期2018年6月Adaptive Angle Tracking Algorithm for Low Altitude Maneuvering TargetLU Peng-wei 1，JIA Fang-xiu 1，WANG Xiao-ming 1，LIU Ming 2（1.Ministerial Key Laboratory of Intelligent Ammunition ，Nanjing University of Science and Technology ，Nanjing 210094，China ；2.Human Settlements and Architectural Engineering ，Xi'an Jiao Tong University ，Xi ’an 710049，China ）Abstract ：In the modern war ，the low-level mode of operations of the maneuvering target seriouslythreatens the combat units on the ground.For the light air defense weapon system ，this paper proposes a tracking algorithm based on the angle.The flight characteristics of the maneuvering target are analyzed.The "current"statistical model and the extended Calman filter are used to realize the angle tracking and prediction ，and adjust the maneuvering frequency adaptively by the least mean square algorithm.The simulation results show that compared with the fixed maneuvering frequency angle tracking algorithm ，the angle tracking algorithm has better tracking performance and less error ，and thealgorithm has the advantages of fast computation speed and good stability.Key words ：angle tracking ，“current ”statistical model ，adaptive filtering ，least mean square Citation format ：LU P W ，JIA F X ，WANG X M ，et al.Adaptive angle tracking algorithm for low altitude maneuvering target ［J ］.Fire Control &Command Control ，2018，43（6）：32-36.32··（总第43-）的防空系统将无法完成对低空、超低空目标的探测与跟踪［1］。
依靠地形地物的遮蔽和地面杂波的干扰，在低空飞行的目标不易被发现，所以低空和超低空突防成为空中单位的主要作战方式。
目前，在轻型防空武器系统超低空防空领域，仅有角测量的条件下，国内外学者提出了大量角跟踪算法。
王勇等人［2］提出了“当前”统计模型下的扩展卡尔曼角跟踪算法，实现了目标的角跟踪，但是跟踪精度不高；Al-Emadi等人［3］提出了一种利用同步信号参考与振荡信号的角度跟踪技术的测量转换器，表现出优异的跟踪性能，但是所利用的系统较为复杂；任波等人［4］提出了非线性预测滤波算法，在选取恰当的加权矩阵条件下，可以实现有效的跟踪，但是，加权矩阵很难解析求到，主要根据经验给定，不能保证最优结果。
本文针对轻型防低空武器系统，以低空机动目标为研究对象，分析其飞行轨迹，提出了一种基于角度的自适应跟踪算法，通过仿真证明了该算法的有效性和优越性，并且可以在简单的系统中实现。
1算法设计空中作战单位在低空或者超低空飞行时，一般采用地形跟随飞行方式，利用地物作掩护，在贴近地面的高度上隐蔽接近攻击目标［5］。
低空机动目标的飞行特性主要分为机动飞行和非机动飞行。
非机动飞行主要表现为巡航飞行；机动飞行包括以下的飞行方式：悬停跳跃、平飞加速、定常爬升、定常转弯、减速转弯、蛇形飞行和高速拉起等。
本文设计的低空机动目标自适应角跟踪算法的基本思路如图1所示。
首先根据低空机动目标的飞行特性，建立合适的模型；然后根据建立的模型对滤波进行初始化设计；在初始化设计后，进行扩展卡尔曼滤波计算；滤波结束后进行模型的自适应修正，以提高接下来跟踪的精度；最后是滤波结果的输出。
图1算法流程图1.1建立目标模型模型的设计对估计精度有很大的影响，在该算法中，将低空机动目标看作质点，观测用的传感器在坐标原点，具体的几何关系如图2所示。
根据图2，建立如下的随机非线性系统：（1）式中，f和h都是非线性方程；w（k）是以Q（k）为协方差矩阵的系统过程噪声，v（k）是量测噪声，其协方差矩阵是R（k）；x（k）是系统的状态向量，z（k）是系统的观测向量，具体关系如下：（2）（3）综合之前的低空机动目标的飞行特性，建立“当前”统计模型，连续系统模型的描述如下［6］：（4）式中，。
将“当前”统计模型离散化，得到离散形式的动力学模型如下：（5）式中：图2目标与观测量的几何关系图鲁鹏威，等：自适应低空机动目标角跟踪算法33··0963（总第43-）火力与指挥控制2018年第6期过程噪声方差Q 的表示如下：（6）式中，，，，，，。
此外，α是“当前”统计模型中的机动频率，σ2x 的计算如下（σ2y 和σ2z 的计算与之类似）：（7）式中，a M 与a -M 是目标在某一方向所能达到的最大正加速度和最大负加速度，a 是目前加速度的预测值。
1.2扩展卡尔曼滤波1.2.1计算预测均值和方差根据前文建立的模型，由于系统的非线性，采用扩展卡尔曼滤波进行状态估计［7］。
首先计算飞行目标的预测均值和方差，计算公式如下：（8）（9）1.2.2计算雅克比矩阵由于观测方程h 是非线性的，所以需要计算h 在处的雅克比矩阵，计算如下：（10）式中，x 、y 、z是中预测的目标位置信息。
1.2.3计算观测预测、协方差矩阵和滤波增益观测预测的计算如下：（11）新息协方差矩阵计算如下：（12）卡尔曼滤波增益计算如下：（13）1.2.4计算后验均值和协方差矩阵后验均值计算如下：（14）协方差矩阵计算如下：（15）1.3机动频率自适应在“当前”统计模型中，机动频率α是根据跟踪目标机动程度来确定的，机动频率越小，目标的机动程度越大。
低空机动目标可以进行机动性小的巡航飞行，也能够实现机动性强的飞行，所以，在该模型中，对机动频率进行自适应调整将会提高跟踪精度。
最小均方（LMS ）算法具有计算简单、易于收敛等特性，可以很好地应用在实时性要求高的跟踪算法中，所以利用LMS 算法对机动频率进行自适应调整［8］。
在LMS 算法中，它的更新方程如下：（16）式中，μ是收敛因子，e （k ）是预期输出与实际输出之差，x （k ）是系统的输入，为了保证算法的收敛性，μ34··0964（总第43-）应该在一定的范围内取值。
将LMS 算法应用到本文的扩展卡尔曼滤波中，将整个滤波系统看作一个网络，机动频率α就是这个网络的权值，e （k ）x （k ）是观测预测与实际观测之差，即，所以，可以得到机动频率α的调整公式：（17）2仿真与分析考虑低空机动目标的飞行特性，利用计算机生成飞行轨迹，仿真的轨迹特征见表1，目标的飞行轨迹见图3，图中每个点之间的时间为0.1s 。
然后计算出目标的高低角和方位角，作为理论数据，将观测误差以正态分布随机加到理论数据上，作为观测数据。
本次仿真的采样时间为0.1s ，观测方位角和高低角的误差方差为1°，方位角理论值和观测值见图4，高低角理论值和观测值见图5。
表1仿真轨迹特征图3飞行轨迹图图4方位角理论值和观测值图5高低角理论值和观测值图6方位角误差图图7高低角误差图在进行仿真时，对固定的机动频率和自适应机动频率分别进行仿真，固定的机动频率取值为0.001和0.01，自适应机动频率的初始取值与其相同，方位角的误差见图6，高低角的误差见图7，计算整个过程的平均均方根误差见下页表2。
