CN102901959B - Ud分解的多目标跟踪方法 - Google Patents
Ud分解的多目标跟踪方法 Download PDFInfo
- Publication number
- CN102901959B CN102901959B CN 201210044945 CN201210044945A CN102901959B CN 102901959 B CN102901959 B CN 102901959B CN 201210044945 CN201210044945 CN 201210044945 CN 201210044945 A CN201210044945 A CN 201210044945A CN 102901959 B CN102901959 B CN 102901959B
- Authority
- CN
- China
- Prior art keywords
- overbar
- target
- lambda
- gamma
- vector
- 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 27
- 239000011159 matrix material Substances 0.000 claims description 26
- 238000005070 sampling Methods 0.000 claims description 6
- 238000009792 diffusion process Methods 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
- 238000005259 measurement Methods 0.000 description 2
- 230000003044 adaptive effect 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
本发明公开了一种UD分解的多目标跟踪方法,用于解决现有的目标跟踪方法数值结构不稳定导致雷达跟踪过程目标失跟的技术问题。技术方案是对估计误差方差阵的UD分解,建立数值稳定结构模型,没有误差估计的方差阵中有两个半正定矩阵相减,在有限字长的处理***中能够保证不会产生含有负特征值的对称矩阵。通过对估计误差方差阵的UD分解,建立了数值稳定的多目标跟踪结构模型,避免了误差估计的方差阵中两个半正定矩阵相减,在有限字长的处理***中不会出现数值发散,从而保证了目标跟踪***的可靠性、避免了雷达跟踪过程目标失跟和整个雷达***性错误。
Description
技术领域
本发明涉及一种雷达跟踪方法,特别涉及一种UD分解的多目标跟踪方法,属于信息技术领域。
背景技术
多目标跟踪技术在军用及民用领域均有广泛的应用,可用于空中目标检测、跟踪与攻击,空中导弹防御,空中交通管制,港口和海洋监视等。近年来,随着战场环境的改变,对抗和反对抗技术的发展,产生了背景强杂波、低信噪比、低检测概率和高虚警率等一系列问题,对多目标跟踪方法的精度和准确性提出了更高的要求。
多目标跟踪的目的是将探测器所接收到的量测对应不同的信息源,形成不同观测集合或轨迹,根据轨迹估计被跟踪目标的数目以及每一目标的运动参数,实现对多个目标的跟踪。用于多目标状态估计的基本滤波方法有α-β滤波、α-β-γ滤波、卡尔曼滤波、扩展卡尔曼滤波、高斯和近似、最优非线性滤波、粒子滤波和自适应滤波等。α-β和α-β-γ滤波器由于结构简单,计算量小,在早期计算机资源短缺时应用很广。卡尔曼滤波是多目标跟踪的一种基本方法,但是需要知道***的精确数学模型,并且只适用于线性***,限制了算法的应用。扩展卡尔曼滤波将卡尔曼滤波理论扩展到非线性领域,用一个高斯分布来近似状态的条件概率分布;而当近似条件不满足时,高斯和滤波器则用一个高斯分布的加权和来近似状态的条件概率分布。最优非线性滤波使用Makov转移概率来描述目标的动力学过程,具有很好的特性,但是计算量较大,因此一直没有得到广泛应用。粒子滤波采用随机采样,由于计算量太大和粒子退化问题,不适合实际应用。为了改进粒子滤波,无迹卡尔曼滤波采用确定性采样,使得采样的粒子点个数减少,避免了粒子滤波中的粒子点退化问题,因此其应用领域很广。自适应滤波方法通过对目标机动的检测,实时调整滤波器参数或增加滤波器的状态,使滤波器实时适应目标运动,特别适合对机动目标的跟踪;目前,在实际雷达跟踪系 统最常用的仍然为JPDA(Joint Probabilistic Data Association,联合概率数据关联)方法(James A.Roecker,A Class of Near Optimal JPDA Algorithms,IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS,1994,VOL.30(2):504-51O),其它方法大多数是对JPDA方法的简化等。然而,JPDA等方法误差估计的方差阵中有两个半正定矩阵相减,在有限字长的处理***中会产生含有正负特征值的对称矩阵,导致雷达跟踪过程目标失跟和整个雷达***性错误。
发明内容
为了解决现有的目标跟踪方法数值结构不稳定导致雷达跟踪过程目标失跟的技术缺陷,本发明提供一种UD分解的多目标跟踪方法,该方法在多目标跟踪的时间和测量更新中,通过对估计误差方差阵的UD分解,建立数值稳定结构模型,没有误差估计的方差阵中有两个半正定矩阵相减,在有限字长的处理***中能够保证不会产生含有负特征值的对称矩阵,可以避免雷达跟踪过程目标失跟和整个雷达***性错误。
本发明解决其技术问题采用的技术方案是,一种UD分解的多目标跟踪方法,其特征包括以下步骤:
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)
Y(k/k-1)=[ΦUi(k-1/k-1) Λ],
Ui(k/k-1)为单位上三角矩阵,Di(k/k-1)为对角矩阵;Dll(k/k-1)为Di(k/k-1)的第l行第l列元素,Ujl(k/k-1)为为Ui(k/k-1)的第j行第l列元素;←表示改写,例如a←b为用b值改写a;初始条件为xi(0/0)和
2、第i个目标观测方程为:zi(k)=gi[xi(k)]+vi(k)
其中:zi(k)为对第i个目标的r维观测向量,gi[xi(k)]为对应的输出,vi(k)表示方差为Ri(k)测量噪声;计算
其中: 为 的第j列向量, 为单位上三角矩阵, 为对角矩阵,zij(k)为雷达对第i个目标的第j(j=1,2,…,m)个回波,xi(k/k)为第i个目标kT时刻的滤波值,λij(k)为权系数,且: ←表示改写,例如a←b为用b值改写a;
3、第i个跟踪估计方法为:
其中:←表示改写,例如c←b为用b值改写c;Ui(k/k)为单位上三角矩阵,Di(k/k)为对角矩阵;Dll(k/k)为Di(k/k)的第l行第l列元素,Ujl(k/k)为为Ui(k/k)的第j行第l列元素; 为第i个目标kT时刻估计误差的方差阵;aj为AT(k)的第j列向量,
Δi,j(k)为第j个候选回波信息向量,
Δi,j(k)=zi,j(k)-gi[xi(k/k-1)]。
本发明的有益结果是:通过对估计误差方差阵的UD分解,建立了数值稳定的多目标跟踪结构模型,避免了误差估计的方差阵中两个半正定矩阵相减,在有限字长的处理***中不会出现数值发散,从而保证了目标跟踪***的可靠性、避免了雷达跟踪过程目标失跟和整个雷达***性错误。
下面结合实例对本发明作详细说明。
具体实施方式
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)
Y(k/k-1)=[ΦUi(k-1/k-1)Λ],
Ui(k/k-1)为单位上三角矩阵,Di(k/k-1)为对角矩阵;Dll(k/k-1)为Di(k/k-1)的第l行第l列元素,Ujl(k/k-1)为为Ui(k/k-1)的第j行第l列元素;←表示改写,例如a←b为用b值改写a;初始条件为xi(0/0)和
2、第i个目标观测方程为:zi(k)=gi[xi(k)]+vi(k)
其中:zi(k)为对第i个目标的r维观测向量,gi[xi(k)]为对应的输出,vi(k)表示方差为Ri(k)测量噪声;计算
其中:←表示改写,例如a←b为用b值改写a; 为 的第j列向量, 为单位上三角矩阵, 为对角矩阵,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个跟踪估计方法为:
其中:←表示改写,例如c←b为用b值改写c;Ui(k/k)为单位上三角矩阵,Di(k/k)为 对角矩阵;Dll(k/k)为Di(k/k)的第l行第l列元素,Ujl(k/k)为为Ui(k/k)的第j行第l列元素; 为第i个目标kT时刻估计误差的方差阵;aj为AT(k)的第j列向量, Gi(k)dT(I-ΩuuT)],
Δi,j(k)为第j个候选回波信息向量,Δi,j(k)=zi,j(k)-gi[xi(k/k-1)]。
Claims (1)
1.一种UD分解的多目标跟踪方法,其特征在于包括以下步骤:
(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)
其中:xi(k/k-1)为第i个目标在kT时刻的一步预测值,为对应的一步预测误差的方差阵,yj为YT(k/k-1)的第j列向量,
Ui(k/k-1)为单位上三角矩阵,Di(k/k-1)为对角矩阵;Dll(k/k-1)为Di(k/k-1)的第l行第l列元素,Ujl(k/k-1)为Ui(k/k-1)的第j行第l列元素;←表示改写,a←b则表示为用b值改写a;初始条件为xi(0/0)和
(2)、第i个目标观测方程为:zi(k)=gi[xi(k)]+vi(k)
其中:zi(k)为对第i个目标的r维观测向量,gi[xi(k)]为对应的输出,vi(k)表示方差为Ri(k)测量噪声;计算
其中:为的第j列向量, 为单位上三角矩阵, 为对角矩阵,zi,j(k)为雷达对第i个目标的第j(j=1,2,…,m)个回波,xi(k/k)为第i个目标kT时刻的滤波值,λi,j(k)为权系数,且: ←表示改写,a←b则表示为用b值改写a;
(3)、第i个目标跟踪估计方法为:
其中:Ui(k/k)为单位上三角矩阵,Di(k/k)为对角矩阵;Dll(k/k)为Di(k/k)的第l行第l列元素,Ujl(k/k)为为Ui(k/k)的第j行第l列元素;为第i个目标kT时刻估计误差的方差阵;aj为AT(k)的第j列向量,
Δ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 |
---|---|---|---|
CN 201210044945 CN102901959B (zh) | 2012-02-27 | 2012-02-27 | Ud分解的多目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210044945 CN102901959B (zh) | 2012-02-27 | 2012-02-27 | Ud分解的多目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102901959A CN102901959A (zh) | 2013-01-30 |
CN102901959B true CN102901959B (zh) | 2013-12-25 |
Family
ID=47574300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201210044945 Expired - Fee Related CN102901959B (zh) | 2012-02-27 | 2012-02-27 | Ud分解的多目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102901959B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103472445B (zh) * | 2013-09-18 | 2015-06-17 | 电子科技大学 | 一种针对多目标场景的检测跟踪一体化方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL1020287C2 (nl) * | 2002-04-02 | 2003-10-03 | Thales Nederland Bv | Werkwijze voor meerdoelendetectie, met name voor toepassing in rondzoekradars met meerbundelvorming in elevatie. |
AU2007201638B2 (en) * | 2006-11-09 | 2010-11-18 | Raytheon Canada Limited | A track quality based multi-target tracker |
PL2452204T3 (pl) * | 2008-07-08 | 2018-05-30 | Thales | Wielokierunkowe przetwarzanie danych dla pasywnych radarów multistatycznych i multikanałowych |
CN101614817A (zh) * | 2009-06-24 | 2009-12-30 | 北京航空航天大学 | 一种基于地面动目标指示雷达***的多目标跟踪方法 |
CN101770024B (zh) * | 2010-01-25 | 2011-08-31 | 上海交通大学 | 多目标跟踪方法 |
-
2012
- 2012-02-27 CN CN 201210044945 patent/CN102901959B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102901959A (zh) | 2013-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Challa | Fundamentals of object tracking | |
CN101975575B (zh) | 基于粒子滤波的被动传感器多目标跟踪方法 | |
CN106372646B (zh) | 基于srck-gmcphd滤波的多目标跟踪方法 | |
CN110361744B (zh) | 基于密度聚类的rbmcda水下多目标跟踪方法 | |
CN105372659A (zh) | 道路交通监控多目标探测跟踪方法及跟踪*** | |
CN116794647A (zh) | 基于轨迹跟踪的导线和塔架分类 | |
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 | |
CN102707279B (zh) | 序列ud分解的多目标跟踪方法 | |
CN102707278B (zh) | 奇异值分解的多目标跟踪方法 | |
CN116047495B (zh) | 一种用于三坐标雷达的状态变换融合滤波跟踪方法 | |
CN109115228A (zh) | 一种基于加权最小二乘容积卡尔曼滤波的目标定位方法 | |
Blackman et al. | Integration of passive ranging with multiple hypothesis tracking (MHT) for application with angle-only measurements | |
CN102901959B (zh) | Ud分解的多目标跟踪方法 | |
CN105353353B (zh) | 多重搜索粒子概率假设密度滤波的多目标跟踪方法 | |
Georgy et al. | Unconstrained underwater multi-target tracking in passive sonar systems using two-stage PF-based technique | |
CN102707280B (zh) | 稳定结构的多目标跟踪方法 | |
CN102707277B (zh) | 信息平方根分解的多目标跟踪方法 | |
CN102590810B (zh) | 三角平方根分解的多目标跟踪方法 | |
CN102901960B (zh) | 协方差分解的多目标跟踪方法 | |
CN102426357B (zh) | 一种具有图像确认的多目标跟踪方法 | |
CN102590809B (zh) | 序列平方根分解的多目标跟踪方法 | |
Liu et al. | Near Neighbor Cheap JPDA IMM based on amplitude information | |
Ikram et al. | A new data association method for 3-D object tracking in automotive applications | |
Moawad et al. | Bearing only tracking in cluttered environment using seismic sensors |
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: 20131225 |
|
CF01 | Termination of patent right due to non-payment of annual fee |