CN109283503A - A kind of radar echo signal Doppler frequency center under uniform scene determines method - Google Patents
A kind of radar echo signal Doppler frequency center under uniform scene determines method Download PDFInfo
- Publication number
- CN109283503A CN109283503A CN201710603597.2A CN201710603597A CN109283503A CN 109283503 A CN109283503 A CN 109283503A CN 201710603597 A CN201710603597 A CN 201710603597A CN 109283503 A CN109283503 A CN 109283503A
- Authority
- CN
- China
- Prior art keywords
- doppler frequency
- data
- data stream
- module
- pulse
- 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/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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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 the radar echo signal Doppler frequency centers under a kind of uniform scene to determine method, is realized by echo data stream generation module, digital quadrature mixing module, pulse compression module and Doppler frequency center determining module.Echo data stream generation module samples radar echo signal, generates echo data stream；Digital quadrature mixing module carries out digital quadrature Frequency mixing processing to echo data stream, obtains baseband inphase and quadrature data stream；Pulse compression module carries out process of pulse-compression to baseband inphase and orthogonal data, obtains data after pulse pressure；Doppler frequency center determining module carries out auto-correlation computation to data after pulse pressure, and determines Doppler frequency center according to operation result.The Doppler frequency center precision that the present invention solves the problems, such as that common Doppler frequency center determines that method determines under uniform scene is lower.
Description
Technical field
The present invention relates to a kind of Doppler frequency centers to determine method, and the radar echo signal under especially a kind of uniform scene is more
General Le mass center determines method.
Background technique
Doppler frequency center is important parameter needed for radar seeker imaging, depends primarily on target seeker and ground appearance
Radial velocity between mark.Doppler frequency center value can be obtained by handling radar echo signal.Uniform scene refers to
Observation scene is not mutated, and contrast is stablized in a certain range.Doppler frequency center determine it is accurate whether directly affect imaging
The quality of quality, because mass center error will lead to signal noise ratio (snr) of image reduction, azimuth ambiguity increases, and makes the target position on image
Set generation offset.Doppler frequency center value can be directly calculated by the flight attitude data of the almanac data of satellite or guided missile,
But since the precision of almanac data and missile flight attitude data is subject to certain restrictions, and many times movement of radar platform
And it is irregular, direct calculated Doppler frequency center value is inaccurate.
Summary of the invention
It is an object of that present invention to provide the radar echo signal Doppler frequency centers under a kind of uniform scene to determine method, solves
Common Doppler frequency center determines the lower problem of the Doppler frequency center precision that method determines under uniform scene.
A kind of radar echo signal Doppler frequency center under uniform scene determines method, the specific steps are that:
The first step builds Doppler frequency center and determines system
Doppler frequency center determines system, comprising: echo data stream generation module, digital quadrature mixing module, pulse compression
Module and Doppler frequency center determining module.
The function of echo data stream generation module are as follows: radar echo signal is sampled, echo data stream is generated；
The function of digital quadrature mixing module are as follows: digital quadrature Frequency mixing processing is carried out to echo data stream, it is same to obtain base band
Phase and quadrature data stream；
The function of pulse compression module are as follows: process of pulse-compression is carried out to baseband inphase and orthogonal data, after obtaining pulse pressure
Data；
The function of Doppler frequency center determining module are as follows: auto-correlation computation is carried out to data after pulse pressure, and according to operation result
Determine Doppler frequency center.
Above-mentioned module is run in fpga chip.
Second step echo data stream generation module samples radar echo signal, generates echo data stream
Echo data stream generation module radar echo signal s (t) within the pulse duration is sampled, and number of echoes is generated
According to stream s (n), sample frequency fs, t is time parameter, and n is time domain point index value, and n=0,1 ..., N-1, N is that data flow is long
Degree.
Sample frequency fsSelection principle: fs≥4(f0+ B/2), wherein f0For the centre frequency of radar emission signal, B is
The bandwidth of radar emission signal, f0+ B/2 is the maximum frequency of radar emission signal.
Third step digital quadrature mixing module carries out digital quadrature Frequency mixing processing to echo data stream
Digital quadrature mixing module is by echo data stream s (n) respectively multiplied by digital local oscillator cos (2 π f0And sin (2 π n)
f0N), low-pass filtering treatment is then carried out respectively:
sI(n)=LPF { s (n) cos (2 π f0n)}
sQ(n)=LPF { s (n) sin (2 π f0n)}
Obtain digital quadrature mixing processing result: baseband inphase data flow sI(n) and base band quadrature data flow sQ(n), wherein
LPF { } indicates low-pass filtering treatment.
4th pace pulse compression module carries out process of pulse-compression to baseband inphase and orthogonal data, obtains data after pulse pressure
Pulse compression module is by baseband inphase data flow sI(n) and base band quadrature data flow sQ(n) it is combined into plural form:
sI(n)+j·sQ(n), and FFT processing is carried out:By sI(n)+j·sQ(n) from when
Domain transforms to frequency domain, obtains frequency domain data S (k)；Construct orientation reference function f (t)=exp { j2 π (- γ t2/ 2) it }, goes forward side by side
Row FFT processing:Orientation reference function f (t) is transformed from the time domain into frequency domain, obtains frequency domain number
According to F (k)；The conjugation of frequency domain data S (k) and frequency domain data F (k) is subjected to complex multiplication, IFFT processing is carried out to result of product:By result of product from frequency-domain transform to time domain, data x (n) after pulse pressure is obtained, wherein
K is frequency domain point indices, k=0,1 ..., N-1；J is indicatede-j2πkn/NFor complex representation form；γ is chirp rate；exp
{ } is a kind of complex exponential representation；F*(k) indicate that plural number takes conjugation.
5th step Doppler frequency center determining module carries out auto-correlation computation to data after pulse pressure, and is determined according to operation result
Doppler frequency center
Doppler frequency center determining module carries out auto-correlation computation to data x (n) after pulse pressure:
Auto-correlation computation result R (m) is obtained, formula is used:
Determine Doppler frequency center fDC, wherein m is translation parameters, and T is pulse width, and angle operation is asked in arg [] expression.
So far the determination of radar echo signal Doppler frequency center value under uniform scene is realized.
This method is determined using the relationship between the power spectrum and correlation function of signal at Fourier transform antithesis in time domain
Radar echo signal Doppler frequency center solves common Doppler frequency center and determines the Doppler that method determines under uniform scene
The lower problem of Centroid accuracy, by repeatedly testing, it is believed that such method is effective, feasible, can be under uniform scene to radar
Echo-signal Doppler frequency center is determined, compared with common Doppler frequency center determines method, identified Doppler frequency center
It is higher to be worth precision.
Specific embodiment
A kind of radar echo signal Doppler frequency center under uniform scene determines method, the specific steps are that:
The first step builds Doppler frequency center and determines system
Doppler frequency center determines system, comprising: echo data stream generation module, digital quadrature mixing module, pulse compression
Module and Doppler frequency center determining module.The function of echo data stream generation module are as follows: radar echo signal is sampled, it is raw
At echo data stream；The function of digital quadrature mixing module are as follows: digital quadrature Frequency mixing processing is carried out to echo data stream, obtains base
Band is the same as phase and quadrature data stream；The function of pulse compression module are as follows: process of pulse-compression is carried out to baseband inphase and orthogonal data,
Obtain data after pulse pressure；The function of Doppler frequency center determining module are as follows: auto-correlation computation is carried out to data after pulse pressure, and according to fortune
It calculates result and determines Doppler frequency center.Above-mentioned module is run in fpga chip.
Second step echo data stream generation module samples radar echo signal, generates echo data stream
Echo data stream generation module radar echo signal s (t) within the pulse duration is sampled, and number of echoes is generated
According to stream s (n), sample frequency fs, t is time parameter, and n is time domain point index value, and n=0,1 ..., N-1, N is that data flow is long
Degree.
Sample frequency fsSelection principle: fs≥4(f0+ B/2), wherein f0For the centre frequency of radar emission signal, B is
The bandwidth of radar emission signal, f0+ B/2 is the maximum frequency of radar emission signal.
Third step digital quadrature mixing module carries out digital quadrature Frequency mixing processing to echo data stream
Digital quadrature mixing module is by echo data stream s (n) respectively multiplied by digital local oscillator cos (2 π f0And sin (2 π n)
f0N), low-pass filtering treatment is then carried out respectively:
sI(n)=LPF { s (n) cos (2 π f0n)}
sQ(n)=LPF { s (n) sin (2 π f0n)}
Obtain digital quadrature mixing processing result: baseband inphase data flow sI(n) and base band quadrature data flow sQ(n), wherein
LPF { } indicates low-pass filtering treatment.
4th pace pulse compression module carries out process of pulse-compression to baseband inphase and orthogonal data, obtains data after pulse pressure
Pulse compression module is by baseband inphase data flow sI(n) and base band quadrature data flow sQ(n) it is combined into plural form:
sI(n)+j·sQ(n), and FFT processing is carried out:By sI(n)+j·sQ(n) from when
Domain transforms to frequency domain, obtains frequency domain data S (k)；Construct orientation reference function f (t)=exp { j2 π (- γ t2/ 2) it }, goes forward side by side
Row FFT processing:Orientation reference function f (t) is transformed from the time domain into frequency domain, obtains frequency domain number
According to F (k)；The conjugation of frequency domain data S (k) and frequency domain data F (k) is subjected to complex multiplication, IFFT processing is carried out to result of product:By result of product from frequency-domain transform to time domain, data x (n) after pulse pressure is obtained,
Middle k is frequency domain point indices, k=0,1 ..., N-1；J is indicatede-j2πkn/NFor complex representation form；γ is chirp rate；
Exp { } is a kind of complex exponential representation；F*(k) indicate that plural number takes conjugation.
5th step Doppler frequency center determining module carries out auto-correlation computation to data after pulse pressure, and is determined according to operation result
Doppler frequency center
Doppler frequency center determining module carries out auto-correlation computation to data x (n) after pulse pressure:
Auto-correlation computation result R (m) is obtained, formula is used:
Determine Doppler frequency center fDC, wherein m is translation parameters, and T is pulse width, and angle operation is asked in arg [] expression.
So far the determination of radar echo signal Doppler frequency center value under uniform scene is realized.
Claims (3)
1. the radar echo signal Doppler frequency center under a kind of uniform scene determines method, the specific steps are that:
The first step builds Doppler frequency center and determines system
Doppler frequency center determines system, comprising: echo data stream generation module, digital quadrature mixing module, pulse compression module
With Doppler frequency center determining module；
The function of echo data stream generation module are as follows: radar echo signal is sampled, echo data stream is generated；
The function of digital quadrature mixing module are as follows: to echo data stream carry out digital quadrature Frequency mixing processing, obtain baseband inphase and
Quadrature data stream；
The function of pulse compression module are as follows: process of pulse-compression is carried out to baseband inphase and orthogonal data, obtains data after pulse pressure；
The function of Doppler frequency center determining module are as follows: auto-correlation computation is carried out to data after pulse pressure, and is determined according to operation result
Doppler frequency center；
Above-mentioned module is run in fpga chip；
Second step echo data stream generation module samples radar echo signal, generates echo data stream
Echo data stream generation module radar echo signal s (t) within the pulse duration is sampled, and echo data stream is generated
S (n), sample frequency fs, t is time parameter, and n is time domain point index value, and n=0,1 ..., N-1, N is data flow length；
Third step digital quadrature mixing module carries out digital quadrature Frequency mixing processing to echo data stream
Digital quadrature mixing module is by echo data stream s (n) respectively multiplied by digital local oscillator cos (2 π f0And sin (2 π f n)0N), so
It is filtered respectively afterwards, obtains digital quadrature mixing processing result: baseband inphase data flow sI(n) and base band quadrature data
Flow sQ(n)；
4th pace pulse compression module carries out process of pulse-compression to baseband inphase and orthogonal data, obtains data after pulse pressure
Pulse compression module is by baseband inphase data flow sI(n) and base band quadrature data flow sQ(n) it is combined into plural form: sI(n)
+j·sQ(n), and FFT processing is carried out:By sI(n)+j·sQ(n) become from time domain
Frequency domain is changed to, frequency domain data S (k) is obtained；Construct orientation reference function f (t)=exp { j2 π (- γ t2/ 2) }, and FFT is carried out
Processing:Orientation reference function f (t) is transformed from the time domain into frequency domain, obtains frequency domain data F
(k)；The conjugation of frequency domain data S (k) and frequency domain data F (k) is subjected to complex multiplication, IFFT processing is carried out to result of product:By result of product from frequency-domain transform to time domain, data x (n) after pulse pressure is obtained, wherein
K is frequency domain point indices, k=0,1 ..., N-1；J is indicatede-j2πkn/NFor complex representation form；γ is chirp rate；exp
{ } is a kind of complex exponential representation；F*(k) indicate that plural number takes conjugation；
5th step Doppler frequency center determining module carries out auto-correlation computation to data after pulse pressure, and more according to operation result determination
General Le mass center
Doppler frequency center determining module carries out auto-correlation computation to data x (n) after pulse pressure:It obtains
Auto-correlation computation result R (m) uses formula:
Determine Doppler frequency center fDC, wherein m is translation parameters, and T is pulse width, and angle operation is asked in arg [] expression；
So far the determination of radar echo signal Doppler frequency center value under uniform scene is realized.
2. the radar echo signal Doppler frequency center under uniformly scene as described in claim 1 determines method, it is characterised in that:
Sample frequency fsSelection principle: fs≥4(f0+ B/2), wherein f0For the centre frequency of radar emission signal, B is radar hair
Penetrate the bandwidth of signal, f0+ B/2 is the maximum frequency of radar emission signal.
3. the radar echo signal Doppler frequency center under uniformly scene as claimed in claim 1 or 2 determines that method, feature exist
In: the digital quadrature mixing module is by echo data stream s (n) respectively multiplied by digital local oscillator cos (2 π f0And sin (2 π f n)0N),
Then low-pass filtering treatment is carried out respectively:
sI(n)=LPF { s (n) cos (2 π f0n)}
sQ(n)=LPF { s (n) sin (2 π f0n)}
Obtain digital quadrature mixing processing result: baseband inphase data flow sI(n) and base band quadrature data flow sQ(n), wherein LPF
{ } indicates low-pass filtering treatment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710603597.2A CN109283503B (en) | 2017-07-23 | 2017-07-23 | Method for determining Doppler centroid of radar echo signal in uniform scene |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710603597.2A CN109283503B (en) | 2017-07-23 | 2017-07-23 | Method for determining Doppler centroid of radar echo signal in uniform scene |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109283503A true CN109283503A (en) | 2019-01-29 |
| CN109283503B CN109283503B (en) | 2022-08-23 |
Family
ID=65185956
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710603597.2A Active CN109283503B (en) | 2017-07-23 | 2017-07-23 | Method for determining Doppler centroid of radar echo signal in uniform scene |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109283503B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102288964A (en) * | 2011-08-19 | 2011-12-21 | 中国资源卫星应用中心 | Imaging processing method for spaceborne high-resolution synthetic aperture radar |
| US20120119943A1 (en) * | 2010-06-28 | 2012-05-17 | Alain Bergeron | Method and apparatus for determining a doppler centroid in a synthetic aperture imaging system |
| CN103869314A (en) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | Dual-flight transfer variation bistatic forward-looking synthetic aperture radar imaging method |
| CN104749570A (en) * | 2015-04-01 | 2015-07-01 | 电子科技大学 | Shift invariant airborne bistatic synthetic aperture radar target positioning method |
| CN104820221A (en) * | 2015-04-27 | 2015-08-05 | 电子科技大学 | Target three-dimensional location method of multi-based synthetic aperture radar |
| CN106353742A (en) * | 2016-08-10 | 2017-01-25 | 北京理工大学 | Quick pulse compression method on basis of sparse inverse Fourier transformation |
-
2017
- 2017-07-23 CN CN201710603597.2A patent/CN109283503B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120119943A1 (en) * | 2010-06-28 | 2012-05-17 | Alain Bergeron | Method and apparatus for determining a doppler centroid in a synthetic aperture imaging system |
| CN102288964A (en) * | 2011-08-19 | 2011-12-21 | 中国资源卫星应用中心 | Imaging processing method for spaceborne high-resolution synthetic aperture radar |
| CN103869314A (en) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | Dual-flight transfer variation bistatic forward-looking synthetic aperture radar imaging method |
| CN104749570A (en) * | 2015-04-01 | 2015-07-01 | 电子科技大学 | Shift invariant airborne bistatic synthetic aperture radar target positioning method |
| CN104820221A (en) * | 2015-04-27 | 2015-08-05 | 电子科技大学 | Target three-dimensional location method of multi-based synthetic aperture radar |
| CN106353742A (en) * | 2016-08-10 | 2017-01-25 | 北京理工大学 | Quick pulse compression method on basis of sparse inverse Fourier transformation |
Non-Patent Citations (1)
| Title |
|---|
| 孙泓波等: "机载雷达成像中的两种改进多普勒质心估计算法", 《兵工学报》, no. 01, 20 March 2002 (2002-03-20) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109283503B (en) | 2022-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106405541B (en) | Full-coherent continuous wave Doppler radar and distance and speed measuring method thereof | |
| US11385348B2 (en) | Wideband sonar receiver and sonar signal processing algorithms | |
| CN104898119B (en) | A kind of moving target parameter estimation method based on correlation function | |
| CA2513883A1 (en) | Technique for non-coherent integration of targets with ambiguous velocities | |
| CN106371083B (en) | Radar monitoring point choosing method and device and echo signal processing method and device | |
| EP2524239A1 (en) | Pulse radar range profile motion compensation | |
| CN107153189B (en) | Signal processing method for linear frequency modulation continuous wave radar ranging | |
| CN109613507B (en) | Detection method for high-order maneuvering target radar echo | |
| FR2669117A1 (en) | METHOD FOR RESOLVING AN AMBIGUUITE DURING DETERMINATION OF ANTENNA ANGLE AND DOPPLER FREQUENCY IN AN OPENING SYNTHESIS RADAR | |
| JP5656505B2 (en) | Radar equipment | |
| CN106646451A (en) | Angle tracking subsystem based on radar signal processing system | |
| CN105738889A (en) | Frequency modulated continuous wave speed measurement and distance measurement method | |
| CN115508820A (en) | Target detection method of linear frequency modulation pulse radar | |
| CN108318879A (en) | ISAR image transverse direction calibrating methods based on IAA Power estimation technologies | |
| Zheng et al. | Radar detection and motion parameters estimation of maneuvering target based on the extended keystone transform (July 2018) | |
| CN109655819A (en) | A kind of clutter recognition three-D imaging method based on real aperture Doppler beam sharpening | |
| CN106772308A (en) | Terahertz wideband radar micro-doppler ambiguity solution method based on arteries and veins internal interference | |
| CN112485772A (en) | Clutter suppression method for inter-pulse frequency agility radar | |
| CN105395214B (en) | Ultrasonic imaging processing method and system based on RF data | |
| Dai et al. | High accuracy velocity measurement based on keystone transform using entropy minimization | |
| CN109061626A (en) | A kind of method that Step Frequency coherent processing detects low signal to noise ratio moving-target | |
| CN109283503A (en) | A kind of radar echo signal Doppler frequency center under uniform scene determines method | |
| CN109001671B (en) | Target detection and parameter estimation method and device for frequency hopping signal | |
| CN106955098B (en) | Blood vessel flow velocity calculation method and device | |
| CN110441749A (en) | A kind of Millimeter Wave Stepped-Frequency High Resolution Radar Target moving parameter estimation 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 |