CN108983189A - A kind of two-dimensional micromotion track estimation method of Vibration Targets - Google Patents
A kind of two-dimensional micromotion track estimation method of Vibration Targets Download PDFInfo
- Publication number
- CN108983189A CN108983189A CN201810796582.7A CN201810796582A CN108983189A CN 108983189 A CN108983189 A CN 108983189A CN 201810796582 A CN201810796582 A CN 201810796582A CN 108983189 A CN108983189 A CN 108983189A
- Authority
- CN
- China
- Prior art keywords
- signal
- time
- antenna
- echo
- receiving antenna
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Abstract
The invention discloses a kind of two-dimensional micromotion track estimation methods of Vibration Targets, which comprises generates radio frequency radar signal, is emitted by transmitting antenna；The radar return of multiple targets is received using three receiving antennas, one of receiving antenna is transmitting antenna；Obtain the echo-signal that each antenna receives；To ICCD processing is carried out after the segmentation of each echo-signal, by treated, signal is spliced to obtain the time-domain signal with the same length of original echoed signals, thus gets the phase of signal；The interferometric phase between reception echo treated the signal of transmitting antenna and reception echo treated the signal of another two antenna is calculated, the Vibration Targets of acquisition are two-dimensional micromotion tracks.The signal wound in time-frequency figure is carried out segment processing by method of the invention, is reapplied ICCD algorithm and is decomposed to vibration signal；This method uses three receiving antennas, then extracts the interferometric phase between each antenna echo signal, can obtain the two-dimensional micromotion track of Vibration Targets.
Description
Technical field
The present invention relates to Radar Signal Processing and micro-Doppler effect fields, are a kind of two-dimensional micromotion rails of Vibration Targets
Mark estimation method.
Background technique
Vibration is a kind of fine motion form, and fine motion causes many concerns in the nearest more than ten years.Doppler caused by fine motion
Effect is referred to as micro-Doppler effect (document [1]: V.C.Chen, " Analysis of radar micro-Doppler
signature with time–frequency transform,”in Proc.IEEE Statistical Signal
Array Process.,2000,pp.463–466.；Document [2]: V.C.Chen, F.Li, S.-S.Ho et al.,
“Analysis of micro-Doppler signatures,”Proc.Inst.Electr.Eng.—Radar Sonar
Navig., vol.150, no.4, pp.271-276, Aug.2003. document [3]: T.Sparr and B.Krane, " Micro-
Doppler analysis of vibrating targets in SAR,”Proc.Inst.Electr.Eng.—Radar
Sonar Navig., vol.150, no.4, pp.277-283, Aug.2003. document [4]: V.C.Chen, F.Li, S.-S.Ho
et al.,“Micro-Doppler effect in radar:Phenomenon,model,and simulation study,”
IEEE Trans.Aerosp.Electron.Syst.,vol.42,no.1,pp.2–21,Jan.2006.).Different target tools
There is respective unique fine motion feature, these fine motion features facilitate the detection and identification of radar target.Currently, the overwhelming majority is ground
Studying carefully all is the individual antenna radar system used, and which can only obtain fine motion information of the target along radar line of sight direction, cannot
Obtain the fine motion track of target.The case where traditional method can not handle fine motion information there are aliasing effects.
Recent years, several methods for extracting fine motion parameter are suggested, including more base station radars, interferometer radar and multiple antennas
Technology.Air force engineering university Zhang Qun et al. by using displaced phase center antenna technology (document [5]: W.Zhang,
C.Tong,Q.Zhang,et al,“Extraction of vibrating features with dual-channel
fixed-receiver bistatic SAR,”IEEE Geoscience&Remote Sensing Letters.,vol.9,
No.3, pp.507-511, Nov.2012), clutter reduction only retains the echo component of Vibration Targets, and extracts fine motion parameter.Guest
Sunset Fa Niya electronic engineering of state university Dustin P.Fairchild et al. extracts vibration using more base station micro-doppler radars
Direction (document [6]: D.P.Fairchild, R.M.Narayanan, " Multistatic of the moving-target with respect to radar line of sight
micro-doppler radar for determining target orientation and activity
classification,”IEEE Transactions on Aerospace&Electronic Systems.,vol.52,
No.1, pp.512-521, Feb.2016), it can also so extract the opposite additional fine motion information of monostatic radar.The Chinese Academy of Sciences
Zhang Yun China of microwave remote sensing technique key lab et al. using Ka wave band interferometer radar carried out train experiment (document [7]:
W.Zhai,Y.Zhang,Q.Yang,X.Shi,“Micro-motion of a moving train observed by a Ka-
band interferometric radar,”Electronics Letters.,vol.52,no.12,pp.1065-1067,
Jun.2016 it), and using interferometric phase is analyzed, it can be observed that the fluctuation of railway carriage, this is also real for the first time in the world
Verifying can analyze fine motion target using interferometric phase.Air force engineering university sieve, which is met et al., obtains sky using multi-antenna technology
Between target two-dimensional imaging (document [8]: Y.Luo, Y.Chen, Y.Sun, Q.Zhang, " Narrowband radar imaging
and scaling for space targets,”IEEE Geoscience&Remote Sensing Letters.,
Vol.14, no.6, pp.946-950, Jun.2017), this method can extract partial target fine motion track, since time-frequency figure exists
Intersect pixel, this method cannot obtain the fine motion track of cross section.
Shanghai Communications University's mechanical system and vibrate National Key Laboratory it is old be skewer et al. propose ICCD algorithm decomposition
Micro-doppler signal (document [9]: S.Chen, X.Dong, G.Xing, Z.Peng, W.Zhang, G.Meng, " Separation
of overlapped non-stationary signals by ridge path regrouping and intrinsic
chirp component decomposition,”IEEE Sensors Journal.,vol.17,no.18,pp.5994-
6005, Sep.2017), ICCD algorithm can decomposed signal time-frequency figure well cross section.
Wrapping phenomena can occur for radar signal, therefore existing ICCD algorithm can not decompose vibration signal, therefore,
Current target fine motion track estimation method can not directly use ICCD algorithm.
Summary of the invention
It is an object of the invention to overcome above-mentioned technological deficiency, provided newly to obtain the two-dimensional micromotion track of Vibration Targets
Technical solution.This method is based on triantennary interferometer radar system, receives echoes using three receiving antennas, and to there is aliasing
The signal of phenomenon carries out segment processing and then splices signal to each segment signal application ICCD algorithm, extracts transmitting antenna
With the interferometric phase between the echo-signal after another two antenna processing, the two-dimensional micromotion track of Vibration Targets is finally obtained.It is imitative
True result demonstrates the validity of this method well.
To achieve the goals above, the present invention proposes a kind of two-dimensional micromotion track estimation method of Vibration Targets, the side
Method includes:
Radio frequency radar signal is generated, is emitted by transmitting antenna；
The radar return of multiple targets is received using three receiving antennas, one of receiving antenna is transmitting antenna；It obtains
The echo-signal for taking each antenna to receive；
To carrying out ICCD decomposition after the segmentation of each echo-signal, the signal after decomposition is spliced to obtain and original echo
The time-domain signal of the same length of signal, thus gets the phase of signal；
Calculate reception echo treated the signal of transmitting antenna and reception echo treated the signal of another two antenna
Between interferometric phase, the Vibration Targets of acquisition are two-dimensional micromotion tracks.
As a kind of improvement of the above method, the line form right angle of the transmitting antenna and other two receiving antennas.
As a kind of improvement of the above method, the method is specifically included:
Step 1) generates radio frequency radar signal, is emitted by transmitting antenna；
The mathematical model of the radar signal is expressed as follows:
T (t)=exp (j2 π fct) (1)
Wherein, t indicates time, fcIndicate the frequency of baseband signal；
T (t) signal is sampled in fixed time period, obtains N number of discrete baseband transmission signal sequence Tn(t),n
=1,2 ... N；
Step 2) emits signal sequence Tn(t), n=1,2 ... N return to receiver, receiving antenna i after being irradiated in target
The signal received is Rn i(t), n=1,2 ... N, i=1,2,3；The backscattering coefficient of k-th of point target is σk, k-th point
The distance between target and receiving antenna i are1≤k≤K, K are the sum of point target, and c is the light velocity, then receiving antenna i is received
The signal arrivedAre as follows:
Step 3) is by three echo-signal Rn i(t) it transforms to time-frequency domain, chooses the data of time-frequency matrix the first row, according to the
The peak position of data line determines the position of signal subsection, and time-domain signal is segmented；ICCD is set according to number of targets to calculate
The parameter of method decomposed signal ingredient；Each section of time-domain signal is decomposed using ICCD algorithm, and by the signal after decomposition into
Row splicing, obtains the time-domain signal with the same length of original echoed signals；Thus the phase of signal is got；
For step 4) using transmitting antenna as receiving antenna 1, other two antennas are receiving antenna 2 and receiving antenna 3, are counted respectively
The interferometric phase between receiving antenna 1 and receiving antenna 2 and receiving antenna 1 and receiving antenna 3 is calculated, is obtained using interferometric phase
The two-dimensional micromotion track of Vibration Targets.
As a kind of improvement of the above method, the step 3) is specifically included:
Step 3-1) by the step 2) obtain discrete reception signal sequence Rn i(t) low-pass filtering treatment is carried out, is passed through
Time-domain signal r after low-pass filtering treatmentn i(t), n=1,2 ... N are as follows:
Step 3-2) to the time-domain signal r after the low-pass filtering of acquisitionn i(t) Short Time Fourier Transform is carried out, time-frequency is obtained
Matrix and time-frequency figure choose the data of time-frequency matrix the first row, determine signal subsection according to the peak position of the first row data
Time-domain signal is segmented by position, so that each segment signal is all without wrapping phenomena, time-domain signal is divided into M sections；
Step 3-3) it chooses and decomposes composition parameter, apply ICCD to calculate respectively the time-domain signal after the step 3-2) segmentation
Method is decomposed, and the time-domain signal after each section of decomposition is stitched together, and spliced time-domain signal has set decomposition
The number of composition parameter, the length of each time-domain signal are N, then the echo that antenna i is received through decomposition after k-th
SubsignalAre as follows:
Wherein,It isAmplitude,It isPhase.
As a kind of improvement of the above method, the step 4) is specifically included:
Step 4-1) using transmitting antenna as receiving antenna 1, the kth of receiving antenna 1 and receiving antenna 2 treated echo
Interferometric phase between a subsignalAre as follows:
Step 4-2) calculate receiving antenna 1 and receiving antenna 3 treated interference phase between k-th of subsignal of echo
PositionAre as follows:
Step 4-3) k-th Scattering Targets point P two-dimensional micromotion track are as follows:
Wherein, τ indicate the time, λ indicate transmitting electromagnetic wave wavelength, d be between transmitting antenna and receiving antenna 1 away from
From.
Present invention has an advantage that
1, the signal wound in time-frequency figure is carried out segment processing by method of the invention, reapplies ICCD algorithm to letter
It number is decomposed, can realize the decomposition of vibration signal well；
2, method of the invention uses three receiving antennas, then extracts the interferometric phase between each antenna echo signal, can
The two-dimensional micromotion track of Vibration Targets is obtained, and traditional single-shot list receipts antenna system can only obtain target along radar line of sight direction
Parameter；
3, method of the invention can decompose interferometric phase when time-frequency figure has intersection by applying ICCD algorithm.
Detailed description of the invention
Fig. 1 is the two-dimensional micromotion track estimation method flow chart of Vibration Targets of the invention；
Fig. 2 is Vibration Targets illustraton of model of the invention；
Fig. 3 is the time-frequency figure before two oscillation point target simulator radar echo signals decompose；
Fig. 4 is the interferometric phase before two oscillation point target simulator radar echo signals decompose；
Fig. 5 is the time-frequency figure for emulating first oscillation point target signal elements after radar echo signal decomposes；
Fig. 6 is the time-frequency figure for emulating second oscillation point target signal elements after radar echo signal decomposes；
Fig. 7 is first vibration point target signal after decomposition in antenna TRX1And RX2Between interferometric phase；
Fig. 8 is second vibration point target signal after decomposition in antenna TRX1And RX2Between interferometric phase；
Fig. 9 is first vibration point target signal after decomposition in antenna TRX1And RX3Between interferometric phase；
Figure 10 is second vibration point target signal after decomposition in antenna TRX1And RX3Between interferometric phase；
Figure 11 is the two-dimensional micromotion track of first vibration point target of reconstruct；
Figure 12 is the two-dimensional micromotion track of second vibration point target of reconstruct.
Specific embodiment
Technical solution for a better understanding of the present invention is made embodiments of the present invention below in conjunction with attached drawing further
Description.
The present invention proposes a kind of based on inherently linear FM signal ingredient breakdown (the Segmental Intrinsic of segmentation
Chirp Component Decomposition, hereinafter referred to as SICCD) algorithm estimation Vibration Targets two-dimensional micromotion track side
Method (vibration parameters refer mainly to Oscillation Amplitude and period here).Different fine motion targets usually has different fine motion features, leads to
The fine motion feature for extracting target is crossed, target preferably can be identified and be classified.Traditional single receiving antenna radar system
Target can only be obtained along the information in radar line of sight direction, target two dimension or three-dimensional fine motion information can not be obtained, in transmitting pulse
When repetition rate is lower, wrapping phenomena but will be generated.Method proposed by the present invention is based on single-shot and receives antenna radar system more, passes through
The two-dimensional micromotion feature of Vibration Targets is effectively extracted to the processing of interferometric phase.This method passes through transmitting antenna first and emits radar
Signal, and using multiple receiving antennas receive echo, echo-signal is then transformed into time-frequency domain, by time-frequency domain signal by when
Domain signal carries out segment processing, then uses ICCD (Intrinsic Chirp Component Decomposition) algorithm
Each section of time-domain signal is decomposed, is finally spliced each section of time-domain signal, and extract interferometric phase, by interfering phase
Position obtains the two-dimensional micromotion track of Vibration Targets.
As shown in Figure 1, the invention proposes a kind of two-dimensional micromotion track estimation method of Vibration Targets, this method is specifically wrapped
It includes:
Step 1): generating radio frequency radar signal, and the mathematical model of the radar signal is expressed as follows:
T (t)=exp (j2 π fct) (1)
Wherein, t indicates time, fcThe frequency for indicating baseband signal, samples T (t) signal in fixed time period,
Obtain N number of discrete baseband transmission signal sequence Tn(t), n=1,2 ... N；
Step 2): transmitting signal Tn(t) it after transmitting, is irradiated in target, and returns to receiver, it is assumed that antenna i connects
The signal received isI=1,2,3；Assuming that sharing K point target, the backscattering coefficient of k-th of point target is σk, kth
The distance between a point target and antenna i areThe light velocity is c, then the signal that antenna i is receivedMathematical model are as follows:
Step 3): the discrete reception signal sequence that the step 2) is obtainedLow-pass filtering treatment is carried out, is passed through
Signal after low-pass filtering treatment are as follows:
Step 4): the time-domain signal after the low-pass filtering that the step 3) is obtained carries out Short Time Fourier Transform, obtains
Time-frequency figure chooses the data of time-frequency matrix the first row, the position of signal subsection is determined according to the peak position of the first row data, will
Time-domain signal is segmented, so that each segment signal is all without wrapping phenomena, it is assumed that time-domain signal is divided into M sections, then respectively
Segment signal length sum should be N；
Segment processing in the step 4) can choose the first row data of time-frequency figure, according to the first row data
Peak value determines the position of segmentation breakpoint.
Step 5): suitable decomposition composition parameter is chosen, the time-domain signal after the step 4) segmentation is applied respectively
ICCD algorithm is decomposed, and the time-domain signal after each section of decomposition is stitched together, and spliced time-domain signal should have set
The number for the decomposition composition parameter set, and the length of each time-domain signal is N, it is assumed that the echo that antenna i is received passes through
K-th of subsignal after decomposition beIts mathematic(al) representation are as follows:
Wherein,It isAmplitude,It isPhase；
Step 6): by the step 5) obtain decomposition after time-domain signal carry out interference processing, and extract it is each at point it
Between interferometric phase, using interferometric phase obtain Vibration Targets two-dimensional micromotion track, and then extract Vibration Targets two dimension throw
Fine motion parameter in shadow plane, what is extracted here is the interferometric phase between transmitting antenna and another two antenna, and interferometric phase mentions
Take process are as follows:
Wherein, conj expression takes conjugate operation, and Pha expression takes phase operation,Indicate antenna i and l treated return
Interferometric phase between k-th of subsignal of wave.
Fig. 2 is Vibration Targets illustraton of model.The model is using scattering point target.(X, Y, Z) is radar fix system, and O is coordinate
It is origin.The model uses three receiving antenna TRX1、RX2、RX3, antenna TRX1Emit radar signal, antenna TRX1、RX2、RX3Together
When receive radar return, antenna TRX1Positioned at coordinate origin O, antenna RX2And RX3Coordinate be (d, 0,0) and (0,0, d) respectively.
Vibrational scattering point P is (x in the coordinate of radar fix system0,y0,z0), the centre of oscillation is O ', O ' be also reference frame (X ', Y ',
Z ') origin, be R at a distance from O and O '1.The angular speed of vibrational scattering point, amplitude and initial phase are used respectivelyA and θ0It indicates.
Point P is α and β at the azimuth of coordinate system (X ', Y ', Z ') direction of vibration and pitch angle.Coordinate of the scattering point P at the τ moment are as follows:
Assuming that scattering point P to antenna TRX1、RX2、RX3Distance use R respectivelyr1(τ)、Rr2(τ)、Rr3(τ) is indicated, then
Antenna TRX1、RX2、RX3The base band echo received are as follows:
In formula (7), σ1、σ2、σ3It is scattering point P and antenna TR respectivelyX1、RX2、RX3Between echo-signal scattering system
Number, λ is carrier wavelength, then antenna TRX1With antenna RX2、RX3Receive the interferometric phase between echo are as follows:
In formula (8),Indicate antenna TRX1With antenna RX2The interferometric phase between echo is received,Indicate antenna
TRX1With antenna RX3Receive the interferometric phase between echo, R1It is antenna TRX1At a distance from the P point centre of oscillation.Assuming that R1It is known
, then the coordinate of scattering point P can be estimated by following formula:
The two-dimension vibration track of Vibration Targets can be reconstructed by formula (9), and then extracts Vibration Targets parameter.
Fig. 3 is the time-frequency figure for emulating the signal before radar signal is decomposed.In the simulation, two vibrational scattering points are false
If their common centres of oscillation are the origin O ' of coordinate system (X ', Y ', Z ').The frequency of first vibrational scattering point, amplitude and
Initial phase is 1Hz, 20m and-π/4rad.The frequency of second vibrational scattering point, amplitude and initial phase are 3Hz, 1.2m and 3
π/4rad.The azimuth angle alpha and pitch angle β of direction of vibration are π/6rad and π/3rad respectively.Antenna TRX1And RX2, RX3Between away from
It is arranged to 1m from d.O ' is (0,50000,0) in the coordinate of coordinate system (X, Y, Z).When pulse recurrence frequency, carrier frequency and accumulation
Between be 3000,10GHz and 1s.The white Gaussian noise of 15dB is added in emulation.From figure 3, it can be seen that emulation signal when
There is wrapping phenomena in frequency figure, and the parameter of target can not be extracted with traditional algorithm.Fig. 4 is before emulation radar signal is decomposed
Interferometric phase, it can be seen that the interferometric phase before decomposition cannot reflect any information.Fig. 5 is the after emulation radar signal is decomposed
The time-frequency figure of one signal, Fig. 6 are the time-frequency figure for emulating second signal after radar signal is decomposed, can from the two figures
Out, SICCD algorithm is so fine that be decomposed to echo-signal, and the cross section of time-frequency figure is also decomposed very well.Fig. 7 is point
First signal is in antenna TR after solutionX1And RX2Between interferometric phase, Fig. 8 be decompose after second signal in antenna TRX1And RX2
Between interferometric phase, Fig. 9 be decompose after first signal in antenna TRX1And RX3Between interferometric phase, Figure 10 be decompose after
Second signal is in antenna TRX1And RX3Between interferometric phase, from this four figure it can be seen that two simulated scatter points interference
Phase is extracted well, by carrying out sine curve fitting to interferometric phase, it is estimated that the frequency of first scattering point
It is 1Hz, the frequency of second scattering point is 3Hz.Figure 11 is first scattering point of reconstruct in the fine motion track of x-z-plane, figure
12 be second scattering point of reconstruct in the fine motion track of x-z-plane.According to the two-dimensional micromotion track reconstructed, it is estimated that
Angle of the direction of vibration of two vibrational scattering points between the projection and x-axis of x-z-plane is 1.11rad, two estimated
The amplitude of scattering point is 19.3m and 1.2m in the projection of x-z-plane.
It should be noted last that the above examples are only used to illustrate the technical scheme of the present invention and are not limiting.Although ginseng
It is described the invention in detail according to embodiment, those skilled in the art should understand that, to technical side of the invention
Case is modified or replaced equivalently, and without departure from the spirit and scope of technical solution of the present invention, should all be covered in the present invention
Scope of the claims in.
Claims (5)
1. a kind of two-dimensional micromotion track estimation method of Vibration Targets, which comprises
Radio frequency radar signal is generated, is emitted by transmitting antenna；
The radar return of multiple targets is received using three receiving antennas, one of receiving antenna is transmitting antenna；It obtains every
The echo-signal that a antenna receives；
To carrying out ICCD decomposition after the segmentation of each echo-signal, the signal after decomposition is spliced to obtain and original echoed signals
The time-domain signal of same length, thus gets the phase of signal；
It calculates between reception echo treated the signal of transmitting antenna and reception echo treated the signal of another two antenna
Interferometric phase, the Vibration Targets of acquisition are two-dimensional micromotion tracks.
2. the two-dimensional micromotion track estimation method of Vibration Targets according to claim 1, which is characterized in that the transmitting day
The line form right angle of line and other two receiving antennas.
3. the two-dimensional micromotion track estimation method of Vibration Targets according to claim 2, which is characterized in that the method tool
Body includes:
Step 1) generates radio frequency radar signal, is emitted by transmitting antenna；
The mathematical model of the radar signal is expressed as follows:
T (t)=exp (j2 π fct) (1)
Wherein, t indicates time, fcIndicate the frequency of baseband signal；
T (t) signal is sampled in fixed time period, obtains N number of discrete baseband transmission signal sequence Tn(t), n=1,
2…N；
Step 2) emits signal sequence Tn(t), n=1,2 ... N return to receiver after being irradiated in target, receiving antenna i is received
The signal arrived is Rn i(t), n=1,2 ... N, i=1,2,3；The backscattering coefficient of k-th of point target is σk, k-th of point target
The distance between receiving antenna i isK is the sum of point target, and c is the light velocity, then receiving antenna i is received
SignalAre as follows:
Step 3) is by three echo-signal Rn i(t) time-frequency domain is transformed to, the data of time-frequency matrix the first row are chosen, according to the first row
The peak position of data determines the position of signal subsection, and time-domain signal is segmented；ICCD algorithm point is set according to number of targets
Solve the parameter of signal component；Each section of time-domain signal is decomposed using ICCD algorithm, and the signal after decomposition is spelled
It connects, obtains the time-domain signal with the same length of original echoed signals；Thus the phase of signal is got；
For step 4) using transmitting antenna as receiving antenna 1, other two antennas are receiving antenna 2 and receiving antenna 3, calculate separately and connect
The interferometric phase between antenna 1 and receiving antenna 2 and receiving antenna 1 and receiving antenna 3 is received, is obtained and is vibrated using interferometric phase
The two-dimensional micromotion track of target.
4. the two-dimensional micromotion track estimation method of Vibration Targets according to claim 3, which is characterized in that the step 3)
It specifically includes:
Step 3-1) by the step 2) obtain discrete reception signal sequence Rn i(t) low-pass filtering treatment is carried out, by low pass
Time-domain signal r after filtering processingn i(t), n=1,2 ... N are as follows:
Step 3-2) to the time-domain signal r after the low-pass filtering of acquisitionn i(t) Short Time Fourier Transform is carried out, time-frequency matrix is obtained
With time-frequency figures, the data of time-frequency matrix the first row are chosen, the position of signal subsection is determined according to the peak position of the first row data,
Time-domain signal is segmented, so that each segment signal is all without wrapping phenomena, time-domain signal is divided into M sections；
Step 3-3) choose decompose composition parameter, by the step 3-2) segmentation after time-domain signal apply respectively ICCD algorithm into
Row decomposes, and the time-domain signal after each section of decomposition is stitched together, and spliced time-domain signal has set decomposition ingredient
The number of parameter, the length of each time-domain signal are N, then k-th son letter of the echo that antenna i is received after decomposing
NumberAre as follows:
Wherein,It isAmplitude,It isPhase.
5. Vibration Targets parameter extracting method according to claim 4, which is characterized in that the step 4) specifically includes:
Step 4-1) using transmitting antenna as receiving antenna 1, k-th of son of receiving antenna 1 and receiving antenna 2 treated echo
Interferometric phase between signalAre as follows:
Step 4-2) calculate receiving antenna 1 and receiving antenna 3 treated interferometric phase between k-th of subsignal of echoAre as follows:
Step 4-3) k-th Scattering Targets point P two-dimensional micromotion track are as follows:
Wherein, τ indicates the time, and λ indicates the wavelength of transmitting electromagnetic wave, and d is the distance between transmitting antenna and receiving antenna 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810796582.7A CN108983189B (en) | 2018-07-19 | 2018-07-19 | Two-dimensional micro-motion track estimation method for vibration target |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810796582.7A CN108983189B (en) | 2018-07-19 | 2018-07-19 | Two-dimensional micro-motion track estimation method for vibration target |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108983189A true CN108983189A (en) | 2018-12-11 |
| CN108983189B CN108983189B (en) | 2020-07-03 |
Family
ID=64548841
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810796582.7A Active CN108983189B (en) | 2018-07-19 | 2018-07-19 | Two-dimensional micro-motion track estimation method for vibration target |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108983189B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114646279A (en) * | 2021-12-22 | 2022-06-21 | 上海交通大学 | Microwave multi-scale vibration and deformation measurement method and system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106569194A (en) * | 2016-10-28 | 2017-04-19 | 中国人民解放军空军工程大学 | Interferometric three-dimensional imaging and micro-motion feature extraction method of broadband radar spatial conical target |
| CN106990398A (en) * | 2016-01-21 | 2017-07-28 | 中国人民解放军空军工程大学 | A kind of body of revolution fine motion feature awareness extracting method |
| CN107085212A (en) * | 2017-04-01 | 2017-08-22 | 中国人民解放军空军工程大学 | A kind of spin target time-varying three-D imaging method based on linearly modulated stepped frequency |
| CN107132536A (en) * | 2017-04-10 | 2017-09-05 | 中国科学院国家空间科学中心 | It is a kind of to eliminate the method that target fine motion is disturbed radar imagery |
| CN107271982A (en) * | 2017-07-31 | 2017-10-20 | 广东工业大学 | A kind of object detection method and device based on millimeter wave |
-
2018
- 2018-07-19 CN CN201810796582.7A patent/CN108983189B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106990398A (en) * | 2016-01-21 | 2017-07-28 | 中国人民解放军空军工程大学 | A kind of body of revolution fine motion feature awareness extracting method |
| CN106569194A (en) * | 2016-10-28 | 2017-04-19 | 中国人民解放军空军工程大学 | Interferometric three-dimensional imaging and micro-motion feature extraction method of broadband radar spatial conical target |
| CN107085212A (en) * | 2017-04-01 | 2017-08-22 | 中国人民解放军空军工程大学 | A kind of spin target time-varying three-D imaging method based on linearly modulated stepped frequency |
| CN107132536A (en) * | 2017-04-10 | 2017-09-05 | 中国科学院国家空间科学中心 | It is a kind of to eliminate the method that target fine motion is disturbed radar imagery |
| CN107271982A (en) * | 2017-07-31 | 2017-10-20 | 广东工业大学 | A kind of object detection method and device based on millimeter wave |
Non-Patent Citations (1)
| Title |
|---|
| SHIQIAN CHEN等: "Nonlinear Chirp Mode Decomposition: A Variational Method", 《IEEE SENSORS JOURNAL》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114646279A (en) * | 2021-12-22 | 2022-06-21 | 上海交通大学 | Microwave multi-scale vibration and deformation measurement method and system |
| CN114646279B (en) * | 2021-12-22 | 2023-03-14 | 上海交通大学 | Microwave multi-scale vibration and deformation measurement method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108983189B (en) | 2020-07-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Imaging of a moving target with rotating parts based on the Hough transform | |
| Malanowski et al. | Detection of moving targets with continuous-wave noise radar: theory and measurements | |
| CN108459308B (en) | Echo simulation method and device based on time-varying RCS data | |
| Fan et al. | Moving target localization behind L-shaped corner with a UWB radar | |
| CN204989459U (en) | Angle deception jamming device to pulse homing guidance radar | |
| Higgins et al. | Passive multistatic radar experiment using WiMAX signals of opportunity. Part 1: Signal processing | |
| CN106680791A (en) | Radar echo simulating method based on broadband sweep frequency data | |
| Garry et al. | Passive ISAR part II: narrowband imaging | |
| CN108983189A (en) | A kind of two-dimensional micromotion track estimation method of Vibration Targets | |
| Jingcheng et al. | UAV detection and identification in the Internet of Things | |
| Abratkiewicz et al. | Time-Frequency Reassigned Micro-Doppler Signature Analysis Using the XY-DemoRad System | |
| Abratkiewicz et al. | Chirp rate estimation and micro-Doppler signatures for pedestrian security radar systems | |
| Webster et al. | Passive multistatic radar experiment using WiMAX signals of opportunity. Part 2: Multistatic velocity backprojection | |
| CN108983190B (en) | Method for acquiring different scattering center micro-motion tracks of rotating target based on interference phase | |
| CN106569195B (en) | A kind of spin fine motion raid cluster resolution method based on the slow time picture of distance- | |
| Liaquat et al. | An End-to-End Modular Framework for Radar Signal Processing: A simulation-based tutorial | |
| Tai et al. | An improved multiplication modulation deception jamming method for countering ISAR | |
| Kang et al. | Micro-Doppler feature extraction for interferometric radar based on Viterbi algorithm and intrinsic chirp component decomposition | |
| Schoenecker et al. | Maximum likelihood probabilistic data association tracker applied to bistatic sonar data sets | |
| Tai et al. | A two-dimensional frequency modulation deception jamming against ISAR | |
| CN111337897A (en) | LFMCW radar rapid target identification method | |
| Burgos-Garcia et al. | Radar signature of a helicopter illuminated by a long LFM signal | |
| Jia et al. | Analysis and simulation of multi-target echo signals from a phased array radar | |
| Dong et al. | Anti-intermittent Sampling Jamming Method Based on Frequency Agile Radar and K-means | |
| Zuo et al. | Detection of UAV target based on Continuous Radon transform and Matched filtering process for Passive Bistatic Radar |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |