CN112630783A - Passive sonar target tracking method - Google Patents
Passive sonar target tracking method Download PDFInfo
- Publication number
- CN112630783A CN112630783A CN202011347206.3A CN202011347206A CN112630783A CN 112630783 A CN112630783 A CN 112630783A CN 202011347206 A CN202011347206 A CN 202011347206A CN 112630783 A CN112630783 A CN 112630783A
- Authority
- CN
- China
- Prior art keywords
- particle
- particles
- tracking
- target tracking
- target
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000002245 particle Substances 0.000 claims abstract description 86
- 238000001514 detection method Methods 0.000 claims abstract description 26
- 238000012952 Resampling Methods 0.000 claims abstract description 8
- 238000010606 normalization Methods 0.000 claims abstract description 7
- 230000004083 survival effect Effects 0.000 claims abstract description 7
- 230000037323 metabolic rate Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000006872 improvement Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 230000005855 radiation Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims 1
- 238000001914 filtration Methods 0.000 abstract description 14
- 230000008569 process Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/66—Sonar tracking systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/537—Counter-measures or counter-counter-measures, e.g. jamming, anti-jamming
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention relates to a sonar target tracking method, in particular to a passive sonar target tracking method, which comprises the following steps: initializing parameters, including determining the number of particles, the metabolic rate and the prior distribution of the particles, and setting the initial positions and states of the particles; calculating a particle space state set, a survival state set and a normalization weight value set; resampling is carried out; and outputting a state estimation result. The passive sonar target tracking method provided by the invention combines particle filtering and a pre-detection tracking algorithm and applies the combined algorithm to the target tracking of the passive sonar, thereby solving the problems of wrong tracking, lost tracking, false alarm and the like of the conventional algorithm in a complex low signal-to-noise ratio, strong interference and multi-target tracking environment, and improving the tracking performance of the equipment.
Description
Technical Field
The invention relates to a sonar target tracking method, in particular to a passive sonar target tracking method.
Background
The passive sonar mainly depends on target radiation signal noise to detect underwater targets on water, and compared with active sonar, the passive sonar has better secrecy and anti-interference performance. However, due to the complexity and specificity of the underwater environment and the specificity of the passive sonar, the sonar is difficult to detect weak targets, and target tracking of the passive sonar becomes a very challenging task.
The DBT method comprises the steps of firstly detecting and then tracking, wherein the DBT method is a traditional target tracking mode, but in an underwater environment with a low signal-to-noise ratio, the detection rate can be greatly reduced, and the traditional DBT method is greatly restricted. And the TBD effectively combines the detection and the tracking without being restricted by a detection result, directly processes the detected previous layer or even more original data, and avoids information loss caused by the data in the post-processing process. On the other hand, "track before detection, TBD" can jointly process multi-frame data, and the influence of signal-to-noise ratio is reduced by using the continuity of target motion in adjacent time periods, so that the tracking performance is improved.
Kalman Filtering (KF) is a tracking filtering method that is currently very widely used. However, this filter is limited by a linear system and is difficult to implement in practical engineering. Extended Kalman Filter (EKF) and traceless Kalman Filter (UKF) improved by Kalman Filter break through the restriction of linear system, but these methods are all premised on gaussian disturbance and do not necessarily meet the actual tracking condition. With the continuous development of computer technology, another well-known nonlinear Filter, Particle Filter (PF), has received high attention from researchers at home and abroad. The Monte Carlo Method (MCM) proposed by MERROPOLIS et al in the 40's of the 20 th century is an rudimentary form of the particle filtering algorithm, which was later modified by Hammersley et al, which has a particle degradation problem. Aiming at the problem of particle degradation, Gordon et al in the last 90 th century proposed a resampling method, which greatly improved the problem of particle degradation. The particle filter is not limited by Gaussian disturbance and a linear system, and is more widely applied.
Based on the background, the invention improves the pre-detection tracking algorithm based on particle filtering and uses the pre-detection tracking algorithm in the target tracking of the passive sonar. The conventional passive sonar target tracking algorithm can better track a low-interference and strong target. However, the conventional algorithm is weak in complex environments such as multiple targets, strong interference, low signal-to-noise ratio and the like. In general, a conventional passive sonar target tracking algorithm detects a target first and then tracks a detection result. However, part of the target information is lost in the detection process, and the subsequent tracking performance is reduced. And the tracking algorithm before detection does not carry out target detection before tracking, thereby reducing target information loss and improving target tracking performance.
Disclosure of Invention
In order to solve the problems, the invention provides a passive sonar target tracking method which combines particle filtering and a pre-detection tracking algorithm and applies the combined algorithm to the target tracking of passive sonar, solves the problems of complex low signal-to-noise ratio, strong interference, frequent false tracking, tracking loss, false alarm and the like of the conventional algorithm in a multi-target tracking environment and has high tracking performance, and the specific technical scheme is as follows:
a passive sonar target tracking method comprises the following steps: initializing parameters, including determining the number of particles, the metabolic rate and the prior distribution of the particles, and setting the initial positions and states of the particles; calculating a particle space state set, a survival state set and a normalization weight value set; resampling is carried out; and outputting a state estimation result.
Further, S1, initialization parameter: setting the number N of particles, a state transition matrix PI and the particle initial rate mu m; according to the priori knowledge, generating N x m particles, and setting the corresponding states as: e0 iTable 1 is survival; the remaining particle states are set as: e0 iDeath is indicated by 0; state transition matrix of particles: PI ═ 1-Pb,Pb;Pd,1-Pd](ii) a Wherein, PnRepresenting the birth rate, PdRepresenting mortality; dividing the observation region into n observation units each having a length of deltaL(ii) a S2, transferring the particle state surviving at the previous moment according to the state transfer matrix PI to form the particle state set at the current moment, aiming at the particle state set at the previous momentParticle state Ek-1 iGenerating a random number mu between (0-1); if E isk-1 i0 and mu<PbWhen E is greaterk i1 is ═ 1; otherwise E k i0; if E isk-1 i1 and mu<PdWhen E is greaterk i0; otherwise Ek i1 is ═ 1; s3, if E k-1 i0 and Ek iIf the particle number is 1, sampling to generate a new particle according to the prior distribution probability; if E isk-1 i1 and Ek iIf the value is 1, a new particle is generated by using a state transfer function F, then the weight value is taken as a weight value according to the likelihood function ratio, and normalization processing is carried out to generate the weight value of the particle: { omega [ [ omega ] ]k iI ═ 1 … N }; determining a likelihood function ratio model: before detection, the target tracking measurement value is the signal energy intensity of the sensor in each direction; let n resolution elements, which may be denoted as j, { j ═ 1 … n }, per observation region, and each resolution element has a length ΔLThen the measurement value received at time k is WhereinThe signal strength of the resolution unit at the moment j is expressed by the expression:
in the formulaTo target at thetakSignal strength of spotted object, noisekCorresponding observation noise; when a target appears in an observation unit, the signal energy of the observation unit is increased and the observation unit is close to the observation unitThe signal energy of the element is also improved, and the improvement degree of the element and the improvement degree of the observation unit present an inverse relation; structure of the device
Where epsilon represents the degree of scattering of the sensor energy. Obtaining a likelihood function:
s4, resampling: generating a new particle space position by taking the space position and the weight of each particle as input values, and setting the weight of all the particles as 1/N;
s5, outputting the result:
furthermore, a plurality of sonars are adopted during tracking, the sonars are arranged into a sonar array along a linear array, the sonar array performs signal preprocessing on the received radiation sound signals, and beam domain data after beam forming is output.
Furthermore, the tracking mode is a tracking mode before detection, and a particle filter nonlinear filter is adopted.
Compared with the prior art, the invention has the following beneficial effects:
the passive sonar target tracking method provided by the invention combines particle filtering and a pre-detection tracking algorithm and applies the combined algorithm to the target tracking of the passive sonar, thereby solving the problems of wrong tracking, lost tracking, false alarm and the like of the conventional algorithm in a complex low signal-to-noise ratio, strong interference and multi-target tracking environment, and improving the tracking performance of the equipment.
Drawings
FIG. 1 is a comparison graph of simulation data of a target real track and three algorithms under different algorithm conditions when SNR is 21;
FIG. 2 is a comparison graph of simulation data of a target real track and three algorithms under different algorithm conditions when the SNR is 9;
FIG. 3 is a comparison graph of simulation data of a target real track and three algorithms under different algorithm conditions when the SNR is 1.5;
FIG. 4 is a first chart of a sea trial bearing;
fig. 5 is a second chart of the azimuth history of a certain sea test.
Detailed Description
The invention will now be further described with reference to the accompanying drawings.
As shown in fig. 1 to 5, a passive sonar target tracking method includes the following steps:
initializing parameters, including determining the number of particles, the metabolic rate and the prior distribution of the particles, and setting the initial positions and states of the particles;
calculating a particle space state set, a survival state set and a normalization weight value set;
resampling is carried out;
and outputting a state estimation result.
S1, initialization parameters:
setting the number N of particles, a state transition matrix PI and the particle initial rate mu m;
according to the priori knowledge, generating N x m particles, and setting the corresponding states as:
E0 itable 1 is survival;
the remaining particle states are set as:
E0 ideath is indicated by 0;
state transition matrix of particles:
PI=[1-Pb,Pb;Pd,1-Pd];
wherein, PbRepresenting the birth rate, PdRepresenting mortality;
dividing the observation region into n observationsUnits each unit having a length of deltaL;
S2, transferring the particle state surviving at the previous moment according to the state transfer matrix PI to form a particle state set at the current moment, and aiming at the particle state E at the previous momentk-1 iGenerating a random number mu between (0-1);
if E isk-1 i0 and mu<PbWhen E is greaterk i1 is ═ 1; otherwise Ek i=0;
If E isk-1 i1 and mu<PdWhen E is greaterk i0; otherwise Ek i=1;
S3, if E k-1 i0 and Ek iIf the particle number is 1, sampling to generate a new particle according to the prior distribution probability;
if E isk-1 i1 and Ek iIf the value is 1, a new particle is generated by using a state transfer function F, then the weight value is taken as a weight value according to the likelihood function ratio, and normalization processing is carried out to generate the weight value of the particle:
{ωk i,i=1…N};
determining a likelihood function ratio model:
before detection, the target tracking measurement value is the signal energy intensity of the sensor in each direction; let n resolution elements, which may be denoted as j, { j ═ 1 … n }, per observation region, and each resolution element has a length ΔLThen the measurement value received at time k isWhereinThe signal strength of the resolution unit at the moment j is expressed by the expression:
in the formulaTo target at thetakSignal strength of spotted object, noisekCorresponding observation noise;
when a target appears in an observation unit, the signal energy of the observation unit is improved, and the signal energy close to the observation unit is also improved, wherein the improvement degree of the observation unit and the improvement degree of the observation unit are in an inverse relation; structure of the device
Where epsilon represents the degree of scattering of the sensor energy. Obtaining a likelihood function:
s4, resampling: generating a new particle space position by taking the space position and the weight of each particle as input values, and setting the weight of all the particles as 1/N;
s5, outputting the result:
the state model of the object motion is:
Xk=F*(Xk-1,Vk-1)
wherein Xk=(θ,Δθ,Ik) Theta is azimuth, Delta theta is azimuth rate of change, IkTarget power spectrum information, F ═ 1, T, 0; 0,1, 0; 0,0,1]T is the time interval, Vk-1Is the system process noise.
During tracking, a plurality of sonars are adopted and arranged into a sonar array along a linear array, and the sonar array performs signal preprocessing on received radiation acoustic signals and outputs beam domain data after beam forming.
The tracking mode is a tracking mode before detection, and a particle filter nonlinear filter is adopted.
1. And (3) algorithm simulation verification:
the simulation parameters are as follows: tracking total time duration Ts 60s, time interval T1 s, number of particles 2000, particle birth rate PbParticle mortality rate P ═ 0.05dThe target moving range angle is 0 to 180, the target initial azimuth is 20, the initial azimuth change speed is 0.45, the change range is 0.2 to 2, the target energy initial size is 0.5, the change range is 0.2 to 1, and the signal-to-noise ratios are respectively set to 1.5db, 9db and 21 db.
As shown in fig. 1 to fig. 3, the results of comparing the tracking performance of the particle filter algorithm and the tracking performance of the kalman filter algorithm and the particle filter algorithm before detection based on the particle filter under different signal-to-noise ratios are shown. When the signal-to-noise ratio is 21db, the three algorithms can track the target more effectively. But from an enlarged view around 30s, it can be seen that the particle filter based pre-detection tracking algorithm (PF-TBD) is closer to the true trajectory. As the signal-to-noise ratio decreases, the KF algorithm fails to track to the target in the initial stage until convergence to the target position around 25s and a large deviation occurs at 48s, as in the 9db condition shown in fig. 2. The PF is basically able to track the upper target, although the error is slightly larger within the initial 20 s. The PF-TBD algorithm accurately tracks the target at all times except when it fails to find the target within the first 3s, and the error of the algorithm is observed to be minimal. As shown in fig. 3, when the signal-to-noise ratio is reduced to 1.5db, neither the PF algorithm nor the KF algorithm can track the target, while the PF-TBD algorithm, although jumping at the initial stage, generally completes the task of target tracking. Therefore, it can be found from fig. 1 that the tracking performance of the three algorithms is reduced with the reduction of the signal-to-noise ratio, but under the condition of low signal-to-noise ratio, the PF-TBD algorithm is obviously superior to the PF algorithm and the KF algorithm. This demonstrates the superiority of PF-TBD at low signal to noise ratios.
2. Verification of test data
As shown in fig. 4 and 5, the diagram is a direction history diagram of a certain sea test, and it can be seen from the diagram that the tracking algorithm before detection based on particle filtering can stably track weak targets and cross targets under strong interference.
The passive sonar target tracking algorithm before detection based on particle filtering is based on particle filtering, and the conventional thought of tracking after detection is abandoned, and beam domain data is directly processed, so that the problems of target information loss and the like in the detection process are avoided. In addition, particle filtering, which is a nonlinear filter that is receiving attention in recent years, breaks through the restriction of using kalman filtering under a linear system, and the application scenarios are wider.
The method is mainly used for target tracking of the passive sonar, and can greatly improve the target tracking performance compared with the conventional algorithm. This algorithm has three advantages:
1. multi-frame data joint processing can be realized, and the target tracking precision under the condition of low signal-to-noise ratio is improved;
2. the detection is not needed, the beam domain data is directly processed, and the information loss caused by multi-layer processing is reduced;
3. compared with Kalman filtering, particle filtering is more widely applied and is not limited by a linear system.
The technical principle of the present invention is described above in connection with specific embodiments. The description is made for the purpose of illustrating the principles of the invention and should not be construed in any way as limiting the scope of the invention. Based on the explanations herein, those skilled in the art will be able to conceive of other embodiments of the present invention without inventive step, which shall fall within the scope of the appended claims.
Claims (4)
1. A passive sonar target tracking method is characterized by comprising the following steps:
initializing parameters, including determining the number of particles, the metabolic rate and the prior distribution of the particles, and setting the initial positions and states of the particles;
calculating a particle space state set, a survival state set and a normalization weight value set;
resampling is carried out;
and outputting a state estimation result.
2. The passive sonar target tracking method according to claim 1,
s1, initialization parameters:
setting the number N of particles, a state transition matrix PI and the particle initial rate mu m;
according to the priori knowledge, generating N x m particles, and setting the corresponding states as:
E0 itable 1 is survival;
the remaining particle states are set as:
E0 ideath is indicated by 0;
state transition matrix of particles:
PI=[1-Pb,Pb;Pd,1-Pd];
wherein, PbRepresenting the birth rate, PdRepresenting mortality;
dividing the observation region into n observation units each having a length of deltaL;
S2, transferring the particle state surviving at the previous moment according to the state transfer matrix PI to form a particle state set at the current moment, and aiming at the particle state E at the previous momentk-1 iGenerating a random number mu between (0-1);
if E isk-1 i0 and mu < PbWhen E is greaterk i1 is ═ 1; otherwise Ek i=0;
If E isk-1 i1 and mu < PdWhen E is greaterk i0; otherwise Ek i=1;
S3, if Ek-1 i0 and Ek iIf the particle number is 1, sampling to generate a new particle according to the prior distribution probability;
if E isk-1 i1 and Ek iIf the value is 1, a new particle is generated by using a state transfer function F, then the weight value is taken as a weight value according to the likelihood function ratio, and normalization processing is carried out to generate the weight value of the particle:
{ωk i,i=1…N};
determining a likelihood function ratio model:
before detection, the target tracking measurement value is the signal energy intensity of the sensor in each direction; let n resolution elements, which may be denoted as j, { j ═ 1 … n }, per observation region, and each resolution element has a length ΔLThen the measurement value received at time k is Zk={Ik jJ ═ 1 … n }, where Ik jThe signal strength of the resolution unit at the moment j is expressed by the expression:
in the formula hk j(θk) To target at thetakSignal strength of spotted object, noisekCorresponding observation noise;
when a target appears in an observation unit, the signal energy of the observation unit is improved, and the signal energy close to the observation unit is also improved, wherein the improvement degree of the observation unit and the improvement degree of the observation unit are in an inverse relation; construction hk j(θk):
Where epsilon represents the degree of scattering of the sensor energy. Obtaining a likelihood function:
s4, resampling: generating a new particle space position by taking the space position and the weight of each particle as input values, and setting the weight of all the particles as 1/N;
s5, outputting the result:
3. the passive sonar target tracking method according to claim 1,
during tracking, a plurality of sonars are adopted and arranged into a sonar array along a linear array, and the sonar array performs signal preprocessing on received radiation acoustic signals and outputs beam domain data after beam forming.
4. The passive sonar target tracking method according to claim 1,
the tracking mode is a tracking mode before detection, and a particle filter nonlinear filter is adopted.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011347206.3A CN112630783A (en) | 2020-11-26 | 2020-11-26 | Passive sonar target tracking method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011347206.3A CN112630783A (en) | 2020-11-26 | 2020-11-26 | Passive sonar target tracking method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112630783A true CN112630783A (en) | 2021-04-09 |
Family
ID=75304147
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011347206.3A Pending CN112630783A (en) | 2020-11-26 | 2020-11-26 | Passive sonar target tracking method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112630783A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113534164A (en) * | 2021-05-24 | 2021-10-22 | 中船海洋探测技术研究院有限公司 | Target path tracking method based on active and passive combined sonar array |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN109061615A (en) * | 2018-10-26 | 2018-12-21 | 海鹰企业集团有限责任公司 | The Target moving parameter estimation method and device of nonlinear system in passive sonar |
CN110673119A (en) * | 2019-10-30 | 2020-01-10 | 海鹰企业集团有限责任公司 | Compressed sensing-based non-regularization orientation estimation method and system |
CN111948657A (en) * | 2020-07-28 | 2020-11-17 | 哈尔滨工程大学 | Maneuvering weak target tracking-before-detection method based on multimode particle filtering |
-
2020
- 2020-11-26 CN CN202011347206.3A patent/CN112630783A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN109061615A (en) * | 2018-10-26 | 2018-12-21 | 海鹰企业集团有限责任公司 | The Target moving parameter estimation method and device of nonlinear system in passive sonar |
CN110673119A (en) * | 2019-10-30 | 2020-01-10 | 海鹰企业集团有限责任公司 | Compressed sensing-based non-regularization orientation estimation method and system |
CN111948657A (en) * | 2020-07-28 | 2020-11-17 | 哈尔滨工程大学 | Maneuvering weak target tracking-before-detection method based on multimode particle filtering |
Non-Patent Citations (2)
Title |
---|
M.G. RUTTEN等: "A comparison of particle filters for recursive track-before-detect", 2005 7TH INTERNATIONAL CONFERENCE ON INFORMATION FUSION * |
徐璐霄: "基于被动阵列的水声弱目标检测前跟踪算法研究", 中国优秀硕士学位论文全文数据库工程科技Ⅱ辑, pages 6 - 8 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113534164A (en) * | 2021-05-24 | 2021-10-22 | 中船海洋探测技术研究院有限公司 | Target path tracking method based on active and passive combined sonar array |
CN113534164B (en) * | 2021-05-24 | 2023-12-12 | 中船海洋探测技术研究院有限公司 | Target path tracking method based on active-passive combined sonar array |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103885057B (en) | Adaptive strain sliding window multi-object tracking method | |
CN103439697B (en) | Target detection method based on dynamic programming | |
CN107202989B (en) | Complex weak target detection and tracking method suitable for passive towed linear array sonar | |
CN107356922B (en) | Multi-frequency radar target tracking method for ambiguity resolution by spread-table-aided method | |
CN106019256B (en) | Radar signal self-adapting detecting method based on autoregression model | |
CN104101876B (en) | Random finite set based multi-target tracking method in outer radiation source radar | |
CN105676181A (en) | Underwater moving target extended Kalman filtering tracking method based on distributed sensor energy ratios | |
CN109901154B (en) | Self-adaptive adjustment method based on recursive RTHT-TBD | |
CN112630783A (en) | Passive sonar target tracking method | |
CN113534133A (en) | Radar multi-target joint detection and tracking method based on expectation maximum iteration algorithm | |
CN113484866B (en) | Multi-target detection tracking method based on passive sonar azimuth history map | |
CN115201799A (en) | Time-varying Kalman filtering tracking method for sonar | |
CN107561506B (en) | Multistage false alarm probability joint recursion target detection method | |
CN116520311A (en) | GLMB-based adaptive track initiation method | |
CN106033120A (en) | Multistation radar asynchronous multiframe joint detection method | |
CN112213706B (en) | Logic TBD detection method for multi-dimensional digital wave gate interframe progressive association | |
CN112114286B (en) | Multi-target tracking method based on line spectrum life cycle and single-vector hydrophone | |
CN109901152B (en) | Asynchronous moving platform multi-base sensor system target detection tracking joint estimation method | |
Blanding et al. | Sequential ML for multistatic sonar tracking | |
CN111722231B (en) | Active sonar target motion information fusion method | |
CN114114170B (en) | Mathematical modeling method for radar signal processor target measurement | |
CN113189550B (en) | Array radar anti-noise interference method and system based on multi-unknown signal extraction | |
CN115469270B (en) | Long baseline positioning method and system based on pre-detection tracking | |
CN111025251B (en) | Multi-target composite detection method based on dynamic programming | |
CN114966710A (en) | Underwater weak target tracking method and device based on multi-Bernoulli RFS algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |