CN109633596A - Multiple target kinematic parameter combined estimation method under pulse Doppler radar system - Google Patents
Multiple target kinematic parameter combined estimation method under pulse Doppler radar system Download PDFInfo
- Publication number
- CN109633596A CN109633596A CN201910064527.3A CN201910064527A CN109633596A CN 109633596 A CN109633596 A CN 109633596A CN 201910064527 A CN201910064527 A CN 201910064527A CN 109633596 A CN109633596 A CN 109633596A
- Authority
- CN
- China
- Prior art keywords
- kinematic parameter
- signal
- transformation
- radar system
- frequency domain
- 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
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
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 present invention provides the multiple target kinematic parameter combined estimation method under a kind of pulse Doppler radar system, detailed process are as follows: step 1: radar return data are modeled, and after carrying out process of pulse-compression to it, fast time frequency domain is transformed to, is obtained apart from frequency domain-orientation time-domain signal；Step 2: based on chockstone transformation to the Range Walk Correction carried out apart from frequency domain-orientation time-domain signal, and to the signal after correction along distance to carrying out inverse Fourier transform；Step 3: the kinematic parameter of Lv Shi transformation estimation target is utilized to the signal obtained after step 2 processing.Multiple target kinematic parameter combined estimation method under pulse Doppler radar system proposed by the present invention, multiple targets can be carried out with chockstone transformation and Lv Shi conversion process simultaneously respectively, therefore Highly precise FFT method can be carried out simultaneously to multiple targets in the same distance unit for being located at same wave beam.
Description
Technical field
The invention belongs to radar surveying technical fields, and in particular to the multiple target fortune under a kind of pulse Doppler radar system
Dynamic parametric joint estimation method.
Background technique
Under the conditions of modern war, often there are multiple targets in radar detection area.When multiple targets are located at same wave
When in beam and in same distance unit, conventional needle can not transport the motion compensation process of single goal to multiple targets simultaneously
Dynamic compensation, therefore accurate multiple target estimates of parameters cannot be obtained.In nineteen ninety " Proceedings of the IEEE
1990 National Aerospace and Electronics Conference " page 210 to page 214, the human hairs such as Wu X
" the Simultaneous imaging of multiple targets in an inverse synthetic of table
One kind is proposed in an aperture radar " text, and parameter Estimation and motion compensation are carried out to multiple target based on maximal possibility estimation
Algorithm, which can provide higher multi-sources distinguishing ability, but need to carry out multi-Dimensional parameters search, computationally intensive.?
Nineteen ninety " IEEE Transactions on Acoustics Speech and Signal Processing " volume 38 the 12nd
Page 2118 to page 2126 of phase, " the Parameter estimation of chirp that Djuric P M et al. is delivered
A kind of algorithm based on solution phase and linear regression is proposed in a signals " text, can be improved computational efficiency, but this calculation
Method is only used for the occasion of simple component signal, and the original signal-to-noise ratio of target echo is required to be higher than 7dB.In " IEEE in 2011
Transactions on Signal Processing " page 3576 to page 3591 of the phase of volume 59 the 8th, Lv X L et al. is delivered
" the text of Lv's distribution:principle, implementation, properties, and performance " one
In propose algorithm-Lv Shi transformation that a kind of pair of Polynomial signal parameter estimated, which is able to solve in time frequency analysis
Contradiction between resolution ratio and cross term, however target across Range cell migration can seriously affect Lv Shi in long integration time
Convert the Parameter Estimation Precision of algorithm.
Summary of the invention
In view of this, the invention proposes the multiple target kinematic parameter Combined estimators under a kind of pulse Doppler radar system
Method, this method, can be simultaneously to the range migrations of multiple targets by carrying out chockstone transformation and Lv Shi conversion process to signal
It is corrected, effectively reduces the computational complexity of parameter search.
Realize that technical scheme is as follows:
A kind of multiple target kinematic parameter combined estimation method under pulse Doppler radar system, detailed process are as follows:
Step 1: modeling radar return data, and after carrying out process of pulse-compression to it, transforms to fast time frequency
Domain is obtained apart from frequency domain-orientation time-domain signal；
Step 2: based on chockstone transformation to the Range Walk Correction carried out apart from frequency domain-orientation time-domain signal, and
To the signal after correction along distance to carry out inverse Fourier transform；
Step 3: the kinematic parameter of Lv Shi transformation estimation target is utilized to the signal obtained after step 2 processing.
Further, step 2 of the present invention is before executing inverse Fourier transform, also according to 1/ (fc+f)≈(1-f/
fc)/fc, wherein f is apart from frequency domain, fcFor carrier frequency, due to chockstone transformation, this is suitable for uniform motion, and the present invention is useful in
In uniformly accelerated motion, it is therefore desirable to do above-mentioned approximation, the result that chockstone is converted makees further transformation, and approximation carries out
After inverse Fourier transform, available linear FM signal can be converted with Lv Shi carry out target component estimation in this way.
Beneficial effect
1) the multiple target kinematic parameter combined estimation method under pulse Doppler radar system proposed by the present invention can be right
Multiple targets carry out chockstone transformation and Lv Shi conversion process simultaneously respectively, therefore can be single to the same distance for being located at same wave beam
Multiple targets in member carry out Highly precise FFT method simultaneously.
2) present invention only passes through quick (inverse) using chockstone transformation and Lv Shi transformation without the prior information for predicting target
Fourier transformation and change mark Fourier transformation, and multi-Dimensional parameters search process is not needed, computational complexity is lower.
Specific embodiment
Below with reference to embodiment, the present invention will be described in detail.
In this example, simulation of Radar System parameter is as follows: carrier frequency fc=1GHz, pulsewidth Tp=4 μ s, bandwidth B=15MHz are adopted
Sample frequency fs=37.5MHz, pulse recurrence frequency PRF=500Hz, correlative accumulation umber of pulse N=2048.
Simulation objectives position and kinematic parameter are as follows: initial radial distance R01=15.3km, R02=15.3km, radial speed
Spend v01=70m/s, v02=60m/s, radial acceleration a01=4m/s2, a02=2m/s2The original for the target echo that radar receives
Beginning signal-to-noise ratio is SNR01=-20dB, SNR02=-25dB.
The embodiment of the present invention provides the multiple target kinematic parameter combined estimation method under a kind of pulse Doppler radar system,
Include the following steps:
Step 1: modeling and process of pulse-compression are carried out to radar return data, the specific method is as follows:
Assuming that in observation scene, there are K point targets.Without loss of generality, the echo model of k-th of target is analyzed.
Assuming that radar emission linear FM signal is
sT(t, τ)=a (τ/Tp)exp(jπγτ2)exp[j2πfc(t+τ)] (1)
Wherein, τ is the fast time, i.e. Distance Time, and t is the slow time, and t=nT (n=0,1 ... N-1), and N is coherent accumulation
The umber of pulse emitted in time, T are the pulse repetition period, and a () is gate function, TpFor pulse width, fcFor carrier frequency, γ is to adjust
Frequency slope.
Target is approximately uniformly accelrated rectilinear motion relative to the radial motion of radar, then between k-th of point target and radar
Instantaneous distance (ignoring the distance change in arteries and veins) can be expressed as
Wherein, R0kFor the initial distance of k-th of point target to radar, v0kFor the radial velocity of k-th of point target, a0kFor
The radial acceleration of k-th of point target.
The base band echo-signal expression formula that radar receives is
Wherein, σ0kFor the reflection coefficient of k-th of point target, λ=c/fcFor wavelength, c is the light velocity.
| v0Under conditions of |/c < < 1, pulse pressure processing is carried out to the base band echo-signal received, and transformed to
Fast time frequency domain can obtain
Wherein B is linear FM signal bandwidth.
(2) formula is substituted into (4) Shi Ke get
Wherein, ξ=1+f/fc。
Step 2: the Range Walk Correction based on chockstone transformation, the specific method is as follows:
To obtain the accurate estimated value of target component, the linear coupling item for the peacekeeping azimuth dimension that needs to adjust the distance is compensated.
Using chockstone transformation adjust the distance after pulse pressure apart from frequency domain --- orientation time-domain signal carries out Range Walk Correction.Chockstone map table
It is up to formula
Bring (6) formula into (5) Shi Ke get
Due to 1/ (fc+f)≈(1-f/fc)/fc, SRk(ta, f) and it can be further represented as
To SRk(ta, f) and it can be obtained along distance to inverse Fourier transform is carried out
From (9), formula be can be seen that, after Range Walk Correction, target tightening in a certain fixed range unit.
Step 3: the kinematic parameter of Lv Shi transformation estimation target, specific side are utilized to the signal obtained after step 2 processing
Method is as follows:
Multiple Target Signals in orientation can be considered that the sum of linear FM signal, formula (9) can be expressed as following form
WhereinG is pulse processing gain.
s′R(ta) parametrization single order symmetrically instantaneous auto-correlation function is defined as
Wherein,For to become the relevant fixed delay of mark operation, and
Wherein,
DefinitionAnd the change mark expression formula is brought intoIt can obtain
As can be seen from the above equation, the coupling between t ' and l has been eliminated, and making FFT to t ' and l respectively can obtain
Fourier transform operation is sought in FFT () expression in above formula.
Then estimates of parameters --- the v of each target can be obtained by threshold judgement0kAnd a0k。
In the present embodiment, radial velocityRadial acceleration
In conclusion the above is merely preferred embodiments of the present invention, being not intended to limit the scope of the present invention.
All within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in of the invention
Within protection scope.
Claims (2)
1. the multiple target kinematic parameter combined estimation method under a kind of pulse Doppler radar system, which is characterized in that specific mistake
Journey are as follows:
Step 1: modeling radar return data, and after carrying out process of pulse-compression to it, transforms to fast time frequency domain,
It obtains apart from frequency domain-orientation time-domain signal；
Step 2: based on chockstone transformation to the Range Walk Correction carried out apart from frequency domain-orientation time-domain signal, and to school
Signal after just is along distance to carrying out inverse Fourier transform；
Step 3: the kinematic parameter of Lv Shi transformation estimation target is utilized to the signal obtained after step 2 processing.
2. the multiple target kinematic parameter combined estimation method under pulse Doppler radar system according to claim 1, special
Sign is that the step 2 is before executing inverse Fourier transform, also according to 1/ (fc+f)≈(1-f/fc)/fc, chockstone is converted
The result arrived makees further transformation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910064527.3A CN109633596A (en) | 2019-01-23 | 2019-01-23 | Multiple target kinematic parameter combined estimation method under pulse Doppler radar system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910064527.3A CN109633596A (en) | 2019-01-23 | 2019-01-23 | Multiple target kinematic parameter combined estimation method under pulse Doppler radar system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109633596A true CN109633596A (en) | 2019-04-16 |
Family
ID=66063349
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910064527.3A Pending CN109633596A (en) | 2019-01-23 | 2019-01-23 | Multiple target kinematic parameter combined estimation method under pulse Doppler radar system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109633596A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111007509A (en) * | 2019-12-17 | 2020-04-14 | 北京理工大学 | Inverse synthetic aperture radar two-dimensional super-resolution imaging method |
| CN111208504A (en) * | 2020-02-28 | 2020-05-29 | 成都汇蓉国科微系统技术有限公司 | PD radar waveform configuration method and device based on DSP |
| CN113740820A (en) * | 2021-09-06 | 2021-12-03 | 西安电子工程研究所 | Mathematical modeling method for pulse Doppler processing of radar signal processor |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7450057B2 (en) * | 2006-10-20 | 2008-11-11 | Northrop Grumman Space & Missions Systems Corp. | Signal processing for accelerating moving targets |
| CN102323575A (en) * | 2011-07-16 | 2012-01-18 | 西安电子科技大学 | Range migration correction method for pulse Doppler (PD) radar in feeble signal detection process |
| CN107450055A (en) * | 2017-07-15 | 2017-12-08 | 西安电子科技大学 | High-speed maneuver object detection method based on Discrete Linear frequency modulation Fourier transform |
| CN108279403A (en) * | 2018-01-04 | 2018-07-13 | 电子科技大学 | Parallel Implementation method is converted based on the Keystone apart from framing |
-
2019
- 2019-01-23 CN CN201910064527.3A patent/CN109633596A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7450057B2 (en) * | 2006-10-20 | 2008-11-11 | Northrop Grumman Space & Missions Systems Corp. | Signal processing for accelerating moving targets |
| CN102323575A (en) * | 2011-07-16 | 2012-01-18 | 西安电子科技大学 | Range migration correction method for pulse Doppler (PD) radar in feeble signal detection process |
| CN107450055A (en) * | 2017-07-15 | 2017-12-08 | 西安电子科技大学 | High-speed maneuver object detection method based on Discrete Linear frequency modulation Fourier transform |
| CN108279403A (en) * | 2018-01-04 | 2018-07-13 | 电子科技大学 | Parallel Implementation method is converted based on the Keystone apart from framing |
Non-Patent Citations (1)
| Title |
|---|
| 田静: "雷达机动目标长时间积累信号处理算法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111007509A (en) * | 2019-12-17 | 2020-04-14 | 北京理工大学 | Inverse synthetic aperture radar two-dimensional super-resolution imaging method |
| CN111208504A (en) * | 2020-02-28 | 2020-05-29 | 成都汇蓉国科微系统技术有限公司 | PD radar waveform configuration method and device based on DSP |
| CN113740820A (en) * | 2021-09-06 | 2021-12-03 | 西安电子工程研究所 | Mathematical modeling method for pulse Doppler processing of radar signal processor |
| CN113740820B (en) * | 2021-09-06 | 2023-07-21 | 西安电子工程研究所 | Mathematical modeling method for pulse Doppler processing of radar signal processor |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Moscardini et al. | Batches algorithm for passive radar: A theoretical analysis | |
| Li et al. | A velocity estimation algorithm of moving targets using single antenna SAR | |
| CN107132534B (en) | Optimization method for high-speed radar target frequency domain detection | |
| Song et al. | Range bias modeling for hyperbolic-frequency-modulated waveforms in target tracking | |
| Chen et al. | A new approach for synthesizing the range profile of moving targets via stepped-frequency waveforms | |
| CN109633596A (en) | Multiple target kinematic parameter combined estimation method under pulse Doppler radar system | |
| CN110850384B (en) | Method for generating broadband deskew echo based on sweep frequency data | |
| CN109613507B (en) | Detection method for high-order maneuvering target radar echo | |
| He et al. | Fast non-searching method for ground moving target refocusing and motion parameters estimation | |
| Wasserzier et al. | On the efficient computation of range and Doppler data in noise radar | |
| CN116559854B (en) | Moving object speed estimation method, device and computer readable storage medium | |
| CN109541540B (en) | Motion single-station pseudo range and velocity combined positioning method | |
| CN112162254A (en) | Method for estimating radial velocity and radial acceleration of target based on ultra-wideband signal | |
| Vu et al. | Direct tracking of multiple targets in MIMO radar | |
| Nagano et al. | Range migration compensation for moving targets with unknown constant velocity in chirp radars | |
| Xia et al. | Doppler rate-of-change estimation method for BPSK signal based on spectrum phase | |
| Garry et al. | Framework and results on passive ISAR | |
| Zou et al. | A sparse representation and GTD model parameter estimation based multiband radar signal coherent compensation method | |
| Moussa et al. | A two-stage sparsity-based method for location and Doppler estimation in bistatic automotive radar | |
| Kang et al. | Micro-Doppler feature extraction for interferometric radar based on Viterbi algorithm and intrinsic chirp component decomposition | |
| Chen et al. | Parameter estimation using nonlinear kalman filters for motion compensation of airborne synthetic aperture radars | |
| CN109633594A (en) | A kind of high-speed moving object parametric joint estimation method based on sub-band division processing | |
| Zhang et al. | Maneuvering target imaging and motion parameter estimation based on improved maximum likelihood estimation | |
| Liu et al. | Fast TDOA and FDOA Estimation for Coherent Pulse Signals | |
| Limin et al. | Synthetic UWB range profile using SFPTs and class MUSIC for high-resolution velocity estimation |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190416 |