CN111220974A - Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals - Google Patents

Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals Download PDF

Info

Publication number
CN111220974A
CN111220974A CN201911254097.8A CN201911254097A CN111220974A CN 111220974 A CN111220974 A CN 111220974A CN 201911254097 A CN201911254097 A CN 201911254097A CN 111220974 A CN111220974 A CN 111220974A
Authority
CN
China
Prior art keywords
pulse
frequency
sub
spectrum
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.)
Granted
Application number
CN201911254097.8A
Other languages
Chinese (zh)
Other versions
CN111220974B (en
Inventor
骆雨桐
齐少宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Ningyuan Electronic And Electrical Technology Co Ltd
Original Assignee
Xi'an Ningyuan Electronic And Electrical Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xi'an Ningyuan Electronic And Electrical Technology Co Ltd filed Critical Xi'an Ningyuan Electronic And Electrical Technology Co Ltd
Priority to CN201911254097.8A priority Critical patent/CN111220974B/en
Publication of CN111220974A publication Critical patent/CN111220974A/en
Application granted granted Critical
Publication of CN111220974B publication Critical patent/CN111220974B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/08Systems for measuring distance only
    • G01S13/10Systems for measuring distance only using transmission of interrupted, pulse modulated waves
    • G01S13/26Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave
    • G01S13/28Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave with time compression of received pulses
    • G01S13/282Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave with time compression of received pulses using a frequency modulated carrier wave
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/66Radar-tracking systems; Analogous systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to a low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals, which adopts a frequency spectrum serial connection mode to perform pulse compression on a total signal spectrum and obtains a high-resolution one-dimensional range profile through inverse Fourier transform. Compared with the existing frequency domain splicing method, the method of segmented up-sampling is adopted, and the sub-pulses are grouped and spliced in pairs, so that the computational complexity in the splicing process is reduced, the processing speed of synthesizing a broadband signal by using a frequency modulation stepping pulse signal is increased, and the real-time imaging of a tracking target by a radar is facilitated.

Description

Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals
Technical Field
The invention belongs to the field of radar signal processing, and particularly relates to a low-complexity frequency domain splicing algorithm based on frequency modulation stepping signals, which can be used for real-time high-resolution imaging and monitoring of a tracking target by a radar.
Background
The method of pulse synthesis is often used to obtain high resolution one-dimensional range profile, and the currently commonly used pulse synthesis methods include: a synthetic distance envelope method, a frequency domain splicing method and a time domain splicing method. The comparison result between the operand and the performance of the 3 algorithms is: the time domain splicing method is inferior to the other two methods in terms of both computational load and performance. The synthetic distance envelope method has the minimum computation amount, but has the problem of energy leakage, and can cause image blurring. The frequency domain splicing algorithm has the best imaging quality, but the operation amount is larger than that of a synthetic envelope method, and the engineering application still has certain limitation. Therefore, if the computation of the frequency domain splicing method can be further reduced and the processing speed of the frequency domain splicing method can be increased, the frequency domain splicing method can have more engineering practicability. Based on the point, the frequency domain splicing method is improved.
The frequency domain stitching algorithm was first proposed in 1998 to perform faster than time domain processing since no resampling of the narrowband pulses is required. However, when the spliced spectrum signal is pulse-compressed, since the spliced spectrum has a certain difference from the ideal broadband signal spectrum, it becomes difficult to construct a suitable matched filter, and the implementation process is inconvenient. In 2013, Nel, W, Tait, J, Lord, R, & Wilkinson and Tao all adopt a method of compressing a sub-pulse and then splicing frequency spectrums, and a matched filter is universal to each sub-pulse, does not need special construction, and is more convenient than a traditional algorithm. And a rapid algorithm is provided on the basis, so that the frequency domain splicing algorithm is more beneficial to engineering realization. However, in order to accurately splice, the position of the center frequency point of the sub-pulse and the position of the center frequency point in the spread spectrum need to be additionally calculated, and a position deviation may be generated in the process, so that the imaging quality is affected.
Disclosure of Invention
Technical problem to be solved
The method aims to solve the problem of difficulty in real-time imaging caused by large data volume under the requirement of high-resolution imaging. The invention provides a low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals, which adopts a frequency spectrum serial connection mode to perform pulse compression on a total signal spectrum and obtains a high-resolution one-dimensional range profile through inverse Fourier transform.
Technical scheme
A low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals is characterized by comprising the following steps:
step 1: echo sub-pulse frequency domain matched filtering
Performing frequency mixing processing on echoes received by a radar to obtain a baseband signal, wherein the sampling rate is set to fs, and the number of sub-pulse sampling points is set to M; and (3) performing inverse Fourier transform on each sub-pulse baseband echo, namely performing M-point FFT to obtain a corresponding signal spectrum:
Figure BDA0002309823150000021
wherein i represents the ith sub-pulse, N represents N sub-pulses in a group of frequency modulation stepping pulse signals, and BsRepresenting the sub-pulse bandwidth, f0Representing the carrier frequency, af the frequency step, tau the target delay,
Figure BDA0002309823150000022
representing the chirp rate, TpIs the sub-pulse time width, t0For the gate delay, rect (-) is a rectangular pulse function,
Figure BDA0002309823150000023
indicates a width of TpThe rectangular pulse of (2);
and performing frequency domain matched filtering on each sub-pulse echo signal, wherein the frequency spectrum of the ideal matched filter is as follows:
Figure BDA0002309823150000024
the result of matched filtering is:
Figure BDA0002309823150000025
step 2: multiplying the matched and filtered signal by a gate delay factor to correct the phase error, wherein the gate delay factor is as follows: exp [ j2 pi (f)0+iΔf)t0];
Figure BDA0002309823150000031
And step 3: zero-filling up-sampling is carried out on the corrected sub-pulse frequency spectrum, so that the number of sampling points is enlarged to 2 times of the original number;
and 4, step 4: carrying out spectrum shifting on the spread spectrum, converting the IFFT of the spread spectrum into a time domain, multiplying the time-shifting factor by the IFFT, shifting the spectrum to a corresponding position, and then carrying out FFT (fast Fourier transform) conversion to a frequency domain;
and 5: adjacent two shifted frequency domains are combined into a group by superposition, and whether only one group of frequency spectrums exists after superposition is judged; if not, repeating the steps 2-5 until the superposition result has only one group of frequency spectrums, namely the required spliced total spectrum;
step 6: and performing IFFT on the finally superposed spliced total spectrum to obtain a high-resolution one-dimensional range profile.
The time shift factor in step (4) is a (i, t) ═ exp (j2 pi m Δ ft),
Figure BDA0002309823150000032
advantageous effects
Compared with the existing frequency domain splicing algorithm, the low-complexity frequency modulation stepping pulse signal-based frequency domain splicing method adopts a segmented up-sampling method to splice sub-pulses in pairs and groups, thereby reducing the computational complexity in the splicing process, improving the processing speed of synthesizing a broadband signal by using the frequency modulation stepping pulse signal, and being beneficial to real-time imaging of a tracking target by a radar.
Drawings
FIG. 1 is a flow chart of the algorithm implementation of the present invention
FIG. 2 is a diagram of a scattering point model of an object
FIG. 3 is a graph of single sub-pulse imaging results
FIG. 4 is a simulation diagram of a conventional frequency domain splicing algorithm
FIG. 5 is a simulation diagram of an improved frequency domain stitching algorithm
FIG. 6 is a complexity contrast plot
Detailed Description
The invention will now be further described with reference to the following examples and drawings:
referring to fig. 1, the specific implementation steps of this example are as follows:
step 1, echo sub-pulse frequency domain matching filtering
The echoes received by the radar can be expressed as:
Figure BDA0002309823150000041
i represents the ith sub-pulse, N represents N sub-pulses in a group of frequency modulated step pulse signals, f0Represents the carrier frequency,. DELTA.f represents the frequency step, τ represents the target time delay, and k represents the chirp rate
Figure BDA0002309823150000042
TpIs the sub-pulse time width, TrRepresenting the sub-pulse repetition time, rect (-) is a rectangular pulse function,
Figure BDA0002309823150000043
mixing the echo and the local oscillator signal to obtain a baseband signal, wherein the expression of the ith sub-pulse after mixing is as follows:
Figure BDA0002309823150000044
performing frequency mixing processing on echoes received by a radar to obtain a baseband signal, wherein the sampling rate is set to fs, and the number of sub-pulse sampling points is set to M (M is 2)q). And (3) performing inverse Fourier transform on each sub-pulse baseband echo, namely performing M-point FFT to obtain a corresponding signal spectrum:
Figure BDA0002309823150000045
i represents the ith sub-pulse, N represents N sub-pulses in a group of frequency modulation stepping pulse signals, BsRepresenting the sub-pulse bandwidth, f0Represents the carrier frequency,. DELTA.f represents the frequency step, τ represents the target time delay, and k represents the chirp rate
Figure BDA0002309823150000046
TpIs the sub-pulse time width, rect (-) is a rectangular pulse function,
Figure BDA0002309823150000047
indicates a width of TpThe square pulse of (2).
And performing frequency domain matched filtering on each sub-pulse echo signal, wherein the frequency spectrum of the ideal matched filter is as follows:
Figure BDA0002309823150000051
the result of matched filtering is:
Figure BDA0002309823150000052
and 2, multiplying the matched and filtered signal by a gate delay factor to correct the phase error. The gate delay factor is: exp [ j2 pi (f)0+iΔf)t0],
Figure BDA0002309823150000053
And 3, carrying out zero filling and up-sampling on the corrected sub-pulse frequency spectrum, so that the number of sampling points is enlarged to 2 times of the original number.
And 4, carrying out spectrum shifting on the expanded spectrum. And after the spectrum after being spread is IFFT converted into a time domain, the time shift factor is multiplied, the spectrum is moved to a corresponding position, and then FFT is converted into a frequency domain. Constructing a time shift factor according to the carrier frequency of the sub-pulse: a (i, t) ═ exp (j2 pi m Δ ft),
Figure BDA0002309823150000054
and 5, overlapping and synthesizing adjacent shifted frequency domains into a group, and judging whether only one group of frequency spectrums exists after overlapping. If not, repeating the steps (2) to (5) until the superposition result has only one group of frequency spectrums, namely the required spliced total spectrum.
And 6, performing IFFT on the finally superposed spliced total spectrum to obtain a high-resolution one-dimensional range profile.
The effects of the invention can be further illustrated by simulation:
1. simulation conditions
The simulation was implemented using MATLAB R2015b software, and the number of sub-pulses in a group of pulse trains was set to 32, and the sub-pulse bandwidth was set to 60M. The sub-pulse time width was 40. mu.s, and the pulse repetition time was 300. mu.s. The longitudinal distances of scattering points on the target relative to the geometric center of the target are set to be three, namely 100m, 100.09m and 103 m. The resolution of the single sub-pulse is 2.5m, and the resolution of the processed composite pulse is 0.078 m.
2. Emulated content and results
The three scattering point models are shown in fig. 2, the imaging result of a single sub-pulse is shown in fig. 3, the simulation result of the traditional frequency domain stitching algorithm is shown in fig. 4, and the simulation result of the improved frequency domain stitching algorithm is shown in fig. 5. As can be seen from a comparison of fig. 2 and fig. 3, the single sub-pulse cannot distinguish between the two targets 100m and 100.09m, and can be distinguished after synthesis. As can be seen from comparison of fig. 4 and 5, the improved method achieves the same high resolution effect as the conventional method. Fig. 6 shows that the improved algorithm has lower computational complexity and higher computational efficiency, so the improved algorithm is superior to the traditional algorithm.

Claims (2)

1. A low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals is characterized by comprising the following steps:
step 1: echo sub-pulse frequency domain matched filtering
Performing frequency mixing processing on echoes received by a radar to obtain a baseband signal, wherein the sampling rate is set to fs, and the number of sub-pulse sampling points is set to M; and (3) performing inverse Fourier transform on each sub-pulse baseband echo, namely performing M-point FFT to obtain a corresponding signal spectrum:
Figure FDA0002309823140000011
wherein i represents the ith sub-pulse, N represents N sub-pulses in a group of frequency modulation stepping pulse signals, and BsRepresenting the sub-pulse bandwidth, f0Represents the carrier frequency,. DELTA.f represents the frequency step size,. tau represents the eyeThe time delay is marked by a time delay mark,
Figure FDA0002309823140000012
representing the chirp rate, TpIs the sub-pulse time width, t0For the gate delay, rect (-) is a rectangular pulse function,
Figure FDA0002309823140000013
indicates a width of TpThe rectangular pulse of (2);
and performing frequency domain matched filtering on each sub-pulse echo signal, wherein the frequency spectrum of the ideal matched filter is as follows:
Figure FDA0002309823140000014
the result of matched filtering is:
Figure FDA0002309823140000015
step 2: multiplying the matched and filtered signal by a gate delay factor to correct the phase error, wherein the gate delay factor is as follows: exp [ j2 pi (f)0+iΔf)t0];
Figure FDA0002309823140000016
And step 3: zero-filling up-sampling is carried out on the corrected sub-pulse frequency spectrum, so that the number of sampling points is enlarged to 2 times of the original number;
and 4, step 4: carrying out spectrum shifting on the spread spectrum, converting the IFFT of the spread spectrum into a time domain, multiplying the time-shifting factor by the IFFT, shifting the spectrum to a corresponding position, and then carrying out FFT (fast Fourier transform) conversion to a frequency domain;
and 5: adjacent two shifted frequency domains are combined into a group by superposition, and whether only one group of frequency spectrums exists after superposition is judged; if not, repeating the steps 2-5 until the superposition result has only one group of frequency spectrums, namely the required spliced total spectrum;
step 6: and performing IFFT on the finally superposed spliced total spectrum to obtain a high-resolution one-dimensional range profile.
2. A low complexity frequency domain concatenation algorithm based on FM stepped pulse signals according to claim 1, wherein said time shift factor in step (4) is A (i, t) ═ exp (j2 π m Δ ft),
Figure FDA0002309823140000021
CN201911254097.8A 2019-12-10 2019-12-10 Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals Active CN111220974B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911254097.8A CN111220974B (en) 2019-12-10 2019-12-10 Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911254097.8A CN111220974B (en) 2019-12-10 2019-12-10 Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals

Publications (2)

Publication Number Publication Date
CN111220974A true CN111220974A (en) 2020-06-02
CN111220974B CN111220974B (en) 2023-03-24

Family

ID=70808283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911254097.8A Active CN111220974B (en) 2019-12-10 2019-12-10 Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals

Country Status (1)

Country Link
CN (1) CN111220974B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112162255A (en) * 2020-09-29 2021-01-01 中国船舶重工集团公司第七二四研究所 Pulse splicing method based on intra-pulse characteristics

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101021561A (en) * 2007-04-06 2007-08-22 清华大学 Wide band rader utilizing multi-transmitting and multi-receiving frequency division signal and imaging method thereof
WO2012151783A1 (en) * 2011-05-12 2012-11-15 苏州全波通信技术有限公司 Double frequency conversion modulation system and frequency conversion method
CN103175893A (en) * 2013-02-04 2013-06-26 中国科学院声学研究所 Defect defecting method based on wide frequency band ultrasonic phase array
US20170111031A1 (en) * 2015-10-15 2017-04-20 Rolls-Royce Plc Method of performing real time decomposition of a signal into components
US20180182410A1 (en) * 2016-12-23 2018-06-28 Synaptics Incorporated Online dereverberation algorithm based on weighted prediction error for noisy time-varying environments
CN109738878A (en) * 2019-02-27 2019-05-10 电子科技大学 Radar range profile's recognition methods based on compressed sensing and frequency stepping waveform

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101021561A (en) * 2007-04-06 2007-08-22 清华大学 Wide band rader utilizing multi-transmitting and multi-receiving frequency division signal and imaging method thereof
WO2012151783A1 (en) * 2011-05-12 2012-11-15 苏州全波通信技术有限公司 Double frequency conversion modulation system and frequency conversion method
CN103175893A (en) * 2013-02-04 2013-06-26 中国科学院声学研究所 Defect defecting method based on wide frequency band ultrasonic phase array
US20170111031A1 (en) * 2015-10-15 2017-04-20 Rolls-Royce Plc Method of performing real time decomposition of a signal into components
US20180182410A1 (en) * 2016-12-23 2018-06-28 Synaptics Incorporated Online dereverberation algorithm based on weighted prediction error for noisy time-varying environments
CN109738878A (en) * 2019-02-27 2019-05-10 电子科技大学 Radar range profile's recognition methods based on compressed sensing and frequency stepping waveform

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
TACHE, M等: "Reliability and performance of optimised Schmitt trigger gates", 《THE JOURNAL OF ENGINEERING》 *
景国彬 等: "一种新的步进频MIMO-SAR带宽合成的处理方法", 《西安电子科技大学学报》 *
纪赛飞: "Wi-Fi外辐射源雷达多频信号合成方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112162255A (en) * 2020-09-29 2021-01-01 中国船舶重工集团公司第七二四研究所 Pulse splicing method based on intra-pulse characteristics

Also Published As

Publication number Publication date
CN111220974B (en) 2023-03-24

Similar Documents

Publication Publication Date Title
AU2010340812B2 (en) Pulse radar range profile motion compensation
CN107247254B (en) A kind of NLFM signal goes oblique processing method
JP4834370B2 (en) Correlation reception processing device
US8022863B1 (en) Spatially variant apodization (SVA) filter
CN102176018B (en) Doppler wave beam sharpening rapid imaging method of mechanical scanning radar
CN106324597A (en) Translational motion compensation and imaging method for PFA-based large-turning-angle ISAR radar
CN107991655B (en) LFM-PC signal and fuzzy function optimization method thereof
CN110764067A (en) Fourier transform LFM-BPSK composite modulation radar signal parameter estimation method
CN111781595B (en) Complex maneuvering group target imaging method based on matching search and Doppler defuzzification
CN110308427A (en) LFM pulse radar frequency-domain impulse compression processing method based on FPGA
CN111551925A (en) Target speed estimation method of pulse-to-pulse frequency agility radar based on FFT
CN111220974B (en) Low-complexity frequency domain splicing method based on frequency modulation stepping pulse signals
CN108008387B (en) Downward-looking three-dimensional imaging method for airborne array antenna
CN111856462A (en) ISAR ship target imaging method based on MPT-MDPT
JP4444057B2 (en) Pulse compression processor
CN103412307A (en) CZT imaging method based on nonlinear frequency modulation signals
CN103954963A (en) Step frequency SAR imaging method based on CS algorithm
Wan et al. Ground moving target imaging based on MSOKT and KT for synthetic aperture radar
CN112698329B (en) Method for suppressing range ambiguity of spaceborne synthetic aperture radar
CN110308448B (en) Method for enhancing two-dimensional image of inverse synthetic aperture radar
CN106950567A (en) Ultra wide band based on high-order sub-aperture CS slides poly- SAR image processing methods
CN111427035A (en) Radar actual measurement data set expansion method based on sub-band micro-Doppler difference
CN112965067A (en) Target speed expanding method suitable for FMCW automobile radar
CN113567983B (en) Radar synthetic aperture sparse imaging method and system using slide rail
CN114509759B (en) Stepping frequency modulation pulse SAR moving target imaging 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