CN106353744A - Multi-parameter combined estimation method based on bi-static FDA-MIMO radars - Google Patents
Multi-parameter combined estimation method based on bi-static FDA-MIMO radars Download PDFInfo
- Publication number
- CN106353744A CN106353744A CN201610962226.9A CN201610962226A CN106353744A CN 106353744 A CN106353744 A CN 106353744A CN 201610962226 A CN201610962226 A CN 201610962226A CN 106353744 A CN106353744 A CN 106353744A
- Authority
- CN
- China
- Prior art keywords
- signal
- distance
- fda
- delta
- sigma
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/581—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets
- G01S13/582—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to a multi-parameter combined estimation method based on bi-static FDA-MIMO radars. The method comprises the following steps: firstly designing a transmitting signal by utilizing characteristics of FDA and MIMO radars; carrying out matched filtering, vectorization and spatial smooth processing on a receiving signal; then estimating a combined steering vector and estimating a DOA and a speed parameter by utilizing an ESPRIT algorithm, and carrying out decoupling and parameter estimation on the DOD and distance information by combining characteristics of a transmitting waveform; and carrying out ambiguity resolution on a distance result estimated by utilizing an ESPRIT algorithm and combining a distance estimated by virtue of a pulse delay estimation algorithm, and carrying out ambiguity resolution on a speed by virtue of an MUSIC algorithm by combining signal characteristics of a large number of pulses. The method provided by the invention has the advantages that the problem that the distance and speed can not be accurately estimated under the condition of a single PRF can be effectively solved, and the estimation of the three-dimensional position and speed of a target can be realized; and a simulation result shows that the method provided by the invention has good estimation accuracy and stability.
Description
Technical field
The invention belongs to MIMO radar technical field, it is based on bistatic fda-mimo radar particularly to a kind of
Multiparameter combined estimation method.
Background technology
So-called mimo radar refers to the waveform using multiple transmitting antennas synchronously transmitting diversity, is connect using multiple simultaneously
Receive antenna and receive echo-signal, and focus on a kind of New Type Radar system of receiving and transmitting signal.Mimo radar and traditional phased array
The maximum difference of radar is that mimo radar can be realized on degree of freedom by changing array structure or the irrelevant waveform of transmitting
Lifting.According to " far and near " of antenna spatial distribution, mimo radar is broadly divided into 2 classes: centralized mimo radar and distributed
Mimo radar.The dual-mode antenna positional distance of centralized mimo radar system, compared with " near ", typically all can not reach incoherent mimo
Spacing size required by radar.Because each antenna can launch different signals, therefore there is good waveform diversity gain,
So as to parameter identification estimated capacity lifting, adaptive technique with and flexible waveform design etc. aspect show good
Effect.The antenna of distributed mimo radar transmit-receive array at a distance of relatively " remote " and is distributed in space diverse location.Each antenna is relatively
Angle in target has obvious difference, therefore shows good space diversity gain.The advantage of the type radar mainly exists
Yu Qineng makes full use of target radar scattering cross-section and amasss (rcs) fluctuation characteristic spatially, thus overcoming target rcs angle scintillations
Improve target detection performance and Parameter Estimation Precision.The present invention is mainly with the bistatic radar in centralized mimo radar as object
Studied.
The method of existing solution prf and contradiction in velocity estimation mainly has multi-pulse repetition period or irregular
The method of repetition.But this method still can produce ambiguous estimation.Because centralized mimo radar can be different by transmitting
Waveform realizes waveform diversity gain.Therefore, it is different from conventional radar, mimo radar can be combined with fda technology, utilize
Fda makes the carrier frequency of transmission signal form a frequency increment with emission array, thus utilizing single prf by waveform diversity
Solve the problems, such as ambiguous estimation.However, this method can produce the coupling of dod and range information on transmitting steering vector, lead to
It is unable to estimate out dod and range information in the case of bistatic mimo radar.A kind of solution is to be divided into emission array
Two submatrixs, the frequency increment of two submatrixs is different;Another kind of method is to make a transmitting pulse medium frequency increment be 0, next
Individual pulse medium frequency increment is not 0, alternate emission, thus using the signal under different pulses, estimate respectively angle information and
Range information.The present invention is carried out decoupling using the first solution.Simultaneously because when Doppler frequency shift and carrier frequency and pulse
Prolong number relevant, in the case of impulse time delay number is less, it can be ignored for the impact of transmitting steering vector, and is counting
Cannot ignore in the case that mesh is larger, therefore can estimate dod and distance in impulse time delay number hour, in impulse time delay number
Solve the problems, such as when larger that velocity estimation obscures.Simultaneously for coherent signal, method main at present is Space planar angle and square
Battle array Reconstruction Method.The former can make the aperture of radar reduce;Thus reducing estimated accuracy, the latter is for the transmission signal of radar
Have certain limitations with parameter estimation, not as the former flexibly, and also can there is the situation of aperture minimizing.
Content of the invention
For overcoming deficiency of the prior art, the present invention provides a kind of connection of the multiparameter based on bistatic fda-mimo radar
Close method of estimation, the problem estimated for bistatic centralized mimo radar parameter, in conjunction with fda technology, solve Sing plus weight
Frequently under (prf) distance estimations and velocity estimation fuzzy problem, realize target 3-dimensional positioning.
According to design provided by the present invention, a kind of multiparameter Combined estimator based on bistatic fda-mimo radar
Method, comprises the steps of:
Step 1, utilize waveform diversity characteristics design transmission signal, according to the frequency increment δ f of transmission signal carrier wave, obtain
The transmitted waveform relevant with dod and range information;
Step 2, the transmitted waveform signal receiving is carried out with matched filtering, vectorization and space smoothing are processed, expired
The signal covariance matrix of order;
Step 3, be based on signal covariance matrix, using esprit algorithm estimate joint steering vector;
Step 4, estimate doa and speed parameter using joint steering vector, and dod and range information are carried out decoupling with
Parameter estimation;
Step 5, information estimated result of adjusting the distance with reference to impulse time delay estimation carry out ambiguity solution process;
Step 6, combine high impulse number laUnder signal characteristic speed parameter estimated result is solved by music algorithm
Fuzzy Processing.
Above-mentioned, step 1 specifically comprises following content: radar emission array is divided into two with array center for reference point
Individual submatrix, the frequency increment of design submatrix 1 is-δ f, and submatrix 2 is that the transmitting steering vector obtaining is:
Wherein, θp
And rpIt is respectively target dod and doa, dtFor the spacing of transmitting antenna, drFor the spacing of reception antenna, λ is signal wavelength, and c is light
Speed, v is target velocity, rpFor target to transmitting terminal and receiving terminal apart from sum.
Above-mentioned, it is as follows that step 2 specifically comprises content:
Step 201: matched filtering process is carried out to receipt signal, obtain l (l=1,2 ..., l) the letter under individual pulse
Number, it is expressed as:
In formula,For receiving steering vector,For launching steering vector, ap(v, l) is multiple for target
Scattering coefficient and Doppler frequency shift, w (l) is noise vector;
Receipt signal under step 202:l pulse is expressed as: x=[x (1), x (2) ..., x (l)], it is carried out to
Quantification treatment obtains signal:
, then signal covariance matrix be expressed as: ry=e (yyh), wherein, ⊙ is that khatri-rao amasss, bl×pV how general () be
Strangle vector, h is target scattering coefficient;
Step 203: signal is carried out with space smoothing process, design (l, n) individual smoothing matrix is as follows:
, then the signal covariance matrix after smoothing is:
, wherein,Amass for kronecker, l0=l-pv+ 1, n0=n-pr+ 1, wherein pvAnd prRepresent respectively and guide to receiving
Vector doppler vector carries out the smooth number of times of space smoothing, h=diag (h) λt,
c0=z11c1, work as pvpr>=p and l0n0During >=p,It is the signal covariance matrix of full rank, wherein p is target total number.
Above-mentioned, it is as follows that step 3 specifically comprises content:
Step 301: feature decomposition is carried out to signal covariance matrix and obtains signal subspaceDue to
span{es}=span { c0, then esMeet es=c0t-1;
Step 302: by esIt is divided into two sub-spaces es1And es2, obtainWherein, usWith
T isCharacteristic vector,For containing the diagonal matrix of dod and velocity information;
Step 303: calculating joint steering vector is
Above-mentioned, it is as follows that described step 4 specifically comprises content:
Step 401: estimate doa and speed using joint steering vector, concrete formula is as follows:
Step 402: using the feature of transmission signal, dod and distance are carried out decoupling, concrete formula is as follows:
Step 403: dod and distance are calculated by formula, specific formula for calculation is as follows:
Above-mentioned, it is as follows that step 5 specifically comprises content: using the actual distance measured by radar pulse is
In formula, kpIt is integer, rut=c/fprfRepresent maximum unambiguous distance,Represent measurement distance;Using fda-mimo thunder
Reaching surveyed actual distance isIn formula, qpIt is integer, ruδf=c/4 δ f represents maximum no mould
Paste distance,Represent estimated distance;Then unambiguous distance can be estimated to draw by following formula:
Above-mentioned, it is as follows that step 6 specifically comprises content: the relation according to speed and frequency increment under big umber of pulse solves no
Fuzzy speed, laReceipt signal under pulse is:
In formula, transmitting steering vector:
It is assumed that true velocity isD is integer, vu=c/2f0T is velocity ambiguity, utilizes
Music Algorithm for Solving d:
In formula, cv1(la) and cv2(la) it is based on two joint steering vectors launching submatrixs.
Beneficial effects of the present invention:
1st, the present invention adopts the method for space smoothing to ensure signal covariance matrix full rank, and doppler vector and reception are led
Carry out space smoothing to vector simultaneously, solve and there is velocity estimation and distance under the Sing plus repetition period in prior art and estimate
Meter is also easy to produce fuzzy problem, improves parameter estimation performance.
2nd, the present invention makes full use of mimo radar waveform diversity feature, by mimo radar and frequency control battle array fda combine so as to
Transmitted waveform does not only have dod information and comprises range information even velocity information, thus realizing dod, doa, distance and speed yet
Combined estimator, solve the velocity estimation brought in traditional phased-array radar and distance by the pulse repetition period (prf)
Estimate the contradiction of no blur estimation, improve the performance of parameter estimation;Achieve the decoupling, thus allowing of dod and range information
Angle and distance Combined estimator is possibly realized, and solves the problems, such as that under single prf, distance estimations obscure;Using Doppler frequency shift and
The relation of carrier frequency, achieves the ambiguity solution of velocity estimation when impulse time delay number is larger using music algorithm, thus realizing mesh
Mark dod, effective estimation of doa, distance and 4 parameters of speed, so that it is determined that the three-dimensional coordinate of target;Utilization space smooths to be calculated
Method, to doppler vector and receive steering vector smoothing processing simultaneously, solve coherent signal and Angle Ambiguity Problem it is ensured that
The full rank condition of covariance matrix.
Brief description:
Fig. 1 is bistatic fda-mimo radar system structural representation;
Fig. 2 is the schematic flow sheet of the present invention;
Fig. 3 is the change curve with snr for the phase place launching steering vector under the present invention;
Fig. 4 is the change curve with snr for the root-mean-square error of lower 4 parameters of the present invention;
Fig. 5 is the change curve with sampled point for the root-mean-square error of lower 4 parameters of the present invention;
Fig. 6 be angle ambiguity under the conditions of 4 parameters root-mean-square error with snr change curve.
Specific embodiment:
The present invention is further detailed explanation with technical scheme below in conjunction with the accompanying drawings, and detailed by preferred embodiment
Describe bright embodiments of the present invention in detail, but embodiments of the present invention are not limited to this.
Embodiment one, referring to shown in Fig. 1~2, a kind of multiparameter Combined estimator side based on bistatic fda-mimo radar
Method, comprises the steps of:
Step 1, utilize waveform diversity characteristics design transmission signal, according to the frequency increment δ f of transmission signal carrier wave, obtain
The transmitted waveform relevant with dod and range information is it is ensured that can realize decoupling to dod and range information in parameter estimation procedure
Close;
Step 2, the transmitted waveform signal receiving is carried out with matched filtering, vectorization and space smoothing are processed, expired
The signal covariance matrix of order;
Step 3, be based on signal covariance matrix, using esprit algorithm estimate joint steering vector;
Step 4, estimate doa and speed parameter using joint steering vector, and dod and range information are carried out decoupling with
Parameter estimation;
Step 5, information estimated result of adjusting the distance with reference to impulse time delay estimation carry out ambiguity solution process;
Step 6, combine high impulse number laUnder signal characteristic speed parameter estimated result is solved by music algorithm
Fuzzy Processing.
The present invention adopts the method for space smoothing to ensure signal covariance matrix full rank, to doppler vector and reception guiding
Vector carries out space smoothing simultaneously, solves and there are velocity estimation and distance estimations under the Sing plus repetition period in prior art
It is also easy to produce fuzzy problem, improve parameter estimation performance.
Embodiment two, referring to shown in Fig. 1~6, a kind of multiparameter Combined estimator side based on bistatic fda-mimo radar
Method, comprises following content:
Step 1, the bistatic mimo radar being constituted based on m transmitting array element, n reception array element, special using waveform diversity
Property design transmission signal, according to the frequency increment δ f of transmission signal carrier wave, obtain the transmitted wave relevant with dod and range information
Shape, in order to solve the problems, such as dod and range information coupling, radar emission array is divided into two with array center for reference point
Submatrix, calculates for convenience, and the frequency increment of design submatrix 1 is-δ f, and submatrix 2 is that the transmitting steering vector obtaining is:
Wherein, θpAnd rpIt is respectively target dod and doa, dtFor the spacing of transmitting antenna, drFor the spacing of reception antenna, λ is signal wave
Long, c is the light velocity, and v is target velocity, rpFor target to transmitting terminal and receiving terminal apart from sum it is ensured that in parameter estimation procedure
Dod and range information can be realized decoupling.
Step 2, the transmitted waveform signal receiving is carried out with matched filtering, vectorization and space smoothing are processed, expired
The signal covariance matrix of order, specifically comprises following content:
Step 201: matched filtering process is carried out to receipt signal, obtain l (l=1,2 ..., l) the letter under individual pulse
Number, it is expressed as:
In formula,For receiving steering vector,For launching steering vector, ap(v, l) is multiple for target
Scattering coefficient and Doppler frequency shift, w (l) is noise vector;
Receipt signal under step 202:l pulse is expressed as: x=[x (1), x (2) ..., x (l)], it is carried out to
Quantification treatment obtains signal:
, then signal covariance matrix be expressed as: ry=e (yyh), wherein,It is that khatri-rao amasss, bl×pV how general () be
Strangle vector, h is target scattering coefficient;
Step 203: space smoothing process is carried out to signal, situation about annexing in order to avoid angle of arrival, designs (l, n)
Individual smoothing matrix is as follows:
, then the signal covariance matrix after smoothing is:
, wherein,Amass for kronecker, l0=l-pv+ 1, n0=n-pr+ 1, wherein pvAnd prRepresent respectively and guide to receiving
Vector doppler vector carries out the smooth number of times of space smoothing, h=diag (h) λt,
c0=z11c1, work as pvpr>=p and l0n0During >=p,It is the signal covariance matrix of full rank.
Step 3, be based on signal covariance matrix, according to the invariable rotary shape of fda-mimo radar signal, using esprit
Algorithm estimates joint steering vector, specifically comprises content as follows:
Step 301: feature decomposition is carried out to signal covariance matrix and obtains signal subspaceMeet es
=c0t-1, wherein, due to span { es}=span { c0, then esMeet es=c0t-1;
Step 302: by esIt is divided into two sub-spaces es1And es2, obtainWherein, usWith
T isCharacteristic vector,For containing the diagonal matrix of dod and velocity information;
Step 303: calculating joint steering vector is
Step 4, estimate doa and speed parameter using joint steering vector, and dod and range information are carried out decoupling with
Parameter estimation, specifically comprises content as follows:
Step 401: estimate doa and speed using joint steering vector, concrete formula is as follows:
Step 402: using the feature of transmission signal, dod and distance are carried out decoupling, concrete formula is as follows:
Step 403: dod and distance are calculated by formula, specific formula for calculation is as follows:
Step 5, information estimated result of adjusting the distance with reference to impulse time delay estimation carry out ambiguity solution process, using radar pulse institute
Measurement actual distance beIn formula, kpIt is integer, rut=c/fprfRepresent maximum unambiguous distance,Represent measurement distance;Using the actual distance that fda-mimo radar is surveyed it isFormula
In, qpIt is integer, ruδf=c/4 δ f represents maximum unambiguous distance,Represent estimated distance;Then unambiguous distance can be by following formula
Estimation draws:
Step 6, combine high impulse number laUnder signal characteristic speed parameter estimated result is solved by music algorithm
Fuzzy Processing, the relation according to speed and frequency increment under big umber of pulse solves no fuzzy speed, laReceipt signal under pulse
For:
In formula, transmitting steering vector:
It is assumed that true velocity isD is integer, vu=c/2f0T is velocity ambiguity, utilizes
Music Algorithm for Solving d:
In formula, cv1(la) and cv2(la) it is based on two joint steering vectors launching submatrixs.
The Cramér-Rao lower bound (cramer rao bound, crb) that doa estimates is unbiased esti-mator deviation theory lower bound.Foundation
The signal model of bistatic fda-mimo radar, the Cramér-Rao lower bound crb of 4 parameters is
Wherein, snr=(e | ξ |2)/(mσ2), l, m, n represent impulse time delay number, transmitting antenna number respectively and receive sky
Line number mesh, κv=2 π f0T/c, κθ=2 π dtCos (θ)/λ, κr=2 π δ f/c,K represents sampling
Count out.
The effect of the present invention can be further illustrated by following emulation experiment.
1) simulated conditions:
Bistatic fda-mimo radar transmit-receive array is even linear array, and parameter setting is as follows:
The parameter setting of the bistatic fda-mimo radar of table 1
2) emulation experiment:
(1) impact to transmitting steering vector for the frequency increment
The impact to transmitting steering vector for the frequency increment is mainly studied in this emulation.The dod of target is 20 °, and distance is
1.5km, speed is 450m/s.Fig. 3 gives the change with pulse for the phase place to launching steering vector for the frequency increment, and dod
With distance with impulse time delay number change curve.By simulation result it is found that due to speed and carrier frequency and impulse time delay number
Mesh is simultaneously relevant, but because its magnitude is relatively small, therefore when umber of pulse less period, the impact to transmitting steering vector is little,
When impulse time delay number is larger, its impact be can not ignore.Using this point, ambiguity solution can be carried out to velocity estimation.
(2) the parameter estimation performance of fda-mimo radar
The parameter estimation performance of fda-mimo radar is mainly studied in this emulation.It is to estimate that performance estimates mean square error first
Change curve with snr.Sampling number is that 256, dod is 30 °, and doa is 35 °, and distance is 150km, and speed is 450m/s, wherein
Distance and speed all can produce ambiguous estimation.By the simulation result of Fig. 4 it is found that with respect to traditional phased-array radar and mimo
Radar, fda-mimo radar can realize the Combined estimator of 4 parameters, and distance and velocity estimation no fuzzy.Simultaneously because
The reason mimo radar virtual aperture, its estimated accuracy increases compared with phased-array radar.
Followed by mean square error is with the change curve of sampling number.Snr=0db, coordinates of targets is constant.Emulation by Fig. 5
It is found that increase with sampling number, the estimated accuracy of radar is also improving constantly result.
(3) the solution angle ambiguity performance of fda-mimo radar
In order to verify the performance to angle ambiguity for the method for the present invention, we by the method for the present invention and are provided without space and put down
Sliding method is compared.There are 2 targets of angle ambiguity in setting, its dod is 30 ° and 45 °, and doa is 35 ° and 36 °, distance
For 150km and 250km, speed is 450m/s and 350m/s.Even if by the simulation result of Fig. 6 it is found that the method for the present invention
Remain to effectively the target that there is angle ambiguity be estimated under conditions of low signal-to-noise ratio.
The present invention is not limited to above-mentioned specific embodiment, and those skilled in the art also can make multiple changes accordingly, but
Any it is equal to the present invention or similar change all should be covered within the scope of the claims.
Claims (7)
1. a kind of multiparameter combined estimation method based on bistatic fda-mimo radar it is characterised in that: comprise the steps of:
Step 1, utilize waveform diversity characteristics design transmission signal, according to the frequency increment δ f of transmission signal carrier wave, obtain with
The dod transmitted waveform relevant with range information;
Step 2, the transmitted waveform signal receiving is carried out with matched filtering, vectorization and space smoothing are processed, and obtain full rank
Signal covariance matrix;
Step 3, be based on signal covariance matrix, using esprit algorithm estimate joint steering vector;
Step 4, utilization joint steering vector are estimated doa and speed parameter, and are carried out decoupling and parameter to dod and range information
Estimate;
Step 5, information estimated result of adjusting the distance with reference to impulse time delay estimation carry out ambiguity solution process;
Step 6, combine high impulse number laUnder signal characteristic by music algorithm, ambiguity solution is carried out to speed parameter estimated result
Process.
2. the multiparameter combined estimation method based on bistatic fda-mimo radar according to claim 1, its feature exists
In: described step 1 specifically comprises following content: radar emission array is divided into two submatrixs with array center for reference point,
The frequency increment of design submatrix 1 is-δ f, and submatrix 2 is that the transmitting steering vector of p-th target obtaining is:
Wherein, θpAnd rp
It is respectively target dod and doa, dtFor the spacing of transmitting antenna, drFor the spacing of reception antenna, λ is signal wavelength, and c is the light velocity, v
For target velocity, rpFor target to transmitting terminal and receiving terminal apart from sum.
3. the multiparameter combined estimation method based on bistatic fda-mimo radar according to claim 2, its feature exists
In: it is as follows that described step 2 specifically comprises content:
Step 201: matched filtering process is carried out to receipt signal, obtain l (l=1,2 ..., l) the signal under individual pulse, table
It is shown as:
In formula,For receiving steering vector,For launching steering vector, ap(v, l) scatters again for target
Coefficient and Doppler frequency shift, w (l) is noise vector;
Receipt signal under step 202:l pulse is expressed as: x=[[x (1), x (2) ..., x (l)]], enters row vector to it
Change processes and obtains signal:
,
Then signal covariance matrix is expressed as: ry=e (yyh), wherein, ⊙ is that khatri-rao amasss, bl×pV () is sweared for Doppler
Amount, h is target scattering coefficient;
Step 203: signal is carried out with space smoothing process, design (l, n) individual smoothing matrix is as follows:
Signal covariance matrix after then smoothing is:
Wherein,Amass for kronecker, l0=l-pv+ 1, n0=n-pr+ 1, wherein pvAnd prRepresent respectively to reception steering vector
Carry out the smooth number of times of space smoothing, h=diag (h) λ with doppler vectort,c0=
z11c1, work as pvpr>=p and l0n0During >=p,It is the signal covariance matrix of full rank, wherein p is target total number.
4. the multiparameter combined estimation method based on bistatic fda-mimo radar according to claim 3, its feature exists
In: it is as follows that described step 3 specifically comprises content:
Step 301: feature decomposition is carried out to signal covariance matrix and obtains signal subspaceDue to span
{es}=span { c0, then esMeet es=c0t-1;
Step 302: by esIt is divided into two sub-spaces es1And es2, obtainWherein, usIt is with tCharacteristic vector,For containing the diagonal matrix of dod and velocity information;
Step 303: calculating joint steering vector is
5. the multiparameter combined estimation method based on bistatic fda-mimo radar according to claim 4, its feature exists
In: it is as follows that described step 4 specifically comprises content:
Step 401: estimate doa and speed using joint steering vector, concrete formula is as follows:
Step 402: using the feature of transmission signal, dod and distance are carried out decoupling, concrete formula is as follows:
Step 403: dod and distance are calculated by formula, specific formula for calculation is as follows:
6. the multiparameter combined estimation method based on bistatic fda-mimo radar according to claim 5 it is characterised in that:
It is as follows that described step 5 specifically comprises content: using the actual distance measured by radar pulse isFormula
In, kpIt is integer, rut=c/fprfRepresent maximum unambiguous distance,Represent measurement distance;Using fda-mimo radar
The actual distance surveyed isIn formula, qpIt is integer, ruδf=c/4 δ f represents maximum no fuzzy
Distance,Represent estimated distance;Then unambiguous distance can be estimated to draw by following formula:
7. the multiparameter combined estimation method based on bistatic fda-mimo radar according to claim 6, its feature exists
In: it is as follows that described step 6 specifically comprises content: the relation according to speed and frequency increment under big umber of pulse solves no fuzzy speed
Degree, laReceipt signal under pulse is:
In formula, transmitting steering vector:
Assume that true velocity isD is integer, vu=c/2f0T is velocity ambiguity, using music algorithm
Solution d:
In formula, cv1(la) and cv2(la) it is based on two joint steering vectors launching submatrixs.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610962226.9A CN106353744B (en) | 2016-10-28 | 2016-10-28 | Multi-parameter combined estimation method based on bistatic FDA-MIMO radars |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610962226.9A CN106353744B (en) | 2016-10-28 | 2016-10-28 | Multi-parameter combined estimation method based on bistatic FDA-MIMO radars |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106353744A true CN106353744A (en) | 2017-01-25 |
CN106353744B CN106353744B (en) | 2018-09-11 |
Family
ID=57864470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610962226.9A Active CN106353744B (en) | 2016-10-28 | 2016-10-28 | Multi-parameter combined estimation method based on bistatic FDA-MIMO radars |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106353744B (en) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106950550A (en) * | 2017-03-31 | 2017-07-14 | 中国人民解放军海军航空工程学院 | Based on the high dynamic deviation On-line Estimation method adjudicated across fuzzy interval under the equal hazy condition of ranging and range rate |
CN107422317A (en) * | 2017-08-03 | 2017-12-01 | 中国人民解放军空军工程大学 | Low angle target arrival direction estimation method based on smoothing matrix collection |
CN107907860A (en) * | 2017-11-15 | 2018-04-13 | 电子科技大学 | A kind of Location fraud method that the anti-time difference measurement of battle array technology is controlled using frequency |
CN108152836A (en) * | 2018-01-03 | 2018-06-12 | 电子科技大学 | A kind of deception measures based on the anti-Samsung frequency difference location method of FDA antennas |
CN108196239A (en) * | 2018-01-18 | 2018-06-22 | 西安电子科技大学 | A kind of frequency diversity MIMO radar without fuzzy parameter method of estimation |
CN108363058A (en) * | 2018-03-06 | 2018-08-03 | 电子科技大学 | Frequency controls the signal parameter design method of battle array imaging radar |
CN108594233A (en) * | 2018-04-24 | 2018-09-28 | 森思泰克河北科技有限公司 | A kind of velocity solution blur method based on MIMO car radars |
CN108594186A (en) * | 2017-08-25 | 2018-09-28 | 西安电子科技大学 | The method that FDA-MIMO radars inhibit main lobe Deceiving interference |
CN108732545A (en) * | 2018-03-06 | 2018-11-02 | 电子信息系统复杂电磁环境效应国家重点实验室 | The frequency control battle array signal acceptance method and structure of matched filtering are mixed based on multichannel |
CN108761414A (en) * | 2018-05-25 | 2018-11-06 | 电子科技大学 | A kind of S-shaped interfering beam test method for controlling battle array based on frequency |
CN108802718A (en) * | 2018-05-30 | 2018-11-13 | 北京理工大学 | Phase decoupling method when based on random exomonental time-division MIMO radar sky |
CN109085541A (en) * | 2018-09-07 | 2018-12-25 | 杭州捍鹰科技有限公司 | MIMO radar array antenna and its signal processing method |
CN109375213A (en) * | 2018-08-08 | 2019-02-22 | 西安电子科技大学 | Frequency diversity array signal processing method based on Subarray partition |
CN109782242A (en) * | 2018-10-29 | 2019-05-21 | 中国传媒大学 | A kind of object localization method based on parallel factor model in monostatic radar systems |
CN110146871A (en) * | 2019-05-21 | 2019-08-20 | 西安电子科技大学 | Target component estimation method based on the inclined FDA-MIMO radar of double frequency |
CN110320489A (en) * | 2019-06-26 | 2019-10-11 | 中国电子科技集团公司第三十八研究所 | A kind of distribution isomorphism face battle array combines the balancing method and system for estimating angular accuracy |
CN110346764A (en) * | 2019-08-21 | 2019-10-18 | 上海无线电设备研究所 | A kind of MIMO radar target range angle decoupling method |
CN110907925A (en) * | 2018-09-14 | 2020-03-24 | 哈尔滨工业大学 | Weight positioning method under high-frequency ground wave radar double-station model |
CN111239719A (en) * | 2020-01-20 | 2020-06-05 | 西安邮电大学 | Stable high-precision low-altitude target height measurement method based on composite reflection coefficient estimation |
WO2021042482A1 (en) * | 2019-09-04 | 2021-03-11 | 南京慧尔视智能科技有限公司 | Method for designing transmitted waveform parameter under mimo system |
CN113287036A (en) * | 2019-09-30 | 2021-08-20 | 华为技术有限公司 | Speed ambiguity resolution method and echo signal processing device |
CN113777599A (en) * | 2021-09-13 | 2021-12-10 | 电子科技大学 | FDA-MIMO radar-based speed ambiguity resolution method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011158430A (en) * | 2010-02-03 | 2011-08-18 | Mitsubishi Electric Corp | Transceiver beam shaping device |
CN104297734A (en) * | 2014-10-20 | 2015-01-21 | 西安电子科技大学 | Deception interference rejection method for MIMO radar based on frequency diversity array |
CN105699945A (en) * | 2016-01-30 | 2016-06-22 | 湖北工业大学 | Waveform optimized design method for frequency controlled array MIMO radar system |
CN105717496A (en) * | 2016-01-30 | 2016-06-29 | 湖北工业大学 | Realization method of FDA (Frequency Diverse Array) MIMO (Multiple-Input Multiple-Output) radar system based on matrix completion |
CN105785328A (en) * | 2016-03-15 | 2016-07-20 | 西安电子科技大学 | Subarray-division-based FDA distance-angle decoupling wave beam formation method |
-
2016
- 2016-10-28 CN CN201610962226.9A patent/CN106353744B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011158430A (en) * | 2010-02-03 | 2011-08-18 | Mitsubishi Electric Corp | Transceiver beam shaping device |
CN104297734A (en) * | 2014-10-20 | 2015-01-21 | 西安电子科技大学 | Deception interference rejection method for MIMO radar based on frequency diversity array |
CN105699945A (en) * | 2016-01-30 | 2016-06-22 | 湖北工业大学 | Waveform optimized design method for frequency controlled array MIMO radar system |
CN105717496A (en) * | 2016-01-30 | 2016-06-29 | 湖北工业大学 | Realization method of FDA (Frequency Diverse Array) MIMO (Multiple-Input Multiple-Output) radar system based on matrix completion |
CN105785328A (en) * | 2016-03-15 | 2016-07-20 | 西安电子科技大学 | Subarray-division-based FDA distance-angle decoupling wave beam formation method |
Non-Patent Citations (1)
Title |
---|
《IET RADAR, SONAR & NAVIGATION》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106950550A (en) * | 2017-03-31 | 2017-07-14 | 中国人民解放军海军航空工程学院 | Based on the high dynamic deviation On-line Estimation method adjudicated across fuzzy interval under the equal hazy condition of ranging and range rate |
CN106950550B (en) * | 2017-03-31 | 2020-03-03 | 中国人民解放军海军航空大学 | High dynamic deviation on-line estimation method based on cross-fuzzy interval judgment under condition of range finding and speed measuring ambiguity |
CN107422317A (en) * | 2017-08-03 | 2017-12-01 | 中国人民解放军空军工程大学 | Low angle target arrival direction estimation method based on smoothing matrix collection |
CN108594186B (en) * | 2017-08-25 | 2020-04-07 | 西安电子科技大学 | Method for inhibiting main lobe deceptive interference by FDA-MIMO radar |
CN108594186A (en) * | 2017-08-25 | 2018-09-28 | 西安电子科技大学 | The method that FDA-MIMO radars inhibit main lobe Deceiving interference |
CN107907860A (en) * | 2017-11-15 | 2018-04-13 | 电子科技大学 | A kind of Location fraud method that the anti-time difference measurement of battle array technology is controlled using frequency |
CN108152836A (en) * | 2018-01-03 | 2018-06-12 | 电子科技大学 | A kind of deception measures based on the anti-Samsung frequency difference location method of FDA antennas |
CN108196239A (en) * | 2018-01-18 | 2018-06-22 | 西安电子科技大学 | A kind of frequency diversity MIMO radar without fuzzy parameter method of estimation |
CN108196239B (en) * | 2018-01-18 | 2022-03-18 | 西安电子科技大学 | Unambiguous parameter estimation method for frequency diversity MIMO radar |
CN108363058A (en) * | 2018-03-06 | 2018-08-03 | 电子科技大学 | Frequency controls the signal parameter design method of battle array imaging radar |
CN108732545A (en) * | 2018-03-06 | 2018-11-02 | 电子信息系统复杂电磁环境效应国家重点实验室 | The frequency control battle array signal acceptance method and structure of matched filtering are mixed based on multichannel |
CN108594233A (en) * | 2018-04-24 | 2018-09-28 | 森思泰克河北科技有限公司 | A kind of velocity solution blur method based on MIMO car radars |
CN108761414A (en) * | 2018-05-25 | 2018-11-06 | 电子科技大学 | A kind of S-shaped interfering beam test method for controlling battle array based on frequency |
CN108761414B (en) * | 2018-05-25 | 2022-05-20 | 电子科技大学 | S-shaped interference wave beam testing method based on frequency control array |
CN108802718A (en) * | 2018-05-30 | 2018-11-13 | 北京理工大学 | Phase decoupling method when based on random exomonental time-division MIMO radar sky |
CN108802718B (en) * | 2018-05-30 | 2021-11-26 | 北京理工大学 | Time division MIMO radar space-time phase decoupling method based on random transmitted pulse |
CN109375213A (en) * | 2018-08-08 | 2019-02-22 | 西安电子科技大学 | Frequency diversity array signal processing method based on Subarray partition |
CN109085541A (en) * | 2018-09-07 | 2018-12-25 | 杭州捍鹰科技有限公司 | MIMO radar array antenna and its signal processing method |
CN110907925B (en) * | 2018-09-14 | 2023-03-24 | 哈尔滨工业大学 | Weight positioning method under high-frequency ground wave radar double-station model |
CN110907925A (en) * | 2018-09-14 | 2020-03-24 | 哈尔滨工业大学 | Weight positioning method under high-frequency ground wave radar double-station model |
CN109782242A (en) * | 2018-10-29 | 2019-05-21 | 中国传媒大学 | A kind of object localization method based on parallel factor model in monostatic radar systems |
CN110146871A (en) * | 2019-05-21 | 2019-08-20 | 西安电子科技大学 | Target component estimation method based on the inclined FDA-MIMO radar of double frequency |
CN110320489A (en) * | 2019-06-26 | 2019-10-11 | 中国电子科技集团公司第三十八研究所 | A kind of distribution isomorphism face battle array combines the balancing method and system for estimating angular accuracy |
CN110346764A (en) * | 2019-08-21 | 2019-10-18 | 上海无线电设备研究所 | A kind of MIMO radar target range angle decoupling method |
WO2021042482A1 (en) * | 2019-09-04 | 2021-03-11 | 南京慧尔视智能科技有限公司 | Method for designing transmitted waveform parameter under mimo system |
CN113287036A (en) * | 2019-09-30 | 2021-08-20 | 华为技术有限公司 | Speed ambiguity resolution method and echo signal processing device |
CN111239719A (en) * | 2020-01-20 | 2020-06-05 | 西安邮电大学 | Stable high-precision low-altitude target height measurement method based on composite reflection coefficient estimation |
CN111239719B (en) * | 2020-01-20 | 2023-08-01 | 西安邮电大学 | Robust high-precision low-altitude target height measurement method based on composite reflection coefficient estimation |
CN113777599A (en) * | 2021-09-13 | 2021-12-10 | 电子科技大学 | FDA-MIMO radar-based speed ambiguity resolution method |
Also Published As
Publication number | Publication date |
---|---|
CN106353744B (en) | 2018-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106353744A (en) | Multi-parameter combined estimation method based on bi-static FDA-MIMO radars | |
CN102156279B (en) | Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output) | |
CN104515971B (en) | Airborne single-station passive positioning method for multiple broadband targets | |
CN101915911B (en) | Space arbitrary configuration distributed SAR moving target parameter estimation method based on cancellation accumulation space-time spectrum | |
CN103383448B (en) | Clutter suppression method suitable for high pulse repetition frequency (HPRF) waveform airborne radar | |
CN102520395B (en) | Clutter suppression method based on bistatic multiple-input and multiple-output radar | |
CN106291540A (en) | A kind of multiple-input and multiple-output GPR backwards projection target imaging method estimated based on DOA | |
CN108693511B (en) | Moving target angle calculation method of time division multiplexing MIMO radar | |
CN109188387B (en) | Target parameter estimation method for distributed coherent radar based on interpolation compensation | |
CN102135617A (en) | Multi-target positioning method of bistatic multi-input multi-output radar | |
CN103901417A (en) | Low-complexity space target two-dimensional angle estimation method of L-shaped array MIMO radar | |
CN103364772B (en) | Target low elevation estimation method based on real number field generalized multiple-signal sorting algorithm | |
CN103901410B (en) | Airborne bistatic MIMO radar clutter suppression method based on sparse recovery | |
CN102707264A (en) | Estimating method of direction of arrival of bistatic MIMO (Multi-Input Multi-Output) radar based on circular array | |
CN103760546B (en) | A kind of radar low target Wave arrival direction estimating method | |
CN103197294B (en) | Elevation angle estimating method of multi-frequency fusion maximum likelihood low-altitude target | |
CN103257344B (en) | Iteration-adaptive-algorithm-based method for detecting coherent MIMO radar target | |
CN103983958A (en) | MIMO radar continuous target angle estimation method based on multi-measurement-vector sparse representation | |
CN103412286B (en) | Transmitting polarization optimizing DOA (direction of arrival) evaluation method based on MIMO (multiple-input multiple-output) radar | |
CN102520399B (en) | Electromagnetic vector array based angle estimation method for metric-wave radar | |
CN104502898B (en) | The maneuvering target method for parameter estimation that modified R FT and amendment MDCFT are combined | |
CN108872971B (en) | Target positioning method and device based on motion single array | |
CN103383449A (en) | ESPRIT algorithm based short-range clutter suppression method for airborne radar | |
CN103760529A (en) | Efficient cascading space-time adaptive processing method based on passive detection | |
CN105044667A (en) | Double-satellite tracking method, device and system for moving target |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |