CN109581299B - Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method - Google Patents

Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method Download PDF

Info

Publication number
CN109581299B
CN109581299B CN201811511900.7A CN201811511900A CN109581299B CN 109581299 B CN109581299 B CN 109581299B CN 201811511900 A CN201811511900 A CN 201811511900A CN 109581299 B CN109581299 B CN 109581299B
Authority
CN
China
Prior art keywords
pulse compression
frequency
fourier transform
inverse fourier
continuous wave
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
Application number
CN201811511900.7A
Other languages
Chinese (zh)
Other versions
CN109581299A (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.)
Huanuo Xingkong Technology Co ltd
Hunan Huanuo Xingkong Electronic Technology Co ltd
Original Assignee
HUNAN NOVASKY ELECTRONIC 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 HUNAN NOVASKY ELECTRONIC TECHNOLOGY CO LTD filed Critical HUNAN NOVASKY ELECTRONIC TECHNOLOGY CO LTD
Priority to CN201811511900.7A priority Critical patent/CN109581299B/en
Publication of CN109581299A publication Critical patent/CN109581299A/en
Application granted granted Critical
Publication of CN109581299B publication Critical patent/CN109581299B/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
    • 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/2813Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays

Landscapes

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

Abstract

The invention discloses a method for suppressing side lobes of ultra-wideband step-frequency continuous wave pulse compression, which comprises the following steps: s1: determining the number of discrete inverse Fourier transform points in pulse compression; s2: calculating discrete inverse Fourier transform of each frequency point of the step frequency continuous wave; s3: calculating the power sum of discrete inverse Fourier transform of each frequency point of the step frequency continuous wave; s4: all frequency points of the step frequency continuous wave are combined, and pulse compression is completed through discrete inverse Fourier transform; s5: calculating the power of the pulse compression of the step frequency continuous wave; dividing the power of pulse compression by the sum of the power of discrete inverse Fourier transform of each frequency point to construct a sidelobe suppression coefficient; and multiplying the sidelobe suppression coefficient by the step frequency continuous wave pulse compression to complete sidelobe suppression. The method has the advantages of small calculated amount, convenience, flexibility, good side lobe suppression effect, no main lobe broadening effect and the like.

Description

Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method
Technical Field
The invention mainly relates to the technical field of radar signal processing, in particular to an ultra wide band step frequency continuous wave pulse compression side lobe suppression method.
Background
The method includes the steps of gradually transmitting a single Frequency point signal by using a Stepped Frequency Continuous Wave (SFCW) according to fixed Frequency steps, obtaining a large bandwidth, and synthesizing an Ultra-Wideband (UWB) signal. The SFCW signal can be considered as a frequency domain signal, and usually an echo signal is transformed into a time domain by using Inverse Fast Fourier Transform (IFFT), so that pulse compression is performed and high range resolution is achieved.
Compared with frequency modulation continuous wave signals, the SFCW signals are more flexible in transmission, can effectively avoid interference frequency points, are discrete in frequency spectrum and are beneficial to digital signal processing. The SFCW signal achieves higher average transmit power and enables coherent processing relative to impulse UWB signals of very short duration.
However, the distance direction obtained by pulse compression of the SFCW signal is high-resolution (approximately represented by a sinc function), and the SFCW signal often has high side lobe level. The commonly used windowing sidelobe suppression method can expand a main lobe while suppressing sidelobes, so that the distance of the SFCW signal is reduced towards the resolution.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: aiming at the technical problems in the prior art, the invention provides the ultra wide band step frequency continuous wave pulse compression side lobe suppression method which is small in calculated amount, convenient and flexible, good in side lobe suppression effect and free of main lobe broadening effect.
In order to solve the technical problems, the invention adopts the following technical scheme:
an ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method comprises the following steps:
step S1: determining the number of discrete inverse Fourier transform points in pulse compression;
step S2: calculating discrete inverse Fourier transform of each frequency point of the step frequency continuous wave;
step S3: calculating the power sum of discrete inverse Fourier transform of each frequency point of the step frequency continuous wave;
step S4: all frequency points of the step frequency continuous wave are combined, and pulse compression is completed through discrete inverse Fourier transform;
step S5: calculating the power of the pulse compression of the step frequency continuous wave; dividing the power of pulse compression by the sum of the power of discrete inverse Fourier transform of each frequency point to construct a sidelobe suppression coefficient; and multiplying the sidelobe suppression coefficient by the step frequency continuous wave pulse compression to complete sidelobe suppression.
As a further improvement of the invention: in step S1, the SFCW signal is composed of K frequency points, and the length of the discrete inverse fourier transform is set to N.
As a further improvement of the invention: in step S2, the discrete inverse fourier transform of each frequency point of the received SFCW signal is successively calculated, and the power sum of the discrete inverse fourier transform of each frequency point is calculated; for the k-th frequency point signal xkThe discrete inverse fourier transform calculation formula is as follows:
Figure BDA0001901005090000021
the power and calculation formula of all frequency point discrete inverse Fourier transform are as follows:
Figure BDA0001901005090000022
wherein XkRepresenting the k-th frequency point signal xkDiscrete inverse Fourier transform of (P)sAnd the power sum of all frequency point discrete inverse Fourier transform is shown.
As a further improvement of the invention: in step S3, all frequency points are combined together according to the following discrete inverse fourier transform formula:
Figure BDA0001901005090000023
realizing pulse compression of all frequency points, wherein RP represents pulse compression results of all frequency points, and j represents unit imaginary number; and calculating the power of the pulse compression according to the following formula:
P(i)=|RP(i)|2,i=1,2,…,N
wherein P represents the power of the pulse compression of all frequency bins.
As a further improvement of the invention: in step S4, the power of pulse compression is divided by the sum of the powers of discrete inverse fourier transform of each frequency point, and a sidelobe suppression coefficient is constructed, where the formula is as follows:
Figure BDA0001901005090000031
where w represents the sidelobe suppression coefficient.
As a further improvement of the invention: in step S5, the pulse compression result is multiplied by a suppression coefficient to complete the side lobe suppression with high resolution in the distance direction, and the formula is as follows:
S(i)=w(i)RP(i),i=1,2,…,N,
where S represents the pulse compression result after sidelobe suppression.
Compared with the prior art, the invention has the advantages that:
the ultra-wideband step frequency continuous wave pulse compression side lobe suppression method is suitable for side lobe suppression without broadening of a main lobe of an ultra-wideband SFCW signal, has the advantages of small calculated amount, convenience, flexibility and the like, and has good side lobe suppression effect and no main lobe broadening effect.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention.
FIG. 2 is a diagram showing the results of a simulation comparison experiment in a specific application example of the present invention.
Fig. 3 is an enlarged schematic view of the main lobe region of the experimental result of fig. 2.
Detailed Description
The invention will be described in further detail below with reference to the drawings and specific examples.
As shown in fig. 1, the ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method of the present invention includes the steps of:
step S1: determining the number of discrete inverse Fourier transform points in pulse compression;
assuming that the SFCW signal consists of K frequency points, setting the length of discrete inverse Fourier transform to be N;
step S2: calculating discrete inverse Fourier transform of each frequency point of the step frequency continuous wave;
and successively calculating the discrete inverse Fourier transform of each frequency point of the received SFCW signal, and calculating the power sum of the discrete inverse Fourier transform of each frequency point. For the k-th frequency point signal xkThe discrete inverse fourier transform calculation formula is as follows:
Figure BDA0001901005090000032
the power and calculation formula of all frequency point discrete inverse Fourier transform are as follows:
Figure BDA0001901005090000041
wherein XkRepresenting the k-th frequency point signal xkDiscrete inverse Fourier transform of (P)sAnd the power sum of all frequency point discrete inverse Fourier transform is shown.
Step S3: calculating the power sum of discrete inverse Fourier transform of each frequency point of the step frequency continuous wave;
all frequency points are combined together according to the following discrete inverse Fourier transform formula:
Figure BDA0001901005090000042
and realizing the pulse compression of all frequency points, wherein RP represents the pulse compression result of all frequency points, and j represents unit imaginary number. And calculating the power of the pulse compression according to the following formula:
P(i)=|RP(i)|2,i=1,2,…,N (4)
wherein P represents the power of the pulse compression of all frequency bins.
Step S4: all frequency points of the step frequency continuous wave are combined, and pulse compression is completed through discrete inverse Fourier transform;
dividing the power of pulse compression by the power sum of discrete inverse Fourier transform of each frequency point to construct a sidelobe suppression coefficient, wherein the formula is as follows:
Figure BDA0001901005090000043
where w represents the sidelobe suppression coefficient.
Step S5: calculating the power of the pulse compression of the step frequency continuous wave; dividing the power of pulse compression by the sum of the power of discrete inverse Fourier transform of each frequency point to construct a sidelobe suppression coefficient; and multiplying the sidelobe suppression coefficient by the step frequency continuous wave pulse compression to complete sidelobe suppression.
And multiplying the pulse compression result by a suppression coefficient to complete the side lobe suppression in the distance direction with high resolution, wherein the formula is as follows:
S(i)=w(i)RP(i),i=1,2,…,N (6)
where S represents the pulse compression result after sidelobe suppression.
In a specific application example, the specific operation flow is as follows:
step S100: assuming that the SFCW signal consists of K frequency bins, the length of the discrete inverse fourier transform is set to N.
In the simulation experiment of the invention, the length of the discrete inverse Fourier transform is set to be N4096. Table 1 below gives the parameters of the SFCW signal in the simulation experiment.
TABLE 1 parameters of SFCW signals
Starting frequency Frequency stepping Number of frequency points K
0.5GHz 2MHz 600
Step S200: and successively calculating the discrete inverse Fourier transform of each frequency point of the received SFCW signal, and calculating the power sum of the discrete inverse Fourier transform of each frequency point. For the k-th frequency point signal xkThe discrete inverse fourier transform calculation formula is as follows:
Figure BDA0001901005090000051
the IFFT power and the calculation formula of all frequency points are as follows;
Figure BDA0001901005090000052
wherein XkRepresenting the k-th frequency point signal xkDiscrete inverse Fourier transform of (P)sAnd the power sum of all frequency point discrete inverse Fourier transform is shown.
Step S300: all frequency points are combined together according to the following discrete inverse Fourier transform formula:
Figure BDA0001901005090000053
and realizing the pulse compression of all frequency points, wherein RP represents the pulse compression result of all frequency points, and j represents unit imaginary number. The power of the pulse compression is calculated according to the following formula:
P(i)=|RP(i)|2,i=1,2,…,N
as can be seen from the simulation experiment results of fig. 2 and fig. 3, the energy level of the side lobe is-13 dB, which is very high, in the original pulse compression result without side lobe suppression.
Step S400: dividing the power of pulse compression by the power sum of discrete inverse Fourier transform of each frequency point to construct a sidelobe suppression coefficient, wherein the formula is as follows:
Figure BDA0001901005090000054
where w represents the sidelobe suppression coefficient.
Step S500: and multiplying the pulse compression result by a suppression coefficient to complete the side lobe suppression in the distance direction with high resolution, wherein the formula is as follows:
S(i)=w(i)RP(i),i=1,2,…,N
where S represents the pulse compression result after sidelobe suppression.
As can be seen from the simulation experiment result of FIG. 2, after the Hamming window is added, the side lobe level of the pulse compression is reduced, but the main lobe is expanded. With the method of the present invention, the main lobe width is not widened but narrowed instead while the level of the side lobe is suppressed, which contributes to the improvement of the distance-to-resolution.
The above is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above-mentioned embodiments, and all technical solutions belonging to the idea of the present invention belong to the protection scope of the present invention. It should be noted that modifications and embellishments within the scope of the invention may be made by those skilled in the art without departing from the principle of the invention.

Claims (5)

1. An ultra-wideband step frequency continuous wave pulse compression sidelobe suppression method is characterized by comprising the following steps:
step S1: determining the number of discrete inverse Fourier transform points in pulse compression;
step S2: calculating discrete inverse Fourier transform of each frequency point of the step frequency continuous wave;
step S3: calculating the power sum of discrete inverse Fourier transform of each frequency point of the step frequency continuous wave;
step S4: all frequency points of the step frequency continuous wave are combined, and pulse compression is completed through discrete inverse Fourier transform;
step S5: calculating the power of the pulse compression of the step frequency continuous wave; dividing the power of pulse compression by the sum of the power of discrete inverse Fourier transform of each frequency point to construct a sidelobe suppression coefficient; the formula is as follows:
Figure FDA0003513230540000011
wherein w represents a sidelobe suppression coefficient; p represents the power of all frequency point pulse compression; psThe power sum of all frequency point discrete inverse Fourier transform is represented; n is the length of discrete inverse Fourier transform; and multiplying the sidelobe suppression coefficient by the step frequency continuous wave pulse compression to complete sidelobe suppression.
2. The method for suppressing sidelobes of ultra-wideband step-frequency continuous wave pulse compression as claimed in claim 1, wherein the step-frequency continuous wave SFCW signal in step S1 is composed of K bins, and the length of the discrete inverse fourier transform is set to N.
3. The ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method according to claim 2, wherein in step S2, discrete inverse fourier transform of each frequency point of the received SFCW signal is successively calculated; for the signal x of the k-th frequency pointkThe discrete inverse fourier transform calculation formula is as follows:
Figure FDA0003513230540000012
and calculating the power sum of the discrete inverse Fourier transform of each frequency point, wherein the power sum calculation formula of the discrete inverse Fourier transform of all the frequency points is as follows:
Figure FDA0003513230540000013
wherein XkRepresenting the k-th frequency point signal xkDiscrete inverse Fourier transform of (P)sAnd the power sum of all frequency point discrete inverse Fourier transform is shown.
4. The ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method according to claim 3, wherein in step S4, all frequency points are combined together according to the following discrete inverse Fourier transform formula:
Figure FDA0003513230540000021
realizing pulse compression of all frequency points, wherein RP represents pulse compression results of all frequency points, and j represents unit imaginary number; and calculating the power of the pulse compression according to the following formula:
P(i)=|RP(i)|2, i=1,2,…,N
wherein P represents the power of the pulse compression of all frequency bins.
5. The ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method according to claim 4, wherein in step S5, the pulse compression result is multiplied by a suppression coefficient to complete sidelobe suppression in a distance direction with high resolution, and the formula is as follows:
S(i)=w(i)RP(i), i=1,2,…,N
where S represents the pulse compression result after sidelobe suppression.
CN201811511900.7A 2018-12-11 2018-12-11 Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method Active CN109581299B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811511900.7A CN109581299B (en) 2018-12-11 2018-12-11 Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811511900.7A CN109581299B (en) 2018-12-11 2018-12-11 Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method

Publications (2)

Publication Number Publication Date
CN109581299A CN109581299A (en) 2019-04-05
CN109581299B true CN109581299B (en) 2022-04-22

Family

ID=65928129

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811511900.7A Active CN109581299B (en) 2018-12-11 2018-12-11 Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method

Country Status (1)

Country Link
CN (1) CN109581299B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007064671A (en) * 2005-08-29 2007-03-15 Toshiba Corp Radar system
CN101881822A (en) * 2010-06-07 2010-11-10 电子科技大学 Method for inhibiting same frequency interference of shared-spectrum radars
CN103138799A (en) * 2013-03-21 2013-06-05 哈尔滨工业大学 Modulation method of low sidelobe random frequency hopping pulse signal
CN103634256A (en) * 2012-08-24 2014-03-12 中兴通讯股份有限公司 Baseband signal sidelobe power suppression and sidelobe power suppression method and device
CN104777461A (en) * 2015-04-29 2015-07-15 中国科学院空间科学与应用研究中心 Method and system for generating broadband chaos radar signals with carrier frequency jumping randomly
CN104898132A (en) * 2015-05-14 2015-09-09 电子科技大学 Navigation anti-interference algorithm combining threshold processing and space-frequency adaptive algorithm
GB201709631D0 (en) * 2017-06-16 2017-08-02 Canon Kk Radar with stepped frequency complementrary code sequences

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007064671A (en) * 2005-08-29 2007-03-15 Toshiba Corp Radar system
CN101881822A (en) * 2010-06-07 2010-11-10 电子科技大学 Method for inhibiting same frequency interference of shared-spectrum radars
CN103634256A (en) * 2012-08-24 2014-03-12 中兴通讯股份有限公司 Baseband signal sidelobe power suppression and sidelobe power suppression method and device
CN103138799A (en) * 2013-03-21 2013-06-05 哈尔滨工业大学 Modulation method of low sidelobe random frequency hopping pulse signal
CN104777461A (en) * 2015-04-29 2015-07-15 中国科学院空间科学与应用研究中心 Method and system for generating broadband chaos radar signals with carrier frequency jumping randomly
CN104898132A (en) * 2015-05-14 2015-09-09 电子科技大学 Navigation anti-interference algorithm combining threshold processing and space-frequency adaptive algorithm
GB201709631D0 (en) * 2017-06-16 2017-08-02 Canon Kk Radar with stepped frequency complementrary code sequences

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Adaptive sidelobes reduction method for stepped frequency continuous waveform radar;Mei Yang 等;《IET International Confrence on Radar Systems(Radar 2012)》;20130408;正文全文 *
用于步进频连续波雷达旁瓣抑制的自适应CLEAN 技术;帅晓飞 等;《电讯技术》;20150228;第55卷(第2期);正文全文 *

Also Published As

Publication number Publication date
CN109581299A (en) 2019-04-05

Similar Documents

Publication Publication Date Title
US8571508B2 (en) Method and wideband antenna system to minimise the influence of interference sources
JP4421487B2 (en) Transmission signal generating device, signal processing device, and antenna device
JP5429699B2 (en) Compression coefficient generator
EP2491617B1 (en) An ultra-wideband radar imaging system using a two-dimensional multiple-input multiple output (MIMO) transducer array
CN113721216B (en) Target detection waveform optimization and processing method of agile coherent radar
CN113376601B (en) Frequency agile radar sidelobe suppression method based on CLEAN algorithm
CN110221268B (en) Phase-amplitude calibration filter and method for improving short pulse non-coherent radar signal processing performance
CN108037494A (en) A kind of radar target method for parameter estimation under impulse noise environment
US20160018512A1 (en) Method for Generating and Compressing Multi-Sweep-Frequency Radar Signals
JP2005085167A (en) Correlation processing device and method, and pulse compression processing device and method
CN113884992A (en) Self-adaptive anti-interference method of frequency agile radar
CN114325633B (en) Coherent processing method and device for dual-band radar
JP2017166918A (en) Signal processor, radar device, signal processing method, and program
CN109581299B (en) Ultra-wideband step-frequency continuous wave pulse compression sidelobe suppression method
CN109358327B (en) Iterative generation method of nonlinear frequency modulation signal, terminal and storage medium
JP4444057B2 (en) Pulse compression processor
CN114609596A (en) Echo characteristic accurate control interference method and system based on intermittent sampling
JP2006284241A (en) Filter apparatus
CN115825953B (en) Forward-looking super-resolution imaging method based on random frequency coding signal
CN112346019A (en) Coherent accumulation processing method for noise radar pulse waveform and low sidelobe pulse compression
Ryabova et al. Processing Multicarrier Phase Coded Signals with OFDM on the USRP Platform for NVIS Sounding of HF Radio Channels
CN115542304A (en) Self-interference signal suppression method based on transfer function fitting and phase matching
CN113238194B (en) Broadband phased array radar anti-decoy interference method based on fractional domain-frequency domain processing
CN110967686B (en) Phase coding signal generation method and system
CN115128563B (en) Amplitude-phase error compensation method for satellite-borne synthetic aperture radar system

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
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: Building B7, Lugu Enterprise Plaza, No. 27 Wenxuan Road, High tech Development Zone, Changsha City, Hunan Province, 410221

Patentee after: Huanuo Xingkong Technology Co.,Ltd.

Address before: Building B7, Lugu Enterprise Plaza, No. 27 Wenxuan Road, High tech Development Zone, Changsha City, Hunan Province, 410221

Patentee before: Hunan Huanuo Xingkong Electronic Technology Co.,Ltd.

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Building B7, Lugu Enterprise Plaza, No. 27 Wenxuan Road, High tech Development Zone, Changsha City, Hunan Province, 410221

Patentee after: Hunan Huanuo Xingkong Electronic Technology Co.,Ltd.

Address before: Building B7, Lugu enterprise Plaza, 27 Wenxuan Road, high tech Zone, Changsha City, Hunan Province, 410205

Patentee before: HUNAN NOVASKY ELECTRONIC TECHNOLOGY Co.,Ltd.