CN108983193A - Quickly non-search ground moving object method for parameter estimation - Google Patents
Quickly non-search ground moving object method for parameter estimation Download PDFInfo
- Publication number
- CN108983193A CN108983193A CN201810847452.1A CN201810847452A CN108983193A CN 108983193 A CN108983193 A CN 108983193A CN 201810847452 A CN201810847452 A CN 201810847452A CN 108983193 A CN108983193 A CN 108983193A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- target
- parameter
- formula
- 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.)
- Pending
Links
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention proposes a kind of quick non-search ground moving object method for parameter estimation, mainly solve the problems, such as that existing method compensation effect is poor, estimated accuracy is low, realize the refocusing of the accurate estimation and target to the parameters of target motion.Its implementation is: carrying out matched filtering to signal is received using azimuth-deramp matched filtering function, eliminates spectrum division；SOKT transformation and phase difference processing uncompensation distance bending and range walk are passed sequentially through again；Doppler ambiguity number and Secondary movement parameter are estimated by the peak coordinate that IFFT and FFT transform obtain；Matched filtering is carried out to signal is received using match adaptation function, then successively carries out KT transformation, IFFT and FFT transform, and first order motion parameter is estimated by peak coordinate.This method compensation effect is good, and computational efficiency is high, and Parameter Estimation Precision is high, can be used for refocusing of the ground moving object in SAR image.
Description
Technical field
The invention belongs to Radar Technology field, in particular to a kind of quick non-search target component estimation method can be used
In refocusing of the ground moving object in SAR image.
Background technique
Traditional synthetic aperture radar SAR can obtain the image of static scene under round-the-clock, all weather conditions, currently
It is widely used in civilian and military field.But when target occurs mobile, it may appear that new Doppler component causes distance to be walked
Dynamic and range curvature causes mobile target to generate dislocation in SAR image or even can not show.Therefore, in order to accurately estimate to transport
Moving-target parameter, needing to adjust the distance to walk about carries out good compensation with range curvature to realize refocusing.
Perry et al. proposes KT transform method, and this method shows good compensation performance under low signal-to-noise ratio environment, can
Compensate the range walk of multiple targets simultaneously in the case where mobile target first order parameter is unknown.But in practical applications, due to
Pulse recurrence frequency is limited, and doppler ambiguity often occurs for the target fast moved, and KT method can not solve Doppler's mould
Paste problem, performance can be remarkably decreased.In order to solve the problems, such as doppler ambiguity, a kind of line by searching moving target is proposed
Property track and speed realize the uncoupled Radon-FFT method in range-azimuth.Hereafter, zheng et al. proposes Deramp-
Keystone transformation, this method can be in the case where not searching for doppler ambiguity number, to target in the case of doppler ambiguity
Range curvature compensates.However practice discovery this method is adjusted the distance, the compensation walked about with range curvature is still not ideal enough, and counts
It is lower to calculate efficiency, it is necessary to carry out further compensation deals.
Summary of the invention
It is an object of the invention in view of the above shortcomings of the prior art, propose a kind of quick non-search ground motion mesh
Method for parameter estimation is marked, to improve computational efficiency, promotes compensation effect.
The technical scheme is that using azimuth-deramp matched filtering-second order Keystone transformation-phase difference
Divide PD technology, uncompensation distance bending and range walk, correct Doppler shift and dopplerbroadening, realizes one, Secondary movement ginseng
The refocusing of several accurate estimation and moving target, implementation step include the following:
(1) to based on SAR platform geometrical model reception signal doing mathematics deformation and abbreviation, obtain about frequency variation
F and slow time tmReception signal S1(f,tm)；
(2) using azimuth-deramp adaptation function to reception signal S1(f,tm) matched filtering is carried out, eliminate Doppler
Spectrum division, the signal S ' after being matched1(f,tm)；
(3) to the signal S ' after matching1(f,tm) second order Keystone transformation is carried out, it realizes range curvature compensation, is closed
In frequency variation f and new slow time ηmSignal S2(f,ηm)；
(4) it adjusts the distance and is bent thermal compensation signal S2(f,ηm) phase difference PD processing is carried out, obtain single order phase signal S3(f,
ηm)；
(5) to single order phase signal S3(f,ηm) execute the domain f IFFT transformation and ηmDomain FFT transform realizes that range walk is mended
It repays, obtains about fast time t and frequency faPeak detection signal S4(t,fa)；
(6) to peak detection signal S4(t,fa) detection is carried out, estimate doppler ambiguity numberWith the two of target
Rank kinematic parameter
(7) using match matched filtering function to reception signal S1(f,tm) carry out matched filtering and single order Keystone change
It changes, to adjust the distance bending and Doppler shift compensates, obtains about frequency variation f and new slow time variable τmWithout frequency
The number of believing one side only S '4(f,τm)；
(8) to no frequency offset signal S '4(f,τm) execute the domain f IFFT transformation and τmDomain FFT transform is obtained comprising first order motion
The peak detection signal S of parameter information5(t,fa)；
(9) to the peak detection signal S comprising first order motion parameter information5(t,fa) detection is carried out, estimate the one of target
Rank kinematic parameter
The invention has the following advantages over the prior art:
The present invention is due to using matched filtering technique, Keystone transformation and phase difference PD technology, thus calculating process
Addition, complex multiplication, IFFT and FFT are only related to, so that computational efficiency is obviously improved；
The present invention is avoided due to using non-search algorithm to doppler ambiguity number and single order, Secondary movement parameter
Search so that the complexity of radar system reduces；
Simulation result shows that the present invention can effectively inhibit the cross term interference under multiple mobile object scene.
Detailed description of the invention
Fig. 1 is usage scenario figure of the invention；
Fig. 2 is implementation flow chart of the invention；
Fig. 3 is with the method for the present invention to the inhibitory effect figure of cross term under multiple mobile object scene；
Fig. 4 is for the present invention to the imaging effect figure of ground moving object under simulated environment；
Fig. 5 is for the present invention to the imaging effect figure of ground moving object under true environment.
Specific embodiment
Referring to Fig.1, usage scenario of the invention includes carried SAR platform and i-th of ground moving object, and wherein SAR is flat
Platform is aircraft, and for detecting the automobile on ground, the pitch angle of aircraft is zero, speed v, and target is in tmIt is moved in time from the position a
The position b is moved, direct motion speed and inter-bank speed are respectively vaiAnd vci。
Referring to Fig. 2, steps are as follows for realization of the invention:
Step 1, to signal S (t, tm) doing mathematics deform and abbreviation, the receives echo-signal S after obtaining abbreviation1(f,tm)。
(1a) radar receiver is in tmEcho-signal S (t, t of the reception to K targetm):
In formula, A0,iFor amplitude, tmFor slow time, Ri(tm) it is the distance between aircraft and i-th of moving target, rect
It (x) is window function that section [- 1/2,1/2] upper value is 1, μ is frequency modulation rate, TpFor the pulse duration for emitting signal, t is fast
Time, TaFor the pulse spacing, c is the light velocity, Ri(tm) be the distance between aircraft and i-th of moving target,Wherein H is podium level, YiIt is aircraft flight path on ground
The distance for projecting to i-th of moving target on face；
(1b) is to Ri(tm) make Taylor series expansion and ignore higher order term:
In formula, R0iFor zero moment aircraft to the distance of i-th of moving target,a1I and a2iRespectively mesh
Target first order motion parameter and Secondary movement parameter, may be expressed as:
(1c) is to a1iIt is decomposed:
In formula, vb,iIndicate base band speed, Mamb,iIndicate doppler ambiguity number, PRF is pulse recurrence frequency；
(1d) obtains the receives echo-signal S after abbreviation in conjunction with (1a), (1b) and (1c)1(f,tm), it is expressed as follows:
In formula, A1,iFor compressed amplitude, fcFor center frequency.
Signal S1(f,tm) in f-tmCoupling terms lead to range walk,Coupling terms lead to range curvature, tmXiang He
Item causes Doppler shift and broadening.
Step 2, to the receives echo-signal S after abbreviation1(f,tm) azimuth-deramp matched filtering is carried out, it is matched
Signal S ' afterwards1(f,tm)。
(2a) selects following azimuth-deramp matched filtering function:
In formula, f is frequency variation, fcFor center frequency, RbIndicate the nearest slant range from aircraft to center line, tm
For the slow time；
(2b) is by azimuth-deramp matched filtering function to receives echo-signal S1(f,tm) matched filtering is carried out,
Signal S ' after being matched1(f,tm):
In formula, A1,iFor the amplitude after matching, rect (x) is the window function that section [- 1/2,1/2] upper value is 1, and μ is frequency modulation
Rate, TpFor the pulse duration for emitting signal, TaFor the pulse spacing, c is the light velocity, R0iFor zero moment SAR platform and i-th of mesh
Target distance, vb,iFor base band speed,Mamb,iFor doppler ambiguity number, PRF is pulse recurrence frequency.
It is handled by above-mentioned azimuth-deramp matched filtering, compensates for the partial distance bending of moving target, reduce
Azimuth of target bandwidth, eliminates orientation angular spectrum separating phenomenon.
Step 3, to the signal S ' after matching1(f,tm) second order Keystone transformation is carried out, obtain no range curvature signal S2
(f,ηm)。
It is as follows that (3a) second order Keystone converts expression formula:
In formula, f is frequency variation, fcFor center frequency, tmFor slow time variable, ηmIndicate new slow time variable；
Second order Ketstone transformation expression formula is brought into the signal S ' after matching by (3b)1(f,tm) in, it obtains about f and new
Slow time variable ηmRange curvature thermal compensation signal S2(f,ηm):
In formula, A1,iFor compressed amplitude, rect (x) is the window function that section [- 1/2,1/2] upper value is 1, and μ is frequency modulation
Rate, TpFor the pulse duration for emitting signal, TaFor the pulse spacing, c is the light velocity, R0iFor zero moment SAR platform and i-th of mesh
Target distance, λ are wavelength, vb,iFor base band speed,PRF is pulse recurrence frequency, Mamb,iIt how general is
Strangle fuzzy number.
It is converted by above-mentioned second order Keystone, compensates for the range curvature of target.
Step 4, it adjusts the distance and is bent thermal compensation signal S2(f,ηm) phase difference PD processing is carried out, obtain single order phase signal S3
(f,ηm)。
The expression formula of (4a) phase difference PD processing is as follows:
In formula, x (ηm) it is linear FM signal,τ0For permanent delay
Time, * are conjugate operation, τ0For difference interval, fcFor center frequency.
(4b), which adjusts the distance, is bent thermal compensation signal S2(f,ηm) phase difference PD processing is carried out, obtain single order phase signal S3(f,
ηm):
In formula, A1,iFor compressed amplitude, rect (x) is the window function that section [- 1/2,1/2] upper value is 1, and f is frequency
Variable quantity, μ are frequency modulation rate, TpFor the pulse duration for emitting signal, ηmFor new slow time, TaFor the pulse spacing, λ is wave
It is long, vb,iFor base band speed,C is the light velocity, Mamb,iFor doppler ambiguity number, PRF is pulse repetition
Frequency, S3,cross(f,ηm) it is cross term,Item is much smaller than π/4, can directly ignore.
It is handled by above-mentioned phase difference PD, reduces the order of signal phase, compensate for the range walk of target.
Step 5, according to single order phase signal S3(f,ηm), obtain the peak detection signal comprising Secondary movement parameter information
S4(t,fa)。
To single order phase signal S3(f,ηm) carry out the domain f IFFT transformation and ηmDomain FFT transform obtains the peak value containing cross term
Rectified signal S4(t,fa):
In formula, t is fast time, faFor frequency, ηmFor new slow time variable, A3,iFor amplitude, λ is wavelength, vb,iFor base band
Speed, τ0For difference interval,B is signal bandwidth, and c is the light velocity, Mamb,iFor doppler ambiguity number, PRF
For pulse recurrence frequency, TaFor pulse spacing, S4,cross(t,fa) it is the domain f IFFT transformation and ηmIntersection after the operation of domain FFT transform
?.
Step 6, to the peak detection signal S comprising Secondary movement parameter information4(t,fa) detection is carried out, estimate target
Secondary movement parameter.
(6a) is to the peak detection signal S comprising Secondary movement parameter information4(t,fa) detection is carried out, learn that wave crest is located atPlace, estimates the parameters of target motion by the wave crest coordinateAnd doppler ambiguity number
In formula,Indicating rounding-off operation, λ is wavelength,For the estimated value of frequency, τ0For difference interval, c is the light velocity,For
The estimated value of fast time t, PRF is pulse recurrence frequency；
(6b) is according to Target moving parameter estimation valueEstimate target Secondary movement parameter
Wherein, RbIndicate that the nearest slant range from radar platform to band center line, v are the flying speed of aircraft.
Step 7, to the receives echo-signal S after abbreviation in step 11(f,tm) carry out match matched filtering and single order
Keystone transformation, obtains no frequency offset signal S '4(f,τm)。
(7a) constructs match matched filtering function by Target moving parameter estimation value:
In formula, f is frequency variation, fcFor center frequency, tmFor the slow time, c is the light velocity,For doppler ambiguity
Number, PRF is pulse recurrence frequency,For the estimated value of target Secondary movement parameter；
It is as follows that (7b) single order Keystone converts expression formula:
(f+fc)tm=fcτm
In formula, τmFor new slow time variable；
(7c) is to the receives echo-signal S after abbreviation1(f,tm) carry out match matched filtering, the signal after being matched
S″1(f,tm):
In formula, A1,iFor compressed amplitude, rect (x) is the window function that section [- 1/2,1/2] upper value is 1, and μ is frequency modulation
Rate, TpFor the pulse duration for emitting signal, TaFor pulse spacing, R0iFor zero moment SAR platform and i-th target away from
From vb,iFor base band speed；
(7d) brings single order Keystone transformation expression formula into signal S " after match matched filtering1(f,tm), obtain nothing
Frequency offset signal S '4(f,τm):
By above-mentioned match matched filtering and single order Keystone conversion process, range walk, the distance of target are compensated for
Bending and Doppler shift and broadening.
Step 8, according to no frequency offset signal S '4(f,τm), obtain the peak detection signal S comprising first order motion parameter information5
(t,fa)。
To no frequency offset signal S '4(f,τm) carry out the domain f IFFT transformation and τmDomain FFT transform is obtained comprising first order motion parameter
The peak detection signal S of information5(t,fa):
In formula, A4,iTo integrate amplitude, t is fast time, faFor frequency, A4,iFor amplitude, B is signal bandwidth, R0iWhen being zero
SAR platform is carved at a distance from i-th of target, TaFor pulse spacing, vb,iFor base band speed, λ is wavelength.
Step 9, estimate the first order motion parameter of target.
To the peak detection signal S comprising first order motion parameter information5(t,fa) detection is carried out, learn that wave crest is located at (fa
=-2vb,i/ λ, t=2R0i/ c) at, the base band speed of moving target is estimated by the wave crest coordinateWith first order motion parameterIt is as follows
In formula,For the estimated value of base band speed, λ is wavelength, and PRF is pulse recurrence frequency,For Doppler's mould
Paste the estimated value of number.
Effect of the invention can be further illustrated by following simulation result:
Emulation experiment 1, inhibitory effect of the verifying the method for the present invention to cross term.
1.1) parameter setting
The corresponding simulation parameter of SAR platform is as follows:
Pulse recurrence frequency is 600Hz, and carrier frequency is 10GHz, and signal bandwidth is 200MHz, and sample frequency is
The flying speed of 400MHz, SAR platform are 120 meter per seconds, and the height of SAR platform is 6000 meters, ground vertical flight tracking range
It is 6000 meters, azimuth anteena size is 0.6 meter.
Assuming that indicated with T1, T2 and T3 respectively in scene there are three moving target, the kinematic parameter of these three targets and wait
It is as shown in table 1 to imitate kinematic parameter.
1 parameters of target motion table of table
1.2) emulation content
It is arranged using above-mentioned parameter, handles receives echo-signal, the change of object observing imaging effect using the method for the present invention
Change situation, as a result as shown in Figure 3, in which:
Fig. 3 (a) be target it is fast when m- slow time diagram, the motion profile of three targets is what curvature was not zero in figure
Oblique line illustrates that three targets all have range walk and range curvature；
Fig. 3 (b) is the doppler spectral of target, since the doppler spectral of T1 and T3 crosses over two PRF, so its orientation angular spectrum
It is divided, and the doppler spectral of T2, in a PRF, orientation angular spectrum does not divide；
Fig. 3 (c) is m- slow time diagram when target is fast after second order Keystone is converted, and the track of three targets is in figure
The oblique line that curvature is zero illustrates that the range curvature of all targets is all compensated；
Fig. 3 (d) is that the peak value figure obtained after phase difference PD is handled forms three due to the continuous accumulation of target energy
A apparent peak value, each peak value represents a target, and remaining peak point is far below at this 3 points, this illustrates that cross term obtains
Good inhibition is arrived.
It can be seen from Fig. 3 (d) under multiple mobile object scene, the method for the present invention can accurately estimate the movement of target
Parameter effectively inhibits cross term and displaying target position.
The kinematic parameter of target is estimated according to peak coordinate, the results are shown in Table 2.
2 Target moving parameter estimation value table of table
Comparison Tables 1 and 2 can be seen that the error of the target component estimated value and true value that obtain using the method for the present invention
Within 0.01, this illustrates that the method for the present invention has good parameter Estimation performance.
Emulation experiment 2, the verifying imaging effect of the method for the present invention to ground moving object under simulated environment.
2.1) parameter setting
The simulation parameter of SAR platform is as follows under simulated environment:
Pulse recurrence frequency is 600Hz, carrier frequency 10GHz, signal bandwidth 60MHz, sample frequency 120MHz,
Platform flying speed is 120 meter per seconds, and the pulse duration is 1 μ s, and nearest slant range is 850 kilometers, in the simulated environment
There are Gaussian noises；
Assuming that indicated with T1, T2 and T3 respectively in scene there are three moving target, the kinematic parameter of these three targets and away from
Tripping contracting signal-to-noise ratio is as shown in table 3.
3 parameters of target motion table of table
| Target 1 | Target 2 | Target 3 | |
| Inter-bank speed (m/s) | vc1=10 | vc2=-15 | vc3=-15 |
| Direct motion speed (m/s) | va1=10 | va2=30 | va3=-30 |
| Range compress signal-to-noise ratio (dB) | SNR1=1 | SNR2=1 | SNR3=2 |
2.2) emulation content
Using the parameter setting in 2.1, under simulated environment using the method for the present invention to the receives echo-signal of target into
Row processing, simulation result are as shown in Figure 4.Wherein:
Fig. 4 (a) be using before this method target it is fast when m- slow time diagram, the motion profile of three targets is equal in the figure
For the oblique line that curvature is not zero, illustrate that three targets all have range walk and range curvature, since there are Gausses for simulated environment
Noise has been filled with Gaussian noise point in entire figure；
Fig. 4 (b) is that the target peak figure for receiving and obtaining after signal is handled using the method for the present invention, occurs three in the figure
Apparent wave crest, the corresponding target of each wave crest；
Fig. 4 (c) be handled using the method for the present invention receive signal after moving target T1 be ultimately imaged result, it can be seen that
The peak value of target is it is obvious that be ultimately imaged result relatively good；
Fig. 4 (d) is that moving target T1 after signal is received using existing improved two-dimensional frequency matched filtering method processing
Imaging results, it can be seen that target does not have specific peak value, and position is spread, and imaging effect is poor；
Comparison diagram 4 (c) and Fig. 4 (d), it can be seen that under the simulated environment containing Gaussian noise, the method for the present invention compared to
Existing method has better focusing performance.
Emulation experiment 3 verifies under true environment the method for the present invention to the imaging effect of ground moving object.
3.1) parameter setting
The simulation parameter of SAR platform is as follows under true environment:
The carried SAR works in x wave band, pulse recurrence frequency 1000Hz, carrier frequency 8.85GHz, signal bandwidth
For 40MHz, sample frequency 60MHz, platform flying speed is 120 meter per seconds, and the pulse duration is 10 μ s, band center away from
From being 9000 meters.
3.2) emulation content
Using the parameter setting in 3.1, under true environment, using the method for the present invention to the receives echo-signal of target into
Row processing, the imaging effect of object observing, and the processing result of two kinds of existing methods is provided as a comparison, simulation result such as Fig. 5
It is shown, in which:
Fig. 5 (a) show to reception signal progress abbreviation and executes the imaging contexts of ground moving object after FFT transform,
It can be seen that moving target is by severe jamming under strong ground clutter background；
Fig. 5 (b) show the imaging contexts of target after inhibition ground clutter, marks one of movement mesh with red circle
Mark, it can be seen that the imaging effect of target is still very poor；
Fig. 5 (c) be labeled target it is fast when m- slow time diagram, since range walk has occurred, which is crossed over
Multiple distance unit；
Fig. 5 (d) is to handle the peak value figure obtained after receives echo-signal using the method for the present invention, and crest location is mesh
Mark；
Fig. 5 (e) is the sectional view tieed up in Doppler of peak value figure shown in Fig. 5 (d), according to peak coordinate estimate target etc.
Effect first order motion parameter is 11.255m/s, and Secondary movement parameter is 0.6563m/s2, inter-bank speed is -7.51m/s；
Fig. 5 (f) is the focusing effect figure that target after receives echo-signal is handled using existing deramp-keystone method,
The position of target is spread in the figure, and focus form is ribbon, and focusing performance is bad；
Fig. 5 (g) is the focusing effect figure that target after receives echo-signal is handled using existing two-dimensional frequency matching matrix,
The position of target is spread in the figure, and focus form is wide ribbon, and focusing performance is bad；
Fig. 5 (h) be target relocate results, it can be seen that since range curvature and range walk has occurred in target,
Its position deviates from highway, according to inter-bank velocity estimation value calculate ranging offset be Δ X=970, by compensation target away from
The refocusing to target is realized from offset.
Comparison diagram 5 (f), Fig. 5 (g) and Fig. 5 (e), it can be seen that deramp-keystone method and two-dimensional frequency matched filtering
Method is very poor to the focusing effect of target, and in contrast, the method for the present invention is under true environment to Ground moving target imaging effect
More preferably, compensation effect is more preferable.
In conclusion emulation experiment demonstrates correctness of the invention, validity and reliability.
Above description is only example of the present invention, does not constitute any limitation of the invention, it is clear that for
It, all may be without departing substantially from the principle of the invention, structure after having understood the content of present invention and principle for one of skill in the art
In the case where, carry out various modifications and change in form and details, but these modifications and variations based on inventive concept
Still within the scope of the claims of the present invention.
Claims (10)
1. a kind of quick non-search ground moving object method for parameter estimation characterized by comprising
(1) to based on SAR platform geometrical model the deformation of reception signal doing mathematics and abbreviation, obtain about frequency variation f and
Slow time tmReception signal S1(f,tm)；
(2) using azimuth-deramp adaptation function to reception signal S1(f,tm) matched filtering is carried out, eliminate doppler spectral point
It splits, the signal S ' after being matched1(f,tm)；
(3) to the signal S ' after matching1(f,tm) second order Keystone transformation is carried out, it realizes range curvature compensation, obtains about frequency
Rate variable quantity f and new slow time ηmSignal S2(f,ηm)；
(4) it adjusts the distance and is bent thermal compensation signal S2(f,ηm) phase difference PD processing is carried out, obtain single order phase signal S3(f,ηm)；
(5) to single order phase signal S3(f,ηm) execute the domain f IFFT transformation and ηmDomain FFT transform is realized range walk compensation, is obtained
To about fast time t and frequency faPeak detection signal S4(t,fa)；
(6) to peak detection signal S4(t,fa) detection is carried out, estimate doppler ambiguity numberWith the Secondary movement of target
Parameter
(7) using match matched filtering function to reception signal S1(f,tm) matched filtering and single order Keystone transformation are carried out,
To adjust the distance bending and Doppler shift compensates, obtain about frequency variation f and new slow time variable τmWithout frequency deviation
Signal S '4(f,τm)；
(8) to no frequency offset signal S '4(f,τm) execute the domain f IFFT transformation and τmDomain FFT transform is obtained comprising first order motion parameter
The peak detection signal S of information5(t,fa)；
(9) to the peak detection signal S comprising first order motion parameter information5(t,fa) detection is carried out, estimate the single order fortune of target
Dynamic parameter
2. the method as described in claim 1, which is characterized in that the reception signal based on SAR platform geometrical model in (1)
Doing mathematics deformation and abbreviation, carry out as follows:
(1a) sets Ri(tm) it is tmMoment SAR platform is at a distance from i-th of target, to Ri(tm) carry out Taylor expansion and retain first three
, it obtains:
In formula, R0iIt is zero moment SAR platform at a distance from i-th of target, a1iFor first order motion parameter, a2iFor Secondary movement ginseng
Number；
(1b) is to first order motion parameter a1iIt is decomposed, is obtained:
In formula, vb,iIndicate base band speed, Mamb,iIndicate doppler ambiguity number；
(1c) combines Ri(tm) Taylor expansion, a1iBreakdown and former receive signal, obtained reception signal S1(f,tm),
It is expressed as follows:
In formula, A1,iFor compressed amplitude, rect (x) is the window function that section [- 1/2,1/2] upper value is 1, and f is frequency variation
Amount, μ are frequency modulation rate, TpFor the pulse duration for emitting signal, tmFor slow time, TaFor pulse spacing, fcFor center frequency,
R0iIt is zero moment SAR platform at a distance from i-th of target, vb,iFor base band speed, a2iFor the Secondary movement parameter of target, Mamb,i
For doppler ambiguity number, c is the light velocity, and PRF is pulse recurrence frequency.
3. method as described in claim 1, which is characterized in that signal S ' after being matched obtained in step (2)1(f,tm), it indicates such as
Under:
In formula, A1,iFor compressed amplitude, rect (x) is the window function that section [- 1/2,1/2] upper value is 1, and f is frequency variation
Amount, μ are frequency modulation rate, TpFor the pulse duration for emitting signal, tmFor slow time, TaFor the pulse spacing, c is the light velocity, fcFor in
Frequency of heart, R0iIt is zero moment SAR platform at a distance from i-th of target, vb,iFor base band speed,
Mamb,iFor doppler ambiguity number, PRF is pulse recurrence frequency.
4. method as described in claim 1, which is characterized in that the range curvature thermal compensation signal S that step (3) obtains2(f,ηm), table
Show as follows:
In formula, f is frequency variation, and μ is frequency modulation rate, ηmFor new slow time, A1,iFor compressed amplitude, rect (x) is area
Between [- 1/2,1/2] upper value be 1 window function, TpFor the pulse duration for emitting signal, TaFor the pulse spacing, c is the light velocity, fc
For center frequency, R0iIt is zero moment SAR platform at a distance from i-th of target, λ is wavelength, vb,iFor base band speed,PRF is pulse recurrence frequency, Mamb,iFor doppler ambiguity number.
5. method as described in claim 1, which is characterized in that single order phase signal S obtained in step (4)3(f,ηm), it indicates
It is as follows:
In formula, A1,iFor compressed amplitude, rect (x) is the window function that section [- 1/2,1/2] upper value is 1, and f is frequency variation
Amount, μ are frequency modulation rate, TpFor the pulse duration for emitting signal, ηmFor new slow time, TaFor the pulse spacing, λ is wavelength,
vb,iFor base band speed, τ0For difference interval,C is the light velocity, Mamb,iFor doppler ambiguity number, PRF is
Pulse recurrence frequency, fcFor center frequency, S3,cross(f,ηm) it is cross term.
6. method as described in claim 1, which is characterized in that peak detection signal S obtained in step (5)4(t,fa), it indicates
It is as follows:
In formula, t is fast time, faFor frequency, A3,iFor amplitude, λ is wavelength, vb,iFor base band speed, τ0For difference interval,B is signal bandwidth, and c is the light velocity, Mamb,iFor doppler ambiguity number, PRF is pulse recurrence frequency, Ta
For pulse spacing, S4,cross(t,fa) it is cross term.
7. the method as described in claim 1, which is characterized in that target Secondary movement parameter obtained in step (6)Table
Show as follows:
In formula,For kinematic parameter aequ,iEstimated value, v be SAR platform speed, RbFor radar platform to band center line
Nearest oblique distance.
8. the method as described in claim 1, which is characterized in that without frequency offset signal S ' obtained in step (7)4(f,τm), it indicates
It is as follows:
In formula, f is frequency variation, τmFor new slow time variable, A1,iFor compressed amplitude, rect (x) is section [- 1/
2,1/2] window function that upper value is 1, μ are frequency modulation rate, TpFor the pulse duration for emitting signal, tmFor slow time, TaFor pulse
Interval, c is the light velocity, fcFor center frequency, R0iIt is zero moment SAR platform at a distance from i-th of target, vb,iFor base band speed.
9. the method as described in claim 1, which is characterized in that obtained in step (8) comprising first order motion parameter information
Peak detection signal S5(t,fa), it is expressed as follows:
In formula, t is fast time, faFor frequency, A4,iFor amplitude, B is signal bandwidth, R0iFor zero moment SAR platform and i-th of mesh
Target distance, TaFor pulse spacing, vb,iFor base band speed, λ is wavelength.
10. the method as described in claim 1, which is characterized in that target first order motion parameter obtained in step (9)Table
Show as follows:
In formula,For the estimated value of base band speed, λ is wavelength, and PRF is pulse recurrence frequency,For doppler ambiguity number
Estimated value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810847452.1A CN108983193A (en) | 2018-07-27 | 2018-07-27 | Quickly non-search ground moving object method for parameter estimation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810847452.1A CN108983193A (en) | 2018-07-27 | 2018-07-27 | Quickly non-search ground moving object method for parameter estimation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108983193A true CN108983193A (en) | 2018-12-11 |
Family
ID=64552060
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810847452.1A Pending CN108983193A (en) | 2018-07-27 | 2018-07-27 | Quickly non-search ground moving object method for parameter estimation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108983193A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110221264A (en) * | 2019-05-21 | 2019-09-10 | 中国人民解放军战略支援部队信息工程大学 | Doppler ambiguity target coherent detection method and device based on period Keystone transformation |
| CN111398913A (en) * | 2020-04-07 | 2020-07-10 | 电子科技大学 | Video SAR moving target echo filtering method based on double PRFs |
| CN112764020A (en) * | 2020-02-28 | 2021-05-07 | 加特兰微电子科技(上海)有限公司 | Method, device and related equipment for resolving speed ambiguity and determining moving speed of object |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1503223A1 (en) * | 2003-07-30 | 2005-02-02 | Raytheon Company | Estimation and correction of phase for focusing search mode SAR images formed by range migration algorithm |
| CN102998674A (en) * | 2012-10-29 | 2013-03-27 | 中国人民解放军空军装备研究院侦察情报装备研究所 | Method and device for detecting multi-channel SAR (synthetic aperture radar) slow ground moving target |
| EP2574957A1 (en) * | 2011-09-30 | 2013-04-03 | Thales | Method for estimating the unambiguous Doppler frequency of a moving target, in particular marine, and radar implementing said method |
| CN106970371A (en) * | 2017-04-28 | 2017-07-21 | 电子科技大学 | A kind of object detection method based on Keystone and matched filtering |
-
2018
- 2018-07-27 CN CN201810847452.1A patent/CN108983193A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1503223A1 (en) * | 2003-07-30 | 2005-02-02 | Raytheon Company | Estimation and correction of phase for focusing search mode SAR images formed by range migration algorithm |
| EP2574957A1 (en) * | 2011-09-30 | 2013-04-03 | Thales | Method for estimating the unambiguous Doppler frequency of a moving target, in particular marine, and radar implementing said method |
| CN102998674A (en) * | 2012-10-29 | 2013-03-27 | 中国人民解放军空军装备研究院侦察情报装备研究所 | Method and device for detecting multi-channel SAR (synthetic aperture radar) slow ground moving target |
| CN106970371A (en) * | 2017-04-28 | 2017-07-21 | 电子科技大学 | A kind of object detection method based on Keystone and matched filtering |
Non-Patent Citations (2)
| Title |
|---|
| XIONGPENGHE等: "Fast non-searching method for ground moving target refocusing and motion parameters estimation", 《DIGITAL SIGNAL PROCESSING》 * |
| 陈一畅等: "一种基于SAR稀疏采样数据的动目标运动参数估计方法", 《电子与信息学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110221264A (en) * | 2019-05-21 | 2019-09-10 | 中国人民解放军战略支援部队信息工程大学 | Doppler ambiguity target coherent detection method and device based on period Keystone transformation |
| CN110221264B (en) * | 2019-05-21 | 2021-04-16 | 中国人民解放军战略支援部队信息工程大学 | Doppler fuzzy target coherent detection method and device based on periodic Keystone transformation |
| CN112764020A (en) * | 2020-02-28 | 2021-05-07 | 加特兰微电子科技(上海)有限公司 | Method, device and related equipment for resolving speed ambiguity and determining moving speed of object |
| CN111398913A (en) * | 2020-04-07 | 2020-07-10 | 电子科技大学 | Video SAR moving target echo filtering method based on double PRFs |
| CN111398913B (en) * | 2020-04-07 | 2022-08-05 | 电子科技大学 | Video SAR moving target echo filtering method based on double PRFs |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107966688B (en) | Broadband radar target speed ambiguity resolving method based on phase interference technology | |
| CN110412558A (en) | The vehicle-mounted fmcw radar velocity ambiguity method of solution based on TDM MIMO | |
| CN1078949C (en) | Apparatus and method for mitigating range-doppler ambiguities in pulse-doppler radars | |
| CN106443671A (en) | SAR radar moving target detecting and imaging method based on FM continuous wave | |
| CN102288963B (en) | Bistatic inverse synthetic aperture radar (ISAR) image fusion method based on sub aperture parameter estimation | |
| CN109507669B (en) | Ground weak moving object parameter estimation method based on coherent accumulation | |
| CN106872974B (en) | High-precision motion target imaging method based on hypersonic platform Two-channels radar | |
| CN102998672B (en) | Step frequency inverse synthetic aperture radar (ISAR) imaging method based on coherent processing | |
| CN108469608B (en) | Method for accurately estimating Doppler centroid of moving platform radar | |
| CN105487060B (en) | A kind of multiple target extracting method of four Slope Modulation of binary channels | |
| CN108427115B (en) | Method for quickly estimating moving target parameters by synthetic aperture radar | |
| CN102707269B (en) | Range walk correction method for airborne radar | |
| CN104749570B (en) | It is a kind of to move constant airborne biradical synthetic aperture radar target localization method | |
| CN112162281A (en) | Multi-channel SAR-GMTI image domain two-step processing method | |
| CN109856635A (en) | CSAR ground moving target refocusing imaging method | |
| CN108983193A (en) | Quickly non-search ground moving object method for parameter estimation | |
| CN107271997B (en) | Airborne multi-channel CSSAR ground moving target motion parameter estimation method | |
| CN101620272B (en) | Target rotate speed estimation method of inverse synthetic aperture radar (ISAR) | |
| CN106597437A (en) | High-frequency airborne FM continuous wave SAR motion compensation and imaging processing method | |
| CN109946659B (en) | Vehicle-mounted millimeter wave radar linear frequency modulation continuous wave motion frequency expansion correction method | |
| CN102121990A (en) | Space-time analysis-based target rotation speed estimating method for inverse synthetic aperture radar | |
| CN111856466A (en) | Efficient ISAR (inverse synthetic aperture radar) translation compensation method for complex moving target | |
| CN106772373B (en) | For the SAR imaging method of any ground moving object | |
| CN112505647B (en) | Moving target azimuth speed estimation method based on sequential sub-image sequence | |
| CN107271996A (en) | A kind of airborne CSSAR Ground moving target imagings 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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181211 |