CN102628937A - Radar detection method based on generalized keystone transformation and non-coherent accumulation - Google Patents
Radar detection method based on generalized keystone transformation and non-coherent accumulation Download PDFInfo
- Publication number
- CN102628937A CN102628937A CN2012101183072A CN201210118307A CN102628937A CN 102628937 A CN102628937 A CN 102628937A CN 2012101183072 A CN2012101183072 A CN 2012101183072A CN 201210118307 A CN201210118307 A CN 201210118307A CN 102628937 A CN102628937 A CN 102628937A
- Authority
- CN
- China
- Prior art keywords
- target
- data
- velocity
- echo
- pulse
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000009825 accumulation Methods 0.000 title claims abstract description 45
- 230000001427 coherent effect Effects 0.000 title claims abstract description 33
- 238000001514 detection method Methods 0.000 title claims abstract description 22
- 230000009466 transformation Effects 0.000 title abstract 3
- 238000000034 method Methods 0.000 claims abstract description 41
- 230000006835 compression Effects 0.000 claims abstract description 17
- 238000007906 compression Methods 0.000 claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 230000001133 acceleration Effects 0.000 claims description 19
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 2
- 238000002592 echocardiography Methods 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 3
- 230000005012 migration Effects 0.000 abstract 1
- 238000013508 migration Methods 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 9
- 238000012937 correction Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005452 bending Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000035485 pulse pressure Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a radar detection method based on generalized keystone transformation and non-coherent accumulation, wherein the method mainly solves the detection problem of a high speed and highly accelerated target in a strong Gaussian white noise background. The method is implemented by the following steps: (1) generating radar target multiple-pulse echoes with Gaussian white noise added, and performing frequency domain pulse compression; (2) correcting range curvature of the pulse-compressed echoes by using the generalized keystone transformation; (3) performing fuzzy compensation on the curvature-corrected echoes by using a fuzzy factor; (4) performing non-coherent accumulation on the fuzzy-compensated echoes, and performing velocity compensation and slow time dimensional quadratic term compensation on the echoes by using a target velocity obtained in the non-coherent accumulation; (5) performing coherent accumulation on the slow time dimensional quadratic term compensated echoes to obtain the target to be detected. The method provided in the invention can eliminate effects of the range curvature and migration of the echoes on energy accumulation, can reduce computation of velocity search, and can be used in the detection of the target in the strong Gaussian white noise background.
Description
Technical field
The invention belongs to the Radar Technology field, relate to object detection method, can be used for efficiently handling in the white Gaussian noise environment based on the detections of radar of coherent pulse echo repeatedly.
Background technology
In the Radar Targets'Detection field; Can run into the situation that target to be detected is a high-speed motion; This class targets has very high speed and acceleration usually; Traditional method is that a plurality of pulse echos are directly carried out slow time dimension FFT, if range walk then can not take place when being low-speed motion in target, thus the purpose that can the energy accumulation of each time target echo got up and to reach target detection come out.But because target is a high-speed motion; Make and bigger range walk can take place at pulse internal object integration time; And the high acceleration of target also can cause the range curvature of echo envelope, if the paired pulses echo directly carries out slow time dimension FFT, can only get up the energy accumulation of same range unit; The signal gross energy is dispersed on a plurality of range units; The peak value of gained pulse accumulation is much lower when range walk and range curvature do not take place, thereby can not obtain good testing result, when serious even can lose objects; Existing research is introduced with the keystone method can get rid of the influence of the range walk of target generation to accumulation; But for reaching good accumulation, the range curvature of the echo envelope that the high acceleration in the time of also need be with target travel causes is removed the influence of accumulation, does not also have good method to solve at present.
Under the movable information condition of unknown of target; If obtain the velocity information of target through the method for accumulation; Usual method is, at first sets very big speed search scope and certain speed step-size in search, then each search speed of calculating of foundation echo envelope that can cause number of walking about; The paired pulses echo carries out the envelope ring shift; Again the pulse echo of carrying out after the envelope ring shift is done non-coherent accumulation, at last that each is the possible result that search speed accumulated compares, and the corresponding speed of maximum accumulation is the velocity estimation value of target.The result who does like this is to be cost with very high search arithmetic amount, therefore, is necessary to look for a kind of method that reduces operand.
Summary of the invention
The objective of the invention is to deficiency to above-mentioned prior art; A kind of radar detecting method based on broad sense keystone conversion and the accumulation of non-coherent is proposed; To reduce the operand of target velocity search; Get rid of high acceleration to the target energy effect of accumulation, thereby help improving radar the high detection performance of quickening target of high speed.
For realizing above-mentioned purpose, radar target detection method of the present invention comprises the steps:
1) approximate model and the noise model with the high accelerated motion target of radar high speed is expressed as respectively
And d (l, t
m);
And paired pulses echo
With the fast time
Be FFT for variable and transform to the frequency of distance territory, obtain frequency domain data z (f, the t of echo apart from dimension
m), wherein f representes frequency of distance;
3) will be apart from frequency domain data z (f, the t of dimension
m) do the frequency domain pulse compression along the distance dimension, the data after the pulse compression are designated as z ' (f, t
m), and to the data z ' after this pulse compression (f, t
m) in slow time variable t
mWith new time variable τ
mReplace, the variable alternate form does
The data that obtain after the variable replacement are designated as z ' (f, τ
m);
4) setting fuzzy factor range is [n
Min, n
Max], n wherein
MinMinimum, n
MaxMaximum, N altogether, and all be positive integer, with data z ' (f, the τ after the fuzzy factor pair variable replacement in this scope
m) carrying out fuzzy compensation, the echo data after the compensation is designated as z '
i(f, τ
m), i=1 ... N;
5) to the echo frequency domain data z ' behind the fuzzy compensation
i(f, τ
m), i=1 ... N carries out interior the inserting of time domain that p is ordered in each range unit, and the data after inserting in the time domain are designated as
I=1 ... N;
6) utilize data after inserting in the time domain, the echo envelope range walk is counted the umber of pulse w in k and the same range unit when calculating target velocity and being v
k
7) according to the target friction speed next time the pitch of waves from number k and the same range unit umber of pulse w of walking about
kDifference, the data after inserting in the time domain are carried out ring shift, and the group of the N after ring shift echo data are carried out non-coherent accumulation, obtain N and organize non-coherent accumulation, be designated as
I=1 ... N, wherein || expression be modulo operation;
8) N is organized non-coherent accumulation
I=1 ... N compares, and getting and wherein accumulating the maximum factor of corresponding fuzzy as a result of peak value is the velocity ambiguity factor of n of target
kEstimated value, be designated as
And the corresponding speed of the accumulation peak value of getting this maximum is fuzzy back target initial velocity v
0Estimated value, be designated as
9) utilize the objective fuzzy obtain after initial velocity estimated value
with the target initial velocity estimated value
that estimation obtains the echo frequency numeric field data
corresponding with fuzzy factor estimated value
carried out velocity compensation then with the estimated value
that fuzzy factor estimated value
is calculated the target initial velocity; Echo data after will compensating again transforms to apart from time domain, is designated as
apart from time domain data after the conversion
10) to conversion later apart from time domain data
Carry out the search of frequency modulation rate with the method for Dechirp, obtain frequency modulation rate γ
aEstimated value, be designated as
And then obtained the estimated value of target travel acceleration
11) with the estimated value
of the aimed acceleration that obtains the data after the velocity compensation
are carried out the compensation of slow time dimension quadratic phase item; And the data after the compensation are done the coherent accumulation of slow time dimension, the peak value of coherent accumulation is the target of detection.
The present invention can effectively suppress noise owing to utilize the method for pulse accumulation to carry out target detection; Simultaneously owing to utilize broad sense keystone conversion correction target envelope range curvature, make the result of non-coherent accumulation after the correction distance bending more near the target echo actual energy; Owing to utilize the method for velocity ambiguity compensation to dwindle the hunting zone of target initial velocity, make that the non-coherent accumulation operand after broad sense keystone conversion reduces greatly in addition.
Below in conjunction with accompanying drawing working of an invention is described in detail.
Description of drawings
Fig. 1 is a realization general flow chart of the present invention;
Fig. 2 is the testing result figure that adopts existing conventional sense method;
Fig. 3 is the testing result figure that adopts existing optimum detection methodology;
Fig. 4 is the testing result figure that adopts the inventive method.
Embodiment
With reference to Fig. 1, concrete performing step of the present invention is following:
Step 1, the radar target multiple-pulse echo of generation plus noise.
1a) the approximate model of radar target multiple-pulse echo:
In the formula,
T
pBe pulsewidth, γ is the frequency modulation rate,
Be the fast time, t is full-time, and m is an integer, T
rBe pulse-recurrence time, f
cBe radar emission signal carrier frequency, t
m=mT
rBe the slow time, A
0Be the intensity of scattering point echo,
Be t
mTarget is to the distance of radar, R constantly
0Be target initial time distance, v is a target velocity, and a is an aimed acceleration, and c is the light velocity;
1b) noise model:
Noise in the continuous N time echo is expressed as with matrix form:
Wherein, l is the range unit label, the noise power of l range unit of u (l) expression, g (l, t
m) be the noise phase of l range unit in m the pulse, column vector g (l) is and u (l) statistical independent M ties up the multiple Gaussian random variable of zero-mean, M representes the transponder pulse number, m=1,2 ..., M;
1c) utilize the approximate model of radar target multiple-pulse echo
And noise
Generate the radar target multiple-pulse echo of plus noise:
Owing to walking about from the unit and the range curvature inconvenience, so paired pulses echo apart from the skip distance of the direct processing target of time domain
Will be with the fast time
Be FFT for variable and transform to the frequency of distance territory, obtain frequency domain data z (f, t
m), what wherein FFT represented is Fast Fourier Transform (FFT), f representes frequency of distance.
3a) frequency domain data z (f, the t through following formula paired pulses echo
m) do the frequency domain pulse compression:
In the formula,
The indicating impulse compression function,
Be the fast time, T
pBe pulsewidth, γ is the frequency modulation rate;
3b) frequency of distance-slow time domain data with the pulse echo after the pulse compression are designated as z ' (f, t
m), these data have comprised target echo and noise information, and this instance only carries out theoretical analysis to target echo, so write out the concrete form of the later target echo signal of pulse pressure is:
In the formula, A ' expression is done the echo strength after the frequency pulse pressure to target echo, and Δ f is the radar emission signal bandwidth,
Be the Doppler frequency of target, λ is a radar wavelength, and target is to the distance of radar
R
0Be target initial time distance, v is a target velocity, and a is an aimed acceleration, t
mBe the slow time, f
cBe radar emission signal carrier frequency, c is the light velocity;
Because target high-speed motion; Make the envelope of target echo signal that the skew that is directly proportional with Doppler frequency can take place; The distance of its skew is
but the distance of this skew is identical for each time echo, does not therefore consider.
Step 4 is utilized broad sense keystone conversion correction target echo range curvature.
4a) because the rapid speed of target, and the repetition frequency of radar transmitted pulse is lower, therefore can cause target generation velocity ambiguity, so the Doppler frequency of target is expressed as:
f
d=f
d0+n
kPRF ,
In the formula, f
D0Be the later Doppler frequency of velocity ambiguity, n
kBe the fuzzy factor, PRF is a pulse repetition rate;
4b) with the expression of target echo after the expression substitution pulse compression of target Doppler frequency, the expression that the target echo that obtains is new is:
In the formula, A ' expression is done the echo strength after the frequency pulse pressure to target echo, and Δ f is the radar emission signal bandwidth,
Be the Doppler frequency of target, λ is a radar wavelength, and target is to the distance of radar
R
0Be target initial time distance, v is a target velocity, and a is an aimed acceleration, γ
aBe doppler frequency rate, t
mBe the slow time, f
cBe radar emission signal carrier frequency, c is the light velocity;
4c) phase place according to the 4th exponential term in the new expression of target echo after the pulse compression is the integral multiple of 2 π, and its value is 1, so after it is removed, obtain the new reduced representation formula s of target echo
3(f, t
m) be:
4c) because the target travel acceleration is very big, so pulse echo is easy to generate range curvature, thereby can influence the result that pulse energy accumulates, for eliminating this influence, introduce broad sense keystone conversion, i.e. paired pulses compression back data z ' (f, t
m) slow time variable t
mIntroduce the variable alternate form
Pulse echo data after the variable replacement are designated as z ' (f, τ
m);
4d) paired pulses compression back target echo signal s
3(f, t
m) carry out broad sense keystone conversion after, the target echo signal s that obtains
4(f, τ
m) expression be:
Step 5 is carried out fuzzy compensation to the echo data after the broad sense keystone conversion.
5a) paired pulses echo data z ' (f, τ
m) in the analysis of target echo signal expression can know that its 3rd exponential term is the velocity ambiguity item, need in this instance to compensate to it;
5b) the fuzzy factor of n of target setting
kScope [n
Min, n
Max], n wherein
MinBe the fuzzy factor estimated value of minimum, n
MaxBe the fuzzy factor estimated value of maximum, total N, and be integer all, utilize the fuzzy factor in this scope, constructed fuction G (f, τ
mn
k) to data z ' (f, τ after the broad sense keystone conversion
m) utilize following formula to carry out fuzzy compensation:
In the formula, G (f, τ
mn
k) expression fuzzy compensation function, n
kThe fuzzy factor of expression, τ
mBe the slow time after variable changes, f representes frequency of distance, f
cBe radar emission signal carrier frequency, PRF is the repetition frequency of radar emission signal, and the pulse echo data behind the fuzzy compensation are designated as z '
i(f, τ
m), i=1 ... N.
5c) with fuzzy compensation function G (f, τ
mn
k) target echo signal s to obtaining after the broad sense keystone conversion
4(f, τ
m) carry out fuzzy compensation, the target echo signal s behind the fuzzy compensation that obtains
5(f, τ
m) expression be:
Step 6 is to the pulse echo data z ' behind the fuzzy compensation
i(f, τ
m), i=1 ... N carries out apart from insert in the time domain, is about to the data z ' of fuzzy compensation afterpulse echo
i(f, τ
m), i=1 ... N carries out interior the inserting of time domain of p point in each range unit, p=1, and 2,3, the data after inserting in the time domain are designated as
I=1 ... N inserts theoretical see for details with method " radar imagery technology " in the time domain.
Step 7, the range walk of echo envelope was counted the umber of pulse w in k and the same range unit when calculating target velocity was v
k
7a) the pulse echo data after the broad sense keystone conversion are after carrying out the velocity ambiguity compensation, and the target velocity scope is 0≤v≤v
MaxSo, only need the speed in this scope is estimated, calculate the fuzzyyest speed of target
According to the fuzzyyest speed v that obtains
Max, calculating target echo envelope range walk number is l=(p+1) * round (2v
Maxf
s/ cPRF) with same range unit in umber of pulse w
l=M/l, wherein PRF representes the radar transmitted pulse repetition frequency, and λ representes radar wavelength, and what round represented is the operation that rounds that rounds up, f
sBe the signal sampling frequency;
7b) count the umber of pulse w in l and the same range unit according to target echo envelope range walk
l, confirm speed search step delta v:
At first, get any two adjacent value l in the several 0~l-1 scopes of range walk
1And l
2, l
1<l
2
Secondly, calculating generation range walk number is l
1The time corresponding target velocity
With the range walk number takes place is l
2The time corresponding target velocity
Then, use generation range walk number to be l
2The time corresponding target velocity v
2With the range walk number takes place is l
1The time corresponding target velocity v
1Difference as speed search step delta v;
Be step-length 7c), obtain 0≤v≤v with Δ v
MaxAll speed v in this scope, and the range walk of computing velocity target echo envelope when being v counts k=(p+1) * round (2v/ (cPRF) Δ fM), c representes the light velocity, Δ f representes the radar emission signal bandwidth;
The range walk of target echo envelope is counted k when 7d) being v according to speed, obtains the umber of pulse w in the same range unit
k=M/k, M indicating impulse number.
Umber of pulse w when at first, being v in the same range unit with the target velocity
kBe the segment length, M pulse is divided into the k group;
Then, the k group data after dividing are carried out successively the ring shift of 0~k-1 position along the distance dimension.
Step 9; Echo data after the ring shift is carried out non-coherent accumulation; Be about to mould value addition
i=1 in all pulse respective distances unit of the data after the N group ring shift ... N; Accomplish the accumulation of non-coherent, N organized non-coherent accumulation and can a peak value occur at a certain range unit this moment.
10a) will own [n
Min, n
Max] N that obtains after the fuzzy compensation factors organizes non-coherent accumulation
I=1 ... N compares, and selects the estimated value of the maximum accumulation corresponding fuzzy factor of peak value as the fuzzy factor of target
And get the corresponding speed of this peak-peak as objective fuzzy after the estimated value of speed
Calculate the estimated value of target initial velocity
In the formula; The estimated value of the Doppler frequency behind
expression target generation velocity ambiguity; λ representes radar wavelength,
estimated value of the fuzzy factor of expression;
Step 11, the echo data
to broad sense keystone conversion and compensating for doppler after fuzzy carries out velocity compensation.
After 11a) also compensating for doppler blurs to the bending of echo data correction distance, by pulse echo data z ' (f, τ
m) in the target echo expression can find out that there is linear coupling in its second exponential term speed with the slow time, for eliminating the influence that this coupling pulsed energy accumulation brings, be similar to
Wherein o () expression high-order is infinitely small, and the desin speed penalty function is:
In the formula,
Expression velocity compensation function,
The fuzzy back of expression target is the estimated value of Doppler frequency just, and f representes frequency of distance, τ
mIt is the slow time after the variable replacement;
11b) fall the exponential term that frequency of distance and doppler velocity in the pulse echo data
are coupled with velocity compensation function
compensation, the pulse echo data after the velocity compensation are designated as
11c) use the velocity compensation function
To target echo signal s
5(f, τ
m) carry out velocity compensation, obtain the target echo signal s after the velocity compensation
6(f, τ
m), be expressed as:
In the formula, the echo strength of A ' expression target, Δ f is the radar emission signal bandwidth, γ
aBe doppler frequency rate,
Be the Doppler frequency of target, λ is a radar wavelength, f
D0Doppler frequency after expression is fuzzy, R
0Be target initial time distance, v is a target velocity, and a is an aimed acceleration, τ
mBe the slow time after the variable replacement, f
cBe radar emission signal carrier frequency, c is the light velocity;
11d) the pulse echo data after the velocity compensation
are distance dimension IFFT; Convert the signal into apart from time domain; Obtain the distance of pulse echo-slow time domain data, be designated as
wherein IFFT represent inverse fast Fourier transform;
11e) to the target echo signal s after the velocity compensation
6(f, τ
m) do the distance-slow time-domain signal of the target echo that distance dimension IFFT obtains
Be expressed as:
" echo strength of target behind the IFFT is done in expression to target echo in the formula, A.
Step 12 is with Dechirp method estimating target motion acceleration.
The concrete realization of this step: the described Dechirp method of list of references " wide and narrow strip Radar Targets'Detection and imaging technique research ";
carries out the estimation of doppler frequency rate to the pulse echo data after the velocity compensation, and the estimated value of doppler frequency rate is designated as
Step 13 utilizes the estimated value
of doppler frequency rate to calculate the estimated value
of target travel acceleration
In the formula; The estimated value of
expression doppler frequency rate, λ representes radar wavelength.
Step 14 compensates the slow time dimension quadratic term of the echo data after the velocity compensation.
14a) distance of the echo signal echo the pulse echo after velocity compensation-slow time domain expression
can be found out; Also there is the quadratic phase item in pulse echo on slow time dimension; For the energy that makes echo can effectively accumulate, carry out the compensation of slow time dimension quadratic term with the data
of following formula after to velocity compensation:
In the formula,
The penalty function of representing slow time dimension quadratic phase item, τ
mBe the slow time after the variable replacement,
The estimated value of expression doppler frequency rate;
14b) the pulse echo data after the slow time dimension quadratic term compensation are designated as
14c) with the function
on the target echo signal
for slow time-dimensional quadratic term compensation target echo signal obtained after
be expressed as:
,
In the formula, A " echo strength of expression target,
Be the fast time, τ
mBe the slow time after the conversion, Δ f is the radar emission signal bandwidth,
Be the Doppler frequency of target, λ is a radar wavelength, f
D0Doppler frequency after expression is fuzzy, R
0Be target initial time distance, v is a target velocity, τ
mBe the slow time after the variable replacement, c is the light velocity.
Step 15, the echo data after the compensation of slow time dimension quadratic term are done the coherent accumulation.
15a) after paired pulses echo data
carries out broad sense keystone conversion, velocity ambiguity compensation, velocity compensation and slow time dimension quadratic term compensation successively; The pulse echo data
that obtain are done the accumulation of multiple-pulse coherent, and the peak value of coherent accumulation is detected target;
15b) target echo signal
is carried out the coherent accumulation, the expression that obtains detected echo signal
is:
In the formula, A " echo strength of expression target,
Be the fast time, f
τBe with τ
mBe the frequency after the variable FFT conversion, FFT representes Fast Fourier Transform (FFT), and γ representes the frequency modulation rate, and Δ f is the radar emission signal bandwidth,
Be the Doppler frequency of target, λ is a radar wavelength, and T is coherent integration time, f
D0Doppler frequency after expression is fuzzy, R
0Be target initial time distance, c is the light velocity.
Effect of the present invention further specifies through following emulation contrast test:
1. experiment scene:
Experimental data comprises with the white Gaussian noise of science software for calculation matlab R2008a emulation generation and the radar return of the high accelerated motion target of high speed; The target number is 1; The initial distance of the relative radar of target is 900 meters, and target radial speed is 3400m/s, and the target radial acceleration is 200m/s
2, the radar emission carrier frequency is 1GHz, and range resolution is 2 meters, and radar pulse repetition frequency is 500Hz, and a coherent comprises 100 pulses in integration time, and signal to noise ratio (S/N ratio) is-10dB.
2. emulation content:
Utilize experimental data; Carry out computer simulation experiment; Provide testing result such as Fig. 2, Fig. 3, shown in Figure 4 of conventional sense method, optimum detection methodology and detection method according to the invention respectively; Wherein Fig. 2 is for adopting the testing result of conventional sense method, and Fig. 3 is for adopting the testing result of optimum detection methodology, and Fig. 4 is for adopting the testing result of the method for the invention.
The conventional sense method is meant that the paired pulses echo data directly is slow time dimension FFT; And optimum detection methodology is meant the paired pulses echo range curvature and the range walk of pulse echo of having carried out compensation of broad sense keystone conversion, velocity ambiguity and velocity compensation post-equalization, but the slow time dimension FFT result of pulse echo who does not carry out slow time dimension quadratic term compensation.
3. analysis of simulation result:
As can beappreciated from fig. 2; When utilizing the conventional sense method to detect,, can cause the range walk and the bending of pulse echo because target is the high accelerated motion of high speed; So the paired pulses echo data can't accumulate out energy peak after directly adopting slow time dimension FFT, thereby can't detect target.
As can beappreciated from fig. 3; When utilizing optimum detection methodology to detect; Though proofreaied and correct the range curvature and the range walk of pulse echo; But do not carry out the compensation of slow time dimension quadratic phase item, thereby make target energy on slow time dimension, still disperse can't energy accumulation still can't be detected target.
As can beappreciated from fig. 4; After utilizing the method for the invention; Be after the paired pulses echo data compensates through broad sense keystone conversion, velocity ambiguity compensation, velocity compensation, slow time dimension quadratic term; The peak value of echo envelope all has been corrected to same range unit, thereby the echoed signal that is corrected to same range unit is carried out can the target echo energy efficient being accumulated when coherent accumulates, thereby has detected target.
In sum, when the high acceleration of high speed target was detected, the inventive method obviously was superior to conventional sense method and optimum detection methodology.
Claims (8)
1. the radar target detection method based on broad sense keystone conversion and the accumulation of non-coherent comprises the steps:
1) approximate model and the noise model with the high accelerated motion target of radar high speed is expressed as respectively
And d (l, t
m);
2), obtain the radar pulse echo according to target approximate model and noise model:
And paired pulses echo
With the fast time
Be FFT for variable and transform to the frequency of distance territory, obtain frequency domain data z (f, the t of echo apart from dimension
m), wherein f representes frequency of distance;
3) will be apart from frequency domain data z (f, the t of dimension
m) do the frequency domain pulse compression along the distance dimension, the data after the pulse compression are designated as z ' (f, t
m), and to the data z ' after this pulse compression (f, t
m) in slow time variable t
mWith new time variable τ
mReplace, the variable alternate form does
The data that obtain after the variable replacement are designated as z ' (f, τ
m);
4) setting fuzzy factor range is [n
Min, n
Max], n wherein
MinMinimum, n
MaxMaximum, N altogether, and all be positive integer, with data z ' (f, the τ after the fuzzy factor pair variable replacement in this scope
m) carrying out fuzzy compensation, the echo data after the compensation is designated as z '
i(f, τ
m), i=1 ... N;
5) to the echo frequency domain data z ' behind the fuzzy compensation
i(f, τ
m), i=1 ... N carries out interior the inserting of time domain that p is ordered in each range unit, and the data after inserting in the time domain are designated as
I=1 ... N;
6) utilize data after inserting in the time domain, the echo envelope range walk is counted the umber of pulse w in k and the same range unit when calculating target velocity and being v
k
7) according to the target friction speed next time the pitch of waves from number k and the same range unit umber of pulse w of walking about
kDifference, the data after inserting in the time domain are carried out ring shift, and the group of the N after ring shift echo data are carried out non-coherent accumulation, obtain N and organize non-coherent accumulation, be designated as
I=1 ... N, wherein || expression be modulo operation;
8) N is organized non-coherent accumulation
I=1 ... N compares, and getting and wherein accumulating the maximum factor of corresponding fuzzy as a result of peak value is the velocity ambiguity factor of n of target
kEstimated value, be designated as
And the corresponding speed of the accumulation peak value of getting this maximum is fuzzy back target initial velocity v
0Estimated value, be designated as
9) utilize the objective fuzzy obtain after initial velocity estimated value
with the target initial velocity estimated value
that estimation obtains the echo frequency numeric field data
corresponding with fuzzy factor estimated value
carried out velocity compensation then with the estimated value
that fuzzy factor estimated value
is calculated the target initial velocity; Echo data after will compensating again transforms to apart from time domain, is designated as
apart from time domain data after the conversion
10) to conversion later apart from time domain data
Carry out the search of frequency modulation rate with the method for Dechirp, obtain frequency modulation rate γ
aEstimated value, be designated as
And then obtained the estimated value of target travel acceleration
11) with the estimated value
of the aimed acceleration that obtains the data after the velocity compensation
are carried out the compensation of slow time dimension quadratic phase item; And the data after the compensation are done the coherent accumulation of slow time dimension, the peak value of coherent accumulation is the target of detection.
2. method according to claim 1, wherein step 3) is described will be apart from frequency domain data z (f, the t of dimension
m) do the frequency domain pulse compression along the distance dimension, carry out through following formula:
3. method according to claim 1, the data after wherein the data after one group of fuzzy factor pair variable replacement of the described usefulness of step 4) are replaced variable are carried out fuzzy compensation, utilize following formula to carry out:
In the formula, n
kThe fuzzy factor of expression, τ
mBe the slow time after variable changes, f representes frequency of distance, f
cBe radar emission signal carrier frequency, PRF is the repetition frequency of radar emission signal.
4. the umber of pulse in the method according to claim 1, wherein step 6) is described when calculating target velocity and being v echo envelope range walk number and same range unit, carry out as follows:
Calculating target echo envelope range walk number is l=(p+1) * round (2v
Maxf
s/ cPRF) with same range unit in umber of pulse w
l=M/l, wherein PRF representes the radar transmitted pulse repetition frequency, and λ representes radar wavelength, and what round represented is the operation that rounds that rounds up, f
sBe signal sampling frequency, p=1,2,3;
6b) count the umber of pulse w in l and the same range unit according to target echo envelope range walk
l, confirm speed search step delta v:
At first, get any two adjacent value l in the several 0~l-1 scopes of range walk
1And l
2, l
1<l
2
Secondly, calculating generation range walk number is l
1The time corresponding target velocity
With the range walk number takes place is l
2The time corresponding target velocity
Then, use generation range walk number to be l
2The time corresponding target velocity v
2With the range walk number takes place is l
1The time corresponding target velocity v
1Difference as speed search step delta v;
Be step-length with Δ v 6c), obtaining the target velocity scope is 0≤v≤v
MaxAll speed v in this scope, and the range walk of computing velocity target echo envelope when being v counts k=(p+1) * round (2v/ (cPRF) Δ fM), c representes the light velocity, Δ f representes the radar emission signal bandwidth;
The range walk of target echo envelope is counted k when 6d) being v according to speed, obtains the umber of pulse w in the same range unit
k=M/k, M indicating impulse number.
5. method according to claim 1; Wherein described in the step 7) according to target friction speed next time the pitch of waves from the difference of the number of walking about; Data to after inserting in the time domain are carried out ring shift; Be to N group data
i=1 after inserting in the time domain ... N, all operate as follows:
Umber of pulse w when at first, being v in the same range unit with the target velocity
kBe the segment length, M pulse is divided into the k group;
Then, the k group data after dividing are carried out successively the ring shift of 0~k-1 position along the distance dimension.
6. method according to claim 1, calculate the estimated value
of target initial velocity and the estimated value
of calculating target travel acceleration in the wherein said step 8) and use computes:
In the formula; The estimated value of first Doppler frequency behind
expression target generation velocity ambiguity; λ representes radar wavelength; The estimated value of the fuzzy factor of
expression, the estimated value of
expression doppler frequency rate.
7. method according to claim 1, wherein the described target initial velocity estimated value pair echo frequency numeric field data corresponding with fuzzy factor estimated value that obtains with estimation of step 8) carried out velocity compensation, utilizes following formula to carry out:
8. according to right 1 described method, wherein the data of the estimated value of the aimed acceleration that obtains of the described usefulness of step 11) after to velocity compensation are carried out the compensation of slow time dimension quadratic phase item, utilize following formula to carry out:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210118307.2A CN102628937B (en) | 2012-04-20 | 2012-04-20 | Radar detection method based on generalized keystone transformation and non-coherent accumulation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210118307.2A CN102628937B (en) | 2012-04-20 | 2012-04-20 | Radar detection method based on generalized keystone transformation and non-coherent accumulation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102628937A true CN102628937A (en) | 2012-08-08 |
CN102628937B CN102628937B (en) | 2014-02-12 |
Family
ID=46587230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210118307.2A Expired - Fee Related CN102628937B (en) | 2012-04-20 | 2012-04-20 | Radar detection method based on generalized keystone transformation and non-coherent accumulation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102628937B (en) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323829A (en) * | 2013-06-04 | 2013-09-25 | 中国人民解放军海军航空工程学院 | Radar moving target long-time phase-coherent accumulation detecting method based on RFRAF |
CN103760555A (en) * | 2014-01-23 | 2014-04-30 | 西安电子科技大学 | Method for improving airborne radar detection and tracking integration precision |
CN104076351A (en) * | 2014-06-30 | 2014-10-01 | 电子科技大学 | Phase-coherent accumulation detection method for high-speed high maneuvering target |
CN104215959A (en) * | 2014-09-22 | 2014-12-17 | 西安电子科技大学 | Method for estimating radial initial speeds and radial acceleration of multiple maneuvering targets |
CN104730498A (en) * | 2015-04-01 | 2015-06-24 | 西安电子科技大学 | Target detection method based on Keystone and weighting rotating FFT |
CN103810325B (en) * | 2014-01-08 | 2016-08-17 | 西安电子科技大学 | The linear thinned array antenna optimization method of low sidelobe based on SQP |
CN106597403A (en) * | 2016-11-29 | 2017-04-26 | 西安电子工程研究所 | High-velocity target coherent accumulation detection method based on piecewise compensation |
CN106646447A (en) * | 2017-01-18 | 2017-05-10 | 武汉雷博合创电子技术有限公司 | Detection method for radar target long-time accumulation based on linear frequency modulation continuous wave |
CN106970371A (en) * | 2017-04-28 | 2017-07-21 | 电子科技大学 | A kind of object detection method based on Keystone and matched filtering |
CN107356908A (en) * | 2017-06-23 | 2017-11-17 | 中国电子科技集团公司第二十研究所 | A kind of frequency agility signal correlative accumulation method |
CN107450055A (en) * | 2017-07-15 | 2017-12-08 | 西安电子科技大学 | High-speed maneuver object detection method based on Discrete Linear frequency modulation Fourier transform |
CN107561508A (en) * | 2017-08-24 | 2018-01-09 | 电子科技大学 | A kind of phase-coherent accumulation detection method for even accelerated motional objects |
CN108398676A (en) * | 2018-05-04 | 2018-08-14 | 电子科技大学 | A kind of external illuminators-based radar Dim moving target detection method |
CN108549064A (en) * | 2018-07-24 | 2018-09-18 | 电子科技大学 | External sort algorithm moving-target detection method based on Doppler frequency fuzzy compensation in arteries and veins |
CN108776332A (en) * | 2018-04-24 | 2018-11-09 | 国家***第海洋研究所 | The method for detecting maneuver targets on sea using high-frequency ground wave radar |
CN108896976A (en) * | 2018-07-05 | 2018-11-27 | 电子科技大学 | A kind of coherent processing method for multichannel external illuminators-based radar |
CN108896971A (en) * | 2018-05-10 | 2018-11-27 | 西安电子科技大学 | A kind of emulation mode of offshore floating Small object echo |
CN109085590A (en) * | 2018-10-16 | 2018-12-25 | 中国人民解放军国防科技大学 | Broadband direct acquisition data ISAR imaging method based on ARP (Address resolution protocol) segmented coherent accumulation |
CN109765540A (en) * | 2019-02-26 | 2019-05-17 | 南京莱斯电子设备有限公司 | A kind of frequency stepping system metre wave radar target extraction method |
CN112327289A (en) * | 2020-10-29 | 2021-02-05 | 内蒙古工业大学 | Method and device for estimating slope distance and speed of moving target |
CN113203998A (en) * | 2021-04-23 | 2021-08-03 | 上海交通大学 | ISAR translation compensation and imaging method, system, medium and device |
CN113238191A (en) * | 2021-05-10 | 2021-08-10 | 电子科技大学 | Double-base MIMO radar echo non-coherent accumulation method based on cuckoo search |
CN115453462A (en) * | 2022-07-28 | 2022-12-09 | 西安电子科技大学 | External radiation source radar echo signal time frequency difference parameter rapid calculation method |
CN115685169A (en) * | 2022-11-09 | 2023-02-03 | 哈尔滨工程大学 | Underwater sound weak moving target detection method based on broadband keystone transformation |
CN116736297A (en) * | 2023-08-09 | 2023-09-12 | 中国科学院空天信息创新研究院 | Heterogeneous multi-frame joint phase-coherent accumulation method |
US12007468B2 (en) | 2021-12-10 | 2024-06-11 | Raytheon Company | Range-doppler keystone processing for direct sampled radar data from targets with long range and high velocity using waveforms with high bandwidth, high duty factor, and long dwell |
CN115453462B (en) * | 2022-07-28 | 2024-07-30 | 西安电子科技大学 | Rapid calculation method for time-frequency difference parameters of radar echo signals of external radiation source |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080174475A1 (en) * | 2006-10-20 | 2008-07-24 | David Charles Clark | Signal processing for accelerating moving targets |
US20100259442A1 (en) * | 2009-04-13 | 2010-10-14 | Abatzoglou Theagenis J | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters |
CN102073043A (en) * | 2010-11-04 | 2011-05-25 | 电子科技大学 | Multi-frame phase coherence accumulation target tracking-before-detecting method |
CN102426354A (en) * | 2011-09-16 | 2012-04-25 | 西安电子科技大学 | Broadband radar detection method based on weighted sequence statistics and multiple-pulse coherence accumulation |
-
2012
- 2012-04-20 CN CN201210118307.2A patent/CN102628937B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080174475A1 (en) * | 2006-10-20 | 2008-07-24 | David Charles Clark | Signal processing for accelerating moving targets |
US20100259442A1 (en) * | 2009-04-13 | 2010-10-14 | Abatzoglou Theagenis J | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters |
CN102073043A (en) * | 2010-11-04 | 2011-05-25 | 电子科技大学 | Multi-frame phase coherence accumulation target tracking-before-detecting method |
CN102426354A (en) * | 2011-09-16 | 2012-04-25 | 西安电子科技大学 | Broadband radar detection method based on weighted sequence statistics and multiple-pulse coherence accumulation |
Non-Patent Citations (1)
Title |
---|
戴奉周 等: "宽带和窄带雷达在噪声中的检测性能比较", 《电子与信息学报》 * |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323829A (en) * | 2013-06-04 | 2013-09-25 | 中国人民解放军海军航空工程学院 | Radar moving target long-time phase-coherent accumulation detecting method based on RFRAF |
CN103323829B (en) * | 2013-06-04 | 2015-09-09 | 中国人民解放军海军航空工程学院 | Based on the radar moving targets long-time phase-coherent accumulation detection method of Radon-fractional order ambiguity function |
CN103810325B (en) * | 2014-01-08 | 2016-08-17 | 西安电子科技大学 | The linear thinned array antenna optimization method of low sidelobe based on SQP |
CN103760555A (en) * | 2014-01-23 | 2014-04-30 | 西安电子科技大学 | Method for improving airborne radar detection and tracking integration precision |
CN103760555B (en) * | 2014-01-23 | 2016-05-18 | 西安电子科技大学 | A kind of method that improves airborne radar detection Tracking Integrative precision |
CN104076351A (en) * | 2014-06-30 | 2014-10-01 | 电子科技大学 | Phase-coherent accumulation detection method for high-speed high maneuvering target |
CN104215959A (en) * | 2014-09-22 | 2014-12-17 | 西安电子科技大学 | Method for estimating radial initial speeds and radial acceleration of multiple maneuvering targets |
CN104730498A (en) * | 2015-04-01 | 2015-06-24 | 西安电子科技大学 | Target detection method based on Keystone and weighting rotating FFT |
CN106597403A (en) * | 2016-11-29 | 2017-04-26 | 西安电子工程研究所 | High-velocity target coherent accumulation detection method based on piecewise compensation |
CN106597403B (en) * | 2016-11-29 | 2018-12-11 | 西安电子工程研究所 | A kind of high-speed target phase-coherent accumulation detection method based on segmented compensation |
CN106646447A (en) * | 2017-01-18 | 2017-05-10 | 武汉雷博合创电子技术有限公司 | Detection method for radar target long-time accumulation based on linear frequency modulation continuous wave |
CN106646447B (en) * | 2017-01-18 | 2019-03-26 | 武汉雷博合创电子科技有限公司 | Radar target long time integration detection method based on linear frequency modulation continuous wave |
CN106970371A (en) * | 2017-04-28 | 2017-07-21 | 电子科技大学 | A kind of object detection method based on Keystone and matched filtering |
CN106970371B (en) * | 2017-04-28 | 2019-05-14 | 电子科技大学 | A kind of object detection method based on Keystone and matched filtering |
CN107356908A (en) * | 2017-06-23 | 2017-11-17 | 中国电子科技集团公司第二十研究所 | A kind of frequency agility signal correlative accumulation method |
CN107356908B (en) * | 2017-06-23 | 2020-07-03 | 中国电子科技集团公司第二十研究所 | Frequency agile signal coherent accumulation method |
CN107450055A (en) * | 2017-07-15 | 2017-12-08 | 西安电子科技大学 | High-speed maneuver object detection method based on Discrete Linear frequency modulation Fourier transform |
CN107450055B (en) * | 2017-07-15 | 2020-04-14 | 西安电子科技大学 | High-speed maneuvering target detection method based on discrete linear frequency modulation Fourier transform |
CN107561508A (en) * | 2017-08-24 | 2018-01-09 | 电子科技大学 | A kind of phase-coherent accumulation detection method for even accelerated motional objects |
CN108776332A (en) * | 2018-04-24 | 2018-11-09 | 国家***第海洋研究所 | The method for detecting maneuver targets on sea using high-frequency ground wave radar |
CN108776332B (en) * | 2018-04-24 | 2021-05-28 | 国家***第一海洋研究所 | Method for detecting marine maneuvering target by using high-frequency ground wave radar |
CN108398676A (en) * | 2018-05-04 | 2018-08-14 | 电子科技大学 | A kind of external illuminators-based radar Dim moving target detection method |
CN108398676B (en) * | 2018-05-04 | 2021-10-26 | 电子科技大学 | External radiation source radar weak moving target detection method |
CN108896971A (en) * | 2018-05-10 | 2018-11-27 | 西安电子科技大学 | A kind of emulation mode of offshore floating Small object echo |
CN108896976A (en) * | 2018-07-05 | 2018-11-27 | 电子科技大学 | A kind of coherent processing method for multichannel external illuminators-based radar |
CN108549064A (en) * | 2018-07-24 | 2018-09-18 | 电子科技大学 | External sort algorithm moving-target detection method based on Doppler frequency fuzzy compensation in arteries and veins |
CN109085590A (en) * | 2018-10-16 | 2018-12-25 | 中国人民解放军国防科技大学 | Broadband direct acquisition data ISAR imaging method based on ARP (Address resolution protocol) segmented coherent accumulation |
CN109765540B (en) * | 2019-02-26 | 2021-03-19 | 南京莱斯电子设备有限公司 | Frequency stepping system meter wave radar target extraction method |
CN109765540A (en) * | 2019-02-26 | 2019-05-17 | 南京莱斯电子设备有限公司 | A kind of frequency stepping system metre wave radar target extraction method |
CN112327289A (en) * | 2020-10-29 | 2021-02-05 | 内蒙古工业大学 | Method and device for estimating slope distance and speed of moving target |
CN113203998A (en) * | 2021-04-23 | 2021-08-03 | 上海交通大学 | ISAR translation compensation and imaging method, system, medium and device |
CN113238191A (en) * | 2021-05-10 | 2021-08-10 | 电子科技大学 | Double-base MIMO radar echo non-coherent accumulation method based on cuckoo search |
CN113238191B (en) * | 2021-05-10 | 2022-06-24 | 电子科技大学 | Double-base MIMO radar echo non-coherent accumulation method based on cuckoo search |
US12007468B2 (en) | 2021-12-10 | 2024-06-11 | Raytheon Company | Range-doppler keystone processing for direct sampled radar data from targets with long range and high velocity using waveforms with high bandwidth, high duty factor, and long dwell |
CN115453462A (en) * | 2022-07-28 | 2022-12-09 | 西安电子科技大学 | External radiation source radar echo signal time frequency difference parameter rapid calculation method |
CN115453462B (en) * | 2022-07-28 | 2024-07-30 | 西安电子科技大学 | Rapid calculation method for time-frequency difference parameters of radar echo signals of external radiation source |
CN115685169A (en) * | 2022-11-09 | 2023-02-03 | 哈尔滨工程大学 | Underwater sound weak moving target detection method based on broadband keystone transformation |
CN116736297A (en) * | 2023-08-09 | 2023-09-12 | 中国科学院空天信息创新研究院 | Heterogeneous multi-frame joint phase-coherent accumulation method |
CN116736297B (en) * | 2023-08-09 | 2023-10-10 | 中国科学院空天信息创新研究院 | Heterogeneous multi-frame joint phase-coherent accumulation method |
Also Published As
Publication number | Publication date |
---|---|
CN102628937B (en) | 2014-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102628937B (en) | Radar detection method based on generalized keystone transformation and non-coherent accumulation | |
CN106872974B (en) | High-precision motion target imaging method based on hypersonic platform Two-channels radar | |
CN106338727B (en) | A kind of vehicle-mounted auxiliary drives the object detection method of radar | |
CN102998672B (en) | Step frequency inverse synthetic aperture radar (ISAR) imaging method based on coherent processing | |
CN104076351B (en) | Phase-coherent accumulation detection method for high-speed high maneuvering target | |
CN101369017B (en) | Shift variant mode double-base synthetic aperture radar imaging method | |
CN102590812B (en) | SAR (synthetic aperture radar) real-time imaging method based on frequency modulated continuous wave | |
CN107132534B (en) | Optimization method for high-speed radar target frequency domain detection | |
CN101738606B (en) | Method for detecting coherent integration of radar target based on generalized Doppler filter bank | |
CN101320086B (en) | Echo signal processing equipment and method of Doppler speed measuring laser radar | |
CN101592733B (en) | Parallel real-time imaging processing method for inverse synthetic aperture radar | |
CN105738889B (en) | A kind of CW with frequency modulation tests the speed distance-finding method | |
CN103399310A (en) | Method for detecting radar weak moving target based on PD (Phase Differentiation) RLVD (Radon-Lv Distribution) | |
CN104360336A (en) | Novel method for extracting radar target micro-motion cycle in self-adaptive mode | |
CN105549002B (en) | A kind of frequency modulated continuous wave radar measuring method based on combined waveform | |
CN102121989B (en) | Method for simultaneously measuring speed and distance of high-speed moving targets | |
CN104849708A (en) | High-speed maneuvering target parameter estimation method based on frequency domain polynomial phase transformation | |
CN101915911A (en) | Space arbitrary configuration distributed SAR moving target parameter estimation method based on cancellation accumulation space-time spectrum | |
Wang et al. | A method of velocity estimation using composite hyperbolic frequency-modulated signals in active sonar | |
CN109613507B (en) | Detection method for high-order maneuvering target radar echo | |
CN110398730A (en) | Based on coordinate rotation and Nonuniform fast Fourier transform maneuvering target coherent detection method | |
CN102778674A (en) | Chirp pulse time delay estimation method for non-uniform sampling | |
CN109375206A (en) | A kind of moving target speed-measuring method based on speed search | |
CN104375128B (en) | Fast high maneuvering target accumulating and detecting method based on cross-correlation functions | |
CN101266294B (en) | Compensation process for inhibiting radar echo signal range migration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140212 |
|
CF01 | Termination of patent right due to non-payment of annual fee |