专利名称：奇异值分解的多目标跟踪方法
技术领域：
本发明涉及一种雷达多目标跟踪方法，特别涉及一种奇异值分解的多目标跟踪方法，属于信息技术领域。
背景技术：
多目标跟踪技术在军用及民用领域均有广泛的应用，可用于空中目标检测、跟踪与攻击，空中导弹防御，空中交通管制，港口和海洋监视等。近年来，随着战场环境的改变，对抗和反对抗技术的发展，产生了背景强杂波、低信噪比、低检测概率和高虚警率等一系列问题，对多目标跟踪方法的精度和准确性提出了更高的要求。
多目标跟踪的目的是将探测器所接收到的量测对应不同的信息源，形成不同观测集合或轨迹，根据轨迹估计被跟踪目标的数目以及每一目标的运动参数，实现对多个目标的跟踪。用于多目标状态估计的基本滤波方法有α-β滤波、α-β-γ滤波、卡尔曼滤波、扩展卡尔曼滤波、高斯和近似、最优非线性滤波、粒子滤波和自适应滤波等。α-β和α-β-Y滤波器由于结构简单，计算量小，在早期计算机资源短缺时应用很广。卡尔曼滤波是多目标跟踪的一种基本方法，但是需要知道系统的精确数学模型，并且只适用于线性系统，限制了算法的应用。扩展卡尔曼滤波将卡尔曼滤波理论扩展到非线性领域，用一个高斯分布来近似状态的条件概率分布；而当近似条件不满足时，高斯和滤波器则用一个高斯分布的加权和来近似状态的条件概率分布。最优非线性滤波使用Makov转移概率来描述目标的动力学过程，具有很好的特性，但是计算量较大，因此一直没有得到广泛应用。粒子滤波采用随机采样，由于计算量太大和粒子退化问题，不适合实际应用。为了改进粒子滤波，无迹卡尔曼滤波采用确定性采样，使得采样的粒子点个数减少，避免了粒子滤波中的粒子点退化问题，因此其应用领域很广。自适应滤波方法通过对目标机动的检测，实时调整滤波器参数或增加滤波器的状态，使滤波器实时适应目标运动，特别适合对机动目标的跟踪；目前，在实际雷达跟踪系统最常用的仍然为JPDA(Joint Probabilistic Data Association,联合概率数据关联)方法(James A. Roecker,A Class of Near Optimal JPDA Algorithms,IEEE TRANSACTIONS 0NAER0SPACE AND ELECTRONIC SYSTEMS, 1994, VOL. 30(2) :504-510),其它方法大多数是对JPDA方法的简化等。然而，JPDA等方法误差估计的方差阵中有两个半正定矩阵相减，在有限字长的处理系统中会产生含有正负特征值的对称矩阵，导致雷达跟踪过程目标失跟和整个雷达系统性错误。
发明内容
为了解决现有目标跟踪方法数值结构不稳定导致雷达跟踪过程目标失跟的技术缺陷，本发明提供一种奇异值分解的多目标跟踪方法，该方法在多目标跟踪的时间和测量更新中，通过对估计误差方差阵的奇异值分解，建立数值稳定结构模型，没有误差估计的方差阵中有两个半正定矩阵相减，在有限字长的处理系统中能够保证不会产生含有负特征值的对称矩阵，可以避免雷达跟踪过程目标失跟和整个雷达系统性错误。
本发明解决其技术问题采用的技术方案是，ー种奇异值分解的多目标跟踪方法，其特征包括以下步骤I、定义N个目标跟踪中第i个目标的离散化模型为Xi (k+Ι) = Φ (k+1, k) Xi (k) + A Qi (k),其中x#) = [x,もも少,夕,见し为状态向量，(X，y, z)为目标在地面參考直角坐标系下的位置坐标，OiQO表示方差为Qi (k)的过程噪声向量，c>(k+l, k) = Φ =
A1 O O _
diagEOp Φ” O1]为状态转移矩阵，A =ダΦ(ん+ 1,ポ(りみ=O A1 O ，Γ⑴为 L O O A1
厂 I JI T -T22
系数矩阵，PrI 0 0I Γ1= [O O UtjO1= O I T ,T
Γ = O Γ1 O，O O I Λ Γ I r3 I r2 r]T
_0 or」L 」2Γ T_，
为采样周期；第i个目标的时间更新为Xi (k/k-Ι) = Φχ (k-l/k-1)通过调用奇异值分解算法得到Vi(k/k-l)Oi(k/k-l)V^(k/k-l) = Y(k/k-l)YT(k/k-l)其中Xi(k/k-l)为第i个目标在kT时亥Ij的ー步预测值，V(k/k-l)D(k/k-l)VT(k/k-l)为对应的一步预测误差的方差阵，
Y(k/k-l) = ^Vt(k-l/k-l)O}(k-l/k-l) AQf (ん-I)]，Vjk/k-l)为正交矩阵，Djk/
k-Ι)为对角矩阵；初始条件为 Xi (0/0)和Vi(OZO)Di(OZO)V;T(0/0) = P;(0/0);2、第 i 个目标观测方程为Zi (k) = gi [Xi (k) ] +Vi (k)其中Zi (k)为对第i个目标的r维观测向量，Lxi (k)]为对应的输出，Vi (k)表示方差为RiGO测量噪声，通过调用奇异值分解算法得到% (k / k)Ot (k / k)\J (klk) = Y(k / ん) τ (k / k)Gi (k) = % (k / k)O1 (k / k)\J (k / ん)(k)x, (klk) = Xi {klk-V) + Gi (ん){文 Alj (ん)zy (k) - gt [x;(ん / ん-1)]}
；=1其中 (ん/ ん)=[v^ (klk- 1)D^ (Jdk-X) (ん)]，R (ん / ん)为正交矩阵，を.(ん / ん)
为对角矩阵；zu(k)为雷达对第i个目标的第j(j = 1,2,…，m)个回波，Xi(k/k)为第i个
目标kT时刻的滤波值，AijGO为权系数，且Σん⑷= i，H#) = gg；[x"(f]；
^ 'ムdx^參,(綱)3、第i个跟踪估计方法为通过调用奇异值分解算法得到 \ (ん / ん)Di (k / k)\J (klk) = A(ん)At (ん)
其中A(ん)=「γτ(k/k-1)D"(k/k-l) G1 (ん)dT(I-Quut) ，
Λ，ι(り O …OI(k)Ω= ° も:(ん)；'° ,u= I，d= Αψ};
OO …々ぶ)」し」[Alm(k)_Δしj (k)为第j个候选回波信息向量，Δ Jj j (k) = Zi, j (k) -gi [Xi (k/k-l)]。本发明的有益结果是通过对估计误差方差阵进行三次奇异值分解，建立了数值稳定的多目标跟踪结构模型，避免了误差估计的方差阵中两个半正定矩阵相减，在有限字长的处理系统中不会出现数值发散，从而保证了多目标跟踪方法的可靠性，避免了雷达跟踪过程目标失跟和整个雷达系统性错误。
·
下面结合实例对本发明作详细说明。
具体实施例方式I、定义N个目标跟踪中第i个目标的离散化模型为Xi (k+Ι) = Φ (k+1, k) Xi (k) + A Qi (k),其中x#) = [x,もも少,夕,见し为状态向量，(X，y, z)为目标在地面參考直角坐标系下的位置坐标，OiQO表示方差为Qi (k)的过程噪声向量，c>(k+l, k) = Φ =
A1 O O _
diagEOp Φ” O1]为状态转移矩阵，Λ =ΓΦ(ん+ 1,ポ(りみ=O A1 O ，Γ⑴为
L O O A1
厂 I J
IT -T2 2
系数矩阵，「ハ 0 0I F1 = [O O UtjO1= O I T ，
Γ = O Γ1 O ，O O I Λ Γ I r3 I r2 r]T
_0 or」L 」2Γ T_，
T为采样周期；第i个目标的时间更新为Xi (k/k-l) = Φχ (k-l/k-1)通过调用奇异值分解算法得到Vi (klk- 1)D! (klk-1) V;T (k I k-l) = Y(k I k-1) Yt (klk-I)其中Xi (k/k-l)为第i个目标在kT时亥Ij的ー步预测值，V (k/k-l) D (k/k-l) Vt (k/k-l)为对应的一步预测误差的方差阵，
Yiklk-^^Xik-llk-l^ik-llk-l) AQf (ん-I)]，Vi(Vk-I)为正交矩阵，Djk/
k-Ι)为对角矩阵；初始条件为 Xi (0/0)和Vi(OZO)Di(OZO)V;T(0/0) = P;(0/0);2、第 i 个目标观测方程为=Zi (k) = gi[Xi (k)]+Vi (k)其中Zi(k)为对第i个目标的r维观测向量，gi[Xi (k)]为对应的输出，ViGO表示方差为Ri (k)测量噪声，通过调用奇异值分解算法得到
权利要求
1.一种奇异值分解的多目标跟踪方法，其特征在于包括以下步骤 (I)、定义N个目标跟踪中第i个目标的离散化模型为
全文摘要
本发明公开了一种奇异值分解的多目标跟踪方法，用于解决现有的目标跟踪方法数值结构不稳定导致雷达跟踪过程目标失跟的技术问题。技术方案是对估计误差方差阵进行奇异值分解，建立数值稳定结构模型，没有误差估计的方差阵中有两个半正定矩阵相减，在有限字长的处理系统中能够保证不会产生含有负特征值的对称矩阵。通过对估计误差方差阵的奇异值分解，建立了数值稳定的多目标跟踪结构模型，避免了误差估计的方差阵中两个半正定矩阵相减，在有限字长的处理系统中不会出现数值发散，从而保证了目标跟踪系统的可靠性、避免了雷达跟踪过程目标失跟和整个雷达系统性错误。
文档编号G01S13/66GK102707278SQ20121004494
公开日2012年10月3日 申请日期2012年2月27日 优先权日2012年2月27日
发明者史忠科 申请人:西北工业大学