CN101975575A - Multi-target tracking method for passive sensor based on particle filtering - Google Patents
Multi-target tracking method for passive sensor based on particle filtering Download PDFInfo
- Publication number
- CN101975575A CN101975575A CN 201010507225 CN201010507225A CN101975575A CN 101975575 A CN101975575 A CN 101975575A CN 201010507225 CN201010507225 CN 201010507225 CN 201010507225 A CN201010507225 A CN 201010507225A CN 101975575 A CN101975575 A CN 101975575A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- target
- sample
- mover
- 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.)
- Granted
Links
- 239000002245 particle Substances 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001914 filtration Methods 0.000 title claims abstract description 10
- 238000005457 optimization Methods 0.000 claims abstract description 27
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims description 34
- 238000012952 Resampling Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 229920001184 polypeptide Polymers 0.000 claims 1
- 102000004196 processed proteins & peptides Human genes 0.000 claims 1
- 108090000765 processed proteins & peptides Proteins 0.000 claims 1
- 238000005070 sampling Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a multi-target tracking method for a passive sensor based on particle filtering, which belongs to the technical field of guidance and mainly solves the problems of easy divergent tracking and inaccurate target state estimation in the traditional multi-target tracking method. The method optimizes distribution of multi-target samples through particle swarm optimization and sample mixing sampling algorithms and tracks the multi-target combined with a joint probability data association algorithm. The method comprises the following steps of: firstly, optimizing the distribution of multi-target joint samples by utilizing the particle swarm optimization algorithm so that the multi-target joint samples are gathered in a high likelihood region with a bigger probability of occurrence of a real target; secondly, calculating an associated probability between the targets and observation and the posterior probability distribution of the targets by utilizing the samples; and finally, decomposing a joint sample weight into the corresponding target sample in a likelihood way according to each target sample in the re-sampling process, and independently re-sampling each target according to the decomposed weight, and further optimizing the distribution of the target sample so as to improve the precision of target tracking.
Description
Technical Field
The invention belongs to the technical field of guidance and relates to target tracking. In particular to a passive sensor multi-target tracking method based on particle swarm optimization and sequential Monte Carlo, which can be used for systems such as infrared guidance and the like.
Background
In multi-target tracking, due to the influence of target missing detection and clutter, uncertainty exists in association between measurement obtained by a sensor and a target, and angle information measured under a passive condition is a nonlinear function of a target state, so that two problems of data association and nonlinear filtering of measurement and the target need to be solved for accurately estimating the target state to realize target tracking.
The traditional multi-target tracking method comprises a nearest neighbor method NN, a joint probability data association JPDA and a multi-hypothesis tracking MHT algorithm, wherein the nearest neighbor method is to directly associate the nearest measurement of an off-target state with a target, when the measurement precision is higher, the tracking performance is better, and when the measurement precision is reduced, the tracking performance is also seriously reduced; multi-hypothesis tracking is exhaustive of all possible correlation events between the target and the measurements, and gradually extends over time, with the disadvantage that the computation time will grow exponentially with the number of targets and measurements; JPDA is one of the most effective methods for solving data association so far, and assigns a certain probability to each pair of target and measured association, and then completes estimation of the posterior probability and state of the target by predicting and updating two steps in combination with Bayesian criterion.
An algorithm SMC based on sequential Monte Carlo is a nonlinear filtering method developed in recent years, a learner combines JPDA with SMC to solve the problem of multi-target tracking, posterior probability distribution of a moving target is fitted by utilizing a certain number of samples and corresponding weights, and the SMC can fit any probability distribution theoretically when the number of sampled samples tends to be infinite. However, in practical applications, considering the comprehensive requirements of tracking accuracy and real-time performance, the number of samples is usually limited, and a sample depletion phenomenon occurs during sampling and resampling, so that the samples lose diversity, and the state estimation is unstable, resulting in tracking divergence.
Disclosure of Invention
Aiming at the problems, the invention provides a passive sensor multi-target tracking method based on particle filtering, so as to keep the diversity of samples and improve the tracking precision of a target.
The key technology for realizing the invention is as follows: the particle swarm optimization algorithm is utilized to optimize the distribution of the multi-target combined samples, so that the multi-target combined samples are gathered to a high-likelihood region of each target state, namely a region with high occurrence probability of a real target, thus the samples for filtering have rich diversity and the importance of each sample is improved; the joint samples are used for calculating the association probability and the target filtering distribution between a target and measurement, in the resampling process, sampling is not carried out according to the weight values of the joint samples generated by multi-target series connection, but the weight values of the joint samples are decomposed into corresponding target samples according to the likelihood of each target sample, the distribution of the target samples is further optimized, and the target tracking precision is improved, and the specific implementation steps comprise the following steps:
(1) extracting target samples according to the initial distribution of each target, and constructing a combined sample:
wherein N represents the combined sample number, i represents the target number, N represents the number of combined samples, c represents the number of targets,the sample of the target i in the nth joint sample at the time point 0 is represented, and the initial weight of each joint sample is taken as
(2) Calculating a prediction joint sample at the time t:
i∈[1,c],n∈[1,N]t is more than or equal to 1, wherein,is the sample of target i in the nth joint sample at time t;
(3) optimizing the particle swarm as follows:
(3a) taking each target sample in the prediction combined sample at the time t as an initial sample for particle swarm optimizationIs a target sampleGiving an initial velocity:
(3b) calculating a target sample at time tThe likelihood of the sensor measurement is expressed asWherein k is 1, L, m is a particle swarm optimization iteration number, and m is more than or equal to 5, which is the set total particle swarm optimization iteration number;
(3c) finding out individual optimal solution of each sample in the target i according to the likelihood of each target sample in the 1 st iteration to the k th iteration
(3d) According to the likelihood of all samples in the ith target, finding out the global optimum solution in all samples of the target
(3e) Obtaining a target sample by utilizing an update equation in a particle swarm optimization algorithmPosition in the (k + 1) th iterationAnd velocity
(3f) Repeating the steps (3b) to (3e) m times to obtain a combined sample after particle swarm optimization:
(4) updating and normalizing the combined sample weight according to the following steps:
(4a) calculating the average value of the target i measured at the time t according to the measurement value corresponding to the optimized target sampleSum varianceSelect out of the satisfactionAll effective measurements ofj∈[1,Mt]Wherein, ytObtained for passive sensorsMeasuring, Epsilon is 9.21 as the set threshold, MtThe number of all effective measurements at the time t;
(4b) enumerating the effective measurementCorrelation event phi with target ii,j;
(4c) Computing effective metricsSample form based association likelihood with target iCalculating an edge correlation event phi in the nth combined sample according to Markov and Bayesian rules of target motioni,jProbability of (c): p (phi)i,j|Yt)nWherein Y istRepresents the set of all valid measurements from time 1 to time t;
(4d) the probability of all the associated events of the nth joint sample is summed to obtain the weight of the nth joint sampleAnd normalizing the weight value to obtain a normalized weight value
(5) From the combined samples and their corresponding weightsEstimating each target state by weighting and summing the combined samples, outputting the result, and simultaneously executing the step (6);
(6) decomposing and resampling the combined sample weight according to the following steps:
(6a) normalizing weight of the nth combined sampleWriting a form of c target sample weight summation:
wherein the weight of the ith target sampleObtaining the likelihood of the ith target sample through likelihood calculation;
(6b) the weight of the ith target sample is taken from the N combined sample weightsBased on these weights, N new samples are sampledWherein the sampleThe corresponding weight is Respectively obtaining the first sample before resampling of a target i at the moment t and a weight value corresponding to the first sample;
(7) and (5) repeating the step (2) and continuing to track the target.
The invention has the following advantages:
(1) according to the sampling particle swarm optimization algorithm, the distribution condition of the target samples is improved, so that the target samples are gathered to a high-likelihood region with high probability of occurrence of the target, the importance of each sample is improved, and higher tracking accuracy can be achieved under the condition of fewer target samples;
(2) the invention considers the mutual influence and coupling condition of the similar targets, and performs mixed sampling on the multi-target samples, namely, in the resampling stage of the target samples, the combined sample weight is decomposed into the corresponding target samples according to the likelihood of each target, and then each target is independently resampled according to the decomposed weight, so that the large weight value sample in each target is copied, the small weight value sample is suppressed, the target sample distribution is further optimized, and the tracking precision is improved.
Drawings
FIG. 1 is an overall flow diagram of the present invention;
FIG. 2 is a schematic diagram of particle swarm optimization particle velocity and position updates used in the present invention;
FIG. 3 is a diagram of the effect of one target tracking with the present invention;
figure 4 is a root mean square error plot for position tracking using the present invention.
Detailed Description
Introduction of basic theory
1. System equation
In a Cartesian coordinate system, the position and the speed of the system state in the x and y directions are taken, and the following nonlinear dynamic system model can be established:
where i is 1, L, c denotes the number of objects, c denotes the total number of objects, respectively representing the coordinates of the object i in the x-direction and the y-direction,respectively representing the speed of the target i in the x-direction and the y-direction, the subscript te N representing the time, the state noiseObey variance ofZero mean gaussian distribution, F, G being the state transition matrix and the input matrix, h being the nonlinear function, and the measurement noise etSubject to a zero mean gaussian distribution with variance R,and etIndependently of each other, ytIs the measured value of the sensor.
In the invention, it is assumed that the passive sensor can only observe the azimuth information of the target, so h is defined as follows:
wherein x iso,yoIs the position of the sensor.
2. Particle swarm optimization
Setting a group X of N particles in a D-dimension search space as X ═ X1,L xn L,xNIn which the N ∈ [1, N ]]The position and velocity of each particle being xn=(xn1,xn2,L xnD) And vn=(vn1,vn2,L,vnD) And the optimal solution of its position is sn=(sn1,sn2,L,snD) And the optimal solution for the whole population position is g ═ g (g)1,g2,L gD) Then, the nth particle position and velocity in the kth particle swarm optimization iteration are updated as follows:
wherein k is 1, L, m is the number of the particle swarm optimization iteration, m is the preset total number of the particle swarm optimization iteration, D is 1, L, D is the number of the particle dimension,d-th dimension data indicating a position of the n-th particle,d-th dimension data representing the velocity of the nth particle,represents an optimal solution of the d-th dimensional data in the position of the n-th particle,representing the optimal solution of the d-th dimension data in the positions of all particles in the whole population, c1And c2Is a learning factor whose classical value is between (0, 2)Normal number, zeta and eta are uniformly distributed pseudo random numbers between (0, 1); due to the fact thatIndicates the difference vector between the current position of the nth particle and its optimal position, so c1The capability of searching the nth particle to the optimal position of the nth particle is characterized; whileThen the difference vector representing the current position of the nth particle and the optimal position of the particle in the whole population, so c2The capability of searching the optimal position of the particle to the whole population is represented; formula 4) includesThe update representing the nth particle velocity is also dependent on its pre-iteration velocity. The updating of the particle position and velocity of the basic particle swarm optimization algorithm is shown in figure 2.
Secondly, the invention relates to a passive sensor multi-target tracking method based on particle filtering
Referring to fig. 1, the specific implementation steps of the present invention include the following:
Let initial time t equal to 0, according to the initial distribution of target iExtracting a target sampleParallel-serial configuration of combined samplesi∈[1,c],n∈[1,N]N is the number of samples extracted, c is the number of targets, wherein, andrespectively representing the coordinates of the ith target sample in the nth combined sample in the x direction and the y direction,andrespectively representing the speeds of the ith target sample in the nth combined sample in the x direction and the y direction, and taking the initial weight of the nth combined sample as
Step 2, calculating prediction combined sample at time t
Target sample according to time t-1And equation of state 1) calculating the predicted sample at time tThese prediction samples are used to construct the joint sample at time t:t is more than or equal to 1, wherein,representing the ith target sample in the nth joint sample at time t.
Step 3, performing particle swarm optimization on the predicted combined sample
(3.1) taking each target sample in the prediction combined samples at the time t as initial samples for particle swarm optimizationSample(s)The initial speeds of (a) are:
Wherein k is 1, L, m is the particle swarm optimization iteration number, m is more than or equal to 5, y is the set total particle swarm optimization iteration number, ytR is a measurement covariance matrix for the measurement values obtained by the sensor,is composed of a target sampleCalculating a measurement value according to the measurement update equation 2);
(3.3) finding the ith target sample in the nth joint sample in the 1 st to kth iterationsMinimum value of (2) byRepresenting that the corresponding sample is taken as the individual optimal solution of the ith target sample in the nth joint sample
(3.4) finding all samples in the ith targetMinimum value of (2) byRepresenting that the corresponding sample is taken as the global optimal solution of all samples in the ith target
(3.5) bondingAnd a sampleSpeed in the k-th iterationUpdating equation 4) according to the particle swarm optimization speedSpeed in the k +1 th iteration
(3.6) bondingAndupdating samples according to position update equation 5) for particle swarm optimization
(3.7) repeating the steps (3.2) - (3.6) m times to obtain an optimized combined sample
Step 4, updating and normalizing the combined sample weight
(4.1) according to the sampleCalculating the mean value of the target i measured at the time tSum covariance
(4.2) use of the mean valueSum covarianceSelecting a set of valid measurements satisfying the condition of equation 11)
Wherein j is 1, L Mt,MtThe total number of effective measurement is shown, and epsilon is 9.21 which is a set threshold value;
(4.5) calculating the edge correlation event phi in the nth combined sample according to Markov property and Bayesian rule of target motioni,jProbability of p (phi)i,j|Yt)n,
Wherein, PfAnd PdRespectively representing false alarm probability and target detection probability, c0Is a correlation event phii,jNumber of undetected targets in, YtRepresents the set of all valid measurements from time 1 to time t;
(4.6) summing the probabilities of all the associated events of the n-th joint sample to obtain the weight of the n-th joint sampleAnd normalizing the weight value to obtain a normalized weight value
Step 5, estimating the target state
Utilizing the combined sample obtained in the step 3And the combined sample weight obtained in step 4As in equation 16) to estimate the target state, output as a result, and simultaneously performs step 6,
(6.1) calculating the likelihood of the ith target sample in the nth joint sample at the time t:
(6.2) calculating the weight of the ith target sample in the nth joint sample according to the likelihood of the ith target sample in the nth joint sample:
the normalized weight of the nth combined sampleThe form of summation of the weights of c target samples can be written:
(6.3) from the N combined sample weights, each weight of the ith target sample is taken out to formBased on these weights, N new samples are sampledWherein the sampleThe corresponding weight is The first sample before target i resampling at time t and the corresponding weight value are respectively.
And 7, repeating the step 2 and continuously tracking the target.
The effect of the invention can be further illustrated by the following experimental simulation:
1. simulation conditions and parameters
Simulation scenario as shown in fig. 3, the real state of each target appearing in the simulation scenario is x ═ x, vx,y,vy]', x, y are the coordinates of each object in the x and y directions of the Cartesian coordinate system, vx,vyRespectively for the speed of each target in the x-direction and the y-direction. The state equation and the measurement equation of the target are respectively shown in the formulas 1) and 2), and each target is subject to a constant speed model:
where T is the sampling time interval, the sensor provides azimuthal information,the simulation parameters are shown in table 1 below,
TABLE 1 Experimental simulation parameters
2. Simulation content and result analysis
In a pure azimuth tracking simulation experiment of three targets under the condition of three sensors, the root mean square error RMSE and the tracking loss rate of the positions of the tracking method of the invention and two multi-target tracking methods of the existing iMC-JPDA and jMC-JPDA are compared in the simulation experiment, and simulation results are respectively shown in FIG. 4 and Table 2, wherein:
FIG. 4(a) is a graph comparing the root mean square error of the positions of the method of the present invention with iMC-JPDA and jMC-JPDA for the sample number N equal to 30;
FIG. 4(b) is a graph comparing the root mean square error of the positions of the method of the present invention with iMC-JPDA and jMC-JPDA at a sample number of N-50;
FIG. 4(c) is a graph comparing the root mean square error of the positions of the method of the present invention with iMC-JPDA and jMC-JPDA at 80 samples N;
FIG. 4(d) is a graph comparing the root mean square error of the positions of the method of the present invention with iMC-JPDA and jMC-JPDA for a sample number of N-100;
as can be seen from FIGS. 4(a) -4 (d), the RMSE for all three tracking methods decreased as the number of samples increased, but the RMSE for the present invention was consistently lower than that for the iMC-JPDA and jMC-JPDA methods.
Table 2 shows the loss of tracking of the inventive method compared to the existing iMC-JPDA and jMC-JPDA,
TABLE 2 comparison of loss of tracking rates of the inventive method with iMC-JPDA and jMC-JPDA
As can be seen from Table 2, the tracking loss rate of the method of the present invention is significantly lower than that of the iMC-JPDA and jMC-JPDA tracking methods under the condition of the same number of samples, and when the number of samples N exceeds 30, the tracking loss of the method of the present invention does not occur.
Claims (2)
1. A passive sensor multi-target tracking method based on particle filtering comprises the following steps:
(1) extracting target samples according to the initial distribution of each target, and constructing a combined sample:
wherein N represents the combined sample number, i represents the target number, N represents the number of combined samples, c represents the number of targets,the sample of the target i in the nth joint sample at the time point 0 is represented, and the initial weight of each joint sample is taken as
(2) Calculating a prediction joint sample at the time t:
i∈[1,c],n∈[1,N]t is more than or equal to 1, wherein,is the sample of target i in the nth joint sample at time t;
(3) optimizing the particle swarm as follows:
(3a) taking each target sample in the prediction combined sample at the time t as an initial sample for particle swarm optimizationIs a target sampleGiving an initial velocity:
(3b) calculating a target sample at time tThe likelihood of the sensor measurement is expressed asWherein k is 1, L, m is a particle swarm optimization iteration number, and m is more than or equal to 5, which is the set total particle swarm optimization iteration number;
(3c) finding out individual optimal solution of each sample in the target i according to the likelihood of each target sample in the 1 st iteration to the k th iteration
(3d) According to the likelihood of all samples in the ith target, finding out the global optimum solution in all samples of the target
(3e) Obtaining a target sample by utilizing an update equation in a particle swarm optimization algorithmPosition in the (k + 1) th iterationAnd velocity
(3f) Repeating the steps (3b) to (3e) m times to obtain a combined sample after particle swarm optimization:
(4) updating and normalizing the combined sample weight according to the following steps:
(4a) calculating the average value of the target i measured at the time t according to the measurement value corresponding to the optimized target sampleSum varianceSelect out of the satisfactionAll effective measurements ofj∈[1,Mt]Wherein, ytFor the measurements obtained by the passive sensor, ∈ 9.21 is the set threshold, MtThe number of all effective measurements at the time t;
(4c) Computing effective metricsSample form based association likelihood with target iCalculating an edge correlation event phi in the nth combined sample according to Markov and Bayesian rules of target motioni,jProbability of (c): p (phi)i,j|Yt)nWherein Y istRepresents the set of all valid measurements from time 1 to time t;
(4d) the probability of all the associated events of the nth joint sample is summed to obtain the weight of the nth joint sampleAnd normalizing the weight value to obtain a normalized weight value
(5) From the combined samples and their corresponding weightsEstimating each target state by weighting and summing the combined samples, outputting the result, and simultaneously executing the step (6);
(6) decomposing and resampling the combined sample weight according to the following steps:
(6a) normalizing weight of the nth combined sampleWriting a form of c target sample weight summation:
wherein the weight of the ith target sampleObtaining the likelihood of the ith target sample through likelihood calculation;
(6b) the weight of the ith target sample is taken from the N combined sample weightsBased on these weights, N new samples are sampledWherein the sampleThe corresponding weight is Respectively obtaining the first sample before resampling of a target i at the moment t and a weight value corresponding to the first sample;
(7) and (5) repeating the step (2) and continuing to track the target.
2. The polypeptide of claim 1The target tracking method, wherein the weight of the ith target sample in step (6a)The likelihood calculation of the ith target sample is used for calculating the likelihood of the ith target sample according to the following steps:
(2.1) calculating the likelihood of the ith target sample in the nth joint sample at the time t:
wherein, ytIn order to obtain a measurement value for the sensor,is a target sampleThe corresponding measured value is measured by the corresponding measuring instrument,the covariance matrix measured at time t for target i,
wherein,
(2.2) calculating the weight of the ith target sample in the nth joint sample according to the likelihood of the ith target sample in the nth joint sample:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105072258A CN101975575B (en) | 2010-10-15 | 2010-10-15 | Multi-target tracking method for passive sensor based on particle filtering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105072258A CN101975575B (en) | 2010-10-15 | 2010-10-15 | Multi-target tracking method for passive sensor based on particle filtering |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101975575A true CN101975575A (en) | 2011-02-16 |
CN101975575B CN101975575B (en) | 2012-03-28 |
Family
ID=43575476
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105072258A Expired - Fee Related CN101975575B (en) | 2010-10-15 | 2010-10-15 | Multi-target tracking method for passive sensor based on particle filtering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101975575B (en) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102338874A (en) * | 2011-06-24 | 2012-02-01 | 浙江大学 | Global probability data correlation method used for passive multi-sensor target tracking |
CN102621542A (en) * | 2012-04-02 | 2012-08-01 | 中国人民解放军海军航空工程学院 | Track method before locomotive weak target detection based on multimode grain filtering and data association |
CN103345258A (en) * | 2013-06-16 | 2013-10-09 | 西安科技大学 | Target tracking method and system of football robot |
CN103401625A (en) * | 2013-08-23 | 2013-11-20 | 西安电子科技大学 | Particle swarm optimization algorithm based cooperative spectrum sensing optimization method |
CN103471591A (en) * | 2013-04-15 | 2013-12-25 | 中国人民解放军海军航空工程学院 | Logical method, global nearest neighbor and target course information based maneuvering multi-target data interconnection algorithm |
CN103759732A (en) * | 2014-01-14 | 2014-04-30 | 北京航空航天大学 | Angle information assisted centralized multi-sensor multi-hypothesis tracking method |
CN103913432A (en) * | 2014-03-25 | 2014-07-09 | 西安交通大学 | Near infrared spectrum wavelength selecting method based on particle swarm optimization |
CN104237845A (en) * | 2014-09-19 | 2014-12-24 | 中国电子科技集团公司第二十八研究所 | Mean and variance analysis method for estimation errors of target positions of passive sensors |
CN104732559A (en) * | 2015-02-02 | 2015-06-24 | 大连民族学院 | Multi-target detecting and tracking method based on RGB-D data |
CN104750998A (en) * | 2015-04-09 | 2015-07-01 | 西安电子科技大学 | Passive multi-sensor target tracking method based on strength filter |
CN104977581A (en) * | 2015-07-15 | 2015-10-14 | 中国电子科技集团公司第三十六研究所 | Multi-moving target situation awareness method, device and system |
EP2975577A4 (en) * | 2013-04-15 | 2016-06-15 | Megachips Corp | State estimation device, program, and integrated circuit |
CN106291530A (en) * | 2016-07-26 | 2017-01-04 | 北京航空航天大学 | A kind of probabilistic data association optimization method based on nearest neighbor method |
CN106767832A (en) * | 2017-01-17 | 2017-05-31 | 哈尔滨工业大学 | A kind of passive multi-source multi-target tracking based on dynamic multidimensional distribution |
CN107202989A (en) * | 2017-05-08 | 2017-09-26 | 电子科技大学 | A kind of complicated Faint target detection and tracking suitable for passive starboard ambiguity of towed linear array sonar |
CN107329131A (en) * | 2017-08-08 | 2017-11-07 | 电子科技大学 | A kind of radar dim target detection tracking of utilization particle filter |
WO2018049602A1 (en) * | 2016-09-14 | 2018-03-22 | 深圳大学 | Multi-target tracking method and tracking system applicable to cluttered environment |
CN107917646A (en) * | 2017-01-10 | 2018-04-17 | 北京航空航天大学 | A kind of anti-interference method of guidance of strong pulsed D based on the prediction of target terminal accessoble region |
WO2018119912A1 (en) * | 2016-12-29 | 2018-07-05 | 深圳大学 | Target tracking method and device based on parallel fuzzy gaussian and particle filter |
CN108320302A (en) * | 2018-01-26 | 2018-07-24 | 西安电子科技大学 | CBMeMBer multi-object tracking methods based on random hypersurface |
CN110738275A (en) * | 2019-10-30 | 2020-01-31 | 中国人民解放军海军航空大学 | UT-PHD-based multi-sensor sequential fusion tracking method |
CN111107070A (en) * | 2019-12-10 | 2020-05-05 | 浙江工业大学 | Target tracking method under DoS attack |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101408983A (en) * | 2008-10-29 | 2009-04-15 | 南京邮电大学 | Multi-object tracking method based on particle filtering and movable contour model |
US20100090823A1 (en) * | 2008-10-09 | 2010-04-15 | Electronics And Telecommunications Research Institute | Hybrid clustering based data aggregation method for multi-target tracking in wireless sensor network |
CN101701826A (en) * | 2009-11-20 | 2010-05-05 | 西安电子科技大学 | Target tracking method of passive multi-sensor based on layered particle filtering |
CN101770024A (en) * | 2010-01-25 | 2010-07-07 | 上海交通大学 | Multi-target tracking method |
CN101783020A (en) * | 2010-03-04 | 2010-07-21 | 湖南大学 | Video multi-target fast tracking method based on joint probability data association |
CN101840579A (en) * | 2010-01-28 | 2010-09-22 | 浙江大学 | Method for realizing multi-target tracking by using video segmentation and particle filter |
-
2010
- 2010-10-15 CN CN2010105072258A patent/CN101975575B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100090823A1 (en) * | 2008-10-09 | 2010-04-15 | Electronics And Telecommunications Research Institute | Hybrid clustering based data aggregation method for multi-target tracking in wireless sensor network |
CN101408983A (en) * | 2008-10-29 | 2009-04-15 | 南京邮电大学 | Multi-object tracking method based on particle filtering and movable contour model |
CN101701826A (en) * | 2009-11-20 | 2010-05-05 | 西安电子科技大学 | Target tracking method of passive multi-sensor based on layered particle filtering |
CN101770024A (en) * | 2010-01-25 | 2010-07-07 | 上海交通大学 | Multi-target tracking method |
CN101840579A (en) * | 2010-01-28 | 2010-09-22 | 浙江大学 | Method for realizing multi-target tracking by using video segmentation and particle filter |
CN101783020A (en) * | 2010-03-04 | 2010-07-21 | 湖南大学 | Video multi-target fast tracking method based on joint probability data association |
Non-Patent Citations (1)
Title |
---|
《系统仿真学报》 20091031 赵欣等 一种基于粒子滤波的被动多传感器多目标跟踪算法 6368-6386,6390 第21卷, 第20期 2 * |
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102338874A (en) * | 2011-06-24 | 2012-02-01 | 浙江大学 | Global probability data correlation method used for passive multi-sensor target tracking |
CN102621542A (en) * | 2012-04-02 | 2012-08-01 | 中国人民解放军海军航空工程学院 | Track method before locomotive weak target detection based on multimode grain filtering and data association |
CN102621542B (en) * | 2012-04-02 | 2014-10-22 | 中国人民解放军海军航空工程学院 | Track method before locomotive weak target detection based on multimode grain filtering and data association |
US9911191B2 (en) | 2013-04-15 | 2018-03-06 | Megachips Corporation | State estimation apparatus, state estimation method, and integrated circuit with calculation of likelihood data and estimation of posterior probability distribution data |
CN103471591A (en) * | 2013-04-15 | 2013-12-25 | 中国人民解放军海军航空工程学院 | Logical method, global nearest neighbor and target course information based maneuvering multi-target data interconnection algorithm |
CN103471591B (en) * | 2013-04-15 | 2017-06-06 | 中国人民解放军海军航空工程学院 | The multiple-moving target data interconnection method of logic-based method, global arest neighbors and bogey heading information |
EP2975577A4 (en) * | 2013-04-15 | 2016-06-15 | Megachips Corp | State estimation device, program, and integrated circuit |
CN103345258A (en) * | 2013-06-16 | 2013-10-09 | 西安科技大学 | Target tracking method and system of football robot |
CN103345258B (en) * | 2013-06-16 | 2016-05-18 | 西安科技大学 | A kind of Soccer robot target tracking method and system |
CN103401625B (en) * | 2013-08-23 | 2015-04-08 | 西安电子科技大学 | Particle swarm optimization algorithm based cooperative spectrum sensing optimization method |
CN103401625A (en) * | 2013-08-23 | 2013-11-20 | 西安电子科技大学 | Particle swarm optimization algorithm based cooperative spectrum sensing optimization method |
CN103759732B (en) * | 2014-01-14 | 2016-06-22 | 北京航空航天大学 | A kind of centralized multisensor multiple hypotheis tracking method of angle information auxiliary |
CN103759732A (en) * | 2014-01-14 | 2014-04-30 | 北京航空航天大学 | Angle information assisted centralized multi-sensor multi-hypothesis tracking method |
CN103913432A (en) * | 2014-03-25 | 2014-07-09 | 西安交通大学 | Near infrared spectrum wavelength selecting method based on particle swarm optimization |
CN104237845A (en) * | 2014-09-19 | 2014-12-24 | 中国电子科技集团公司第二十八研究所 | Mean and variance analysis method for estimation errors of target positions of passive sensors |
CN104732559A (en) * | 2015-02-02 | 2015-06-24 | 大连民族学院 | Multi-target detecting and tracking method based on RGB-D data |
CN104732559B (en) * | 2015-02-02 | 2017-08-01 | 大连民族学院 | A kind of multi-target detection and tracking method based on RGB D data |
CN104750998B (en) * | 2015-04-09 | 2017-08-25 | 西安电子科技大学 | Target tracking method of passive multi-sensor based on density filter |
CN104750998A (en) * | 2015-04-09 | 2015-07-01 | 西安电子科技大学 | Passive multi-sensor target tracking method based on strength filter |
CN104977581A (en) * | 2015-07-15 | 2015-10-14 | 中国电子科技集团公司第三十六研究所 | Multi-moving target situation awareness method, device and system |
CN106291530B (en) * | 2016-07-26 | 2018-12-18 | 北京航空航天大学 | A kind of probabilistic data association optimization method based on nearest neighbor method |
CN106291530A (en) * | 2016-07-26 | 2017-01-04 | 北京航空航天大学 | A kind of probabilistic data association optimization method based on nearest neighbor method |
US10935653B2 (en) | 2016-09-14 | 2021-03-02 | Shenzhen University | Multi-target tracking method and tracking system applicable to clutter environment |
WO2018049602A1 (en) * | 2016-09-14 | 2018-03-22 | 深圳大学 | Multi-target tracking method and tracking system applicable to cluttered environment |
WO2018119912A1 (en) * | 2016-12-29 | 2018-07-05 | 深圳大学 | Target tracking method and device based on parallel fuzzy gaussian and particle filter |
CN107917646B (en) * | 2017-01-10 | 2020-11-24 | 北京航空航天大学 | Infrared air-to-air missile anti-interference guidance method based on target terminal reachable area prediction |
CN107917646A (en) * | 2017-01-10 | 2018-04-17 | 北京航空航天大学 | A kind of anti-interference method of guidance of strong pulsed D based on the prediction of target terminal accessoble region |
CN106767832B (en) * | 2017-01-17 | 2020-11-13 | 哈尔滨工业大学 | Passive multi-source multi-target tracking method based on dynamic multi-dimensional distribution |
CN106767832A (en) * | 2017-01-17 | 2017-05-31 | 哈尔滨工业大学 | A kind of passive multi-source multi-target tracking based on dynamic multidimensional distribution |
CN107202989A (en) * | 2017-05-08 | 2017-09-26 | 电子科技大学 | A kind of complicated Faint target detection and tracking suitable for passive starboard ambiguity of towed linear array sonar |
CN107329131B (en) * | 2017-08-08 | 2020-08-28 | 电子科技大学 | Radar weak target detection tracking method using particle filtering |
CN107329131A (en) * | 2017-08-08 | 2017-11-07 | 电子科技大学 | A kind of radar dim target detection tracking of utilization particle filter |
CN108320302A (en) * | 2018-01-26 | 2018-07-24 | 西安电子科技大学 | CBMeMBer multi-object tracking methods based on random hypersurface |
CN108320302B (en) * | 2018-01-26 | 2019-10-11 | 西安电子科技大学 | CBMeMBer multi-object tracking method based on random hypersurface |
CN110738275A (en) * | 2019-10-30 | 2020-01-31 | 中国人民解放军海军航空大学 | UT-PHD-based multi-sensor sequential fusion tracking method |
CN110738275B (en) * | 2019-10-30 | 2022-03-25 | 中国人民解放军海军航空大学 | UT-PHD-based multi-sensor sequential fusion tracking method |
CN111107070A (en) * | 2019-12-10 | 2020-05-05 | 浙江工业大学 | Target tracking method under DoS attack |
Also Published As
Publication number | Publication date |
---|---|
CN101975575B (en) | 2012-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101975575B (en) | Multi-target tracking method for passive sensor based on particle filtering | |
CN110503071B (en) | Multi-target tracking method based on variational Bayesian label multi-Bernoulli superposition model | |
CN111985093A (en) | Adaptive unscented Kalman filtering state estimation method with noise estimator | |
CN103644903B (en) | Synchronous superposition method based on the tasteless particle filter of distributed edge | |
CN103345577B (en) | Variation Bayesian probability assumed density multi-object tracking method | |
CN111722214B (en) | Method for realizing radar multi-target tracking PHD | |
CN111711432B (en) | Target tracking algorithm based on UKF and PF hybrid filtering | |
CN106772524A (en) | A kind of agricultural robot integrated navigation information fusion method based on order filtering | |
WO2018098926A1 (en) | Multi-target tracking method and system applicable to flicker noise | |
CN111693984B (en) | Improved EKF-UKF moving target tracking method | |
Havangi | Target tracking based on improved unscented particle filter with Markov chain Monte Carlo | |
CN103247057A (en) | Road target multi-hypothesis tracking algorithm under target-echo-road network data association | |
CN106054167A (en) | Intensity filter-based multi-extended target tracking method | |
CN108871365B (en) | State estimation method and system under course constraint | |
CN111259332B (en) | Fuzzy data association method and multi-target tracking method in clutter environment | |
CN116047495B (en) | State transformation fusion filtering tracking method for three-coordinate radar | |
CN107329131A (en) | A kind of radar dim target detection tracking of utilization particle filter | |
CN113030945B (en) | Phased array radar target tracking method based on linear sequential filtering | |
Hu et al. | Hybrid sampling-based particle filtering with temporal constraints | |
CN104467742A (en) | Sensor network distribution type consistency particle filter based on Gaussian mixture model | |
Wang et al. | Particle Filter with Hybrid Proposal Distribution for Nonlinear State Estimation. | |
Wang et al. | Improving particle filter with better proposal distribution for nonlinear filtering problems | |
CN108981707B (en) | Passive tracking multi-target method based on time difference measurement box particle PHD | |
CN102707278B (en) | Multi-target tracking method for singular value decomposition | |
CN114577213B (en) | Distributed multi-platform underwater multi-target association and passive positioning method |
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: 20120328 Termination date: 20151015 |
|
EXPY | Termination of patent right or utility model |