CN105068072B - A kind of speed compensation method of the one-dimensional range profile of moving target - Google Patents
A kind of speed compensation method of the one-dimensional range profile of moving target Download PDFInfo
- Publication number
- CN105068072B CN105068072B CN201510434889.9A CN201510434889A CN105068072B CN 105068072 B CN105068072 B CN 105068072B CN 201510434889 A CN201510434889 A CN 201510434889A CN 105068072 B CN105068072 B CN 105068072B
- Authority
- CN
- China
- Prior art keywords
- range profile
- dimensional range
- buffer
- dimensional
- velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
-
- 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/40—Means for monitoring or calibrating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9064—Inverse SAR [ISAR]
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of speed compensation method of the one-dimensional range profile of moving target, the method includes:The one-dimensional range profile sequence that radar is continuously obtained within a preset time interval is cached in buffer, and one-dimensional range profile sequence to being cached in buffer is pre-processed, and obtains the analysis window of one-dimensional range profile, and removes noise;Then, all one-dimensional range profile sequences within analysis window in buffer are carried out into cross correlation process;Estimate target envelope amount of movement, and plus reference signal time delay correspondence position, obtain the actual range of target;The velocity estimation value of one-dimensional range profile is obtained by the method for a curve matching;Velocity estimation value according to one-dimensional range profile, velocity compensation is carried out to one-dimensional range profile;Variance evaluation is carried out to carrying out each one-dimensional range profile after velocity compensation, it is determined that the umber of pulse as velocity compensation meets error requirements.By using method provided by the present invention, the adverse effect that high-speed motion brings can be effectively eliminated.
Description
Technical field
The present invention relates to radar target feature identification technique, more particularly to a kind of speed of the one-dimensional range profile of moving target
Compensation method.
Background technology
ISAR (ISAR) imaging technique grows up in Synthetic Aperture Radar Technique.ISAR is imaged
For radar is static, the imaging system of target motion, different from conventional radar, ISAR is a kind of high-resolution imaging radar to system,
The exact image of non-cooperative moving targets can at a distance be obtained.
ISAR imaging techniques as a kind of valuable detection means in the past few decades in obtained extensive attention.
A large amount of in running order in-orbit spacecrafts are dispersed with space, space junk, aerolite in addition with enormous amount and
Failure spacecraft etc., is the important content of space exploration to the imaging of these targets.But the research of ISAR technologies so far
It is concentrated mainly on the relatively low aerial target of speed (such as aircraft), the ISAR imagings of the Space Object in Orbit higher to speed are ground
Study carefully also less.
Extraterrestrial target is target speed very high with the significant difference of tradition ISAR imaging objects, than the fortune of conventional target
Dynamic speed improves an order of magnitude so that target distance image is produced fuzzy and distorted, and then causes ISAR images radially to defocus,
Therefore the influence that high-speed motion brings must take into consideration to extraterrestrial target imaging.
When measured target radial motion speed is relatively low, the echo after De-chirp method is approximately sinusoidal signal, by arteries and veins
Punching press contracting can just obtain a Range Profile of target;But when object of observation is the high-speed moving objects such as space orbit spacecraft,
Because target is very high relative to the radial velocity of radar so that receive scattering point echo for linear frequency modulation (LFM) signal (is no longer
Simple signal under low-speed motion), the Range Profile Severe distortion of the high-speed moving object obtained after process of pulse-compression.
The content of the invention
In view of this, the present invention provides a kind of speed compensation method of the one-dimensional range profile of moving target, such that it is able to have
Effect ground eliminates the influence that high-speed motion brings.
What technical scheme was specifically realized in:
A kind of speed compensation method of the one-dimensional range profile of moving target, the method includes:
Step A, the one-dimensional range profile sequence for continuously obtaining radar within a preset time interval are cached in buffer, and
One-dimensional range profile sequence to being cached in buffer is pre-processed, and obtains the analysis window of one-dimensional range profile, and removes noise;
Step B, all one-dimensional range profile sequences within analysis window in buffer are carried out into cross correlation process;Estimate
Target envelope amount of movement, and plus reference signal time delay correspondence position, obtain the actual range of target;
Step C, the velocity estimation value that one-dimensional range profile is obtained by the method for a curve matching;
Step D, the velocity estimation value according to one-dimensional range profile, velocity compensation is carried out to one-dimensional range profile;
Step E, variance evaluation is carried out to carrying out each one-dimensional range profile after velocity compensation, it is determined that as velocity compensation
Umber of pulse meets error requirements.
Optionally, the one-dimensional range profile sequence to caching in buffer carries out pretreatment includes:
One-dimensional range profile course figure according to the one-dimensional range profile sequence cached in buffer is into gray level image;
The edge of the gray level image is extracted using image processing method, the starting of the analysis window of one-dimensional range profile is obtained
Position and final position, and it is removed the operation of noise.
Optionally, cross correlation process is carried out to the one-dimensional range profile sequence within analysis window using formula below:
Wherein, r(k)(nk) it is cross correlation process result, nkIt is cross correlation process variable, represents to ξ 'H (k)Carry out nkIt is individual away from
From unit movement, JmaxRepresent the size of buffer;ξ′H (k-1) iRepresent the one-dimensional range profile after normalization.
Optionally, the actual range p of target is obtained using formula belowk:
pk=pk-1+cnk/ 2 Δ f, k=1,2 ..., Jmax-1;
Wherein, p0It is the target range of initial time, c is the light velocity, and Δ f is radar swept bandwidth.
Optionally, the velocity estimation value according to one-dimensional range profile, carrying out velocity compensation to one-dimensional range profile includes:
Velocity estimation value according to one-dimensional range profile, translation compensation treatment is carried out to each one-dimensional range profile;
Make FFT to the one-dimensional range profile after translation compensation treatment to estimate, obtain the phase difference ΔΦ of adjacent echoes;
Preset the first threshold value Z*, as ΔΦ > Z*When, return and perform step B;And work as ΔΦ≤Z*When, perform step
Rapid E.
Optionally, the step E includes:
Variance evaluation is carried out to carrying out each one-dimensional range profile after velocity compensation, the variance of each one-dimensional range profile is obtained
Estimate;
Estimate of variance to all one-dimensional range profiles is averaged, and obtains mean of variance;
Pre-set the second threshold value NN*, whenWhen, increase buffer size and return to execution step A;WhenWhen, terminate whole flow process.
Optionally, variance evaluation is carried out by formula below:
Wherein,Represent that i-th one-dimensional range profile is J in sizemaxBuffer in scattering strength it is flat
Average.
Optionally, mean of variance is obtained by formula below:
As above it is visible, in the speed compensation method of the one-dimensional range profile of moving target of the invention, due to radar being existed
The continuous one-dimensional range profile sequence for obtaining is cached in buffer in prefixed time interval, and to cached in buffer it is one-dimensional away from
Pre-processed from as sequence, obtained the analysis window of one-dimensional range profile, and removed noise;Then, by buffer it is all
One-dimensional range profile sequence within analysis window carries out cross correlation process;Estimate target envelope amount of movement, and plus reference signal
Time delay correspondence position, obtains the actual range of target;The speed for obtaining one-dimensional range profile by the method for a curve matching is estimated
Evaluation;Velocity estimation value according to one-dimensional range profile, velocity compensation is carried out to one-dimensional range profile;It is every after velocity compensation to carrying out
Individual one-dimensional range profile carries out variance evaluation, it is determined that the umber of pulse as velocity compensation meets error requirements, therefore can be to echo
Signal carries out velocity compensation, effectively eliminates the adverse effect that high-speed motion brings.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of the speed compensation method of the one-dimensional range profile of the moving target in the embodiment of the present invention.
Fig. 2 is the speed compensation method effect diagram of the one-dimensional range profile of the moving target in the embodiment of the present invention.
Specific embodiment
To make the objects, technical solutions and advantages of the present invention become more apparent, develop simultaneously embodiment referring to the drawings, right
The present invention is further described.
In the inventive solutions, the target echo of broadband signal is one-dimensional range profile, therefore, it is difficult to as arrowband feelings
Condition like that, is tracked with range gate.So, in the inventive solutions, adjacent twice return (i.e. can be used
Dimension Range Profile) " slide related " be tracked." sliding related " of twice return, similar to " auto-correlation " said in reception,
Except that auto-correlation is made time delay slip only with echo itself, with time delay slippage as variable, two are calculated
Coefficient correlation, the position where representing target with the peak value of coefficient correlation, and estimate the speed of target.
Motion compensation is the key technology in ISAR imaging techniques, generally includes envelope alignment and phase compensation.In this hair
In bright technical scheme, the experience of the envelope alignment in ISAR imaging techniques, envelope alignment is generally related using sliding, and
Maximum as alignment criterion using coefficient correlation, its alignment error (is regarded with the time interval of two echoes (i.e. one-dimensional range profile)
Angle change) size it is relevant.When adjacent echoes visual angle change is minimum, Delay Estima-tion error is very little, can be compared
± 1/4 Range resolution unit is much smaller.But, this refers to signal to noise ratio situation very high.When to Long Range Target Tracking, noise
Than much lower, Delay Estima-tion error now essentially from the influence of noise, cause rather than changing due to visual angle away from
From as changing, signal to noise ratio is lower, then evaluated error is also bigger, about ± 1/2 Range resolution unit.Due to the bag that noise causes
The error of network alignment is entirely random, in the absence of correlation, averagely can reduce it by multiple pulses,.
In the time of submicrosecond level, radial velocity substantially constant, its radial acceleration can be ignored extraterrestrial target, this
When, the range difference of two neighboring echo time delay should be constant, thus multiecho can one after the other be made adjacent slip correlation, obtain
A series of delay inequality, is then averaging.If each time noise is statistical iteration, the Delay Estima-tion error that it causes also not phase
Close, average with n times, the variance of error can be reduced to 1/N.If the time is more long, in addition to the error that noise causes, visual angle change
The Delay Estima-tion error for causing will also account for certain ratio, and the error that it causes has correlation, and average effect is less better.
A kind of speed compensation method of the one-dimensional range profile of moving target is present embodiments provided, such that it is able to effectively disappear
Except the influence that high-speed motion brings.
Fig. 1 is the schematic flow sheet of the speed compensation method of the one-dimensional range profile of the moving target in the embodiment of the present invention.
As shown in figure 1, the speed compensation method of the one-dimensional range profile of moving target in the embodiment of the present invention can include it is as described below
The step of:
Step 101, the one-dimensional range profile sequence that radar is continuously obtained within a preset time interval is cached in buffer,
And the one-dimensional range profile sequence to being cached in buffer is pre-processed, the analysis window of one-dimensional range profile is obtained, and removal is made an uproar
Sound.
In the inventive solutions, radar can continuously obtain multiple one-dimensional range profiles within a preset time interval,
Wherein, single one-dimensional range profile can be expressed as:
ξ={ ξ1, ξ2..., ξN} (1)
Multiple one-dimensional range profiles that above-mentioned radar is obtained constitute an one-dimensional range profile sequence and are cached in buffer,
Then one-dimensional range profile sequence can be expressed as in buffer:
ξ(k)={ ξ1 (k), ξ2 (k)..., ξN (k), k=0,1 ..., Jmax-1 (2)
Wherein, JmaxRepresent the size (maximum number of i.e. storable one-dimensional range profile) of buffer.
In addition, one-dimensional range profile sequence is also denoted as in buffer:
ξ(tk)={ ξ1(tk), ξ2(tk) ..., ξN(tk), tk=t0+k·dt (3)
Wherein, t0It is initial time, dt is time interval, k=0,1 ..., Jmax-1。
After one-dimensional range profile sequence is cached in buffer, you can to the one-dimensional range profile sequence cached in buffer
Pre-processed, to obtain the analysis window of one-dimensional range profile, and removed noise.
In the inventive solutions, it is possible to use various specific implementations carry out above-mentioned pretreatment.Hereinafter,
Technical scheme will be described in detail by taking a kind of specific implementation therein as an example.
For example, in a preferred embodiment of the invention, the one-dimensional range profile sequence to being cached in buffer is entered
Row pretreatment can include:
One-dimensional range profile course figure according to the one-dimensional range profile sequence cached in buffer is into gray level image;
The edge of the gray level image is extracted using image processing method, the starting of the analysis window of one-dimensional range profile is obtained
Position and final position, and it is removed the operation of noise.
Wherein, if the original position and final position of the analysis window of one-dimensional range profile are respectively NSAnd NE, then this it is one-dimensional away from
Analysis window one from picture has Nr=(NE-NS+ 1) individual range cell.
After the operation for being removed noise, the one-dimensional range profile sequence after pretreatment can be expressed as:
ξ′(k)={ ξ 'i (k), i=0,1 ..., Nr- 1, k=0,1 ..., Jmax-1 (4)
In the inventive solutions, can also further in the one-dimensional range profile sequence after above-mentioned pretreatment
Each one-dimensional range profile is normalized, the one-dimensional range profile ξ ' after can normalizingH (k) i:
Wherein, ξ 'i (k)Represent the scattering strength of each range cell in one-dimensional range profile;ξ′max (k)Represent that each is one-dimensional
The maximum of all range cell scattering strengths in Range Profile.
Step 102, cross correlation process is carried out by all one-dimensional range profile sequences within analysis window in buffer;Estimate
Meter target envelope amount of movement, and plus reference signal time delay correspondence position, obtain the actual range of target.
For example, in a preferred embodiment of the invention, it is possible to use formula as described below to analysis window it
Interior one-dimensional range profile sequence carries out cross correlation process:
Wherein, r(k)(nk) it is cross correlation process result, nkIt is cross correlation process variable, represents to ξ 'H (k)Carry out nkIt is individual away from
From unit movement, JmaxRepresent the size of buffer.
Assuming that two neighboring envelope (one-dimensional range profile after normalizing) ξ 'H (k-1)With ξ 'H (k)Between side-play amount be Δ
rk, then n is worked ask=Δ rkWhen, r(k)(nk) reach maximum.
In addition, in a preferred embodiment of the invention, it is possible to use formula as described below obtains the reality of target
Apart from pk:
pk=pk-1+cnk/ 2 Δ f, k=1,2 ..., Jmax-1 (7)
Wherein, p0It is the target range of initial time, c is the light velocity, and Δ f is radar swept bandwidth.
Step 103, the velocity estimation value of one-dimensional range profile is obtained by the method for a curve matching.
In the inventive solutions, it can be assumed that target range p is a curve on time independent variable t, can
To be expressed as:
P=VTt+p0+e (8)
Wherein, e is observation error, VTIt is the radial velocity of target, the variances sigma of e2=E [e2], if had k times to the target
Observation, k is the number of times of observation, and observation equation is respectively:
Write above-mentioned observation equation as vector form:P=[p1, p2... pk]T, e=[e1, e2..., ek]T, V=[VT, p0
]T,Then above-mentioned observation equation can be expressed as:
P=TV+e (10)
The cost function of the least-squares estimation of V is:
The estimate of VCost function should be made to reach minimum, that is, meet condition as described below:
Therefore, the estimate of V(i.e. the velocity estimation value of one-dimensional range profile) is:
Step 104, the velocity estimation value according to one-dimensional range profile, velocity compensation is carried out to one-dimensional range profile.
In the inventive solutions, it is possible to use various specific implementations carry out above-mentioned step 104.With
Under, technical scheme will be described in detail by taking a kind of specific implementation therein as an example.
For example, in a preferred embodiment of the invention, the velocity estimation value according to one-dimensional range profile, to one-dimensional
Range Profile carries out velocity compensation can be included:
Velocity estimation value according to one-dimensional range profile, translation compensation treatment is carried out to each one-dimensional range profile;
After being processed by translation compensation, the one-dimensional range profile of each adjacent echoes basic marshalling, therefore can be again
Make FFT to one-dimensional range profile after translation compensation treatment to estimate, obtain the phase difference ΔΦ of adjacent echoes;
Preset the first threshold value Z*, as ΔΦ > Z*When, return and perform step 102;And work as ΔΦ≤Z*When, judge
Velocity compensation terminates, and performs follow-up step 105.
Step 105, variance evaluation is carried out to carrying out each one-dimensional range profile after velocity compensation, it is determined that being used as velocity compensation
Umber of pulse meet error requirements.
In the inventive solutions, it is possible to use various specific implementations carry out above-mentioned step 104.With
Under, technical scheme will be described in detail by taking a kind of specific implementation therein as an example.
For example, in a preferred embodiment of the invention, the step 105 can include:.
Variance evaluation is carried out to carrying out each one-dimensional range profile after velocity compensation, the variance of each one-dimensional range profile is obtained
Estimate;
For example, variance evaluation can be carried out by formula as described below:
Wherein,Represent that i-th one-dimensional range profile is J in sizemaxBuffer in scattering strength it is flat
Average.
Then, then to the estimate of variance of all one-dimensional range profiles average, obtain mean of variance.For example, can be with
Mean of variance is obtained by formula as described below:
Finally, the second threshold value NN is pre-set*, whenWhen, then increase buffer size and return and perform step
Rapid 101;And work asWhen, it is determined that the umber of pulse as velocity compensation has met error requirements, terminates whole stream
Journey, so as to the group velocity for completing one-dimensional range profile is estimated and compensation process.
In the inventive solutions, the motion provided in the present invention can be checked by the method for simulation analysis
The technique effect of the speed compensation method of the one-dimensional range profile of target.
For example, in the inventive solutions, can be emulated in the case of known to radar parameter to point target
Analysis.
For Ka band radars, if centre frequency is 37.5GHz, target radial speed is in its maximum measuring speed scope
Maximum measuring speed VT=12km/s, actual radial velocity is generally less than the maximum measuring speed, target range radar
300km.When pulse width is Tp=50 μ s, Range Profile will be made to produce the translation of 36m, it is single equivalent to 480 high-resolution distances
Unit, after target velocity is estimated, Range Profile translation can be corrected for.
Fig. 2 is the speed compensation method effect diagram of the one-dimensional range profile of the moving target in the embodiment of the present invention.Such as
Shown in Fig. 2, in this simulation analysis, radar parameter as described below, the one-dimensional range profile shown in the left figure of Fig. 2 have been used
It is an one-dimensional range profile for scattering point when extraterrestrial target is 8000m/s with the relative velocity of radar.Shown in the right figure of Fig. 2
Be by the one-dimensional range profile after a velocity compensation.It can be seen from contrast according to two figures, original one-dimensional range profile is entered
After row velocity compensation, the translation of one-dimensional range profile and broadening are obtained for obvious correction, so as to effectively eliminate at a high speed
The adverse effect that motion brings.
In summary, in the speed compensation method of the one-dimensional range profile of moving target of the invention, due to radar being existed
The continuous one-dimensional range profile sequence for obtaining is cached in buffer in prefixed time interval, and to cached in buffer it is one-dimensional away from
Pre-processed from as sequence, obtained the analysis window of one-dimensional range profile, and removed noise;Then, by buffer it is all
One-dimensional range profile sequence within analysis window carries out cross correlation process;Estimate target envelope amount of movement, and plus reference signal
Time delay correspondence position, obtains the actual range of target;The speed for obtaining one-dimensional range profile by the method for a curve matching is estimated
Evaluation;Velocity estimation value according to one-dimensional range profile, velocity compensation is carried out to one-dimensional range profile;It is every after velocity compensation to carrying out
Individual one-dimensional range profile carries out variance evaluation, it is determined that the umber of pulse as velocity compensation meets error requirements, therefore can be to echo
Signal carries out velocity compensation, effectively eliminates the adverse effect that high-speed motion brings.
Presently preferred embodiments of the present invention is the foregoing is only, is not intended to limit the invention, it is all in essence of the invention
Within god and principle, any modification, equivalent substitution and improvements done etc. should be included within the scope of protection of the invention.
Claims (7)
1. a kind of speed compensation method of the one-dimensional range profile of moving target, it is characterised in that the method includes:
Step A, the one-dimensional range profile sequence for continuously obtaining radar within a preset time interval are cached in buffer, and to slow
The one-dimensional range profile sequence cached in storage is pre-processed, and obtains the analysis window of one-dimensional range profile, and removes noise;
Step B, all one-dimensional range profile sequences within analysis window in buffer are carried out into cross correlation process;Estimate target
Envelope amount of movement, and plus reference signal time delay correspondence position, obtain the actual range of target;
Step C, the velocity estimation value that one-dimensional range profile is obtained by the method for a curve matching;
Step D, the velocity estimation value according to one-dimensional range profile, velocity compensation is carried out to one-dimensional range profile, is specifically included:
Velocity estimation value according to one-dimensional range profile, translation compensation treatment is carried out to each one-dimensional range profile;
Make FFT to the one-dimensional range profile after translation compensation treatment to estimate, obtain the phase difference ΔΦ of adjacent echoes;
Preset the first threshold value Z*, work as ΔΦ>Z*When, return and perform step B;And work as ΔΦ≤Z*When, perform step E;
Step E, variance evaluation is carried out to carrying out each one-dimensional range profile after velocity compensation, it is determined that the pulse as velocity compensation
Number meets error requirements.
2. method according to claim 1, it is characterised in that the one-dimensional range profile sequence to being cached in buffer is entered
Row pretreatment includes:
One-dimensional range profile course figure according to the one-dimensional range profile sequence cached in buffer is into gray level image;
The edge of the gray level image is extracted using image processing method, the original position of the analysis window of one-dimensional range profile is obtained
And final position, and it is removed the operation of noise.
3. method according to claim 2, it is characterised in that using formula below to one-dimensional within analysis window
Range Profile Sequence carries out cross correlation process:
Wherein, r(k)(nk) it is cross correlation process result, nkIt is cross correlation process variable, it is right to representCarry out nkIndividual range cell
It is mobile, JmaxRepresent the size of buffer;Represent the one-dimensional range profile after normalization.
4. method according to claim 3, it is characterised in that the actual range p of target is obtained using formula belowk:
pk=pk-1+cnk/ 2 Δ f, k=1,2 ..., Jmax-1;
Wherein, p0It is the target range of initial time, c is the light velocity, and Δ f is radar swept bandwidth.
5. method according to claim 1, it is characterised in that the step E includes:
Variance evaluation is carried out to carrying out each one-dimensional range profile after velocity compensation, the variance evaluation of each one-dimensional range profile is obtained
Value;
Estimate of variance to all one-dimensional range profiles is averaged, and obtains mean of variance;
Pre-set the second threshold value NN*, whenWhen, increase buffer size and return to execution step A;WhenWhen, terminate whole flow process.
6. method according to claim 5, it is characterised in that variance evaluation is carried out by formula below:
Wherein,Represent that i-th one-dimensional range profile is J in sizemaxBuffer in scattering strength average value.
7. method according to claim 6, it is characterised in that mean of variance is obtained by formula below:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510434889.9A CN105068072B (en) | 2015-07-22 | 2015-07-22 | A kind of speed compensation method of the one-dimensional range profile of moving target |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510434889.9A CN105068072B (en) | 2015-07-22 | 2015-07-22 | A kind of speed compensation method of the one-dimensional range profile of moving target |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105068072A CN105068072A (en) | 2015-11-18 |
CN105068072B true CN105068072B (en) | 2017-07-07 |
Family
ID=54497487
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510434889.9A Active CN105068072B (en) | 2015-07-22 | 2015-07-22 | A kind of speed compensation method of the one-dimensional range profile of moving target |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105068072B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106872974B (en) * | 2017-01-23 | 2019-07-09 | 西安电子科技大学 | High-precision motion target imaging method based on hypersonic platform Two-channels radar |
CN109444842B (en) * | 2019-01-04 | 2020-09-25 | 北京环境特性研究所 | Target electromagnetic scattering characteristic data reconstruction method and device |
CN109741278B (en) * | 2019-01-04 | 2020-09-11 | 北京环境特性研究所 | Image denoising method |
CN112630778B (en) * | 2020-11-16 | 2022-08-16 | 南京理工大学 | Ultrahigh sound speed target imaging method |
US20220373674A1 (en) * | 2021-05-21 | 2022-11-24 | Kabushiki Kaisha Toshiba | Method and device for improved range tracking for inverse synthetic aperture radar |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09211121A (en) * | 1996-02-02 | 1997-08-15 | Osaka Gas Co Ltd | Searching method and device |
JPH09243743A (en) * | 1996-03-06 | 1997-09-19 | Osaka Gas Co Ltd | Searching method and device |
CN103529444A (en) * | 2013-09-27 | 2014-01-22 | 安徽师范大学 | Vehicle-mounted millimeter-wave radar moving target recognizer and recognition method |
CN103558596A (en) * | 2013-11-14 | 2014-02-05 | 上海电机学院 | Stepped frequency radar fuze velocity compensation method |
-
2015
- 2015-07-22 CN CN201510434889.9A patent/CN105068072B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09211121A (en) * | 1996-02-02 | 1997-08-15 | Osaka Gas Co Ltd | Searching method and device |
JPH09243743A (en) * | 1996-03-06 | 1997-09-19 | Osaka Gas Co Ltd | Searching method and device |
CN103529444A (en) * | 2013-09-27 | 2014-01-22 | 安徽师范大学 | Vehicle-mounted millimeter-wave radar moving target recognizer and recognition method |
CN103558596A (en) * | 2013-11-14 | 2014-02-05 | 上海电机学院 | Stepped frequency radar fuze velocity compensation method |
Non-Patent Citations (1)
Title |
---|
《频率步进雷达高速目标成像中的运动补偿方法研究》;左衍琴;《中国优秀硕士学位论文全文数据库 信息科技辑》;20111215;第22-32页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105068072A (en) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105068072B (en) | A kind of speed compensation method of the one-dimensional range profile of moving target | |
CN107966688B (en) | Broadband radar target speed ambiguity resolving method based on phase interference technology | |
CN103760558B (en) | Terahertz radar ISAR imaging method | |
CN108051809B (en) | Moving target imaging method and device based on Radon transformation and electronic equipment | |
CN105807267B (en) | A kind of MIMO radar extends mesh object detection method | |
CN101509972B (en) | Wideband radar detecting method for correcting correlation matrix based on high resolution target distance image | |
CN107271993B (en) | Scanning radar angle super-resolution imaging method based on maximum posterior | |
CN103163523B (en) | Low level wind shear velocity estimation method based on compressed sensing | |
CN107132534B (en) | Optimization method for high-speed radar target frequency domain detection | |
CN101738606A (en) | Method for detecting coherent integration of radar target based on generalized Doppler filter bank | |
CN110954885B (en) | Adaptive target reconstruction method for frequency agile radar based on SBL | |
CN103744068A (en) | Moving target detection imaging method of dual-channel frequency modulation continuous wave SAR system | |
CN104730496B (en) | A kind of Squint SAR ground moving target extracting method | |
CN103777178B (en) | A kind of synchronous error compensation method, equipment and system | |
CN104950305A (en) | Real beam scanning radar angle super-resolution imaging method based on sparse constraint | |
CN104502912A (en) | Imaging method for inverse synthetic aperture radar of high-speed moving targets | |
CN104849708A (en) | High-speed maneuvering target parameter estimation method based on frequency domain polynomial phase transformation | |
Li et al. | ISAR imaging of maneuvering target with complex motions based on ACCF–LVD | |
CN110568434A (en) | multi-channel uniform acceleration trajectory SAR moving target two-dimensional speed estimation method | |
CN109031299B (en) | ISAR (inverse synthetic aperture radar) translation compensation method based on phase difference under low signal-to-noise ratio condition | |
CN106054154B (en) | A kind of step frequency signal high-resolution imaging method of the maneuvering target based on ICPF | |
Popović et al. | Autofocusing of SAR images based on parameters estimated from the PHAF | |
CN104155653B (en) | SAR back projection imaging method based on feature distance subspace | |
Ma et al. | Fast sidelobe suppression based on two-dimensional joint iterative adaptive filtering | |
CN106814352A (en) | A kind of multi-target detection method based on Golay complementary waveform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |