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 PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
- G01S13/26—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave
- G01S13/28—Systems 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/282—Systems 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
-
- 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/66—Radar-tracking systems; Analogous systems
-
- 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- 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/28—Details of pulse systems
- G01S7/285—Receivers
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
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:
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,representing the chirp rate, TpIs the sub-pulse time width, t0For the gate delay, rect (-) is a rectangular pulse function,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:
the result of matched filtering is:
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];
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.
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:
The echoes received by the radar can be expressed as:
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 rateTpIs the sub-pulse time width, TrRepresenting the sub-pulse repetition time, rect (-) is a rectangular pulse function,
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:
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:
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 rateTpIs the sub-pulse time width, rect (-) is a rectangular pulse function,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:
the result of matched filtering is:
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],
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),
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:
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,representing the chirp rate, TpIs the sub-pulse time width, t0For the gate delay, rect (-) is a rectangular pulse function,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:
the result of matched filtering is:
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];
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.
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)
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)
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 |
-
2019
- 2019-12-10 CN CN201911254097.8A patent/CN111220974B/en active Active
Patent Citations (6)
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)
Title |
---|
TACHE, M等: "Reliability and performance of optimised Schmitt trigger gates", 《THE JOURNAL OF ENGINEERING》 * |
景国彬 等: "一种新的步进频MIMO-SAR带宽合成的处理方法", 《西安电子科技大学学报》 * |
纪赛飞: "Wi-Fi外辐射源雷达多频信号合成方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (1)
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 |