单模型机动目标跟踪算法的仿真研究

《单模型机动目标跟踪算法的仿真研究》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、单模型机动目标跟踪算法的仿真研究摘耍：分析研究了用CV模型，CA模型，Singer模型对单机动Fl标进行跟踪算法的matlab仿真。比较了各模型的滤波性能，得出了一些有意义的结论。关键词：1=1标跟踪；卡尔曼滤波；CV模型；CA模型；Singer模型；一.问题描述本文研究的实例是一个二维平面的雷达这个问题属于单目标跟踪问题，一般来说，如果目标做匀速直线运动时，跟踪问题十分容易；但当目标做机动时，由于无法准确预知目标下一时刻的运动状态，使得跟踪变得很怵1难。这就需要发展合适的目标运动模型，现在的各种模型大致分为单模型和多模型方法，出于多模型较为复杂，这里我们仅对单模型方法进行讨论。常用的单模型有

2、匀速模型（CV）、匀加速模型（CA）、Signer模型。现实世界中的大部分运动冃标都存在各种机动，冃标做匀速直线飞行的概率很小，采用CV模型一般是不可取的，只有当H标做匀速直线飞行或者近似匀速直线飞行时才能取得很好的效果。机动强度不大时，可以采用CA模型或者Singer模型雷达对口标的量测并不真实准确，而是存在一定的随机噪声干扰，一般假设噪声符合高斯分布。由于量测数据大多含有噪声和杂波，为了提高冃标状态（位置、速度等）估计精度，通常要对量测数据进行预处理以提高数据的准确度和精度。情景想定假定有一二座标雷达对一平面上运动的冃标进行观测，冃标在0-400秒沿轴作恒速直线运动，运动速度为15米/秒，

3、冃标的起始点为（-10000米,2000米），在=ty=t400-600秒向x轴方向做的慢转弯，加速度为ax=-0.075米/秒，ay=0.075米/秒，完成慢转弯。雷达扫描周期02=T2秒，和独立地进行观测，观测噪声的标准差均为100米。试建立雷达对目标的跟踪算法，并进行仿真分析，给出仿真分析结果。二.模型建立1・考虑随机干扰情况。当冃标无机动，即冃标作匀速或匀加速直线运动时,可分别采用常速CV模型或三阶常加速CA模型。其CV模型可以表示为：旳）範）叩）其CA模型可以表示为:01o'0x(0二001+0W)x(00001式中x,x,x分別为运动目标的位置、速度和加速度分量，w⑴是均值为零，方

4、差为/的高斯白噪声。从上式可以看出，CV和CA模型都是线性模型，这给目标跟踪算法的事先带来了方便，简化了计算。。当H标处于机动状态即H标的加速度向量发生变化时,采用以上模型会引起较大的误差，这时需要全面考虑H标的机动状态采用其他模型，如下面介绍的忖标模型。时间相关模型(Singer模型)Singer首次假设机动加速度a(t)服从一阶时间相关过程，其时间相关函数Ra(t)为指数衰减形式,即:R(T)=E{a[t)a(t+T)}=crya^(a>0)式中，为目标加速度方差：a为机动频率。假定机动加速度均值为零，机动加速度的概率密度函数近似服从均匀分布,方差曲近似服从均匀分布概率密度计算得來，即为：

5、c侶牛[1+4也7]式屮，Amax为最大机动加速度；Pmax为其发生概率；P。为非机动发生概率。对时间相关函数Ra(t)应用Wiener一Kolmogorov口化方法后，机动加速度a(t)可用输入为白噪声的一阶时间相关模型来表示，即：d(r)=)+w(/)则这时机动U标模型可表示为:x(t)x(0■0°1迩)+»0一a掩)1■■w(/)式屮，W(t)是均值为零，方差为2d，的白噪声。一.具体实现1.CV(恒速)模型取状态变量为：XX=yy状态方程为:观测方程为:其屮：X(k+l)=QX(k)Z(k)=HX(k)+V(k)z=1000H=0100Q=[1，TO0;0丄(W)O1,T;0,001]

6、对目标位置和速度的最佳滤波和最佳预测如下:预测：X(k/k-l)=QX(k-l/k-l)预测误差协方差：P(k/k-l)=QP(k-l/k-l)QT卡尔曼增益：K(k)=P(k/k-l)HT[HP(k/k-l)HT4-R]_1滤波:AAAX(k/k)=X(k/k-l)+K(k)[Z(k)-HX(k/k-1)]滤波协方差：其屮：P(k/k)=[I-K(k)H]P(k/k-l)R=滤波的初始化在实际小，我们通常无法得知H标的初始状态，这时我们可以利用前几个观测值建立状态的起始估计。由于只考虑H标位置和速度，这里利用前两个观测值建立起始估计，即两点起始法：AX(2/2)=[zx(l);vx;zy(l

7、);vy]P(2/2)=varvai7T,2varA2/(TA2),0,0;0,0,varA2,varA2/T;0,0,varA2/T,2varA2/(TA2)<7；=(tJ=10000仿真分析采用MATLAB编写仿真程序，利用蒙特卡罗方法对跟踪滤波器进行仿真分析，次数为50次。以下给出仿真图和结果分析。(1)・假设为匀速直线运动(0-400s)图一图一显示了在匀速直线运动下CV模型的滤波效果。图