CN105242259B - Time series packet filtering and Magnitude Difference approach direction-finding method gradually - Google Patents
Time series packet filtering and Magnitude Difference approach direction-finding method gradually Download PDFInfo
- Publication number
- CN105242259B CN105242259B CN201510566904.5A CN201510566904A CN105242259B CN 105242259 B CN105242259 B CN 105242259B CN 201510566904 A CN201510566904 A CN 201510566904A CN 105242259 B CN105242259 B CN 105242259B
- Authority
- CN
- China
- Prior art keywords
- arteries
- msub
- mrow
- decibels
- veins group
- 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
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Direction-finding method is approached gradually the invention discloses time series packet filtering and Magnitude Difference, including long short pulse is grouped, selects filter coefficient, determines current goal place arteries and veins group and determine the step of target bearing four.By the single arteries and veins group in CPI, packet filtering processing is limited the direction-finding method based on time series packet filtering and Magnitude Difference approximation Strategy gradually that the present invention is provided to reduce CPI width again again in prior art basis, can improve the azimuth resolution of target detection and the bearing accuracy of detection target.
Description
Technical field
The present invention relates to a kind of radar target direction-finding method, more particularly to a kind of time series packet filtering and Magnitude Difference
Gradually direction-finding method is approached.
Background technology
Radar based on CPI (coherent processing interval) tupe is mostly handled using MTD (moving-target detection) mode, but
Because MTD processing target bearing precision is limited by CPI width, the azimuth resolution of target detection is than relatively low.Therefore,
Generally using two in the CPI adjacent filtered results of arteries and veins group by interpolation method, centroid algorithm or antenna beam correlation method come
Improve the bearing accuracy of detection target.Nonetheless, the precision of this several algorithm is still limited by the limitation of CPI width, target inspection
The azimuth resolution of survey is not fully up to expectations.
The content of the invention
In order to solve the above problems, survey is approached gradually the invention provides a kind of time series packet filtering and Magnitude Difference
To method.By single arteries and veins group in CPI, packet filtering processing reduces CPI's to this method again again in prior art basis
Width is limited, so as to improve the azimuth resolution of target detection and the bearing accuracy of detection target.
Above-mentioned purpose is achieved by the following technical solution：
A kind of direction-finding method based on time series packet filtering and Magnitude Difference approximation Strategy gradually, it is characterised in that bag
Include：
Step one, long short pulse packet：According to radar horizon, long short pulse in each CPI is grouped to form multiple arteries and veins
Group, each arteries and veins group uses the wave filter of corresponding exponent number；
Step 2, selects filter coefficient：Clutter map is set up according to the pulse echo of different arteries and veins groups, and it is strong according to clutter map
The degree selection arteries and veins group uses the coefficient of wave filter；
Step 3, determines arteries and veins group where current goal：After the exponent number and coefficient that determine each arteries and veins group wave filter, by CPI two
The pulse echo of individual adjacent arteries and veins group carries out MTD filtering process respectively, and compares wave filter output amplitude value, and current goal falls defeated
Go out in the larger arteries and veins group of range value；
Step 4, determines target bearing：Arteries and veins group where the current goal that step 3 is determined is further divided into left and right two
By wave filter output amplitude value after individual half arteries and veins group, each half arteries and veins group filtering process, and according to below equation, current goal is calculated
Orientation：
Wherein, θtargetFor current goal orientation；θ1、θ2For the center hold of the arteries and veins group of left and right two and half；A1、A2For left and right two
Individual half arteries and veins group wave filter output amplitude value；K value calculation formula are as follows：
Wherein, θ3dBFor antenna 3dB beam angles.
Further, the method for selection filter coefficient is：Clutter map intensity is divided into 4 grades, clutter map intensity is small
Weak noise filter coefficient is used when 15 decibels, using the filter of moderate clutter when clutter map intensity is more than 15 decibels less than 35 decibels
Ripple device coefficient, clutter map intensity uses strong noise filter coefficient when being more than 35 decibels less than 55 decibels, and clutter map intensity is more than
Superpower noise filter coefficient is used at 55 decibels.
Further, the frequency response curve of the weak noise filter coefficient is 35 points in the maximum notch depth of zero-frequency
Shellfish；The frequency response curve of the middle noise filter coefficient is 40 decibels in the maximum notch depth of zero-frequency；The strong clutter
The frequency response curve of filter coefficient is 60 decibels in the maximum notch depth of zero-frequency；The superpower noise filter coefficient frequency
Rate response curve is 75 decibels in the maximum notch depth of zero-frequency.
Further, long short pulse is divided into short pulse sequence, a short long pulse intervening sequence, long pulse in each CPI
Rush three arteries and veins groups of sequence.
Beneficial effects of the present invention：
The direction-finding method based on time series packet filtering and Magnitude Difference approximation Strategy gradually that the present invention is provided is existing
Have in technical foundation by the single arteries and veins group in CPI again again packet filtering processing come reduce CPI width limit, can improve
The azimuth resolution of target detection and the bearing accuracy of detection target.
Brief description of the drawings
Fig. 1：Long short pulse is grouped schematic diagram；
Fig. 2：Different noise intensity filter response curve figures；
Fig. 3：Half arteries and veins group filtering process schematic diagram.
Embodiment
Embodiment of the present invention is described in further detail below in conjunction with embodiment and accompanying drawing.
The invention provides a kind of direction-finding method based on time series packet filtering and Magnitude Difference approximation Strategy gradually,
Including：
Step one, long short pulse packet：According to radar horizon, long short pulse in each CPI is grouped to be formed it is multiple not
With arteries and veins group, different arteries and veins groups use the wave filter of different rank；As shown in figure 1, Fig. 1 middle arteries group 1 is short pulse sequence, tranmitting frequency
F1；Arteries and veins group 2 uses a short long pulse intervening sequence, short pulse tranmitting frequency F2, long pulse tranmitting frequency F1；Arteries and veins group 3 is length
Pulse train, tranmitting frequency F2.Wherein short pulse is blind near region benefit, and frequency changes the detection probability for improving target.Arteries and veins
Group 1 uses the rank wave filter of 10 pulse 12, and arteries and veins group 2 and arteries and veins group 3 use the rank wave filter of 5 pulse 7.The exponent number of wave filter is by radar day
Line swing circle, antenna main lobe width and pulse repetition period determine jointly.When the antenna rotation period and antenna main lobe of radar
After width (i.e. antenna beam residence time) is determined, because bursts period is short, the cycle of occupancy is more, can use higher order filter
Handled, and the long pulse cycle is long, the cycle of occupancy is few, therefore is handled using lower order filter.
Step 2, selects filter coefficient：Clutter map is set up according to the pulse echo of different arteries and veins groups, and it is strong according to clutter map
The degree selection arteries and veins group uses the coefficient of wave filter；As shown in Fig. 2 clutter map intensity is divided into 4 grades, clutter map intensity is small
Weak noise filter coefficient is used when 15 decibels, using the filter of moderate clutter when clutter map intensity is more than 15 decibels less than 35 decibels
Ripple device coefficient, clutter map intensity uses strong noise filter coefficient when being more than 35 decibels less than 55 decibels, and clutter map intensity is more than
Superpower noise filter coefficient is used at 55 decibels；Maximum of the frequency response curve of the weak noise filter coefficient in zero-frequency
Notch depth is 35 decibels；The frequency response curve of the middle noise filter coefficient is 40 points in the maximum notch depth of zero-frequency
Shellfish；The frequency response curve of the strong noise filter coefficient is 60 decibels in the maximum notch depth of zero-frequency；It is described superpower miscellaneous
Wave filter coefficient frequencies response curve is 75 decibels in the maximum notch depth of zero-frequency；
Step 3, determines arteries and veins group where current goal：After the exponent number and coefficient that determine each arteries and veins group wave filter, by CPI two
The pulse echo of individual adjacent arteries and veins group carries out MTD filtering process respectively, and compares wave filter output amplitude value, and current goal falls defeated
Go out in the larger arteries and veins group of range value；
Step 4, determines target bearing：Arteries and veins group where the current goal that step 3 is determined is further divided into left and right two
By wave filter output amplitude value after individual half arteries and veins group (as shown in Figure 3), each half arteries and veins group filtering process, and according to below equation,
Calculate current goal orientation：
Wherein, θtargetFor current goal orientation；θ1、θ2For the center hold of the arteries and veins group of left and right two and half；A1、A2For left and right two
Individual half arteries and veins group wave filter output amplitude value；K value calculation formula are as follows：
Wherein, θ3dBFor antenna 3dB beam angles.
The direction-finding method based on time series packet filtering and Magnitude Difference approximation Strategy gradually that the present invention is provided is existing
Have in technical foundation by the single arteries and veins group in CPI again again packet filtering processing come reduce CPI width limit, can improve
The azimuth resolution of target detection and the bearing accuracy of detection target.
The effect of above-described embodiment indicates that the essentiality content of the present invention, but the protection of the present invention is not limited with this
Scope.It will be understood by those within the art that, technical scheme can be modified or equivalent substitution,
Without departing from the essence and protection domain of technical solution of the present invention.
Claims (4)
1. a kind of time series packet filtering and Magnitude Difference approach direction-finding method gradually, it is characterised in that including：
Step one, long short pulse packet：According to radar horizon, long short pulse in each CPI is grouped to form multiple arteries and veins groups, respectively
Arteries and veins group uses the wave filter of corresponding exponent number；
Step 2, selects filter coefficient：Clutter map is set up according to the pulse echo of different arteries and veins groups, and selected according to clutter map intensity
Select the coefficient that the arteries and veins group uses wave filter；
Step 3, determines arteries and veins group where current goal：After the exponent number and coefficient that determine each arteries and veins group wave filter, by two phases in CPI
The pulse echo of adjacent arteries and veins group carries out MTD filtering process respectively, and compares wave filter output amplitude value, and current goal falls in output width
In the larger arteries and veins group of angle value；
Step 4, determines target bearing：Arteries and veins group where the current goal that step 3 is determined is further divided into left and right two and half
By wave filter output amplitude value after arteries and veins group, each half arteries and veins group filtering process, and according to below equation, current goal side is calculated
Position：
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mo>&times;</mo>
<mi>L</mi>
<mi>n</mi>
<mfrac>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
Wherein, θtargetFor current goal orientation；θ1、θ2For the center hold of the arteries and veins group of left and right two and half；A1、A2For left and right two and half
Arteries and veins group wave filter output amplitude value；K value calculation formula are as follows：
<mrow>
<mi>K</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
<mrow>
<mn>3</mn>
<mi>d</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
Wherein, θ3dBFor antenna 3dB beam angles；
CPI refers to coherent processing interval, and MTD refers to moving-target detection.
2. time series packet filtering according to claim 1 and Magnitude Difference approach direction-finding method gradually, its feature exists
In the method for selection filter coefficient is：Clutter map intensity is divided into 4 grades, clutter map intensity is used when being less than 15 decibels
Weak noise filter coefficient, clutter map intensity uses moderate noise filter coefficient, clutter when being more than 15 decibels less than 35 decibels
Figure intensity uses strong noise filter coefficient when being more than 35 decibels less than 55 decibels, clutter map intensity uses super when being more than 55 decibels
Strong noise filter coefficient.
3. time series packet filtering according to claim 2 and Magnitude Difference approach direction-finding method gradually, its feature exists
In：The frequency response curve of the weak noise filter coefficient is 35 decibels in the maximum notch depth of zero-frequency；The moderate is miscellaneous
The frequency response curve of wave filter coefficient is 40 decibels in the maximum notch depth of zero-frequency；The strong noise filter coefficient
Frequency response curve is 60 decibels in the maximum notch depth of zero-frequency；The superpower noise filter coefficient frequencies response curve exists
The maximum notch depth of zero-frequency is 75 decibels.
4. the time series packet filtering and Magnitude Difference according to claim 1,2 or 3 approach direction-finding method gradually, it is special
Levy and be：Long short pulse is divided into short pulse sequence, a short long pulse intervening sequence, long pulse sequence three in each CPI
Arteries and veins group.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510566904.5A CN105242259B (en) | 2015-09-08 | 2015-09-08 | Time series packet filtering and Magnitude Difference approach direction-finding method gradually |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510566904.5A CN105242259B (en) | 2015-09-08 | 2015-09-08 | Time series packet filtering and Magnitude Difference approach direction-finding method gradually |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105242259A CN105242259A (en) | 2016-01-13 |
| CN105242259B true CN105242259B (en) | 2017-08-29 |
Family
ID=55039965
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510566904.5A Active CN105242259B (en) | 2015-09-08 | 2015-09-08 | Time series packet filtering and Magnitude Difference approach direction-finding method gradually |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105242259B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101470202A (en) * | 2007-12-26 | 2009-07-01 | 清华大学 | Pulse Doppler radar system and its signal processing method |
| CN102565763A (en) * | 2011-12-12 | 2012-07-11 | 中国船舶重工集团公司第七二四研究所 | Adaptive clutter suppression moving-target signal processing technology and realizing method |
| CN104155631A (en) * | 2014-08-22 | 2014-11-19 | 西安电子科技大学 | Self-adaptive pulse number distribution method based on airborne radar clutter spectral width |
| CN104330782A (en) * | 2014-11-04 | 2015-02-04 | 西安电子科技大学 | Time domain and modulation domain parameter combined measuring method of triangular frequency-modulation pulse signals |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5640583B2 (en) * | 2010-09-06 | 2014-12-17 | 日本電気株式会社 | Target detection system, detection method, and detection information processing program |
-
2015
- 2015-09-08 CN CN201510566904.5A patent/CN105242259B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101470202A (en) * | 2007-12-26 | 2009-07-01 | 清华大学 | Pulse Doppler radar system and its signal processing method |
| CN102565763A (en) * | 2011-12-12 | 2012-07-11 | 中国船舶重工集团公司第七二四研究所 | Adaptive clutter suppression moving-target signal processing technology and realizing method |
| CN104155631A (en) * | 2014-08-22 | 2014-11-19 | 西安电子科技大学 | Self-adaptive pulse number distribution method based on airborne radar clutter spectral width |
| CN104330782A (en) * | 2014-11-04 | 2015-02-04 | 西安电子科技大学 | Time domain and modulation domain parameter combined measuring method of triangular frequency-modulation pulse signals |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105242259A (en) | 2016-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103885053B (en) | A kind of radar data moving-target based on tracking filter detects processing method | |
| CN104215951B (en) | System and method for detecting low-speed small target under sea cluster background | |
| CN101226236B (en) | Sky wave over-the-horizon radar self-adaption interference rejection method based on sidelobe constraint | |
| CN102520395B (en) | Clutter suppression method based on bistatic multiple-input and multiple-output radar | |
| CN103018727A (en) | Sample-training-based non-stationary clutter suppression method of vehicle-mounted radar | |
| CN102890272A (en) | Method for processing millimeter wave cloud radar signal | |
| CN104076350B (en) | Based on the offshore floating object detection method of normalization doppler spectral | |
| CN105403864B (en) | Based on the two-dimentional boat-carrying high-frequency ground wave radar ocean clutter cancellation method for improving oblique projection | |
| CN104237866A (en) | Improved constant false alarm rate detecting method for ship-borne linear frequency modulation continuous wave radar | |
| CN104168232B (en) | The assay method of multidiameter delay and Doppler frequency shift in a kind of underwater acoustic channel | |
| CN104656074A (en) | Target detecting method of weighted robust broadband beam forming | |
| CN103954941B (en) | The multipulse cognitive clutter suppression method of a kind of airborne phased array radar two dimension | |
| CN104316914A (en) | Radar target self-adaptation detection method depending on shape parameters | |
| CN105137437A (en) | Target detection method based on spatial domain phase variance weighting | |
| CN108804736A (en) | A kind of method and apparatus multiple degrees of freedom FM signal design and optimized | |
| CN103064077A (en) | Amplitude-weighting-based target detection method and equipment | |
| CN104407337B (en) | Towards the data filtering methods and system of turning maneuvering target | |
| CN106546949A (en) | A kind of double array element sinusoidal signal arrival bearing's methods of estimation based on frequency estimation meter | |
| CN109581362A (en) | Signal processing method of the synthetic aperture radar altimeter under variable pulse cluster mode | |
| CN106199549A (en) | A kind of method using spectrum-subtraction to promote LFMCW radar signal to noise ratio | |
| CN107526077A (en) | A kind of high inversion method of portable high frequency groundwave radar wave | |
| CN102788970B (en) | High-efficiency self-adapting suppression method for short-range clutters of airborne digital array radar | |
| CN103792523B (en) | UHF wave band Multichannel radar radial velocity detection method based on tensor product | |
| CN105242259B (en) | Time series packet filtering and Magnitude Difference approach direction-finding method gradually | |
| CN105548987B (en) | A kind of continuous wave radar aimed acceleration blind estimating method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| CB03 | Change of inventor or designer information |
Inventor after: Geng Zhaoqian Inventor after: Yang Hui Inventor after: Li Xuming Inventor after: Ding Zhiqing Inventor after: Zhang Hui Inventor before: Geng Zhaoqian |
|
| COR | Change of bibliographic data | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder |
Address after: 210038 No.9 Hengyi Road, Nanjing Economic and Technological Development Zone, Jiangsu Province Patentee after: NANJING CHANGJIANG ELECTRONICS GROUP Co.,Ltd. Address before: No. 9, Hengyi Road, Xingang economic and Technological Development Zone, Nanjing, Jiangsu 210038 Patentee before: NANJING CHANGJIANG ELECTRONICS GROUP Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder |