CN102707278B - 奇异值分解的多目标跟踪方法 - Google Patents
奇异值分解的多目标跟踪方法 Download PDFInfo
- Publication number
- CN102707278B CN102707278B CN201210044944.XA CN201210044944A CN102707278B CN 102707278 B CN102707278 B CN 102707278B CN 201210044944 A CN201210044944 A CN 201210044944A CN 102707278 B CN102707278 B CN 102707278B
- Authority
- CN
- China
- Prior art keywords
- target
- target tracking
- singular value
- value decomposition
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 238000005070 sampling Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract 1
- 230000000087 stabilizing effect Effects 0.000 abstract 1
- 238000001914 filtration Methods 0.000 description 10
- 239000002245 particle Substances 0.000 description 7
- 238000001514 detection method Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000008485 antagonism Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种奇异值分解的多目标跟踪方法,用于解决现有的目标跟踪方法数值结构不稳定导致雷达跟踪过程目标失跟的技术问题。技术方案是对估计误差方差阵进行奇异值分解,建立数值稳定结构模型,没有误差估计的方差阵中有两个半正定矩阵相减,在有限字长的处理***中能够保证不会产生含有负特征值的对称矩阵。通过对估计误差方差阵的奇异值分解,建立了数值稳定的多目标跟踪结构模型,避免了误差估计的方差阵中两个半正定矩阵相减,在有限字长的处理***中不会出现数值发散,从而保证了目标跟踪***的可靠性、避免了雷达跟踪过程目标失跟和整个雷达***性错误。
Description
技术领域
本发明涉及一种雷达多目标跟踪方法,特别涉及一种奇异值分解的多目标跟踪方法,属于信息技术领域。
背景技术
多目标跟踪技术在军用及民用领域均有广泛的应用,可用于空中目标检测、跟踪与攻击,空中导弹防御,空中交通管制,港口和海洋监视等。近年来,随着战场环境的改变,对抗和反对抗技术的发展,产生了背景强杂波、低信噪比、低检测概率和高虚警率等一系列问题,对多目标跟踪方法的精度和准确性提出了更高的要求。
多目标跟踪的目的是将探测器所接收到的量测对应不同的信息源,形成不同观测集合或轨迹,根据轨迹估计被跟踪目标的数目以及每一目标的运动参数,实现对多个目标的跟踪。用于多目标状态估计的基本滤波方法有α-β滤波、α-β-γ滤波、卡尔曼滤波、扩展卡尔曼滤波、高斯和近似、最优非线性滤波、粒子滤波和自适应滤波等。α-β和α-β-γ滤波器由于结构简单,计算量小,在早期计算机资源短缺时应用很广。卡尔曼滤波是多目标跟踪的一种基本方法,但是需要知道***的精确数学模型,并且只适用于线性***,限制了算法的应用。扩展卡尔曼滤波将卡尔曼滤波理论扩展到非线性领域,用一个高斯分布来近似状态的条件概率分布;而当近似条件不满足时,高斯和滤波器则用一个高斯分布的加权和来近似状态的条件概率分布。最优非线性滤波使用Makov转移概率来描述目标的动力学过程,具有很好的特性,但是计算量较大,因此一直没有得到广泛应用。粒子滤波采用随机采样,由于计算量太大和粒子退化问题,不适合实际应用。为了改进粒子滤波,无迹卡尔曼滤波采用确定性采样,使得采样的粒子点个数减少,避免了粒子滤波中的粒子点退化问题,因此其应用领域很广。自适应滤波方法通过对目标机动的检测,实时调整滤波器参数或增加滤波器的状态,使滤波器实时适应目标运动,特别适合对机动目标的跟踪;目前,在实际雷达跟踪***最常用的仍然为JPDA(Joint Probabilistic Data Association,联合概率数据关联)方法(JamesA.Roecker,A Class of Near Optimal JPDA Algorithms,IEEE TRANSACTIONS ONAEROSPACEAND ELECTRONIC SYSTEMS,1994,VOL.30(2):504-51O),其它方法大多数是对JPDA方法的简化等。然而,JPDA等方法误差估计的方差阵中有两个半正定矩阵相减,在有限字长的处理***中会产生含有正负特征值的对称矩阵,导致雷达跟踪过程目标失跟和整个雷达***性错误。
发明内容
为了解决现有目标跟踪方法数值结构不稳定导致雷达跟踪过程目标失跟的技术缺陷,本发明提供一种奇异值分解的多目标跟踪方法,该方法在多目标跟踪的时间和测量更新中,通过对估计误差方差阵的奇异值分解,建立数值稳定结构模型,没有误差估计的方差阵中有两个半正定矩阵相减,在有限字长的处理***中能够保证不会产生含有负特征值的对称矩阵,可以避免雷达跟踪过程目标失跟和整个雷达***性错误。
本发明解决其技术问题采用的技术方案是,一种奇异值分解的多目标跟踪方法,其特征包括以下步骤:
1、定义N个目标跟踪中第i个目标的离散化模型为
xi(k+1)=Φ(k+1,k)xi(k)+Λωi(k),
其中:为状态向量,(x,y,z)为目标在地面参考直角坐标系下的位置坐标,ωi(k)表示方差为Qi(k)的过程噪声向量,Φ(k+1,k)=Φ=diag[Φ1,Φ1,Φ1]为状态转移矩阵, Γ(t)为系数矩阵, Γ1=[0 0 1]T, T为采样周期;
第i个目标的时间更新为:
xi(k/k-1)=Φxi(k-1/k-1)
通过调用奇异值分解算法得到
其中:xi(k/k-1)为第i个目标在kT时刻的一步预测值,V(k/k-1)D(k/k-1)VT(k/k-1)为对应的一步预测误差的方差阵, Vi(k/k-1)为正交矩阵,Di(k/k-1)为对角矩阵;初始条件为xi(0/0)和
2、第i个目标观测方程为:zi(k)=gi[xi(k)]+vi(k)
其中:zi(k)为对第i个目标的r维观测向量,gi[xi(k)]为对应的输出,vi(k)表示方差为Ri(k)测量噪声,通过调用奇异值分解算法得到
3、第i个跟踪估计方法为:通过调用奇异值分解算法得到
其中:
Δi,j(k)为第j个候选回波信息向量,
Δi,j(k)=zi,j(k)-gi[xi(k/k-1)]。
本发明的有益结果是:通过对估计误差方差阵进行三次奇异值分解,建立了数值稳定的多目标跟踪结构模型,避免了误差估计的方差阵中两个半正定矩阵相减,在有限字长的处理***中不会出现数值发散,从而保证了多目标跟踪方法的可靠性,避免了雷达跟踪过程目标失跟和整个雷达***性错误。
下面结合实例对本发明作详细说明。
具体实施方式
1、定义N个目标跟踪中第i个目标的离散化模型为
xi(k+1)=Φ(k+1,k)xi(k)+Λωi(k),
其中:为状态向量,(x,y,z)为目标在地面参考直角坐标系下的位置坐标,ωi(k)表示方差为Qi(k)的过程噪声向量,Φ(k+1,k)=Φ=diag[Φ1,Φ1,Φ1]为状态转移矩阵, Γ(t)为系数矩阵, Γ1=[0 0 1]T, T为采样周期;
第i个目标的时间更新为:
xi(k/k-1)=Φxi(k-1/k-1)
通过调用奇异值分解算法得到
其中:xi(k/k-1)为第i个目标在kT时刻的一步预测值,V(k/k-1)D(k/k-1)VT(k/k-1)为对应的一步预测误差的方差阵, Vi(k/k-1)为正交矩阵,Di(k/k-1)为对角矩阵;初始条件为xi(0/0)和
2、第i个目标观测方程为:zi(k)=gi[xi(k)]+vi(k)
其中:zi(k)为对第i个目标的r维观测向量,gi[xi(k)]为对应的输出,vi(k)表示方差为Ri(k)测量噪声,通过调用奇异值分解算法得到
其中: 为正交矩阵,为对角矩阵;zij(k)为雷达对第i个目标的第j(j=1,2,…,m)个回波,xi(k/k)为第i个目标kT时刻的滤波值,λij(k)为权系数,且:例如取gi[xi(k)]=[ri(k) αi(k) βi(k)]T,ri为雷达能测量斜距、αi为高低角、βi方位角,且
3、第i个跟踪估计方法为:通过调用奇异值分解算法得到
其中:
Δi,j(k)为第j个候选回波信息向量,Δi,j(k)=zi,j(k)-gi[xi(k/k-1)]。
Claims (1)
1.一种奇异值分解的多目标跟踪方法,其特征在于包括以下步骤:
(1)、定义N个目标跟踪中第i个目标的离散化模型为
xi(k+1)=Φ(k+1,k)xi(k)+Λωi(k),
其中:为状态向量,(x,y,z)为目标在地面参考直角坐标系下的位置坐标,ωi(k)表示方差为Qi(k)的过程噪声向量,Φ(k+1,k)=Φ=diag[Φ1,Φ1,Φ1]为状态转移矩阵,Γ为系数矩阵, T为采样周期;
第i个目标的时间更新为:
xi(k/k-1)=Φxi(k-1/k-1)
通过调用奇异值分解算法得到
Vi(k/k-1)Di (k/k-1)Vi T(k/k-1)=Y(k/k-1)YT(k/k-1)
其中:xi(k/k-1)为第i个目标在kT时刻的一步预测值,V(k/k-1)D(k/k-1)VT(k/k-1)为对应的一步预测误差的方差阵,Vi(k/k-1)为正交矩阵,Di(k/k-1)为对角矩阵;初始条件为xi(0/0)和Vi(0/0)Di (0/0)Vi T(0/0)=Pi(0/0);
(2)、第i个目标观测方程为:zi(k)=gi[xi(k)]+vi(k)
其中:zi(k)为对第i个目标的r维观测向量,gi[xi (k)]为对应的输出,vi(k)表示方差为Ri(k)的测量噪声,通过调用奇异值分解算法得到
(3)、第i个目标跟踪估计方法为:通过调用奇异值分解算法得到
Vi(k/k)Di(k/k)Vi T(k/k)=A(k)AT(k)
△i,j(k)为第j个候选回波信息向量,
△i,j(k)=zi,j(k)-gi[xi(k/k-1)]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210044944.XA CN102707278B (zh) | 2012-02-27 | 2012-02-27 | 奇异值分解的多目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210044944.XA CN102707278B (zh) | 2012-02-27 | 2012-02-27 | 奇异值分解的多目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102707278A CN102707278A (zh) | 2012-10-03 |
CN102707278B true CN102707278B (zh) | 2014-02-05 |
Family
ID=46900208
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210044944.XA Expired - Fee Related CN102707278B (zh) | 2012-02-27 | 2012-02-27 | 奇异值分解的多目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102707278B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104793201B (zh) * | 2015-05-04 | 2018-04-24 | 哈尔滨工业大学 | 一种跟踪临近空间高超声速目标的修正变结构网格交互多模型滤波方法 |
CN106872970B (zh) * | 2017-03-07 | 2019-10-18 | 中国电子科技集团公司第三十八研究所 | 一种基于差分进化的多目标变数据率跟踪装置及其方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5798942A (en) * | 1994-04-05 | 1998-08-25 | Trw Inc. | N-best feasible hypotheses multitarget tracking system for space-based early warning systems |
CN101697007A (zh) * | 2008-11-28 | 2010-04-21 | 北京航空航天大学 | 一种基于雷达图像的飞鸟目标识别与跟踪方法 |
CN101770024A (zh) * | 2010-01-25 | 2010-07-07 | 上海交通大学 | 多目标跟踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3179355B2 (ja) * | 1996-11-07 | 2001-06-25 | 三菱電機株式会社 | 多目標追尾方法及び多目標追尾装置 |
-
2012
- 2012-02-27 CN CN201210044944.XA patent/CN102707278B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5798942A (en) * | 1994-04-05 | 1998-08-25 | Trw Inc. | N-best feasible hypotheses multitarget tracking system for space-based early warning systems |
CN101697007A (zh) * | 2008-11-28 | 2010-04-21 | 北京航空航天大学 | 一种基于雷达图像的飞鸟目标识别与跟踪方法 |
CN101770024A (zh) * | 2010-01-25 | 2010-07-07 | 上海交通大学 | 多目标跟踪方法 |
Non-Patent Citations (1)
Title |
---|
JP特开平10-142325A 1998.05.29 |
Also Published As
Publication number | Publication date |
---|---|
CN102707278A (zh) | 2012-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Challa | Fundamentals of object tracking | |
CN106372646B (zh) | 基于srck-gmcphd滤波的多目标跟踪方法 | |
CN109791200B (zh) | 基于轨迹跟踪的导线和塔架分类 | |
WO2008038889A1 (en) | Radar tracking device and method thereof | |
Gongguo et al. | Sensor scheduling for ground maneuvering target tracking in presence of detection blind zone | |
Nordlund et al. | Probabilistic noncooperative near mid-air collision avoidance | |
CN102707278B (zh) | 奇异值分解的多目标跟踪方法 | |
CN116047495B (zh) | 一种用于三坐标雷达的状态变换融合滤波跟踪方法 | |
Blackman et al. | Integration of passive ranging with multiple hypothesis tracking (MHT) for application with angle-only measurements | |
CN102707279B (zh) | 序列ud分解的多目标跟踪方法 | |
CN102707280B (zh) | 稳定结构的多目标跟踪方法 | |
CN116679296A (zh) | 一种基于主辅复合滤波的远程高速目标雷达跟踪方法 | |
CN102426357B (zh) | 一种具有图像确认的多目标跟踪方法 | |
Coetzee et al. | Multifunction radar resource management using tracking optimisation | |
CN102901959B (zh) | Ud分解的多目标跟踪方法 | |
CN102901960B (zh) | 协方差分解的多目标跟踪方法 | |
CN102707277B (zh) | 信息平方根分解的多目标跟踪方法 | |
CN102590810B (zh) | 三角平方根分解的多目标跟踪方法 | |
CN102590809B (zh) | 序列平方根分解的多目标跟踪方法 | |
Wang et al. | Robust CPHD algorithm for maneuvering targets tracking via airborne pulsed Doppler radar | |
Liu et al. | Near Neighbor Cheap JPDA IMM based on amplitude information | |
Sun et al. | Research and Implementation on Multi-beacon Aided AUV Integrated Navigation Algorithm Based on UKF | |
Moawad et al. | Bearing only tracking in cluttered environment using seismic sensors | |
Samuel | Tracking Filter Design for a Full-time Off-road Ground Maneuvering Target | |
Fong | Multisensor data fusion with estimated weights |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140205 |
|
CF01 | Termination of patent right due to non-payment of annual fee |