CN112684454B - Track cross target association method based on sub-frequency bands - Google Patents
Track cross target association method based on sub-frequency bands Download PDFInfo
- Publication number
- CN112684454B CN112684454B CN202011402804.6A CN202011402804A CN112684454B CN 112684454 B CN112684454 B CN 112684454B CN 202011402804 A CN202011402804 A CN 202011402804A CN 112684454 B CN112684454 B CN 112684454B
- Authority
- CN
- China
- Prior art keywords
- sub
- target
- band
- frequency
- bands
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 25
- 238000012216 screening Methods 0.000 claims abstract description 13
- 238000001228 spectrum Methods 0.000 claims abstract description 10
- 238000001514 detection method Methods 0.000 claims description 8
- 230000005574 cross-species transmission Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a track cross target association method based on sub-bands, which comprises the steps of firstly decomposing the whole bandwidth of a broadband into a plurality of sub-bands, then respectively carrying out beam forming processing on each sub-band, then screening out target characteristic sub-bands by utilizing acoustic energy information and continuous spectrum information of different sub-bands of each target, and carrying out track association before and after target cross by utilizing the characteristic sub-band information. The invention has the beneficial effects that: the method aims to comprehensively utilize angle information, energy information and spectrum characteristics of the target to solve the problem that the target loses or is mistakenly followed after track crossing. And screening the characteristic sub-frequency bands of the target by using the acoustic energy information and the continuous spectrum information of the target, and performing track association before and after cross tracking by using the characteristic sub-frequency band information, thereby solving the problem of cross tracking loss of the multi-target tracking track.
Description
Technical Field
The invention relates to the field of passive target cross tracking in sonar array signal processing, in particular to a track cross target association method based on a sub-frequency band.
Background
If the target is out of tracking or mistracking, the realization of the functions of sonar system such as target positioning, movement element settlement, feature analysis and identification, listening and the like will not be mentioned.
The key point of effectively tracking the track crossed target is mainly the correlation method of the tracks before and after crossing and the signal to noise ratio of the weak target in the track is improved as much as possible. If the weak target can be successfully detected and the signal-to-noise ratio of the weak target is higher after detection, the probability of successful association with the track is high; if the trajectories can be correctly correlated, the target is also tracked smoothly. Therefore, how to improve the signal-to-noise ratio of weak targets in the crossed targets and how to improve the correlation accuracy before and after the track crossing are particularly important. For passive sonar, a method with strong target detection capability is found, and the accuracy of association is improved, so that part or all of the problems can be solved.
The conventional detection tracking processing of the existing passive sonar mainly has the following two problems:
conventional passive detection processing: the entire frequency band is typically divided into several frequency bands. Once the design is complete, the bin width is fixed. In practice, however, there are often limited frequency bands (but not necessarily narrow bands) for weak targets. If the processing frequency band is not consistent with the target frequency band, either the signal energy is not fully utilized (the processing frequency band is smaller than the target frequency band) or more noise is introduced to cause the detection capability to be reduced (the processing frequency band is larger than the target frequency band).
Conventional passive tracking: when multi-target tracking is carried out, strong targets interfere weak targets to influence tracking of the weak targets; when the target tracks are crossed, the conventional algorithm only uses the angle of the target for tracking, and the track is often associated with a measurement error, so that the target is missed or mistakenly tracked.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a track cross target association method based on a sub-frequency band.
The purpose of the invention is achieved by the following technical scheme. A track cross target association method based on sub-bands comprises the steps of firstly decomposing the whole bandwidth of a broadband into a plurality of sub-bands, then respectively carrying out beam forming processing on each sub-band, then screening out target characteristic sub-bands by utilizing acoustic energy information and continuous spectrum information of different sub-bands of each target, and carrying out track association before and after target cross by utilizing characteristic sub-band information, so that the problem of multi-target tracking track cross tracking loss is solved.
The beneficial effects of the invention are as follows: the method aims to comprehensively utilize angle information, energy information and spectrum characteristics of the target to solve the problem that the target loses or is mistakenly followed after track crossing. And screening the characteristic sub-frequency bands of the target by using the acoustic energy information and the continuous spectrum information of the target, and performing track association before and after cross tracking by using the characteristic sub-frequency band information, thereby solving the problem of cross tracking loss of the multi-target tracking track.
Drawings
FIG. 1 is a schematic cross-track diagram;
FIG. 2 is an azimuth history map;
FIG. 3 is a schematic diagram of target sub-band energy spillover;
the target of fig. 4 is the energy fraction in the different sub-bands.
Fig. 5 shows a azimuth history of sub-band 4.
Detailed Description
The invention will be described in detail below with reference to the following drawings:
the invention discloses a track crossing target association method based on sub-bands, which comprises the steps of firstly decomposing the whole bandwidth of a broadband into a plurality of sub-bands, then respectively carrying out beam forming processing on each sub-band, then screening out target characteristic sub-bands by utilizing sound energy information and continuous spectrum information of different sub-bands of each target, and carrying out track association before and after target crossing by utilizing the characteristic sub-band information, thereby solving the problem of multi-target tracking track crossing and tracking loss. The algorithm principle is as follows:
(1) Sub-band spectral formation
Suppose that the time snapshot received by the sonar array element is x (n) = [ x = 1 (n),…,x M (n)]FFT conversion thereof to frequency snapshot X (f) = [ X = [) 1 (f),…,X M (f)]Wherein n is a time snapshot point, f is an FFT frequency point, and M is the number of array elements. The beam output obtained after frequency domain beamforming is B (f, θ) i ):
B(f,θ)=w(f,θ)·X(f) (1)
Where θ is the beam number and w (f, θ) is its corresponding response vector.
Suppose the entire processing band is f min f max ]Is divided into L sub-bands, the first sub-band having a frequency range of [ f [ ] l.min f l.max ]. Performing beam forming on each sub-band to obtain beam output of each sub-band, i.e.
For each B l (f, theta) square detection is carried out to finally obtain the output P of the first sub-band at the time t l (θ,t)。
All sub-frequency bands at the time t are synthesized to obtain the output P of the whole processing frequency band FULL (θ, t), i.e.
(2) Characteristic sub-band screening
In the passive target detection and tracking process, the frequency spectrum characteristics of the target are not continuously stable in the whole processing frequency band. But the spectral characteristics of the target in certain sub-bands, which are defined as characteristic sub-bands, may be considered to remain stationary for a limited period of time. In a typical sub-frequency, the target may achieve a greater signal-to-noise ratio than the entire processing bin. The characteristic sub-band screening method comprises the following steps:
normalizing each sub-frequency band and the synthesized result respectively to obtain
P l '(θ,t)=P l (θ,t)/max(P l (θ,t)) (4)
P' FULL (θ,t)=P FULL (θ,t)/max(P FULL (θ,t)) (5)
Defining the energy overflow of the ith target in the ith sub-band relative to the whole band as
Indicating an energy spill over for a certain period of time,energy spillover representing accumulation time, where the kth time period is denoted as [ t ] i,start t i,end ]Length of N t The time length can be set according to the actual condition, theta i Indicates that the i-th object is at t i The beam number of the time instant. They measure the relative strength of a certain sub-band of the target in the whole processing band, and a larger value indicates that the energy of the target in the sub-band is higher in the whole processing band.
the energy ratio of each sub-band wave beam theta at the time t in the whole processing frequency band is
Calculating the average energy ratio of the target in each sub-frequency band in a certain time period as
Wherein,represents the average energy ratio of the i number target in the kth sub-band in the kth time period. The cumulative energy ratio of the target number i in the ith sub-band is expressed as
Defining the energy intensity difference of the target I and the target j in the l sub-band asNamely, it is
When a certain sub-band satisfies the following condition
This sub-band is called the characteristic sub-band, where E 0 、Δ 0 And (4) selecting. Assuming that L' characteristic sub-bands are obtained by sub-band screening, the target cross-tracking process is performed in these characteristic sub-bands.
(3) Sub-band target cross-tracking processing
Assume that the time period before the intersection of two target tracks is k before The time period after the track crossing is k after The average energy ratio of the two target tracking front and back in the characteristic sub-frequency band is respectively Calculating the cross-correlation coefficient of the target before and after tracking:
wherein-1 is not more than rho ij ≤1,ρ ii =1。
Respectively calculating rho according to the cross track schematic diagram 13 、ρ 24 、ρ 14 、ρ 23 Defining the correlation coefficient as eta
η=(ρ 13 +ρ 24 )-(ρ 14 +ρ 23 ) (17)
When in use
Wherein eta 0 Greater than or equal to 0, general eta 0 And taking zero. The greater the value of η, the greater the confidence of the association. And finishing track cross tracking association.
An even linear array is assumed, the number of array elements is 128, the spacing between the array elements is 0.1m, the working frequency range is 1000Hz-4000Hz, and the sampling frequency is 15kHz. For 3 existing targets, the frequency of a target 1 is 1000Hz-3800Hz, the signal-to-noise ratio is-5 dB, the initial angle is 16 degrees, and the end angle is-11 degrees; the frequency of the target 2 is 1000Hz-3800Hz, the signal-to-noise ratio is-16 dB, the starting angle is 7 degrees, and the ending angle is 7 degrees; the frequency of target 3 is 1000Hz-4000Hz, the signal-to-noise ratio is 0dB, the starting angle is-45 degrees, and the ending angle is-63 degrees. Where target 1 and target 2 begin to intersect at t =23 and end at t =47, the bearing history is shown in fig. 2.
The process divides the whole band into 8 sub-bands and the cumulative energy spill over for each sub-band for target 1 and target 2 is shown in figure 3.
The different sub-band energies before the target 1 and target 2 trajectories cross are shown in fig. 4.
According to fig. 3 and 4, the characteristic sub-band is selected as [1 3 45 ]. The azimuth course of the target in the 4 th sub-band is shown in FIG. 5
Fig. 5 compares with fig. 2, and the signal-to-noise ratio of the weak target 2 in the sub-band 4 is significantly higher than that in the full band. Firstly, calculating the correlation coefficient before and after the track intersection of the target 1 and the target 2 of 8 sub-frequency bands:
and performing association processing by using the characteristic sub-frequency band to obtain an association coefficient:
it can be seen from the above that, the cross processing performed by using 8 sub-bands and the characteristic sub-band can be associated successfully, but the association reliability can be improved by using the characteristic sub-band.
The calculation process is as follows:
(1) Performing beam forming on the L sub-bands at the time t to obtain a sub-band spectrum P l (theta, t) and full band beam output P FULL (θ,t);
(2) Calculating the target before track crossingAndscreening to obtain a characteristic sub-frequency band;
(3) And calculating the correlation coefficient eta before and after the target track is crossed by utilizing the characteristic sub-frequency bands to carry out correlation judgment.
It should be understood that equivalent substitutions and changes to the technical solution and the inventive concept of the present invention should be made by those skilled in the art to the protection scope of the appended claims.
Claims (1)
1. A track cross target association method based on sub-frequency bands is characterized in that: firstly, decomposing the whole bandwidth of a broadband into a plurality of sub-frequency bands, then respectively carrying out beam forming processing on each sub-frequency band, then screening out target characteristic sub-frequency bands by utilizing acoustic energy information and continuous spectrum information of different sub-frequency bands of each target, and carrying out track association before and after target crossing by utilizing the characteristic sub-frequency band information so as to solve the problem of multi-target tracking track crossing and tracking loss; the method specifically comprises the following steps:
(1) Sub-band spectral formation
Suppose that the time snapshot received by the sonar array element is x (n) = [ x = 1 (n),…,x M (n)]FFT conversion thereof to frequency snapshot X (f) = [ X = [) 1 (f),…,X M (f)]Wherein n is a time snapshot point, f is an FFT frequency point, M is the number of array elements, and the output of the wave beam obtained after the wave beam of the frequency domain is formed is B (f, theta) i ):
B(f,θ)=w(f,θ)·X(f) (1)
Wherein, theta is the beam number, and w (f, theta) is the corresponding response vector;
suppose the entire processing band is f min f max ]Is divided into L sub-bands, the first sub-band having a frequency range of [ f [ ] l.min f l.max ]Performing beam forming on each sub-band to obtain beam output of each sub-band, i.e. beam forming
For each B l (f, theta) square detection is carried out to finally obtain the output P of the first sub-band at the time t l (θ,t);
All sub-frequency bands at the time t are synthesized to obtain the output P of the whole processing frequency band FULL (θ, t), i.e.
(2) Characteristic sub-band screening
The characteristic sub-band screening method comprises the following steps:
normalizing each sub-frequency band and the synthesized result respectively to obtain
P l '(θ,t)=P l (θ,t)/max(P l (θ,t)) (4)
P' FULL (θ,t)=P FULL (θ,t)/max(P FULL (θ,t)) (5)
Defining the energy overflow of the ith target in the ith sub-band relative to the whole band as
Indicating an energy spill over for a certain period of time,energy spillover representing accumulation time, where the kth time period is denoted as [ t ] i,start t i,end ]Length of N t The time length can be set according to the actual condition, theta i Indicates that the i-th object is at t i A beam number of a time;
the energy ratio of each sub-band wave beam theta at the time t in the whole processing frequency band is
Calculating the average energy ratio of the target in each sub-frequency band in a certain time period as
Wherein,the average energy ratio of the i number of targets in the ith sub-band in the kth time period is represented, and the accumulated energy ratio of the i number of targets in the ith sub-band is represented as
defining the energy intensity difference of the target I and the target j in the l sub-band asNamely, it is
When a certain sub-band satisfies the following condition
This sub-band is called the characteristic sub-band, where E 0 、Δ 0 Selecting; assuming that L' characteristic sub-bands are obtained through sub-band screening, the target cross-tracking processing is carried out on the characteristic sub-bandsCarrying out the steps of (1);
(3) Sub-band target cross-tracking processing
Assume that the time period before the intersection of two target tracks is k before The time period after the track crossing is k after The average energy ratio of the two target tracking front and back in the characteristic sub-frequency band is respectively Calculating the cross-correlation coefficient of the target before and after tracking:
wherein-1 ≦ ρ ij ≤1,ρ ii =1;
Respectively calculating rho according to the cross track schematic diagram 13 、ρ 24 、ρ 14 、ρ 23 Defining the correlation coefficient as eta
η=(ρ 13 +ρ 24 )-(ρ 14 +ρ 23 ) (17)
When in use
Wherein eta 0 And the greater the eta value is, the greater the association confidence coefficient is, and the track cross tracking association is completed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011402804.6A CN112684454B (en) | 2020-12-04 | 2020-12-04 | Track cross target association method based on sub-frequency bands |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011402804.6A CN112684454B (en) | 2020-12-04 | 2020-12-04 | Track cross target association method based on sub-frequency bands |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112684454A CN112684454A (en) | 2021-04-20 |
CN112684454B true CN112684454B (en) | 2022-12-06 |
Family
ID=75447261
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011402804.6A Active CN112684454B (en) | 2020-12-04 | 2020-12-04 | Track cross target association method based on sub-frequency bands |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112684454B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113239775B (en) * | 2021-05-09 | 2023-05-02 | 西北工业大学 | Method for detecting and extracting tracks in azimuth lineage diagram based on hierarchical attention depth convolution neural network |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5119432A (en) * | 1990-11-09 | 1992-06-02 | Visidyne, Inc. | Frequency division, energy comparison signal processing system |
US6104345A (en) * | 1998-04-09 | 2000-08-15 | State Of Israel, Ministry Of Defense Armament Development Authority-Rafael | Direction of arrival tracking of multiple targets |
CN109188443A (en) * | 2018-06-29 | 2019-01-11 | 中国船舶重工集团公司第七〇五研究所 | A kind of passive target tracking method based on Interactive Multiple-Model |
CN109298413A (en) * | 2018-09-01 | 2019-02-01 | 哈尔滨工程大学 | A kind of method that specific aim solves the problems, such as the multiple target tracking data correlation under complex electromagnetic environment |
CN109444897A (en) * | 2018-09-13 | 2019-03-08 | 中国船舶重工集团公司第七〇五研究所 | A kind of more gusts of Data Associations based on multiple features |
JP2020034363A (en) * | 2018-08-29 | 2020-03-05 | 沖電気工業株式会社 | Signal tracking device, signal tracking system, signal tracking method, and program |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6819230B2 (en) * | 2002-08-08 | 2004-11-16 | The United States Of America As Represented By The Secretary Of The Navy | Target track crossing prediction/detection |
-
2020
- 2020-12-04 CN CN202011402804.6A patent/CN112684454B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5119432A (en) * | 1990-11-09 | 1992-06-02 | Visidyne, Inc. | Frequency division, energy comparison signal processing system |
US6104345A (en) * | 1998-04-09 | 2000-08-15 | State Of Israel, Ministry Of Defense Armament Development Authority-Rafael | Direction of arrival tracking of multiple targets |
CN109188443A (en) * | 2018-06-29 | 2019-01-11 | 中国船舶重工集团公司第七〇五研究所 | A kind of passive target tracking method based on Interactive Multiple-Model |
JP2020034363A (en) * | 2018-08-29 | 2020-03-05 | 沖電気工業株式会社 | Signal tracking device, signal tracking system, signal tracking method, and program |
CN109298413A (en) * | 2018-09-01 | 2019-02-01 | 哈尔滨工程大学 | A kind of method that specific aim solves the problems, such as the multiple target tracking data correlation under complex electromagnetic environment |
CN109444897A (en) * | 2018-09-13 | 2019-03-08 | 中国船舶重工集团公司第七〇五研究所 | A kind of more gusts of Data Associations based on multiple features |
Non-Patent Citations (5)
Title |
---|
A modified Eckart post-beamformer filter for improved detection using broadband features;Sanjay K. Mehta等;《1996 IEEE International Conference on Acoustics, Speech, and Signal Processing Conference Proceedings》;20020806;第3045-3048页 * |
一种基于调制频段补偿的宽带解调方法;饶柏斌 等;《声学技术》;20190430;第38卷(第2期);第346-348页 * |
不同类型水听器布阵阵元域融合可行性分析;陈卓;《电子世界》;20190331;第107-109页 * |
基于Hough变换的目标交叉跟踪算法;楼万翔 等;《舰船电子工程》;20160930;第36卷(第9期);第35-38页 * |
被动声呐多目标自动检测与跟踪算法研究;高运 等;《声学与电子工程》;20180430(第4期);第15-17页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112684454A (en) | 2021-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108680910B (en) | Frequency modulation broadcast external radiation source radar target detection method based on waveform cognition | |
CN108828566B (en) | Underwater pulse signal identification method based on towed linear array | |
CN101825707B (en) | Monopulse angular measurement method based on Keystone transformation and coherent integration | |
CN104267379A (en) | Active and passive radar cooperative anti-interference method based on waveform design | |
CN104777450B (en) | A kind of two-stage MUSIC microphone array direction-finding method | |
CN104811886B (en) | Microphone array direction-finding method based on phase difference measurement | |
CN112684454B (en) | Track cross target association method based on sub-frequency bands | |
CN111751792A (en) | Passive radar signal sorting PDW preprocessing method | |
CN111751797B (en) | Method for determining first-order and second-order echo spectrum boundary of high-frequency ground wave radar based on azimuth angle | |
Yang et al. | Improved CFAR detection and direction finding on time–frequency plane with high-frequency radar | |
CN113835068A (en) | Blind source separation real-time main lobe interference resisting method based on independent component analysis | |
CN113589255B (en) | Arrival angle estimation method based on multi-frequency joint sparse Bayesian learning | |
CN115494497A (en) | Phase coding pulse radar signal direction finding method | |
CN115267660A (en) | Underwater narrowband signal anti-interference orientation and automatic ambiguity resolution method | |
CN111796241A (en) | Side lobe suppression method for amplitude comparison between multi-beam array beams | |
CN115436896A (en) | Rapid radar single-snapshot MUSIC angle measurement method | |
CN109655794B (en) | Detection and identification method for suppressing interference by narrow-band self-defense noise | |
CN105353340A (en) | Double-layer cylindrical array underwater passive target detection method | |
CN110764066A (en) | Target detection method based on real signal subspace under existence of error | |
CN104391288B (en) | The false-alarm elimination method of matching interferometric phase | |
CN109001678B (en) | Thunder detection and positioning method based on three-dimensional microphone array | |
CN109814065B (en) | Beam forming method based on phase factor weighting | |
CN110221289B (en) | Target detection method for three-coordinate active phased array radar | |
CN105652256B (en) | A kind of high-frequency ground wave radar TBD methods based on polarization information | |
CN109061597B (en) | Ionospheric clutter suppression method based on blind source separation and time-frequency ridge domain filtering |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |