CN104808190A - Steady waveform design method for improving worst parameter estimation performance of MIMO (multiple-input and multiple-output) radar - Google Patents
Steady waveform design method for improving worst parameter estimation performance of MIMO (multiple-input and multiple-output) radar Download PDFInfo
- Publication number
- CN104808190A CN104808190A CN201510164914.6A CN201510164914A CN104808190A CN 104808190 A CN104808190 A CN 104808190A CN 201510164914 A CN201510164914 A CN 201510164914A CN 104808190 A CN104808190 A CN 104808190A
- Authority
- CN
- China
- Prior art keywords
- theta
- circletimes
- centerdot
- delta
- beta
- 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/42—Diversity systems specially adapted for radar
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)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
A steady waveform design method for improving the worst parameter estimation performance of an MIMO radar belongs to the field of signal processing.
Description
Technical field
The invention belongs to signal transacting field, further relate to the sane waveform design method of MIMO radar the poorest parameter estimation performance under the improvement clutter environment of Waveform Design technical field.
Background technology
In the last few years, the optimization of MIMO radar waveform was subject to the attention of increasing scholar and slip-stick artist.According to the object module used in waveform optimization problem, current waveform optimization method can be divided into following two classes: (1) is based on the waveform optimization of point target (point target); (2) based on the waveform optimization of Extended target (extended target).Based on the Waveform Design of point target, the object of optimization is waveform Correlation Matrix (WCM, waveform covariance matrix) or radar ambiguity function (radar ambiguityfunction).Waveform optimization method based on WCM only designs the spatial domain of transmitted waveform instead of the overall feature of transmitted waveform.Specifically, the people such as D.R.Fuhrmann and G.S.Antonio designs to realize specific energy spatial domain to WCM and distributes.And the people such as S.Peter have not only paid close attention to energy spatial domain and distribute, and have also contemplated that the spatial domain cross-correlation between different target, namely minimize spatial domain cross-correlation between different azimuth to improve the detection estimated performance of system.
Under the hypothesis that Received signal strength is not polluted by the clutter depending on transmitted waveform, the people such as J.Li propose the waveform optimization criterion of a few class based on CRB to optimize WCM thus the Parameter Estimation Precision of raising point target.Received signal strength is by under clutter pollution condition, and the people such as H.Y Wang consider target prior imformation and know under condition based on the MIMO radar waveform of CRB and the combined optimization problem of biased estimator.It should be noted that in these methods, solving of waveform optimization problem all needs parameter to know.But in Practical Project, these parameters must be obtained by estimation, thus inevitably there is evaluated error.Thus, optimizing the parameter estimation performance that obtains of waveform to evaluated error and uncertainty based on estimated parameter is than more sensitive.
Summary of the invention
The object of the invention under being to overcome clutter conditions traditional waveform optimization method to initial parameter evaluated error sensitive issue, propose a kind of sane waveform design method improving MIMO radar the poorest parameter estimation performance, the method comprises Parameter uncertainties convex set, based on the alternative manner solving-optimizing problem of DL technology, to alleviate parameter estimating error or the uncertain system sensitivity problem brought, thus improve the MIMO radar waveform optimize parameter estimate performance under worst case.
The basic ideas of the inventive method are: first MIMO radar Received signal strength model under structure clutter conditions, based on the CRB of this model inference parameter to be estimated, then set up the explicit sane waveform optimization model comprising parameter uncertainty, for solving this nonlinear optimal problem, propose a kind of iterative algorithm based on DL technology, each step of iteration can relax as Semidefinite Programming thus obtain Efficient Solution.After obtaining an optimum intermediate solution based on iterative algorithm, optimum waveform covariance matrix can reconstruct under least square meaning.
The present invention improves the sane waveform design method of MIMO radar the poorest parameter estimation performance, and it comprises the steps:
Step one, structure MIMO radar Received signal strength model
Suppose that MIMO radar Received signal strength is:
Wherein,
for being proportional to the complex magnitude of target RCS (radar cross section),
for target location parameter, K is target numbers, and ρ (θ) is for being in the reflection coefficient of θ position clutter block, and W represents interference noise, and often row are separate and with distribution circle symmetric complex random vector, have zero-mean, its covariance B the unknown,
for the matrix that transmits, a (θ
k) and v (θ
k) represent reception respectively, launch steering vector, be specifically expressed as:
In formula, f
0for carrier frequency, τ
m(θ
k), m=1,2 ... M
rwith
for the transmission time, a
c(θ) and v
c(θ) θ is represented respectively
kthe reception of place's target and transmitting steering vector;
If rang ring is divided into N
c(N
c> > NML) individual resolution element, MIMO radar Received signal strength model is rewritten as
Wherein,
represent clutter transport function, ρ (θ
i) be θ
ithe reflection coefficient of place's clutter block, N
c(N
c> > M
tm
r) be clutter spatial sampling quantity, a
c(θ
i) and v
c(θ
i) represent θ respectively
ithe reception of place's clutter block, transmitting steering vector; Vec (H
c) be the same multiple Gaussian random vector distributed, its average is zero, and covariance is
Step 2, build sane waveform optimization model based on CRB
Consider unknown parameter θ=[θ
1, θ
2..., θ
k]
t,
cRB under condition, through deriving, this CRB can be expressed as follows:
Wherein,
In formula, R
s=SS
h,
For positive semidefinite Hermitian matrix,
Only consider direction of arrival angle, namely θ evaluated error is on the impact of system performance, can be as follows to a kth destination channel matrix modeling:
Wherein,
h
kbe respectively a kth destination channel matrix that is actual and hypothesis, δ
kfor
error, belong to following convex set:
And,
wherein,
be respectively h
kreal and suppose vector reciprocal,
for
error, belong to following convex set:
Based on foregoing, improving the sane waveform optimization problem of worst condition parameter estimation performance under clutter conditions can be expressed as: about under the constraint of WCM, based on Parameter uncertainties convex set
optimize WCM to minimize the CRB under worst case; Under Trace-opt criterion, optimization problem can be described as:
tr(R
S)=LP
Wherein, P represents total emissive power; In formula, the 3rd constraint establishment is because each transmitter unit emissive power can not be less than zero;
Solving of step 3, sane waveform internal layer optimization problem
Internal layer optimization problem solve based on following lemma 1:
Lemma 1. supposes that A is the positive semidefinite hermitian matrix of a M × M, then the inequality below is set up:
time and if only if A is diagonal matrix, equation is set up; According to lemma 1, internal layer optimization problem can relax and be:
Based on CRB, above formula can be rewritten as:
Real part operational character Re{} is got in deletion, is because in above formula, each and item are real numbers;
From above formula, and in formula, the denominator of kth item only depends on δ
kwith
two, the problem equivalent therefore in above formula in, under corresponding constraint, maximize and each in formula, can be expressed as:
For solving above formula, to R
sapplication diagonal angle loading technique, that is:
Wherein, ε < < λ
max(R
s) be load factor, λ
max() representing matrix eigenvalue of maximum, selects ε=λ
max(R
s)/1000; Use respectively
replace the R in sane optimization problem
s, can obtain
with
respectively for
and δ
kconvex;
Thus, above formula can be rewritten as:
Above formula can be torn open and be write as two independently minimization problems below:
s.t. ||δ
k||
F≤ζ
k
Above-mentioned two minimization problems can be solved by lemma 2 below:
Lemma 2, suppose hermitian matrix
's
then and if only if
time, Z
wherein, Δ C=A-B
hc
-1b is that the Schur of C in Z mends;
Lemma 2 by reference, above-mentioned two minimization problems can be converted into following SDP problem 1:
Wherein, t is auxiliary variable;
By what obtain from above two formulas
with
bring in sane optimization problem, consider outside optimization problem;
Solving of step 4, the sane outer optimization problem of waveform
Following proposition is utilized to solve outside optimization problem
Proposition: utilize matrix manipulation, the constraint in sane optimization problem can be equivalent to following LMI:
Wherein
Use lemma 2 in conjunction with above-mentioned proposition, outer optimization problem can be expressed as following SDP problem:
Wherein, X is an auxiliary variable;
After obtaining optimum E, under least square meaning, R
sby following model construction:
s.t. tr(R
S)=LP
Use lemma 2 in conjunction with above-mentioned proposition, above formula can be equivalent to following SDP problem 2:
tr(R
S)=LP
Step 5, employing alternative manner solve sane waveform optimization problem
Step 5.1, given waveform covariance matrix initial value;
Step 5.2, solve above-mentioned SDP problem 1 to obtain optimum δ
k,
Step 5.3, solve SDP problem 2 to obtain optimum E;
Step 5.4, return step 5.2 iteration again, until CRB no longer significantly reduces.
Step 6, based on least square method, the waveform covariance matrix that reconstruct is optimum, can obtain R
s.
The invention has the beneficial effects as follows: the method can be used for releiving traditional waveform optimization method to parameter estimating error and uncertain sensitive issue.First set up MIMO radar Received signal strength model under clutter scene, characterize the lower bound-Cramér-Rao lower bound (CRB) of Parameter Estimation Precision to be estimated based on this model inference, then Parameter uncertainties convex set explicitly is comprised in traditional waveform optimization problem; For solving this nonlinear optimal problem, the present invention proposes a kind of alternative manner loading (DL) technology based on diagonal angle, each step in iteration all can be converted into Semidefinite Programming (SDP) problem, thus can Efficient Solution be obtained, to realize the sane waveform optimization of MIMO radar under worst case, and then the parameter estimation performance under worst case is promoted, compared with irrelevant waveform, the method has obvious lifting to the parameter estimation performance under worst case.
Accompanying drawing explanation
Fig. 1 is the process flow diagram that the present invention realizes;
Fig. 2 is the process flow diagram of iterative algorithm of the present invention;
Fig. 3 is the optimum transmit beam direction figure of the present invention when initial angle exists evaluated error and array signal to noise ratio (S/N ratio) is 10dB;
Fig. 4 under there is evaluated error situation in initial angle, the present invention carry the CRB changed with ANSR that algorithm and irrelevant waveform obtain.
Fig. 5 is the optimum transmit beam direction figure of the present invention when array calibration exists evaluated error and array signal to noise ratio (S/N ratio) is 10dB;
Fig. 6 under there is evaluated error situation in array calibration, the present invention put forward the worst case that algorithm and irrelevant waveform obtain under the CRB that changes with ANSR.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail.
As shown in Figures 1 to 6, to improve the implementation procedure of the sane waveform design method of MIMO radar the poorest parameter estimation performance as follows in the present invention:
1, sane waveform optimization problem model is set up
1) MIMO radar signal model is built
Suppose that MIMO radar Received signal strength is:
Wherein,
for being proportional to the complex magnitude of target RCS (radar cross section),
for target location parameter, both need to estimate.K is target numbers, and ρ (θ) is for being in the reflection coefficient of θ position clutter block, and W represents interference noise, and often row are separate and with distribution circle symmetric complex random vector, have zero-mean, its covariance B the unknown,
for the matrix that transmits.A (θ
k) and v (θ
k) represent reception respectively, launch steering vector, be specifically expressed as:
In formula, f
0for carrier frequency, τ
m(θ
k), m=1,2 ... M
rwith
for the transmission time, a
c(θ) and v
c(θ) θ is represented respectively
kthe reception of place's target and transmitting steering vector.
If rang ring is divided into N
c(N
c> > NML) individual resolution element, Received signal strength model can be rewritten as
Wherein,
represent clutter transport function, ρ (θ
i) be θ
ithe reflection coefficient of place's clutter block, N
c(N
c> > M
tm
r) be clutter spatial sampling quantity, a
c(θ
i) and v
c(θ
i) represent θ respectively
ithe reception of place's clutter block, transmitting steering vector.Vec (H
c) may be thought of as the same multiple Gaussian random vector distributed, its average is zero, and covariance is
R
hccan also be expressed as further:
wherein,
2) the sane waveform optimization model based on CRB is built
Consider unknown parameter θ=[θ
1, θ
2..., θ
k]
t,
cRB under condition, through deriving, this CRB can be expressed as follows:
Wherein,
In formula, R
s=SS
h,
For positive semidefinite Hermitian matrix,
Significantly, CRB be about θ,
h
c, W function, and these parameters need by estimate obtain, thus inevitably there is evaluated error.Thus, the systematic parameter estimated performance utilizing the CRB optimization waveform based on certain group estimates of parameters to obtain is worse than other one group of more rational estimates of parameters possibly.And then, in engineer applied, the tender subject of system performance to initial parameter evaluated error must be considered.
In the present invention, only consider direction of arrival angle, namely θ evaluated error is on the impact of system performance.Thus, can be as follows to a kth destination channel matrix modeling:
Wherein,
h
kbe respectively a kth destination channel matrix that is actual and hypothesis, δ
kfor
error, belong to following convex set:
And,
wherein,
be respectively h
kreal and suppose vector reciprocal,
for
error, belong to following convex set:
Based on above-mentioned discussion, improving the sane waveform optimization problem of worst condition parameter estimation performance under clutter conditions can be expressed as: about under the constraint of WCM, based on Parameter uncertainties convex set
optimize WCM to minimize the CRB under worst case.Under Trace-opt criterion, optimization problem can be described as:
tr(R
S)=LP
Wherein, P represents total emissive power; In formula, the 3rd constraint establishment is because each transmitter unit emissive power can not be less than zero.
Clearly, CRB matrix trace, the i.e. objective function of above formula, be one about R
sand δ
k,
k=1,2 ..., the very complicated nonlinear function of K, thus utilizes the classic methods such as such as convex optimization to be very difficult to solve.
2. the solving of sane waveform optimization problem
1) the solving of internal layer optimization problem
As mentioned above, the objective function of optimization problem is very complicated nonlinear function, is difficult to utilize traditional Optimization Method.For solving this problem, first consider internal layer optimization problem.Internal layer optimization problem solve based on following lemma 1:
Lemma 1. supposes that A is the positive semidefinite hermitian matrix of a M × M, then the inequality below is set up:
time and if only if A is diagonal matrix, equation is set up.According to lemma 1, internal layer optimization problem can relax and be:
Based on CRB, above formula can be rewritten as:
Real part operational character Re{} is got in deletion, is because in above formula, each and item are real numbers.
From above formula, and in formula, the denominator of kth item only depends on δ
kwith
two, the problem equivalent therefore in above formula in, under corresponding constraint, maximize and each in formula, can be expressed as:
It should be noted that due to
known
for indefinite matrix, therefore above formula is difficult to solve, for solving this problem, to R
sapplication diagonal angle loading technique, that is:
Wherein, ε < < λ
max(R
s) be load factor, λ
max() representing matrix eigenvalue of maximum, below in l-G simulation test, select ε=λ
max(R
s)/1000.Use respectively
replace the R in sane optimization problem
s, can obtain
clearly,
with
respectively for
and δ
kconvex.
Thus, above formula can be rewritten as:
Similar, above formula can be write as two independently minimization problems below:
s.t. ||δ
k||
F≤ζ
k
Above two problems can be solved by lemma 2 below:
Hermitian matrix is supposed in lemma 2.
's
then and if only if
time,
wherein, Δ C=A-B
hc
-1b is that the Schur of C in Z mends.
Lemma 2 by reference, above two problems can clearly be converted into following SDP problem:
Wherein, t is auxiliary variable.
By what obtain from above two formulas
with
bring in sane optimization problem, consider outside optimization problem.
2) the solving of outside optimization problem
The present invention utilizes following proposition for solving outside optimization problem
Proposition: utilize matrix manipulation, the constraint in sane optimization problem can be equivalent to following LMI:
Wherein
Use lemma 2 and in conjunction with above-mentioned proposition, outer optimization problem can be expressed as following SDP problem:
Wherein, X is an auxiliary variable.
After obtaining optimum E, under least square meaning, R
sby following model construction:
s.t. tr(R
S)=LP
Be similar to above-mentioned discussion, above formula can be equivalent to following SDP problem:
tr(R
S)=LP
3) iterative algorithm
Given WCM initial value, δ
k,
and R
sbe optimized by following steps:
1. solve internal layer SDP problem and obtain optimum δ
k,
2. solve outer SDP problem and obtain E;
Repeat step 1. 2., until CRB no longer significantly reduces.After this, solve constructed model R
s.
Effect of the present invention further illustrates by following emulation:
Simulated conditions:
MIMO radar is 33 and receives, utilize two MIMO radar system, its antenna configuration is respectively: MIMO radar (0.5,0.5), MIMO radar (1.5,0.5) the numeral transmitter, here in bracket and the array element distance (in units of wavelength) in receiver.It is 256 that systematic sampling is counted.Array signal to noise ratio (S/N ratio) is defined as
span is-10dB to 30dB.Wherein, P refers to total emissive power,
for additive white Gaussian noise variance.Modeling clutter is discrete sampling, and its RCS is modeled as independent identically distributed gaussian random variable vector, and average is zero, and variance is
and suppose coherent processing inteval internal fixtion.Clutter signal to noise ratio (S/N ratio) is defined as
equal 30dB.-5
°there is a strong jamming in direction, and signal to noise ratio (S/N ratio) is 60dB.Only in θ=20, ° place has a reflection coefficient to be the point target of 1.In following emulation, suppose two kinds of situations, the first only considers that initial angle estimates to there is error; It two is the correction errors only considering to exist in transmitting-receiving array.
Emulation content:
A: initial angle estimates to there is uncertain situation
Suppose that the uncertainty that initial angle is estimated is Δ θ=[-3 °, 3 °], namely
wherein
for the estimation of θ, through calculating, obtain data: MIMO (0.5,0.5) is ζ=5.4382, σ=7.6593, MIMO (1.5,0.5) is ζ=27.6329, σ=29.6754.
Fig. 3 is optimum transmit beam direction figure under ASNR=10dB condition.Can observe, the peak value of the beam pattern that transmits is positioned at around target location, this means, under this convex uncertain worst case, systematic parameter estimated performance can improve.In addition, due to sparse emission array, MIMO radar (1.5,0.5) there will be graing lobe situation, as shown in Fig. 3 (b).
The CRB of Fig. 4 for being changed with ASNR by proposed algorithm and uncorrelated waveform gained.Clearly, CRB reduces with the increase of ASNR.In addition, can observe, under the worst case that institute's extracting method obtains, parameter estimation performance is better than uncorrelated waveform.And along with the increase of ASNR, institute extracting method gained CRB is progressive in uncorrelated waveform.In addition, the CRB of MIMO radar (1.5,0.5) shown in Fig. 4 (b) is starkly lower than the CRB of MIMO radar (0.5,0.5) shown in Fig. 4 (a).
B: transmitting-receiving array exists the situation of correction error
In this case, be no matter that transmitting and receiving array is assumed that to have correction error (amplitude of sensor and phase error and site error).Each element of transmitting and receiving array steering vector disturb by the disturbance variable, this disturbance variable is the multiple Gaussian random variable of Cyclic Symmetry of zero-mean, and variance is
after calculating, obtain ζ=13.4764 of MIMO (0.5,0.5), σ=14.5712, ζ=29.8362 of MIMO (1.5,0.5), σ=32.6573.
Fig. 5 features the optimum transmit beam pattern that ASNR=10dB obtains.From Fig. 5, the conclusion similar in appearance to Fig. 3 can be drawn.Carry CRB under the worst case that algorithm and uncorrelated waveform gained change with ASNR as shown in Figure 6, the conclusion obtained from Fig. 6 is similar to Fig. 4.
In sum, under the present invention is directed to clutter conditions, waveform optimization method is to initial parameter evaluated error sensitive issue, propose the sane waveform optimization method based on the convex uncertain collection of parameter, and propose a kind of iterative method loaded based on diagonal angle for this complex nonlinear optimization problem.For the sane performance of parameter estimation of MIMO radar system under raising clutter conditions, first there is error condition for direction of arrival angle and carry out modeling in the present invention, and this parameter estimating error convex set explicitly is comprised in afferent echo shape optimization problem, for solving this nonlinear optimal problem, the present invention proposes a kind of alternative manner loaded based on diagonal angle and carries out alternative optimization to transmitted waveform and parameter estimating error, to obtain optimum transmitted waveform covariance matrix.Each step of iteration all can load based on diagonal angle and relax as Semidefinite Programming, thus can obtain Efficient Solution.Known based on discussing above, institute of the present invention extracting method be can be the sane performance estimated by design transmitted waveform raising radar parameter in engineer applied and provides solid theory and realize foundation.
Claims (1)
1. improve the sane waveform design method of MIMO radar the poorest parameter estimation performance, it is characterized in that, the method comprises the steps:
Step one, structure MIMO radar Received signal strength model
Suppose that MIMO radar Received signal strength is:
Wherein,
for being proportional to the complex magnitude of target RCS (radar cross section),
for target location parameter, K is target numbers, and ρ (θ) is for being in the reflection coefficient of θ position clutter block, and W represents interference noise, and often row are separate and with distribution circle symmetric complex random vector, have zero-mean, its covariance B the unknown,
for the matrix that transmits, a (θ
k) and v (θ
k) represent reception respectively, launch steering vector, be specifically expressed as:
In formula, f
0for carrier frequency, τ
m(θ
k), m=1,2 ... M
rwith
n=1,2 ... M
tfor the transmission time, a
c(θ) and v
c(θ) θ is represented respectively
kthe reception of place's target and transmitting steering vector;
If rang ring is divided into N
c(N
c>>NML) individual resolution element, MIMO radar Received signal strength model is rewritten as
Wherein,
represent clutter transport function, ρ (θ
i) be θ
ithe reflection coefficient of place's clutter block, N
c(N
c>>M
tm
r) be clutter spatial sampling quantity, a
c(θ
i) and v
c(θ
i) represent θ respectively
ithe reception of place's clutter block, transmitting steering vector; Vec (H
c) be the same multiple Gaussian random vector distributed, its average is zero, and covariance is
Step 2, build sane waveform optimization model based on CRB
Consider unknown parameter θ=[θ
1, θ
2..., θ
k]
t,
cRB under condition, through deriving, this CRB can be expressed as follows:
Wherein,
In formula, R
s=SS
h,
for positive semidefinite Hermitian matrix,
Only consider direction of arrival angle, namely θ evaluated error is on the impact of system performance, can be as follows to a kth destination channel matrix modeling:
Wherein,
be respectively a kth destination channel matrix that is actual and hypothesis, δ
kfor
error, belong to following convex set:
And,
wherein,
be respectively h
kreal and suppose vector reciprocal,
for
error, belong to following convex set:
Based on foregoing, improving the sane waveform optimization problem of worst condition parameter estimation performance under clutter conditions can be expressed as: about under the constraint of WCM, based on Parameter uncertainties convex set
, optimize WCM to minimize the CRB under worst case; Under Trace-opt criterion, optimization problem can be described as:
Wherein, P represents total emissive power; In formula, the 3rd constraint establishment is because each transmitter unit emissive power can not be less than zero;
Solving of step 3, sane waveform internal layer optimization problem
Internal layer optimization problem solve based on following lemma 1:
Lemma 1. supposes that A is the positive semidefinite hermitian matrix of a M × M, then the inequality below is set up:
time and if only if A is diagonal matrix, equation is set up; According to lemma 1, internal layer optimization problem can relax and be:
Based on CRB, above formula can be rewritten as:
real part operational character Re{} is got in deletion, is because in above formula, each and item are real numbers;
From above formula, and in formula, the denominator of kth item only depends on δ
kwith
two, the problem equivalent therefore in above formula in, under corresponding constraint, maximize and each in formula, can be expressed as:
Wherein, ε << λ
max(R
s) be load factor, λ
max() representing matrix eigenvalue of maximum, selects ε=λ
max(R
s)/1000; Use respectively
replace the R in sane optimization problem
s, can obtain
with
respectively for
and δ
kconvex;
Thus, above formula can be rewritten as:
Above formula can be torn open and be write as two independently minimization problems below:
Above-mentioned two minimization problems can be solved by lemma 2 below:
Lemma 2, suppose hermitian matrix
's
, then and if only if
time,
, wherein, Δ C=A-B
hc
-1b is that the Schur of C in Z mends;
Lemma 2 by reference, above-mentioned two minimization problems can be converted into following SDP problem 1:
Wherein, t is auxiliary variable;
By what obtain from above two formulas
with
bring in sane optimization problem, consider outside optimization problem;
Solving of step 4, the sane outer optimization problem of waveform
Following proposition is utilized to solve outside optimization problem
Proposition: utilize matrix manipulation, the constraint in sane optimization problem can be equivalent to following LMI:
Wherein
Use lemma 2 in conjunction with above-mentioned proposition, outer optimization problem can be expressed as following SDP problem:
Wherein, X is an auxiliary variable;
After obtaining optimum E, under least square meaning, R
sby following model construction:
Use lemma 2 in conjunction with above-mentioned proposition, above formula can be equivalent to following SDP problem 2:
Step 5, employing alternative manner solve sane waveform optimization problem
Step 5.1, given waveform covariance matrix initial value;
Step 5.2, solve above-mentioned SDP problem 1 to obtain optimum
Step 5.3, solve SDP problem 2 to obtain optimum E;
Step 5.4, return step 5.2 iteration again, until CRB no longer significantly reduces.
Step 6, based on least square method, the waveform covariance matrix that reconstruct is optimum, can obtain R
s.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510164914.6A CN104808190B (en) | 2015-04-09 | 2015-04-09 | Improve the sane waveform design method of the worst parameter Estimation performance of MIMO radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510164914.6A CN104808190B (en) | 2015-04-09 | 2015-04-09 | Improve the sane waveform design method of the worst parameter Estimation performance of MIMO radar |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104808190A true CN104808190A (en) | 2015-07-29 |
CN104808190B CN104808190B (en) | 2018-02-23 |
Family
ID=53693178
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510164914.6A Expired - Fee Related CN104808190B (en) | 2015-04-09 | 2015-04-09 | Improve the sane waveform design method of the worst parameter Estimation performance of MIMO radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104808190B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974391A (en) * | 2016-04-28 | 2016-09-28 | 大连大学 | MIMO (multiple-input multiple-output) radar robust waveform design method with target prior knowledge unknown |
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 |
CN110493811A (en) * | 2019-08-15 | 2019-11-22 | 南京航空航天大学 | Based on the stealthy steady resource allocation methods of radar-communication integration system of radio frequency |
CN111399395A (en) * | 2020-03-23 | 2020-07-10 | 武汉科技大学 | Implementation method of F-M II state space model based on radar target prediction system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103852750A (en) * | 2014-01-28 | 2014-06-11 | 大连大学 | Robust MIMO radar waveform optimization method for improving worst estimated performance |
CN103852749A (en) * | 2014-01-28 | 2014-06-11 | 大连大学 | Robust waveform optimization method for improving MIMO-STAP detection performance |
CN104375121A (en) * | 2014-01-28 | 2015-02-25 | 大连大学 | Combined optimizing method of MIMO radar waveform and biased estimator based on prior information |
-
2015
- 2015-04-09 CN CN201510164914.6A patent/CN104808190B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103852750A (en) * | 2014-01-28 | 2014-06-11 | 大连大学 | Robust MIMO radar waveform optimization method for improving worst estimated performance |
CN103852749A (en) * | 2014-01-28 | 2014-06-11 | 大连大学 | Robust waveform optimization method for improving MIMO-STAP detection performance |
CN104375121A (en) * | 2014-01-28 | 2015-02-25 | 大连大学 | Combined optimizing method of MIMO radar waveform and biased estimator based on prior information |
Non-Patent Citations (6)
Title |
---|
HONGYAN WANG等: "Joint optimization of MIMO radar waveform and biased estimator with prior information in the presence of clutter", 《JOURNAL ON ADVANCES IN SIGNAL PROCESSING》 * |
HONGYAN WANG等: "Robust Waveform Design for MIMO Radar to Improve the Worst-Case Estimation Performance in the Presence of Clutter", 《PROCEEDING OF GENERAL ASSEMBLY AND SCIENTIFIC SYMPOSIUM,2014 XXXITH URSI》 * |
HONGYAN WANG等: "Waveform Optimization for MIMO Radar Based on Cramer-Rao Bound in the Presence of Clutter", 《IEICE TRANS. COMMUN.》 * |
张裔东: "杂波环境下MIMO雷达波形优化研究", 《中国优秀硕士学位论文全文数据库 I136-471》 * |
王洪雁: "MIMO雷达波形优化", 《中国博士学位论文全文数据库 I136-71》 * |
王洪雁等: "提高MIMO雷达估计性能的稳健波形优化设计", 《电子信息对抗技术》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974391A (en) * | 2016-04-28 | 2016-09-28 | 大连大学 | MIMO (multiple-input multiple-output) radar robust waveform design method with target prior knowledge unknown |
CN105974391B (en) * | 2016-04-28 | 2018-09-25 | 大连大学 | The non-steady waveform design method of MIMO radar under the conditions of knowing target priori |
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 |
CN110320489B (en) * | 2019-06-26 | 2022-02-08 | 中国电子科技集团公司第三十八研究所 | Measurement method and system for distributed isomorphic area array joint angle estimation precision |
CN110493811A (en) * | 2019-08-15 | 2019-11-22 | 南京航空航天大学 | Based on the stealthy steady resource allocation methods of radar-communication integration system of radio frequency |
CN111399395A (en) * | 2020-03-23 | 2020-07-10 | 武汉科技大学 | Implementation method of F-M II state space model based on radar target prediction system |
CN111399395B (en) * | 2020-03-23 | 2022-11-25 | 武汉科技大学 | Implementation method of F-M II state space model based on radar target prediction system |
Also Published As
Publication number | Publication date |
---|---|
CN104808190B (en) | 2018-02-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103364764B (en) | Airborne radar non-stationary clutter suppression method | |
CN104808180A (en) | Steady waveform optimizing method for MIMO radar in clutter background | |
CN104267379B (en) | A kind of active radar and passive radar based on Waveform Design works in coordination with anti-interference method | |
CN102156279B (en) | Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output) | |
CN103760519B (en) | High resolution DOA estimation method and system | |
CN103245956B (en) | A kind of GPS anti-multipath method based on robust ada-ptive beamformer algorithm | |
CN107167785A (en) | A kind of sane big array MIMO radar target transmitting-receiving angle combined estimation method | |
CN105182313A (en) | MIMO-STAP steady waveform design method based on incomplete clutter prior knowledge | |
CN104833959A (en) | MIMO radar waveform optimization method based on target prior information | |
CN107290730A (en) | Bistatic MIMO radar angle evaluation method under array mutual-coupling condition | |
CN104808190A (en) | Steady waveform design method for improving worst parameter estimation performance of MIMO (multiple-input and multiple-output) radar | |
CN103926573B (en) | Mono-static MIMO radar distribution type target angle estimation method based on fourth-order cumulant | |
CN101770022A (en) | Multiple input multiple output (MIMO) radar array position error self-correcting method based on genetic algorithm | |
CN105807275A (en) | MIMO-OFDM-STAP steady waveform design method based on partial clutter priori knowledge | |
CN102135617A (en) | Multi-target positioning method of bistatic multi-input multi-output radar | |
CN104808179A (en) | Cramer-rao bound based waveform optimizing method for MIMO radar in clutter background | |
CN104991232A (en) | Signal-level fusion networking radar anti-cheating interference method under object signal correlation | |
CN109633591A (en) | External illuminators-based radar is biradical away from localization method under a kind of observation station location error | |
CN103852749A (en) | Robust waveform optimization method for improving MIMO-STAP detection performance | |
CN102928828B (en) | Phase difference estimation method of distributed radar based on orthogonal waveforms | |
CN105022040A (en) | Array element error estimation method based on clutter data combined fitting | |
CN104391279A (en) | Ionosphere propagation characteristic based phase diameter disturbance suppression method | |
CN104950296A (en) | Robustness nonhomogeneity detecting method based on heavily weighted adaptive power residue | |
CN104991233A (en) | Networking radar anti-cheating interference method based on signal level fusion | |
CN108226893A (en) | A kind of low complex degree receiver design method for MIMO radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180223 Termination date: 20190409 |
|
CF01 | Termination of patent right due to non-payment of annual fee |