CN111999725B - Wideband polynomial phase signal declivity method and device under guidance of narrowband signal - Google Patents
Wideband polynomial phase signal declivity method and device under guidance of narrowband signal Download PDFInfo
- Publication number
- CN111999725B CN111999725B CN202010904054.6A CN202010904054A CN111999725B CN 111999725 B CN111999725 B CN 111999725B CN 202010904054 A CN202010904054 A CN 202010904054A CN 111999725 B CN111999725 B CN 111999725B
- Authority
- CN
- China
- Prior art keywords
- signal
- narrowband
- lfm
- coefficient
- term coefficient
- 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
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/0209—Systems with very large relative bandwidth, i.e. larger than 10 %, e.g. baseband, pulse, carrier-free, ultrawideband
-
- 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
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
-
- 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
- G01S7/32—Shaping echo pulse signals; Deriving non-pulse signals from echo pulse signals
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
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 wideband polynomial phase signal declivity method and a device under the guidance of a narrowband signal, wherein the method comprises the following steps: transmitting a plurality of continuous narrowband LFM signals; acquiring pulse echoes of a plurality of continuous narrowband LFM signals; performing delay correlation processing on each pulse echo to obtain a plurality of narrow-band echo delay signals; performing fractional Fourier transform to obtain a quadratic term coefficient and a cubic term coefficient of each narrowband echo delay signal estimation; compensating pulse echo of each narrowband LFM signal, and carrying out Fourier transform on the compensated result to obtain a first term coefficient and a zero term coefficient of each narrowband echo delay signal estimation; acquiring a cubic term coefficient of a broadband signal; obtaining k-1 times of coefficients of the broadband signal, and performing declivity processing on the broadband signal by using the declivity signal; the invention has the advantages that: the novel wideband polynomial phase signal declivity method is provided, and the declivity effect is good.
Description
Technical Field
The invention relates to the technical field of microwaves, in particular to a wideband polynomial phase signal declivity method and device under the guidance of narrowband signals.
Background
With the continuous development of microwave technology, the radar can transmit and receive signals with wider bandwidth than the prior art, and the radar can perform clearer imaging processing by means of high distance resolution of the signals with wider bandwidth and Doppler generated by rotation of a space target. The large bandwidth signal encounters a spatial target that moves at high speed and with acceleration, while the transmitted signal is a chirp signal, the received signal becomes a polynomial phase signal. Because the in-band fluctuation of the receiving and transmitting response of the antenna with large bandwidth is not ideal linear frequency modulation signal, the signal needs to be subjected to declivity treatment, and the inherent declivity signal is loaded on the broadband signal for declivity treatment in the traditional declivity treatment.
The invention discloses a method for generating broadband deskew echo based on sweep frequency data, belonging to the technical field of radar signal processing, wherein the method is disclosed in Chinese invention bulletin No. CN 110850384B. First, a pulsed radar transmits a chirp signal to a detection target containing Q scattering centers. And each scattering center in the detection target is fed back to the echo signal of each radar to obtain broadband echo signals of the Q scattering centers. Setting a reference signal, and performing declining treatment on broadband echo signals of each scattering center to obtain a frequency response containing the characteristic of a target RCS; carrying out quick Fourier inverse transformation on the frequency response containing the RCS characteristic of the target, wherein each scattering center inversely transforms to obtain a complex amplitude value, and Q complex amplitude values form a one-dimensional range profile of the detection target; and multiplying the one-dimensional range profile of the detection target with the deskewing echo signals of each scattering center, accumulating and summing to finally obtain the whole broadband deskewing echo of the detection target. The invention greatly reduces the operation amount and simultaneously ensures the same precision and accuracy as convolution operation. But the main emphasis is not the deskewing process, and the deskewing of the wideband polynomial phase signal under the guidance of the narrowband signal cannot be well performed.
Disclosure of Invention
The invention aims to solve the technical problem that the wideband polynomial phase signal under the guidance of the narrowband signal in the prior art has poor declivity effect.
The invention solves the technical problems by the following technical means: a method for wideband polynomial phase signal deskewing under narrowband signal steering, the method comprising:
step one: transmitting a plurality of continuous narrowband LFM signals;
step two: acquiring pulse echoes of a plurality of continuous narrowband LFM signals;
step three: performing delay correlation processing on each pulse echo to obtain a plurality of narrow-band echo delay signals;
step four: traversing all the narrow-band echo delay signals, taking values at fixed time, and performing fractional Fourier transform to obtain a quadratic term coefficient and a cubic term coefficient estimated by each narrow-band echo delay signal;
step five: according to the quadratic term coefficient and the cubic term coefficient of each narrowband echo delay signal estimation, compensating the pulse echo of each narrowband LFM signal, and carrying out Fourier transform on the compensated result to obtain the primary term coefficient and the zero-order term coefficient of each narrowband echo delay signal estimation;
step six: obtaining the third term coefficient of the broadband signal according to the third term coefficient and the frequency modulation slope of the ith narrowband echo delay signal and the frequency modulation slope of the broadband signal;
step seven: obtaining k-1 times term coefficients of the broadband signals according to the k-1 times term coefficients of the ith and jth narrowband echo delay signals, the frequency modulation slope and the frequency modulation slope of the broadband signals, wherein the value of k is 1,2, 3 or 4;
step eight: and obtaining the frequency modulation slope of the declivity signal according to the quadratic term coefficient and the cubic term coefficient of the broadband signal, and performing declivity processing on the broadband signal by using the declivity signal.
The method is based on a plurality of continuous narrowband LFM signals with different waveforms, realizes the extraction of waveform parameters of a space target broadband signal, does not need to add extra hardware cost, does not change the existing radar imaging signal processing system, has small calculation load, obtains the frequency modulation slope of the declivity signal through the quadratic term coefficient and the cubic term coefficient of the broadband signal, obtains the declivity signal through the frequency modulation slope of the declivity signal, performs declivity processing on the broadband signal by using the declivity signal, realizes accurate declivity real-time declivity aiming at the broadband signal, has good declivity effect, and meanwhile, the waveform parameters of the broadband signal have important application value to the radar accurate imaging field and the follow-up target identification.
Further, the ith narrowband LFM signal in the step one isWherein i=1, 2,..n, N is the total number of narrowband LFM signals, S i (t) is the ith narrowband LFM signal, j is the imaginary symbol, e () As an exponential function based on e, f 0 For the common working frequency point of all narrowband LFM signals, the narrowband LFM signals and the broadband signals share the same working frequency point, k i Frequency modulation slope of ith narrowband LFM signal, B i Is the amplitude of the i-th narrowband LFM signal.
Further, the pulse echo in the second step is expressed asWherein S is r i (t) is the ith pulse echo, alpha 1i 、α 2i 、α 3i And alpha 4i The zero order term coefficient, the first order term coefficient, the second order term coefficient and the third order term coefficient of the ith narrowband LFM signal respectively, A i Is the echo amplitude of the i-th narrowband LFM signal.
Further, the third step includes: delaying the ith pulse echo by tau components to obtain an ith narrow-band echo delay signal
Wherein (alpha) 3i +3α 4i t) is the instantaneous frequency of the ith narrowband echo delay signal.
Further, the fourth step includes: by the formulaObtaining a first intermediate quantity->And a second intermediate amount->Time values t1 and t2 can obtain the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantity
The third order term coefficient of the narrow-band echo delay signal estimation can be obtained according to the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantityAnd quadratic term coefficient->Respectively is
Further, the fifth step includes: according to the quadratic term coefficient and the cubic term coefficient estimated by each narrow-band echo delay signal, compensating the pulse echo of each narrow-band LFM signal to obtain a compensated result
For the compensated result, the result is calculated by the formulaPerforming Fourier transform to obtain a first order coefficient of each narrowband echo delay signal estimation>And zero order term coefficient->
Further, the sixth step includes: by the formulaAcquiring a cubic coefficient of a broadband signal, wherein alpha 4b Is the third order coefficient of the wideband signal, beta br Is the frequency modulation slope of the wideband signal,/->Is the third order coefficient, beta, of the estimated ith narrowband LFM signal i Is the chirp rate of the i-th narrowband LFM signal.
Further, the step seven includes: by the formulaObtaining the linear relation between the frequency modulation slope of the narrowband LFM signal and the k-1 th order coefficient, wherein beta 1i Frequency modulation slope beta of ith narrowband LFM signal 2j For the frequency modulation slope of the jth narrowband LFM signal, pa1 and pa2 are the first and second parameters, respectively For the k-1 th order term coefficients of the estimated i-th narrowband LFM signal,k-1 th order term coefficients for the estimated jth narrowband LFM signal;
by the formula alpha kb =β br Pa1+pa2 obtaining the k-1 th order coefficients of the wideband signal, where α kb Is the k-1 th order term coefficient of the wideband signal.
Further, the step eight includes:
obtaining expression Sigmatch of broadband signal according to zero order coefficient, primary order coefficient, secondary order coefficient and tertiary order coefficient of broadband signal b (t)=exp[j2π((α 1 +α 2b t+α 3b t 2 +α 4b t 3 ))]
Obtaining the frequency modulation slope alpha of the declassified signal according to the expression of the broadband signal 3b +α 4b And t, obtaining a declivity signal conjugated with the broadband signal according to the frequency modulation slope of the declivity signal, and performing declivity processing on the broadband signal by using the secondary term and the tertiary term of the declivity signal for counteracting the broadband signal.
The invention also provides a wideband polynomial phase signal declivity device under the guidance of a narrowband signal, which comprises:
the transmitting module is used for transmitting a plurality of continuous narrowband LFM signals;
the receiving module is used for acquiring pulse echoes of a plurality of continuous narrowband LFM signals;
the delay correlation processing module is used for carrying out delay correlation processing on each pulse echo to obtain a plurality of narrow-band echo delay signals;
the first Fourier transform module is used for traversing all the narrow-band echo delay signals, taking the value at fixed time, and carrying out fractional Fourier transform to obtain a quadratic term coefficient and a cubic term coefficient estimated by each narrow-band echo delay signal;
the second Fourier transform module is used for compensating the pulse echo of each narrowband LFM signal according to the quadratic term coefficient and the cubic term coefficient estimated by each narrowband echo delay signal, and carrying out Fourier transform on the compensated result to obtain the quadratic term coefficient and the zero-order term coefficient estimated by each narrowband echo delay signal;
the first broadband signal coefficient acquisition module is used for acquiring the cubic term coefficient of the broadband signal according to the cubic term coefficient and the frequency modulation slope of the ith narrowband echo delay signal and the frequency modulation slope of the broadband signal;
the second broadband signal coefficient acquisition module is used for acquiring k-1 times of coefficients of the broadband signals according to the k-1 times of coefficients of the ith and jth narrowband echo delay signals, the frequency modulation slope and the frequency modulation slope of the broadband signals, wherein the value of k is 1,2 or 3;
and the declivity module is used for obtaining the frequency modulation slope of the declivity signal according to the quadratic term coefficient and the cubic term coefficient of the broadband signal, and declivity processing is carried out on the broadband signal by utilizing the declivity signal.
Further, the ith narrowband LFM signal in the transmitting module isWherein i=1, 2,..n, N is the total number of narrowband LFM signals, S i (t) is the ith narrowband LFM signal, j is the imaginary symbol, e () As an exponential function based on e, f 0 For the common working frequency point of all narrowband LFM signals, the narrowband LFM signals and the broadband signals share the same working frequency point, k i Frequency modulation slope of ith narrowband LFM signal, B i Is the amplitude of the i-th narrowband LFM signal.
Further, the pulse echo in the receiving module is expressed asWherein S is r i (t) is the ith pulse echo, alpha 1i 、α 2i 、α 3i And alpha 4i The zero order term coefficient, the first order term coefficient, the second order term coefficient and the third order term coefficient of the ith narrowband LFM signal respectively, A i Is the echo amplitude of the i-th narrowband LFM signal.
Further, the delay related processing module is further configured to: delaying the ith pulse echo by tau components to obtain an ith narrow-band echo delay signal
Wherein (alpha) 3i +3α 4i t) is the instantaneous frequency of the ith narrowband echo delay signal.
Further, the first fourier transform module is further configured to: by the formulaObtaining a first intermediate quantity->And a second intermediate amount->Time values t1 and t2 can obtain the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantity
The third order term coefficient of the narrow-band echo delay signal estimation can be obtained according to the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantityAnd quadratic term coefficient->Respectively is
Further, the second fourier transform module is further configured to: according to the quadratic term coefficient and the cubic term coefficient estimated by each narrow-band echo delay signal, compensating the pulse echo of each narrow-band LFM signal to obtain a compensated result
For the compensated result, the result is calculated by the formulaPerforming Fourier transform to obtain a first order coefficient of each narrowband echo delay signal estimation>And zero order term coefficient->
Further, the first wideband signal coefficient acquisition module is further configured to: by the formulaAcquiring a cubic coefficient of a broadband signal, wherein alpha 4b Is the third order coefficient of the wideband signal, beta br Is the frequency modulation slope of the wideband signal,/->Is the third order coefficient, beta, of the estimated ith narrowband LFM signal i Is the chirp rate of the i-th narrowband LFM signal.
Further, the second wideband signal coefficient acquisition module is further configured to: by the formulaObtaining the linear relation between the frequency modulation slope of the narrowband LFM signal and the k-1 th order coefficient, wherein beta 1i Frequency modulation slope beta of ith narrowband LFM signal 2j For the frequency modulation slope of the jth narrowband LFM signal, pa1 and pa2 are the first and second parameters, respectively For the k-1 th order term coefficients of the estimated i-th narrowband LFM signal,k-1 th order term coefficients for the estimated jth narrowband LFM signal;
by the formula alpha kb =β br Pa1+pa2 obtaining the k-1 th order coefficients of the wideband signal, where α kb Is the k-1 th order term coefficient of the wideband signal.
Further, the deskewing module is further configured to:
obtaining expression Sigmatch of broadband signal according to zero order coefficient, primary order coefficient, secondary order coefficient and tertiary order coefficient of broadband signal b (t)=exp[j2π((α 1 +α 2b t+α 3b t 2 +α 4b t 3 ))]
Obtaining the frequency modulation slope alpha of the declassified signal according to the expression of the broadband signal 3b +α 4b And t, obtaining a declivity signal conjugated with the broadband signal according to the frequency modulation slope of the declivity signal, and performing declivity processing on the broadband signal by using the secondary term and the tertiary term of the declivity signal for counteracting the broadband signal.
The invention has the advantages that: the method is based on a plurality of continuous narrowband LFM signals with different waveforms, realizes the extraction of waveform parameters of a space target broadband signal, does not need to add extra hardware cost, does not change the existing radar imaging signal processing system, has small calculation load, obtains the frequency modulation slope of the declivity signal through the quadratic term coefficient and the cubic term coefficient of the broadband signal, obtains the declivity signal through the frequency modulation slope of the declivity signal, performs declivity processing on the broadband signal by using the declivity signal, realizes accurate declivity real-time declivity aiming at the broadband signal, has good declivity effect, and meanwhile, the waveform parameters of the broadband signal have important application value to the radar accurate imaging field and the follow-up target identification.
Drawings
FIG. 1 is a flow chart of a wideband polynomial phase signal deskewing method under narrowband signal guidance provided by an embodiment of the present invention;
fig. 2 is a schematic timing diagram of a narrowband LFM signal and a wideband signal in a wideband polynomial phase signal deskewing method under narrowband signal guidance according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions in the embodiments of the present invention will be clearly and completely described in the following in conjunction with the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
As shown in fig. 1, a wideband polynomial phase signal declivity method under narrowband signal guidance, the method comprising:
step S1: transmitting a plurality of continuous narrowband LFM signals; specific: the i-th narrowband LFM signal isWherein i=1, 2,..n, N is the total number of narrowband LFM signals, S i (t) is the ith narrowband LFM signal, j is the imaginary symbol, e () As an exponential function based on e, f 0 For the common working frequency point of all narrowband LFM signals, the narrowband LFM signals and the broadband signals share the same working frequency point, k i Frequency modulation slope of ith narrowband LFM signal, B i Is the amplitude of the i-th narrowband LFM signal.
Step S2: acquiring pulse echoes of a plurality of continuous narrowband LFM signals; specific: the pulse echo is expressed asWherein S is r i (t) is the ith pulse echo, alpha 1i 、α 2i 、α 3i And alpha 4i The zero order term coefficient, the first order term coefficient, the second order term coefficient and the third order term coefficient of the ith narrowband LFM signal respectively, A i Is the echo amplitude of the i-th narrowband LFM signal.
Step S3: performing delay correlation processing on each pulse echo to obtain a plurality of narrow-band echo delay signals; the specific process is as follows: delaying the ith pulse echo by tau components to obtain an ith narrow-band echo delay signal
Wherein (alpha) 3i +3α 4i t) is the instantaneous frequency of the ith narrowband echo delay signal.
Step S4: traversing all the narrow-band echo delay signals, taking values at fixed time, and performing fractional Fourier transform to obtain a quadratic term coefficient and a cubic term coefficient estimated by each narrow-band echo delay signal; the specific process is as follows: by the formulaObtaining a first intermediate quantity->And a second intermediate amount->Time values t1 and t2 can obtain the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantity
The third order term coefficient of the narrow-band echo delay signal estimation can be obtained according to the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantityAnd quadratic term coefficient->Respectively is
Step S5: according to the quadratic term coefficient and the cubic term coefficient estimated by each narrow-band echo delay signal, compensating the pulse echo of each narrow-band LFM signal, and carrying out Fourier transform on the compensated result to obtain each narrow-bandA first order term coefficient and a zero order term coefficient of echo delay signal estimation; the specific process is as follows: according to the quadratic term coefficient and the cubic term coefficient estimated by each narrow-band echo delay signal, compensating the pulse echo of each narrow-band LFM signal to obtain a compensated result
For the compensated result, the result is calculated by the formulaPerforming Fourier transform to obtain a first order coefficient of each narrowband echo delay signal estimation>And zero order term coefficient->
Step S6: obtaining the third term coefficient of the broadband signal according to the third term coefficient and the frequency modulation slope of the ith narrowband echo delay signal and the frequency modulation slope of the broadband signal; the specific process is as follows: by the formulaAcquiring a cubic coefficient of a broadband signal, wherein alpha 4b Is the third order coefficient of the wideband signal, beta br Is the frequency modulation slope of the wideband signal,/->Is the third order coefficient, beta, of the estimated ith narrowband LFM signal i Is the chirp rate of the i-th narrowband LFM signal.
Step S7: obtaining k-1 times term coefficients of the broadband signals according to the k-1 times term coefficients of the ith and jth narrowband echo delay signals, the frequency modulation slope and the frequency modulation slope of the broadband signals, wherein the value of k is 1,2, 3 or 4; the specific process is as follows: by the formulaObtaining the linear relation between the frequency modulation slope of the narrowband LFM signal and the k-1 th order coefficient, wherein beta 1i Frequency modulation slope beta of ith narrowband LFM signal 2j For the frequency modulation slope of the jth narrowband LFM signal, pa1 and pa2 are the first and second parameters, respectively, and +.> For the k-1 th order coefficient of the estimated i-th narrowband LFM signal,/I>K-1 th order term coefficients for the estimated jth narrowband LFM signal;
by the formula alpha kb =β br Pa1+pa2 obtaining the k-1 th order coefficients of the wideband signal, where α kb Is the k-1 th order term coefficient of the wideband signal.
Step S8: the frequency modulation slope of the declivity signal is obtained according to the quadratic term coefficient and the cubic term coefficient of the broadband signal, and the declivity signal is utilized to declivity the broadband signal, and the specific process is as follows:
obtaining expression Sigmatch of broadband signal according to zero order coefficient, primary order coefficient, secondary order coefficient and tertiary order coefficient of broadband signal b (t)=exp[j2π((α 1 +α 2b t+α 3b t 2 +α 4b t 3 ))]
Obtaining the frequency modulation slope alpha of the declassified signal according to the expression of the broadband signal 3b +α 4b And t, obtaining a declivity signal conjugated with the broadband signal according to the frequency modulation slope of the declivity signal, and performing declivity processing on the broadband signal by using the secondary term and the tertiary term of the declivity signal for counteracting the broadband signal.
In order to intuitively understand the scheme of the present invention, as shown in fig. 2, a schematic timing diagram of a wideband pulse and a narrowband pulse is shown, where the wideband pulse is a wideband signal of the present invention, and the narrowband pulse is a narrowband LFM signal of the present invention, and the timing sequence of the wideband pulse and the narrowband pulse is formed by a plurality of wideband-narrowband pulse matching pulse sequences, where each wideband-narrowband pulse matching pulse sequence includes N narrowband LFM signals with different frequency modulation slopes and one wideband signal.
According to the technical scheme, the wideband polynomial phase signal declining method under the guidance of the narrowband signals, provided by the invention, is based on a plurality of continuous narrowband LFM signals with different waveforms, so that the extraction of waveform parameters of a space target wideband signal is realized, no extra hardware cost is required to be added, an existing radar imaging signal processing system is not changed, the calculation load is small, the declining signal frequency modulation slope is obtained through the quadratic term coefficient and the cubic term coefficient of the wideband signal, the declining signal is obtained through the declining signal frequency modulation slope, the declining signal is utilized to declinize the wideband signal, the precise declining real-time declining aiming at the wideband signal is realized, the declining effect is good, and meanwhile, the waveform parameters of the wideband signal have important application value to the radar precise imaging field and the subsequent target identification.
Example 2
Corresponding to embodiment 1 of the present invention, embodiment 2 of the present invention further provides a wideband polynomial phase signal deskewing apparatus under narrowband signal guidance, the apparatus including:
the transmitting module is used for transmitting a plurality of continuous narrowband LFM signals;
the receiving module is used for acquiring pulse echoes of a plurality of continuous narrowband LFM signals;
the delay correlation processing module is used for carrying out delay correlation processing on each pulse echo to obtain a plurality of narrow-band echo delay signals;
the first Fourier transform module is used for traversing all the narrow-band echo delay signals, taking the value at fixed time, and carrying out fractional Fourier transform to obtain a quadratic term coefficient and a cubic term coefficient estimated by each narrow-band echo delay signal;
the second Fourier transform module is used for compensating the pulse echo of each narrowband LFM signal according to the quadratic term coefficient and the cubic term coefficient estimated by each narrowband echo delay signal, and carrying out Fourier transform on the compensated result to obtain the quadratic term coefficient and the zero-order term coefficient estimated by each narrowband echo delay signal;
the first broadband signal coefficient acquisition module is used for acquiring the cubic term coefficient of the broadband signal according to the cubic term coefficient and the frequency modulation slope of the ith narrowband echo delay signal and the frequency modulation slope of the broadband signal;
the second broadband signal coefficient acquisition module is used for acquiring k-1 times of coefficients of the broadband signals according to the k-1 times of coefficients of the ith and jth narrowband echo delay signals, the frequency modulation slope and the frequency modulation slope of the broadband signals, wherein the value of k is 1,2 or 3;
and the declivity module is used for obtaining the frequency modulation slope of the declivity signal according to the quadratic term coefficient and the cubic term coefficient of the broadband signal, and declivity processing is carried out on the broadband signal by utilizing the declivity signal.
Specifically, the ith narrowband LFM signal in the transmitting module isWherein i=1, 2,..n, N is the total number of narrowband LFM signals, S i (t) is the ith narrowband LFM signal, j is the imaginary symbol, e () As an exponential function based on e, f 0 For the common working frequency point of all narrowband LFM signals, the narrowband LFM signals and the broadband signals share the same working frequency point, k i Frequency modulation slope of ith narrowband LFM signal, B i Is the amplitude of the i-th narrowband LFM signal.
Specifically, the expression of the pulse echo in the receiving module is thatWherein S is r i (t) is the ith pulse echo, alpha 1i 、α 2i 、α 3i And alpha 4i The zero order term coefficient, the first order term coefficient, the second order term coefficient and the third order term coefficient of the ith narrowband LFM signal respectively, A i Is the echo amplitude of the i-th narrowband LFM signal.
Specifically, the delay related processing module is further configured to: delaying the ith pulse echo by tau components to obtain an ith narrow-band echo delay signal
Wherein (alpha) 3i +3α 4i t) is the instantaneous frequency of the ith narrowband echo delay signal.
Specifically, the first fourier transform module is further configured to: by the formulaObtaining a first intermediate quantity->And a second intermediate amount->Time values t1 and t2 can obtain the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantity
The third order term coefficient of the narrow-band echo delay signal estimation can be obtained according to the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantityAnd quadratic term coefficient->Respectively is
Specifically, the second fourier transform module is further configured to: compensating pulse echo of each narrowband LFM signal according to the quadratic term coefficient and the cubic term coefficient estimated by each narrowband echo delay signal to obtain compensationResults after
For the compensated result, the result is calculated by the formulaPerforming Fourier transform to obtain a first order coefficient of each narrowband echo delay signal estimation>And zero order term coefficient->
Specifically, the first wideband signal coefficient acquisition module is further configured to: by the formulaAcquiring a cubic coefficient of a broadband signal, wherein alpha 4b Is the third order coefficient of the wideband signal, beta br Is the frequency modulation slope of the wideband signal,/->Is the third order coefficient, beta, of the estimated ith narrowband LFM signal i Is the chirp rate of the i-th narrowband LFM signal.
Specifically, the second wideband signal coefficient acquisition module is further configured to: by the formulaObtaining the linear relation between the frequency modulation slope of the narrowband LFM signal and the k-1 th order coefficient, wherein beta 1i Frequency modulation slope beta of ith narrowband LFM signal 2j For the frequency modulation slope of the jth narrowband LFM signal, pa1 and pa2 are the first and second parameters, respectively For the k-1 th order term coefficients of the estimated i-th narrowband LFM signal,k-1 th order term coefficients for the estimated jth narrowband LFM signal;
by the formula alpha kb =β br Pa1+pa2 obtaining the k-1 th order coefficients of the wideband signal, where α kb Is the k-1 th order term coefficient of the wideband signal.
Specifically, the declivity module is further configured to:
obtaining expression Sigmatch of broadband signal according to zero order coefficient, primary order coefficient, secondary order coefficient and tertiary order coefficient of broadband signal b (t)=exp[j2π((α 1 +α 2b t+α 3b t 2 +α 4b t 3 ))]
Obtaining the frequency modulation slope alpha of the declassified signal according to the expression of the broadband signal 3b +α 4b And t, obtaining a declivity signal conjugated with the broadband signal according to the frequency modulation slope of the declivity signal, and performing declivity processing on the broadband signal by using the secondary term and the tertiary term of the declivity signal for counteracting the broadband signal.
The above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for wideband polynomial phase signal deskewing under narrowband signal steering, the method comprising:
step one: transmitting a plurality of continuous narrowband LFM signals;
step two: acquiring pulse echoes of a plurality of continuous narrowband LFM signals;
step three: performing delay correlation processing on each pulse echo to obtain a plurality of narrow-band echo delay signals;
step four: traversing all the narrow-band echo delay signals, taking values at fixed time, and performing fractional Fourier transform to obtain a quadratic term coefficient and a cubic term coefficient estimated by each narrow-band echo delay signal;
step five: according to the quadratic term coefficient and the cubic term coefficient of each narrowband echo delay signal estimation, compensating the pulse echo of each narrowband LFM signal, and carrying out Fourier transform on the compensated result to obtain the primary term coefficient and the zero-order term coefficient of each narrowband echo delay signal estimation;
step six: obtaining the third term coefficient of the broadband signal according to the third term coefficient and the frequency modulation slope of the ith narrowband echo delay signal and the frequency modulation slope of the broadband signal;
step seven: obtaining k-1 times term coefficients of the broadband signals according to the k-1 times term coefficients of the ith and jth narrowband echo delay signals, the frequency modulation slope and the frequency modulation slope of the broadband signals, wherein the value of k is 1,2, 3 or 4;
step eight: and obtaining the frequency modulation slope of the declivity signal according to the quadratic term coefficient and the cubic term coefficient of the broadband signal, and performing declivity processing on the broadband signal by using the declivity signal.
2. The wideband polynomial phase signal deskewing method according to claim 1, wherein the i-th narrowband LFM signal in the step one isWherein i=1, 2,..n, N is the total number of narrowband LFM signals, S i (t) is the ith narrowband LFM signal, j is the imaginary symbol, e () As an exponential function based on e, f 0 For the common working frequency point of all narrowband LFM signals, the narrowband LFM signals and the broadband signals share the same working frequencyPoints k i Frequency modulation slope of ith narrowband LFM signal, B i Is the amplitude of the i-th narrowband LFM signal.
3. The wideband polynomial phase signal deskewing method according to claim 2, wherein the pulse echo in step two is expressed asWherein S is r i (t) is the ith pulse echo, alpha 1i 、α 2i 、α 3i And alpha 4i The zero order term coefficient, the first order term coefficient, the second order term coefficient and the third order term coefficient of the ith narrowband LFM signal respectively, A i Is the echo amplitude of the i-th narrowband LFM signal.
4. A wideband polynomial phase signal deskewing method under narrowband signal steering as claimed in claim 3 wherein step three includes: delaying the ith pulse echo by tau components to obtain an ith narrow-band echo delay signal
Wherein (alpha) 3i +3α 4i t) is the instantaneous frequency of the ith narrowband echo delay signal.
5. The wideband polynomial phase signal deskewing method under narrowband signal steering of claim 4, wherein step four comprises: by the formulaObtaining a first intermediate quantity->And a second intermediate amount->Time values t1 and t2 can obtain the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantity
The third order term coefficient of the narrow-band echo delay signal estimation can be obtained according to the relation between the instantaneous frequency of the narrow-band echo delay signal and the first intermediate quantity and the second intermediate quantityAnd quadratic term coefficient->Respectively is
6. The wideband polynomial phase signal deskewing method under narrowband signal steering of claim 5, wherein step five comprises: according to the quadratic term coefficient and the cubic term coefficient estimated by each narrow-band echo delay signal, compensating the pulse echo of each narrow-band LFM signal to obtain a compensated result
7. The wideband polynomial phase signal deskewing method under narrowband signal steering of claim 1, wherein step six includes: by the formulaAcquiring a cubic coefficient of a broadband signal, wherein alpha 4b Is the third order coefficient of the wideband signal, beta br Is the frequency modulation slope of the wideband signal,/->Is the third order coefficient, beta, of the estimated ith narrowband LFM signal i Is the chirp rate of the i-th narrowband LFM signal.
8. The wideband polynomial phase signal deskewing method under narrowband signal steering of claim 7, wherein step seven comprises: by the formulaObtaining the linear relation between the frequency modulation slope of the narrowband LFM signal and the k-1 th order coefficient, wherein beta 1i Frequency modulation slope beta of ith narrowband LFM signal 2j For the frequency modulation slope of the jth narrowband LFM signal, pa1 and pa2 are the first and second parameters, respectively, and +.> For the k-1 th order coefficient of the estimated i-th narrowband LFM signal,/I>K-1 th order term coefficients for the estimated jth narrowband LFM signal;
by the formula alpha kb =β br Pa1+pa2 obtaining the k-1 th order coefficients of the wideband signal, where α kb Is the k-1 th order term coefficient of the wideband signal.
9. The wideband polynomial phase signal deskewing method under narrowband signal steering of claim 8, wherein step eight comprises:
obtaining expression Sigmatch of broadband signal according to zero order coefficient, primary order coefficient, secondary order coefficient and tertiary order coefficient of broadband signal b (t)=exp[j2π((α 1 +α 2b t+α 3b t 2 +α 4b t 3 ))]
Obtaining the frequency modulation slope alpha of the declassified signal according to the expression of the broadband signal 3b +α 4b And t, obtaining a declivity signal conjugated with the broadband signal according to the frequency modulation slope of the declivity signal, and performing declivity processing on the broadband signal by using the secondary term and the tertiary term of the declivity signal for counteracting the broadband signal.
10. A wideband polynomial phase signal deskewing apparatus under narrowband signal steering, the apparatus comprising:
the transmitting module is used for transmitting a plurality of continuous narrowband LFM signals;
the receiving module is used for acquiring pulse echoes of a plurality of continuous narrowband LFM signals;
the delay correlation processing module is used for carrying out delay correlation processing on each pulse echo to obtain a plurality of narrow-band echo delay signals;
the first Fourier transform module is used for traversing all the narrow-band echo delay signals, taking the value at fixed time, and carrying out fractional Fourier transform to obtain a quadratic term coefficient and a cubic term coefficient estimated by each narrow-band echo delay signal;
the second Fourier transform module is used for compensating the pulse echo of each narrowband LFM signal according to the quadratic term coefficient and the cubic term coefficient estimated by each narrowband echo delay signal, and carrying out Fourier transform on the compensated result to obtain the quadratic term coefficient and the zero-order term coefficient estimated by each narrowband echo delay signal;
the first broadband signal coefficient acquisition module is used for acquiring the cubic term coefficient of the broadband signal according to the cubic term coefficient and the frequency modulation slope of the ith narrowband echo delay signal and the frequency modulation slope of the broadband signal;
the second broadband signal coefficient acquisition module is used for acquiring k-1 times of coefficients of the broadband signals according to the k-1 times of coefficients of the ith and jth narrowband echo delay signals, the frequency modulation slope and the frequency modulation slope of the broadband signals, wherein the value of k is 1,2 or 3;
and the declivity module is used for obtaining the frequency modulation slope of the declivity signal according to the quadratic term coefficient and the cubic term coefficient of the broadband signal, and declivity processing is carried out on the broadband signal by utilizing the declivity signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010904054.6A CN111999725B (en) | 2020-09-01 | 2020-09-01 | Wideband polynomial phase signal declivity method and device under guidance of narrowband signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010904054.6A CN111999725B (en) | 2020-09-01 | 2020-09-01 | Wideband polynomial phase signal declivity method and device under guidance of narrowband signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111999725A CN111999725A (en) | 2020-11-27 |
CN111999725B true CN111999725B (en) | 2023-07-11 |
Family
ID=73466211
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010904054.6A Active CN111999725B (en) | 2020-09-01 | 2020-09-01 | Wideband polynomial phase signal declivity method and device under guidance of narrowband signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111999725B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391298A (en) * | 2014-11-28 | 2015-03-04 | 中国电子科技集团公司第三十八研究所 | Digital delaying compensation device and method and achieving device and method |
CN105301590A (en) * | 2015-11-03 | 2016-02-03 | 南京信息职业技术学院 | Maneuvering target frequency modulation stepping inverse synthetic aperture imaging method |
CN108919221A (en) * | 2018-07-17 | 2018-11-30 | 武汉大学 | A kind of phase-coherent accumulation detection method for variable accelerated motion target |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5418892A (en) * | 1993-01-12 | 1995-05-23 | Board Of Trustees Of The Leland Stanford Junior University | Subspace-based line detection |
US7054504B2 (en) * | 1999-02-25 | 2006-05-30 | Ludwig Lester F | Relative optical path phase reconstruction in the correction of misfocused images using fractional powers of the fourier transform |
CN101833082B (en) * | 2010-04-20 | 2012-06-27 | 中国科学院空间科学与应用研究中心 | Wideband frequency-modulation stepping signal processing method based on full deskew |
CN101963662B (en) * | 2010-09-20 | 2013-09-25 | 北京理工大学 | Self-focusing preprocessing method based on short-time fractional order Fourier domain filter |
EP2976624A1 (en) * | 2013-03-21 | 2016-01-27 | The Technical University of Denmark | Refractive index based measurements |
CN106597440B (en) * | 2016-12-12 | 2018-11-27 | 南京信息职业技术学院 | Low signal-to-noise ratio imaging method for frequency modulation stepping radar |
KR101879516B1 (en) * | 2016-12-22 | 2018-07-17 | 한국항공우주연구원 | Apparatus and method for digital dechirp processing |
CN106872974B (en) * | 2017-01-23 | 2019-07-09 | 西安电子科技大学 | High-precision motion target imaging method based on hypersonic platform Two-channels radar |
CN108107430B (en) * | 2017-11-09 | 2023-03-03 | 北京理工大学 | Ship target ISAR imaging method based on fractional Fourier transform |
CN108008369B (en) * | 2017-11-30 | 2019-09-24 | 中国科学院国家空间科学中心 | A kind of NLFM signal lack sampling processing method |
CN108152820B (en) * | 2017-12-20 | 2021-10-08 | 西安电子科技大学 | A Bistatic Radar Imaging Method Based on Tomography |
CN108318879B (en) * | 2017-12-29 | 2021-08-31 | 哈尔滨工业大学 | ISAR Image Lateral Calibration Method Based on IAA Spectral Estimation Technology |
CN110146889B (en) * | 2019-06-17 | 2020-07-14 | 中国人民解放军国防科技大学 | Large-rotation-angle ISAR imaging method based on optimal echo sub-region selection |
CN110988874B (en) * | 2019-11-08 | 2022-03-22 | 西北大学 | An ISAR Imaging Method for Complex Moving Targets |
CN111505640B (en) * | 2020-06-24 | 2021-04-13 | 深圳大学 | Method and system for deskew one-bit acquisition of broadband radar echo |
-
2020
- 2020-09-01 CN CN202010904054.6A patent/CN111999725B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391298A (en) * | 2014-11-28 | 2015-03-04 | 中国电子科技集团公司第三十八研究所 | Digital delaying compensation device and method and achieving device and method |
CN105301590A (en) * | 2015-11-03 | 2016-02-03 | 南京信息职业技术学院 | Maneuvering target frequency modulation stepping inverse synthetic aperture imaging method |
CN108919221A (en) * | 2018-07-17 | 2018-11-30 | 武汉大学 | A kind of phase-coherent accumulation detection method for variable accelerated motion target |
Non-Patent Citations (1)
Title |
---|
基于FRFT的一维距离像高速运动补偿;袁斌;唐鹏飞;徐世友;陈曾平;;雷达科学与技术(第03期);第84-89,97页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111999725A (en) | 2020-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104898119B (en) | A kind of moving target parameter estimation method based on correlation function | |
CN106842166A (en) | A kind of solution velocity ambiguity method suitable for LFMCW radar system | |
US5237329A (en) | Method of correcting range migration in image generation in synthetic aperture radar | |
Shah et al. | Step-frequency radar with compressive sampling (SFR-CS) | |
WO2011083300A1 (en) | Pulse radar range profile motion compensation | |
CN104237886A (en) | High-precision synthetic aperture radar imaging method | |
CN107356908A (en) | A kind of frequency agility signal correlative accumulation method | |
CN111999725B (en) | Wideband polynomial phase signal declivity method and device under guidance of narrowband signal | |
Dai et al. | High accuracy velocity measurement based on keystone transform using entropy minimization | |
CN112630740A (en) | Dense multi-target distinguishing method and device | |
CN108363047B (en) | Four-phase code Doppler compensation method based on sliding window interpolation method | |
CN113296064B (en) | SAR radar receiving channel time delay calibration method and system based on Frank code | |
CN101738615B (en) | Method for estimating rotating speed of target of inverse synthetic aperture radar | |
US20230324509A1 (en) | Automotive radar with range migration mitigation capability | |
CN106093907B (en) | A kind of broadband LFMCW radar channel-equalization methods based on special aobvious point | |
Liu et al. | Joint estimation of range difference and range rate difference for unknown frequency hopping signals in passive localization | |
CN113567950B (en) | Millimeter wave radar distance and velocity spectrum estimation method and system | |
CN119179077B (en) | Azimuth envelope aligned high-resolution moving ship SAR fine focusing method and system | |
Zeng et al. | Doppler ambiguity correction based on keystone transform for wideband radar | |
CN114488023B (en) | Dual matched filtering method and device and electronic equipment | |
CN112578381B (en) | A target three-dimensional rotation vector estimation method based on interferometric inverse synthetic aperture radar | |
Peng et al. | Study on transmitting mode and imaging algorithm of MIMO-SAR | |
EP4249939A1 (en) | Radar receiver | |
Zhang et al. | A novel method for estimating the baseband Doppler centroid of conventional synthetic aperture radar | |
Tian et al. | Efficient motion parameter estimation for maneuvering target based on segment processing |
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 |