CN108535704A - A kind of signal Pre-sorting method based on self-adaption two-dimensional cluster - Google Patents
A kind of signal Pre-sorting method based on self-adaption two-dimensional cluster Download PDFInfo
- Publication number
- CN108535704A CN108535704A CN201810313620.9A CN201810313620A CN108535704A CN 108535704 A CN108535704 A CN 108535704A CN 201810313620 A CN201810313620 A CN 201810313620A CN 108535704 A CN108535704 A CN 108535704A
- Authority
- CN
- China
- Prior art keywords
- distance
- bearing
- sorting
- threshold value
- dimension
- 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
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of signal Pre-sorting methods based on self-adaption two-dimensional cluster, can be used for the signal pre-identification of new system radar under complex environment, this approach includes the following steps：Using aerial array to target two-dimensional localization, distance and bearing information is obtained；The threshold value initial value of distance and bearing is set；Calculate the similarity of target range peacekeeping azimuth dimension；According to the similarity and threshold value of calculating, the two dimension cluster of distance and bearing is obtained；The threshold value of adaptive adjustment distance and bearing, the similarity mean value after Statistical Clustering Analysis apart from peacekeeping azimuth dimension complete signal Pre-sorting until when mean value is less than the threshold value in distance or orientation.The method that the present invention uses Pre-sorting after first positioning, breach the limitation of conventional radar signal sorting, using the two-dimentional clustering processing of distance and bearing, increasing batch and the leakage batch phenomenon of getting worse in signal sorting are efficiently solved, the validity of radar signal sorting is improved.
Description
Technical field
The present invention relates to Radar Signal Processing Technology fields in electronic countermeasure, especially a kind of to be clustered based on self-adaption two-dimensional
Signal Pre-sorting method.
Background technology
Electronic countermeasure is the important means of warfare of attack and defense in modernized war, and radar signal sorting is then multiple
It is the important component of electronic warfare system to the radar emission signal sorting of unknown parameters under strays magnetic environment.
Warp-wise has been multi-functional for modern radar, the direction of multiduty New System, new technology is developed, and a radar may have more
Kind working condition and a variety of systems, and various complicated wave forms are often devised, this just destroys signal sorting and identification institute profit
Regularity makes reconnaissance system for radar greatly be challenged, and the increase of various electronic counter-measures equipment quantity so that radar
Reconnaissance system is under complexity, intensive electromagnetic signal environment, and the complexity and operand of signal processing are concentrated mainly on signal
In sorting process, the level of signal sorting has been to weigh the important indicator of countermeasures set advance.
For complex electromagnetic environment, radar emitter signal can usually overlap on time domain, spatial domain, frequency domain, tradition
Radar Signal Sorting Method there are serious " increase batch " and " leakage batch " phenomenon, the high density of signal environment, but also tradition is divided
Select computational processing big, these can all lead to the failure of signal sorting, therefore there is an urgent need to explore new signal sorting method.This
Invention is reduced the operand of follow-up signal sorting process, it is correct can be reached signal using the method for Pre-sorting after first positioning
The purpose of sorting.
Invention content
The goal of the invention of the present invention is：In view of the above problems, it is small to provide a kind of operand, sorting can be improved
The method of correctness.
A kind of signal Pre-sorting method based on self-adaption two-dimensional cluster proposed by the present invention, this method include：
1) it utilizes aerial array to target two-dimensional localization, obtains distance and bearing (R, θ) information sequence；
2) threshold value (Δ R, Δ θ) initial value of distance and bearing is set；
3) similarity (d of target range peacekeeping azimuth dimension is calculatedR,dθ)；
4) according to the similarity of calculating and threshold value, the two dimension cluster of distance and bearing is obtained；
5) the adaptive threshold value of adjustment distance and bearing, the similarity mean value after Statistical Clustering Analysis apart from peacekeeping azimuth dimension
When mean value is less than the threshold value of distance or orientation, signal Pre-sorting is completed.
In the above method：The array antenna is specially to target two-dimensional localization：Using heavy caliber thinned array, when R≤
2D2Target handles array near-field region when/λ, and target range and azimuth information, R can be obtained simultaneously based on array spherical wave model
For target range, D is array bore, and λ is target wavelength；The non-sorting signals of wideband radar are first implemented to position, different
It is scanned in distance R and azimuth angle theta, spatial spectrum is calculated to each step frequency point k：
Wherein, N is the covariance matrix of noise,For the array flow vector under different frequency point, then broadband thunder
It is up to signal space spectrum：
K is stepping frequency points in above formula；Search radar signal space spectrum peak value, the corresponding scanning distance of this peak value and
Orientation is exactly the distance and bearing information of target.
In the step 2), according to thinned array to the precision of object ranging and direction finding, distance and bearing threshold value is set
Initial value.
More specifically, in the step 3), the similarity apart from peacekeeping azimuth dimension is calculated using Euclidean distance：
Wherein dRFor the Euclidean distance of distance dimension, M is pulse train number, RrefFor reference
Distance, it is reference distance generally to take the distance measured for the first time；
Wherein dθFor the Euclidean distance of azimuth dimension, θrefFor reference azimuth, generally take for the first time
The orientation measured is reference azimuth.
More specifically, in the step 4), when pulse train is meeting d simultaneously apart from peacekeeping azimuth dimension similarityR≤
Δ R and dθWhen≤Δ θ, sequence is gathered for one kind, and a new class is otherwise rebuild, and so on, until all pulse trains
Complete two dimension cluster.
More specifically, in the step 5), mean value of the pulse train apart from peacekeeping azimuth dimension similarity is calculated, with as follows
Formula：
With
Meet when simultaneouslyWithWhen, adjust the threshold value of distance and bearing, repeating said steps 3), 4) and
5) until no longer needing to adjustment threshold value, then signal Pre-sorting is completed.
The invention has the advantages that the present invention breaches conventional radar signal using the method for Pre-sorting after first positioning
The limitation of sorting efficiently solves the increasing batch of getting worse in signal sorting using the two-dimentional clustering processing of distance and bearing
With leakage batch phenomenon, the validity of radar signal sorting is improved.
Description of the drawings
Fig. 1 is the flow chart of the signal Pre-sorting method proposed by the present invention based on self-adaption two-dimensional cluster；
Fig. 2 is frequency Pre-sorting result schematic diagram in emulation experiment of the present invention；
Fig. 3 is in emulation experiment of the present invention apart from Pre-sorting result schematic diagram；
Fig. 4 is orientation Pre-sorting result schematic diagram in emulation experiment of the present invention.
Specific implementation mode
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with the attached drawing in the present invention, to this
Technical solution in invention is clearly and completely described.Based on the method in the present invention, those of ordinary skill in the art are not having
There is the every other method obtained under the premise of making creative work, shall fall within the protection scope of the present invention.
Fig. 1 is the flow chart of the signal Pre-sorting method proposed by the present invention based on self-adaption two-dimensional cluster, as shown,
This method comprises the following steps：
Step S1 obtains distance and bearing (R, θ) information sequence using aerial array to target two-dimensional localization；
Described is to utilize heavy caliber thinned array to target two-dimensional localization, as R≤2D2Target processing array near field region when/λ
Domain can obtain target range and azimuth information simultaneously based on array spherical wave model, and R is target range, and D is array bore, and λ is
Target wavelength.The non-sorting signals of wideband radar are first implemented to position, are scanned on different distance R and azimuth angle theta, to each step
Spatial spectrum is calculated into Frequency point k：
Wherein, N is the covariance matrix of noise, αfk(R, θ) is the array flow vector under different frequency point, then wideband radar
Signal space is composed：
K is stepping frequency points in above formula.Search radar signal space spectrum peak value, the corresponding scanning distance of this peak value and
Orientation is exactly the distance and bearing information of target.
Threshold value (Δ R, Δ θ) initial value of distance and bearing is arranged in step S2；
The threshold value is the precision to object ranging and direction finding according to thinned array, and the initial of distance and bearing threshold value is arranged
Value.
Step S3 calculates the similarity of target range peacekeeping azimuth dimension；
Similarity apart from peacekeeping azimuth dimension is calculated using Euclidean distance：
Wherein dRFor the Euclidean distance of distance dimension, M is pulse train number, RrefFor reference
Distance, it is reference distance generally to take the distance measured for the first time.
Wherein dθFor the Euclidean distance of azimuth dimension, θrefFor reference azimuth, generally take for the first time
The orientation measured is reference azimuth.
Step S4 obtains the two dimension cluster of distance and bearing according to the similarity and threshold value of calculating；
The two dimension cluster is apart from peacekeeping azimuth dimension similarity while to meet d using pulse trainR≤ Δ R and dθ≤
When Δ θ, sequence is gathered for one kind, and a new class is otherwise rebuild, and so on, until to complete two dimension poly- for all pulse trains
Class.
Step S5, the adaptive threshold value for adjusting distance and bearing, the similarity after Statistical Clustering Analysis apart from peacekeeping azimuth dimension are equal
ValueUntil when mean value is less than the threshold value in distance or orientation, signal Pre-sorting is completed.
Mean value of the pulse train apart from peacekeeping azimuth dimension similarity is calculated, with following formula：
With
Meet when simultaneouslyWithWhen, adjust the threshold value of distance and bearing, repeating said steps S3, S4 and
S5 then completes signal Pre-sorting until no longer needing to adjustment threshold value.
Illustrate the attainable technique effect of institute of the invention below by emulation experiment.
Needle simultaneous four broadbands frequency-agile radar radiation source, range-azimuth information in spatial domain are respectively
(10200m, 1 °), (10400m, 2 °), (10600m, 3 °) and (10800m, 4 °), phase between target one (3260MHz, 3290MHz)
Every 10MHz frequency agilities, target two (3300MHz, 3340MHz) is separated by 40MHz frequency agilities, target three (3340MHz, 3370MHz)
It is separated by 10MHz frequency agilities, target four (3220MHz, 3250MHz) is separated by 10MHz and 20MHz frequency agilities.If pulse number M=
100, the initial threshold of distance and bearing is Δ R=100m and Δ θ=0.5 °, and frequency Pre-sorting, distance are pre- during two dimension clusters
Sorting and orientation Pre-sorting result difference are as shown in Figure 2, Figure 3 and Figure 4, reach under complex electromagnetic environment correct Pre-sorting sky
Domain, time domain and the overlapped multi-section radar signal of frequency domain purpose.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
All any modification, equivalent and improvement etc., should all be included in the protection scope of the present invention made by within refreshing and principle.
Claims (6)
1. a kind of signal Pre-sorting method based on self-adaption two-dimensional cluster, which is characterized in that this method includes：
1) it utilizes aerial array to target two-dimensional localization, obtains distance and bearing information (R, θ) sequence；
2) threshold value (Δ R, Δ θ) initial value of distance and bearing is set；
3) similarity (d of target range peacekeeping azimuth dimension is calculatedR,dθ)；
4) according to the similarity of calculating and threshold value, the two dimension cluster of distance and bearing is obtained；
5) the adaptive threshold value of adjustment distance and bearing, the similarity mean value after Statistical Clustering Analysis apart from peacekeeping azimuth dimension
Until when mean value is less than the threshold value in distance or orientation, signal Pre-sorting is completed.
2. according to the method described in claim 1, it is characterized in that, institute's array antenna is specially to target two-dimensional localization：
Using heavy caliber thinned array, as R≤2D2Target handles array near-field region when/λ, can be same based on array spherical wave model
When obtain target range and azimuth information, R is target range, and D is array bore, and λ is target wavelength；Wideband radar is not divided
It selects signal first to implement to position, is scanned on different distance R and azimuth angle theta, spatial spectrum is calculated to each step frequency point k：
Wherein, N is the covariance matrix of noise,For the array flow vector under different frequency point, then wideband-radar signal
Spatial spectrum is：
K is stepping frequency points in above formula；The peak value of search radar signal space spectrum, the corresponding scanning distance of this peak value and orientation
It is exactly the distance and bearing information of target.
3. according to the method described in claim 2, it is characterized in that, in the step (2), foundation thinned array is to object ranging
With the precision of direction finding, the initial value of distance and bearing threshold value is set.
4. according to the method described in claim 3, it is characterized in that, in the step (3), distance dimension is calculated using Euclidean distance
With the similarity of azimuth dimension：
Wherein dRFor the Euclidean distance of distance dimension, M is pulse train number, RrefFor reference distance,
The general distance for taking first time to measure is reference distance；
Wherein dθFor the Euclidean distance of azimuth dimension, θrefFor reference azimuth, generally takes and measure for the first time
Orientation be reference azimuth.
5. according to the method described in claim 4, it is characterized in that, in the step (4), when pulse train is apart from peacekeeping side
Position dimension similarity meets d simultaneouslyR≤ Δ R and dθWhen≤Δ θ, sequence is gathered for one kind, a new class is otherwise rebuild, with this
Analogize, until all pulse trains complete two dimension cluster.
6. according to the method described in claim 5, it is characterized in that, in the step (5),
Mean value of the pulse train apart from peacekeeping azimuth dimension similarity is calculated, with following formula：
With
Meet when simultaneouslyWithWhen, adjust the threshold value of distance and bearing, repeating said steps (3), (4) and (5)
Until no longer needing to adjustment threshold value, then signal Pre-sorting is completed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810313620.9A CN108535704B (en) | 2018-04-10 | 2018-04-10 | Signal pre-sorting method based on self-adaptive two-dimensional clustering |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810313620.9A CN108535704B (en) | 2018-04-10 | 2018-04-10 | Signal pre-sorting method based on self-adaptive two-dimensional clustering |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108535704A true CN108535704A (en) | 2018-09-14 |
| CN108535704B CN108535704B (en) | 2022-01-21 |
Family
ID=63479672
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810313620.9A Active CN108535704B (en) | 2018-04-10 | 2018-04-10 | Signal pre-sorting method based on self-adaptive two-dimensional clustering |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108535704B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110018461A (en) * | 2019-04-16 | 2019-07-16 | 西安电子工程研究所 | Group of Targets Recognition Method based on High Range Resolution and Monopulse estimation |
| CN113447907A (en) * | 2021-09-01 | 2021-09-28 | 湖南艾科诺维科技有限公司 | Radar sorting system control method and radar sorting system |
| CN114296443A (en) * | 2021-11-24 | 2022-04-08 | 贵州理工学院 | Unmanned modular combine harvester |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105891815A (en) * | 2016-05-03 | 2016-08-24 | 西南大学 | Combined estimation algorithm based on broadcast signal passive positioning |
| CN106056098A (en) * | 2016-06-23 | 2016-10-26 | 哈尔滨工业大学 | Pulse signal cluster sorting method based on class merging |
| CN107238817A (en) * | 2017-07-04 | 2017-10-10 | 中国人民解放军海军航空工程学院 | A kind of parameter adaptive setting and the radar emitter signal method for separating of adjust automatically |
-
2018
- 2018-04-10 CN CN201810313620.9A patent/CN108535704B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105891815A (en) * | 2016-05-03 | 2016-08-24 | 西南大学 | Combined estimation algorithm based on broadcast signal passive positioning |
| CN106056098A (en) * | 2016-06-23 | 2016-10-26 | 哈尔滨工业大学 | Pulse signal cluster sorting method based on class merging |
| CN107238817A (en) * | 2017-07-04 | 2017-10-10 | 中国人民解放军海军航空工程学院 | A kind of parameter adaptive setting and the radar emitter signal method for separating of adjust automatically |
Non-Patent Citations (1)
| Title |
|---|
| 乐明军: ""复杂信号环境下雷达信号分选技术方法研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110018461A (en) * | 2019-04-16 | 2019-07-16 | 西安电子工程研究所 | Group of Targets Recognition Method based on High Range Resolution and Monopulse estimation |
| CN110018461B (en) * | 2019-04-16 | 2023-03-24 | 西安电子工程研究所 | Group target identification method based on high-resolution range profile and monopulse angle measurement |
| CN113447907A (en) * | 2021-09-01 | 2021-09-28 | 湖南艾科诺维科技有限公司 | Radar sorting system control method and radar sorting system |
| CN113447907B (en) * | 2021-09-01 | 2021-11-12 | 湖南艾科诺维科技有限公司 | Radar sorting system control method and radar sorting system |
| CN114296443A (en) * | 2021-11-24 | 2022-04-08 | 贵州理工学院 | Unmanned modular combine harvester |
| CN114296443B (en) * | 2021-11-24 | 2023-09-12 | 贵州理工学院 | Unmanned modularized combine harvester |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108535704B (en) | 2022-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Luo et al. | Three-dimensional precession feature extraction of space targets | |
| CN109946665A (en) | The method of acquisition real goal based on array radar | |
| CN108732543B (en) | Airborne networking radar radiation parameter joint optimization method based on radio frequency stealth | |
| US10001550B2 (en) | Method for automatic classification of radar objects | |
| CN111079859B (en) | Passive multi-station multi-target direction finding cross positioning and false point removing method | |
| CN104502911B (en) | Wall parameter estimation method of through-wall imaging radar | |
| CN108535704A (en) | A kind of signal Pre-sorting method based on self-adaption two-dimensional cluster | |
| CN107656238A (en) | A kind of high-frequency ground wave radar direction finding new method based on double LUT Methods | |
| CN104535972A (en) | Coherent transferring type jamming inhibition method for airborne radar | |
| CN103675808B (en) | A kind of monopulse radar seeker can not differentiate multi-target detection method | |
| CN105787081A (en) | Radiation platform correlation method based on radiation source space positions | |
| Wan et al. | Accurate estimation the scanning cycle of the reconnaissance radar based on a single unmanned aerial vehicle | |
| CN110261828B (en) | Distance-angle error two-dimensional clustering-based transverse interference judgment method | |
| CN103064084A (en) | Ambiguity solving method based on distance frequency domain | |
| Liu et al. | Multiple-element retrodirective cross-eye jamming against amplitude-comparison monopulse radar | |
| Pengcheng et al. | Radar group target recognition based on HRRPs and weighted mean shift clustering | |
| CN103728615A (en) | Method and system for detecting a plurality of targets of phased array secondary radars | |
| RU2665032C2 (en) | Device for recognition of aerospace objects in two-radio radar complexes with active phased antenna arrays (apaa) | |
| CN107037398A (en) | A kind of two-dimentional MUSIC algorithms estimate the parallel calculating method of direction of arrival | |
| CN112418275B (en) | Millimeter wave fuse foil strip interference identification method based on range profile feature extraction | |
| RU2304789C1 (en) | Method of radar tracking of object's trajectory | |
| Bingbing et al. | An indoor positioning algorithm and its experiment research based on RFID | |
| CN106353751A (en) | High-resolution fuze and target detecting method | |
| Wang et al. | GLNet for target detection in millimeter wave images | |
| Li et al. | Radar high resolution range profile recognition via multi-SV method |
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 |