CN116559853B - Low repetition frequency solution speed blurring method - Google Patents
Low repetition frequency solution speed blurring method Download PDFInfo
- Publication number
- CN116559853B CN116559853B CN202310833476.2A CN202310833476A CN116559853B CN 116559853 B CN116559853 B CN 116559853B CN 202310833476 A CN202310833476 A CN 202310833476A CN 116559853 B CN116559853 B CN 116559853B
- Authority
- CN
- China
- Prior art keywords
- speed
- target
- measurement result
- pulse
- low repetition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 8
- 238000005259 measurement Methods 0.000 claims description 20
- 238000005070 sampling Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000001427 coherent effect Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 3
- 230000009466 transformation Effects 0.000 description 1
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/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/581—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets
-
- 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
Abstract
A low repetition frequency solution speed blurring method adopts a technical means of roughly estimating the target speed, matching different filter coefficients for different speeds and accurately solving the speed in a slow time dimension, solves the technical effects that in a traditional multi-pulse working mode, aiming at low repetition frequency radar speed blurring, in order to obtain the speed information of a target, priori knowledge or multiple groups of solution blurring is needed to be adopted, so that the application scene is limited or radar time is wasted.
Description
Technical Field
The invention belongs to the technical field of radar target speed measurement, and particularly relates to a low-repetition-frequency de-blurring technology.
Background
The pulse Doppler radar can obtain the speed and distance information of a target at the same time, and is widely applied to modern radar design. The common pulse Doppler radar is divided into three types of low repetition frequency, medium repetition frequency and high repetition frequency according to the pulse repetition frequency.
The low-repetition-frequency radar has the characteristics of low repetition frequency and no blurring of ranging, and is applied to the fields of aerospace, meteorological detection, early warning reconnaissance and the like. For high speed targets, the doppler frequency shift is greater than the radar repetition frequency, creating overlap in the speed dimension, and thus severe speed ambiguity. Radial velocity is becoming an important characteristic parameter as the information basis for target identification, true-false resolution and collaborative detection of numerous tasks.
The following references are referred to:
xu Min research on speed measurement technology of monopulse measurement radar, modern radar 2005
Chen Chao, wang Yufei, qiangli, hu Shiyou research on new method for resolving speed ambiguity of Low-repetition-frequency pulse Doppler radar- "tactical missile technique" 2012
Hu Lang, gaoqingsong, chen Chun, wang Jianjiang method for resolving speed ambiguity of pulse Doppler radar- "sky early warning research report" 2022
Low-priority Zhi Ling, zhang Daizhong, zhang Xiaoju "solution speed fuzzy algorithm based on table lookup" -radar science and technology "2011
Wang Likun, zhou Lin and Guo Doppler solution ambiguity solving method based on multi-group heavy frequency coherent accumulation-modern radar 2020
In the current stage, the method for solving the speed ambiguity is divided into two main types: the speed range of the target is roughly determined by using priori knowledge, and the speed range is solved by methods such as difference, filtering and the like, and the method only can be used for tracking the working state due to the fact that the priori knowledge is required and has certain requirements on the precision of the priori knowledge, and the method cannot be used for searching the working mode without the priori knowledge; the other type of the method is to perform the combined defuzzification through a plurality of pulse groups with different repetition frequencies by using methods such as Chinese remainder theorem or table lookup, and the like, and the method can be used for searching without prior knowledge, but needs a plurality of pulse groups, thereby wasting precious time resources.
Disclosure of Invention
Aiming at low-repetition-frequency radar speed blurring in a traditional multi-pulse working mode, the invention solves the technical problems that the application scene is limited or radar time is wasted because priori knowledge or multiple groups of solutions are required to be utilized for obtaining the speed information of a target; a technical means of roughly estimating the target speed, matching different filter coefficients for different speeds, and accurately solving the speed in a slow time dimension is adopted; the method has the technical effects of realizing the deblurring of the target speed by using only one group of pulses without priori knowledge, saving radar detection time, avoiding the problem of filter mismatch caused by Doppler speed and improving the target signal-to-noise ratio. In order to achieve the above object, the present invention adopts the following technical scheme.
Step one: generating a group of pulse signals S ti (N), where n=1, 2 … N, i=1, 2 … K.
Further, each transmit waveform is made identical and coherent, i.e. S t1 （n）= S t2 （n）= …=S tK （n）=S t （n）。
Further, S with spike blur function feature is selected ti (n) such as single-point frequency signals and Costas code signals, thereby improving the accuracy of speed measurement.
Step two: will S ti (n) up-conversion, transmitting by antenna, returning to target, receiving by antenna, for phased array radar, considering aperture transit effect, mixing to intermediate frequency, sampling pulse by pulse, because of adopting low repetition frequency, distance has no ambiguity, K transmitting pulses correspond to K receiving pulses, i.e. sampling signal is S ri (M), wherein m=1, 2 … M.
Step three: dividing the object to be searched into Q sections along the speed interval, wherein the average speed of each section is V 1 、V 2 …V Q 。
Step four: based on the average speed of each segment, the formula is usedCalculating Q matched filter coefficients MF q (n), where pi is the circumference ratio, V q Is the center speed of the Q-th speed segment, Δt is the sampling interval, λ is the wavelength, q=1, 2 … Q.
Further, the second, third and fourth steps can be performed simultaneously.
Step five: using Q matched filter coefficients, using the formulaThe first pulse is matched filtered, where p=1, 2 … M-n+1, representing the conjugate.
Step six: taking A q (p) absolute values forming a two-dimensional matrixFor the pair ofCFAR detection is carried out on the matrix, if no target exists, detection is continued, and if the target exists, the Q < th > corresponding to the target is recorded T Line N T Column, obtain coarse speed measurement result V QT 。
Step seven: the pulse repetition period divided by the wavelength by 2 times gives the blur velocity V ams According to the rough speed measurement result, use the Q T Matched filter coefficient MF corresponding to line QT (n) matching and filtering other pulses, and measuring the speed in the FFT in the slow time dimension to obtain a refined speed measurement result V JT And adding the rough speed measurement result and the refined speed measurement result to obtain the target speed.
Further, a rounding function round (), and the rough speed measurement result and the refined speed measurement result are substituted into a formulaAnd calculating, updating the target speed and improving the accuracy.
Detailed Description
Assuming that a moving target is at 90km, the speed is 405m/s, the radar selects a C wave band, a frequency point is 6GHz, a single carrier frequency signal, a pulse width is 300us, a pulse repetition period is 1000us,8 pulses are accumulated, and the sampling rate is 0.1MHz, the technical scheme of the invention is specifically described:
selecting single-point frequency signal to construct intermediate frequency transmitting signal (n=1, 2,3 …) S t （n）=1（n=1，2…30）；
Up-converting the signal to 6GHz, transmitting back to target, receiving via antenna, superimposing noise, mixing to intermediate frequency, sampling pulse by pulse, and sampling signal due to DopplerWherein n is mi Is complex Gaussian white noise with standard deviation of 1, the speed is 400m/s, the sampling interval is 10us, and the wavelength is 0.05m; at the same time, the segments are segmented at 20m/s intervals, covering a range of + -500 m/s, i.e. V 1 =-500m/s，V 2 =-480m/s…V 51 =500 m/s, generating 51 matched filter coefficients;
the first pulse is matched filtered with the 51 matched filter coefficients;
constructing a two-dimensional matrix, drawing an amplitude matrix by matlab, and obtaining a matrix with a matrix width of V 46 The amplitude corresponding to =400 m/s is the largest, i.e. the target coarse speed is 400m/s;
dividing the wavelength by 0.05m and the pulse repetition period by 1000us and then by 2 gives a blur speed=25 m/s;
the 46 th filter coefficient is used for matching and filtering other pulses, FFT transformation is carried out in a slow time dimension, and the target accurate velocity is obtained to be 5.02m/s;
substituting the rough speed measurement result and the refined speed measurement result into a formulaAnd calculating to obtain the target speed of 405.02m/s, wherein the target speed is within an error allowable range compared with the actual speed.
The foregoing is illustrative of the present invention and is not to be construed as limiting thereof, but rather as being included within the spirit and scope of the present invention.
Claims (4)
1. A low repetition rate solution speed ambiguity method, comprising:
step one: generating a group of pulse signals S ti (N), where n=1, 2 … N, i=1, 2 … K, makes each transmit waveform identical and coherent, i.e. S t1 （n）= S t2 （n）= …=S tK （n）=S t （n）；
Step two: will S ti (n) up-conversion, transmitting by antenna, returning when meeting target, receiving by antenna, mixing to intermediate frequency, sampling pulse by pulse to obtain sampling signal as S ri (M), wherein M = 1,2 … M;
step three: dividing the object to be searched into Q sections along the speed interval, and setting the average speed of each section to be V 1 、V 2 …V Q ；
Step four: based on the average speed of each segment, the formula is usedCalculate Q matchesMatched filter coefficient MF q (n), where pi is the circumference ratio, V q Is the center speed of the Q-th speed segment, Δt is the sampling interval, λ is the wavelength, q=1, 2 … Q;
step five: using Q matched filter coefficients, using the formulaRespectively carrying out matched filtering on the first pulse, wherein p=1, 2 … M-N+1, and the sign represents taking conjugate;
step six: taking A q (p) absolute values forming a two-dimensional matrixCFAR detection is carried out on the matrix, if no target exists, detection is continued, and if the target exists, the Q < th > corresponding to the target is recorded T Line N T Column, obtain coarse speed measurement result V QT ；
Step seven: the pulse repetition period divided by the wavelength by 2 times gives the blur velocity V ams According to the rough speed measurement result, use the Q T Matched filter coefficient MF corresponding to line QT (n) matching and filtering other pulses, and measuring the speed in the FFT in the slow time dimension to obtain a refined speed measurement result V JT And adding the rough speed measurement result and the refined speed measurement result to obtain the target speed.
2. The low repetition rate resolution speed ambiguity resolution method as in claim 1, wherein said step one, comprises: selecting S with spike blur function feature ti （n）。
3. The low repetition rate resolution and blur method of claim 1 wherein steps two, three and four are performed simultaneously.
4. The low repetition rate resolution speed ambiguity resolution method of claim 1 further comprising the step eight: rounding function round (), substituting the rough speed measurement result and the refined speed measurement result into a formulaAnd calculating, updating the target speed and improving the accuracy.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310833476.2A CN116559853B (en) | 2023-07-10 | 2023-07-10 | Low repetition frequency solution speed blurring method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310833476.2A CN116559853B (en) | 2023-07-10 | 2023-07-10 | Low repetition frequency solution speed blurring method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116559853A CN116559853A (en) | 2023-08-08 |
| CN116559853B true CN116559853B (en) | 2023-09-12 |
Family
ID=87486495
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310833476.2A Active CN116559853B (en) | 2023-07-10 | 2023-07-10 | Low repetition frequency solution speed blurring method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116559853B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2047813A1 (en) * | 1990-08-29 | 1992-03-01 | Kapriel V. Krikorian | Acceleration compensation by matched filtering |
| US8044846B1 (en) * | 2007-11-29 | 2011-10-25 | Lockheed Martin Corporation | Method for deblurring radar range-doppler images |
| CN104239683A (en) * | 2014-07-23 | 2014-12-24 | 西安空间无线电技术研究所 | Resolving distance-velocity ambiguity method on basis of changing signal rate |
| CN106646446A (en) * | 2017-03-15 | 2017-05-10 | 中国人民解放军国防科学技术大学 | Detection method for moving target of pulse compression frequency-agile radar |
| CN110208786A (en) * | 2019-05-29 | 2019-09-06 | 西安空间无线电技术研究所 | A kind of two repetition ambiguity solution method of space based radar |
| CN110673105A (en) * | 2019-09-25 | 2020-01-10 | 武汉滨湖电子有限责任公司 | Method for resolving velocity ambiguity of pulse Doppler radar |
| CN112834992A (en) * | 2021-01-14 | 2021-05-25 | 中山大学 | Signal processing method and device of pulse Doppler radar and storage medium |
-
2023
- 2023-07-10 CN CN202310833476.2A patent/CN116559853B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2047813A1 (en) * | 1990-08-29 | 1992-03-01 | Kapriel V. Krikorian | Acceleration compensation by matched filtering |
| US5113194A (en) * | 1990-08-29 | 1992-05-12 | Hughes Aircraft Company | Acceleration compensation by matched filtering |
| US8044846B1 (en) * | 2007-11-29 | 2011-10-25 | Lockheed Martin Corporation | Method for deblurring radar range-doppler images |
| CN104239683A (en) * | 2014-07-23 | 2014-12-24 | 西安空间无线电技术研究所 | Resolving distance-velocity ambiguity method on basis of changing signal rate |
| CN106646446A (en) * | 2017-03-15 | 2017-05-10 | 中国人民解放军国防科学技术大学 | Detection method for moving target of pulse compression frequency-agile radar |
| CN110208786A (en) * | 2019-05-29 | 2019-09-06 | 西安空间无线电技术研究所 | A kind of two repetition ambiguity solution method of space based radar |
| CN110673105A (en) * | 2019-09-25 | 2020-01-10 | 武汉滨湖电子有限责任公司 | Method for resolving velocity ambiguity of pulse Doppler radar |
| CN112834992A (en) * | 2021-01-14 | 2021-05-25 | 中山大学 | Signal processing method and device of pulse Doppler radar and storage medium |
Non-Patent Citations (1)
| Title |
|---|
| 低重频雷达回波准最大似然速度估计;金胜等;《雷达科学与技术》;第129-136页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116559853A (en) | 2023-08-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107861117B (en) | Multi-target parameter measuring method suitable for continuous wave perimeter surveillance radar | |
| CN110412558B (en) | Method for resolving speed ambiguity of vehicle-mounted FMCW radar based on TDM MIMO | |
| CN106443615B (en) | A kind of bistatic across range gate method for locating speed measurement of MIMO radar high-speed target | |
| CN107561508B (en) | Coherent accumulation detection method for uniformly accelerated moving target | |
| CN108415010B (en) | Radar multi-target detection method based on trapezoidal LFMCW modulation | |
| CN111693983B (en) | Complementary waveform construction method and module, cognitive radar system and waveform emission method | |
| CN105629211B (en) | The signal processing method and automobile lane change auxiliary system of the combined waveform automobile lane change auxiliary system of multi-target detection | |
| CN109799488B (en) | Nonparametric search radar maneuvering target long-time coherent accumulation method | |
| CN111505618B (en) | Decoupling correction ranging method based on frequency estimation and suitable for field of vehicle-mounted millimeter wave radar | |
| CN109633622B (en) | Radar maneuvering target cross-unit accumulation method based on time reversal second-order Keystone transformation | |
| CN112882006B (en) | Millimeter wave short-range target detection method based on composite frequency modulation continuous wave | |
| CN104239683B (en) | Resolving distance-velocity ambiguity method on basis of changing signal rate | |
| CN107356908B (en) | Frequency agile signal coherent accumulation method | |
| CN108919249A (en) | A kind of radar target Joint estimation method based on two-dimentional local interpolation | |
| CN109613507B (en) | Detection method for high-order maneuvering target radar echo | |
| CN113093120A (en) | Method for estimating PRI agile radar target parameters based on capon algorithm | |
| CN111856466A (en) | Efficient ISAR (inverse synthetic aperture radar) translation compensation method for complex moving target | |
| CN110824439B (en) | Radar target rapid long-time coherent accumulation method | |
| CN111239721A (en) | Entropy-solving and speed-ambiguity-solving method for vehicle-mounted MIMO radar | |
| Zheng et al. | Radar detection and motion parameters estimation of maneuvering target based on the extended keystone transform (July 2018) | |
| CN116559853B (en) | Low repetition frequency solution speed blurring method | |
| CN108490425A (en) | A kind of angle-measuring method of bistatic MIMO radar | |
| CN115407279A (en) | Pulse compression radar range-Doppler coupling error compensation method | |
| CN115877381A (en) | Bistatic radar collaborative imaging method based on complementary random waveform | |
| CN112748434B (en) | Doppler weather radar speed de-ambiguity method based on polynomial fitting |
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 |