CN102707277B - Multi-target tracking method for information square root factorization - Google Patents
Multi-target tracking method for information square root factorization Download PDFInfo
- Publication number
- CN102707277B CN102707277B CN201210044943.5A CN201210044943A CN102707277B CN 102707277 B CN102707277 B CN 102707277B CN 201210044943 A CN201210044943 A CN 201210044943A CN 102707277 B CN102707277 B CN 102707277B
- Authority
- CN
- China
- Prior art keywords
- target
- matrix
- target tracking
- centerdot
- variance
- 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 35
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000005070 sampling Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 abstract 1
- 238000001914 filtration Methods 0.000 description 10
- 239000002245 particle Substances 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 2
- 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
Abstract
The invention discloses a multi-target tracking method for information square root factorization. The multi-target tracking method is used for solving the technical problem of failure in target tracking by radars due to instability of numerical structures of existing target tracking methods. The technical scheme includes that the multi-target tracking method is characterized in that an estimation error variance array is subjected to information square root factorization to set up a numerical stable structural model. Subtraction of two positive semidefinite matrixes in the estimation error variance array is avoided, and symmetric matrixes with negative characteristic values cannot be generated in a processing system with limited word length. By means of information square root factorization of the estimation error variance array, the multi-target tracking structural model is set up, subtraction of two positive semidefinite matrixes in the estimation error variance array is avoided, numerical diffusion in the processing system with limited word length is avoided, and accordingly reliability of a target tracking system is guaranteed while failure in target tracking by radars and faults of a whole radar system are avoided.
Description
Technical field
The present invention relates to a kind of Radar Multi Target tracking, the multi-object tracking method that particularly a kind of information square root decomposes, belongs to areas of information technology.
Background technology
Multitarget Tracking is all widely used at military and civil area, can be used for aerial target and detects, follows the tracks of and attack, Air Missile defence, air traffic control, harbour and marine surveillance etc.In recent years, along with the change of battlefield surroundings, the development of antagonism and anti-countermeasure techniques, has produced the series of problems such as the strong clutter of background, low signal-to-noise ratio, low detection probability and high false alarm rate, and the precision of multi-object tracking method and accuracy are had higher requirement.
The object of multiple target tracking is by received information source corresponding to measurement of detector, forms different observation set or track, according to track, estimates the number of tracked target and the kinematic parameter of each target, realizes the tracking to a plurality of targets.For the basic filtering method of multiple goal state estimation have that alpha-beta filtering, alpha-beta-γ filtering, Kalman filtering, EKF, gaussian sum are approximate, optimal nonlinear filtering, particle filter and auto adapted filtering etc.Alpha-beta and alpha-beta-γ wave filter are due to simple in structure, and calculated amount is little, and when computer resource is short in early days, application is very wide.Kalman filtering is a kind of basic skills of multiple target tracking, but need to know the mathematical models of system, and is only applicable to linear system, has limited the application of algorithm.EKF expands to non-linear field by kalman filtering theory, is similar to the conditional probability distribution of state by a Gaussian distribution; And when approximate condition does not meet, Gaussian sum filter device is similar to the conditional probability distribution of state by the weighted sum of a Gaussian distribution.Optimal nonlinear filtering is described the dynamic process of target with Makov transition probability, have good characteristic, but calculated amount is larger, is therefore never used widely.Particle filter adopts stochastic sampling, because calculated amount is too large and particle degenerate problem, is not suitable for practical application.In order to improve particle filter, Unscented kalman filtering adopts deterministic sampling, and the particle point number of sampling is reduced, and avoided the particle point degenerate problem in particle filter, so its application is very wide.Adaptive filter method, by the detection to target maneuver, is adjusted the state of filter parameter or increase wave filter in real time, makes wave filter adapt in real time target travel, is particularly suitable for the tracking to maneuvering target; At present, actual radar tracking system the most frequently used be still JPDA(Joint Probabilistic Data Association, JPDA) method (James A.Roecker, A Class of Near Optimal JPDA Algorithms, IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS, 1994, VOL.30(2): 504-51O), other method great majority are to simplification of JPDA method etc.Yet, in the variance battle array that the method errors such as JPDA are estimated, there are two positive semidefinite matrixs to subtract each other, in the disposal system of limited wordlength, can produce the symmetric matrix that contains positive and negative eigenwert, cause radar tracking enabling objective to lose and follow and whole radar system mistake.
Summary of the invention
In order to solve the unstable technological deficiency that causes radar tracking enabling objective to lose and follow of existing method for tracking target value structure, the invention provides the multi-object tracking method that a kind of information square root decomposes, the method is in the time and measurement renewal of multiple target tracking, by variance of estimaion error battle array being carried out to the decomposition of information square root, set up numerical stability structural model, not having has two positive semidefinite matrixs to subtract each other in the variance battle array of estimation of error, in the disposal system of limited wordlength, can guarantee can not produce the symmetric matrix that contains negative eigenwert, can avoid radar tracking enabling objective to lose follows and whole radar system mistake.
The present invention solves the technical scheme that its technical matters adopts, the multi-object tracking method that a kind of information square root decomposes, and its feature comprises the following steps:
1, in N target following of definition, the discretization model of i target is
x
i(k+1)=Φ(k+1,k)x
i(k)+Λω
i(k),
Wherein:
for state vector, (x, y, z) is the position coordinates of target under ground reference rectangular coordinate system, ω
i(k) represent that variance is Q
i(k) process noise vector, Φ (k+1, k)=Φ=diag[Φ
1, Φ
1, Φ
1] be state-transition matrix,
Γ (t) is matrix of coefficients,
Γ
1=[0?0?1]
T,
T is the sampling period;
The time of i target is updated to:
x
i(k+1/k)=S
i(k+1/k)b
i(k+1/k)
Wherein: n=9, x
i(k+1/k) be i target to (k+1) T one-step prediction value constantly,
for the variance battle array of corresponding one-step prediction error, S
i(k+1/k) be upper triangular matrix, T
1for orthogonal transform matrix, A
i(k+1), B
i(k+1), c
i(k+1) matrix producing for computation process; Starting condition is x
iand S (0/0)
i(0/0);
2, i target observation equation is: z
i(k)=g
i[x
i(k)]+v
i(k)
Wherein: z
i(k) be the r dimension observation vector to i target, g
i[x
i(k)] be corresponding output, v
i(k) represent that variance is R
i(k) measure noise; Calculate
Wherein:
E
i(k+1/k+1) matrix producing for computation process,
for upper triangular matrix, T
2for orthogonal transform matrix, z
ij(k+1) be the j(j=1 of radar to i target, 2 ..., m) individual echo, x
i(k+1/k+1) be i target (k+1) T filter value constantly, λ
ij(k+1) be weight coefficient, and:
3, i tracking method of estimation is:
Wherein: n=9, y
jfor Y
t(k+1) j column vector,
S
ljfor
the capable j column element of l, S
i(k+1/k+1) be upper triangular matrix,
it is i target (k+1) T variance of estimaion error battle array constantly;
Δ
i, j(k+1) be j candidate's echo information vector,
Δ
i,j(k+1)=z
i,j(k+1`)-g
i[x
i(k+1/k)]。
Useful result of the present invention is: by variance of estimaion error battle array being carried out to the decomposition of information square root, set up the multiple target tracking structural model of numerical stability, avoided two positive semidefinite matrixs in the variance battle array of estimation of error to subtract each other, in the disposal system of limited wordlength, there will not be numerical value to disperse, thus guaranteed multi-object tracking method reliability, avoided radar tracking enabling objective lose with and whole radar system mistake.
Below in conjunction with example, the present invention is elaborated.
Embodiment
1, in N target following of definition, the discretization model of i target is
x
i(k+1)=Φ(k+1,k)x
i(k)+Λω
i(k),
Wherein:
for state vector, (x, y, z) is the position coordinates of target under ground reference rectangular coordinate system, ω
i(k) represent that variance is Q
i(k) process noise vector, Φ (k+1, k)=Φ=diag[Φ
1, Φ
1, Φ
1] be state-transition matrix,
Γ (t) is matrix of coefficients,
Γ
1=[0?0?1]
T,
T is the sampling period;
The time of i target is updated to:
x
i(k+1/k)=S
i(k+1/k)b
i(k+1/k)
Wherein: n=9, x
i(k+1/k) be i target to (k+1) T one-step prediction value constantly,
for the variance battle array of corresponding one-step prediction error, S
i(k+1/k) be upper triangular matrix, T
1for orthogonal transform matrix, A
i(k+1), B
i(k+1), c
i(k+1) matrix producing for computation process; Starting condition is x
iand S (0/0)
i(0/0);
2, i target observation equation is: z
i(k)=g
i[x
i(k)]+v
i(k)
Wherein: z
i(k) be the r dimension observation vector to i target, g
i[x
i(k)] be corresponding output, v
i(k) represent that variance is R
i(k) measure noise; Calculate
Wherein:
E
i(k+1/k+1) matrix producing for computation process,
for upper triangular matrix, T
2for orthogonal transform matrix, z
ij(k+1) be the j(j=1 of radar to i target, 2 ..., m) individual echo, x
i(k+1/k+1) be i target (k+1) T filter value constantly, λ
ij(k+1) be weight coefficient, and:
for example get g
i[x
i(k+1)]=[r
i(k+1) α
i(k+1) β
i(k+1)]
t, r
ifor radar can be measured oblique distance, α
ifor angular altitude, β
iposition angle, and
3, i tracking method of estimation is:
Wherein: n=9, s
ljfor
the capable j column element of l, S
i(k+1/k+1) be upper triangular matrix,
be i target (k+1) T variance of estimaion error battle array constantly, y
jfor Y
t(k+1) j column vector,
Δ
i, j(k+1) be j candidate's echo information vector,
Δ
i,j(k+1)=z
i,j(k+1`)-g
i[x
i(k+1/k)]。
3, i tracking method of estimation is:
Wherein: n=9, s
ljfor
the capable j column element of l, S
i(k+1/k+1) be upper triangular matrix,
be i target (k+1) T variance of estimaion error battle array constantly, y
jfor Y
t(k+1) j column vector,
Δ
i, j(k+1) be j candidate's echo information vector,
Δ
i,j(k+1)=z
i,j(k+1`)-g
i[x
i(k+1/k)]。
Claims (1)
1. the multi-object tracking method that information square root decomposes, is characterized in that comprising the following steps:
(1), in N target following of definition, the discretization model of i target is
x
i(k+1)=Φ(k+1,k)x
i(k)+Λω
i(k),
Wherein:
for state vector, (x, y, z) is the position coordinates of target under ground reference rectangular coordinate system, ω
i(k) represent that variance is Q
i(k) process noise vector, Φ (k+1, k)=Φ=diag[Φ
1, Φ
1, Φ
1] be state-transition matrix,
Γ is matrix of coefficients,
t is the sampling period;
The time of i target is updated to:
x
i(k+1/k)=S
i(k+1/k)b
i(k+1/k)
Wherein: n=9, x
i(k+1/k) be i target to (k+1) T one-step prediction value constantly,
for the variance battle array of corresponding one-step prediction error, S
i(k+1/k) be upper triangular matrix, T
1for orthogonal transform matrix, A
i(k+1), B
i(k+1), c
i(k+1) matrix producing for computation process; Starting condition is x
iand S (0/0)
i(0/0);
(2), i target observation equation is: z
i(k)=g
i[x
i(k)]+v
i(k)
Wherein: z
i(k) be the r dimension observation vector to i target, g
i[x
i(k)] be corresponding output, v
i(k) represent that variance is R
i(k) measurement noise; Calculate
Wherein:
e
i(k+1/k+1) matrix producing for computation process,
for upper triangular matrix, T
2for orthogonal transform matrix, z
i,j(k+1) be the j(j=1 of radar to i target, 2 ..., m) individual echo, x
i(k+1/k+1) be i target (k+1) T filter value constantly, λ
i,j(k+1) be weight coefficient, and:
(3), i target following method of estimation is:
Wherein: n=9, b
jfor Y
t(k+1) j column vector,
s
ljfor
the capable j column element of l, S
i(k+1/k+1) be upper triangular matrix,
it is i target (k+1) T variance of estimaion error battle array constantly;
△
i,j(k+1) be j candidate's echo information vector,
△
i,j(k+1)=z
i,j(k+1)-g
i[x
i(k+1/k)]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210044943.5A CN102707277B (en) | 2012-02-27 | 2012-02-27 | Multi-target tracking method for information square root factorization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210044943.5A CN102707277B (en) | 2012-02-27 | 2012-02-27 | Multi-target tracking method for information square root factorization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102707277A CN102707277A (en) | 2012-10-03 |
CN102707277B true CN102707277B (en) | 2014-02-05 |
Family
ID=46900207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210044943.5A Expired - Fee Related CN102707277B (en) | 2012-02-27 | 2012-02-27 | Multi-target tracking method for information square root factorization |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102707277B (en) |
Family Cites Families (4)
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 |
JP3179355B2 (en) * | 1996-11-07 | 2001-06-25 | 三菱電機株式会社 | Multi-target tracking method and multi-target tracking device |
CN101697007B (en) * | 2008-11-28 | 2012-05-16 | 北京航空航天大学 | Radar image-based flyer target identifying and tracking method |
CN101770024B (en) * | 2010-01-25 | 2011-08-31 | 上海交通大学 | Multi-target tracking method |
-
2012
- 2012-02-27 CN CN201210044943.5A patent/CN102707277B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102707277A (en) | 2012-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106372646B (en) | Multi-target tracking method based on SRCK-GMCPHD filtering | |
CN101975575B (en) | Multi-target tracking method for passive sensor based on particle filtering | |
CN110361744B (en) | RBMCDA underwater multi-target tracking method based on density clustering | |
CN105301584A (en) | IPPHDF maneuvering multi-target tracking method of simultaneously solving range ambiguities | |
CN110738275A (en) | UT-PHD-based multi-sensor sequential fusion tracking method | |
CN111428369A (en) | Method for calculating confidence of space target collision early warning result | |
Nordlund et al. | Probabilistic noncooperative near mid-air collision avoidance | |
CN102707278B (en) | Multi-target tracking method for singular value decomposition | |
CN102707279B (en) | Multi-target tracking method for sequence UD decomposition | |
Blackman et al. | Integration of passive ranging with multiple hypothesis tracking (MHT) for application with angle-only measurements | |
CN102707280B (en) | Structurally stable multi-target tracking method | |
CN102901959B (en) | Multi-target tracking method by means of UD factorization | |
CN116679296A (en) | Remote high-speed target radar tracking method based on main and auxiliary composite filtering | |
CN102707277B (en) | Multi-target tracking method for information square root factorization | |
CN102426357B (en) | Multi-target tracking method with function of image confirmation | |
CN102590810B (en) | Multiple target tracking method for triangular square root decomposition | |
CN102901960B (en) | Multi-target tracking method by means of covariance factorization | |
CN102590809B (en) | Multiple target tracking method for sequential square root decomposition | |
CN104931932A (en) | Improved debiased coordinate transform Kalman filtering method | |
CN116047495B (en) | State transformation fusion filtering tracking method for three-coordinate radar | |
CN104330772A (en) | Bi-station location method based on multidirection optimizing all tracking type UKF filtering algorithm | |
Liu et al. | Near neighbor cheap JPDA IMM based on amplitude information | |
Xie et al. | Research for underwater target tracking by using multi-sonar | |
Kaplan et al. | Comparison of proposed target tracking algorithm, GRNNa, to Kalman Filter in 3D environment | |
Gu et al. | Research on Asynchronous Information Fusion Technology of the Distributed Active/Passive Imaging Detection System |
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 |