CN103592642A - Method for designing MIMO radar waveforms - Google Patents
Method for designing MIMO radar waveforms Download PDFInfo
- Publication number
- CN103592642A CN103592642A CN201310499194.XA CN201310499194A CN103592642A CN 103592642 A CN103592642 A CN 103592642A CN 201310499194 A CN201310499194 A CN 201310499194A CN 103592642 A CN103592642 A CN 103592642A
- Authority
- CN
- China
- Prior art keywords
- waveform
- signal
- energy function
- abs
- frequency interval
- 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
Images
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/28—Details of pulse systems
- G01S7/282—Transmitters
-
- 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/46—Indirect determination of position data
- G01S13/48—Indirect determination of position data using multiple beams at emission or reception
-
- 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/28—Details of pulse systems
- G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
-
- 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)
Abstract
The invention discloses a method for designing MIMO radar waveforms. The method mainly solves the problem in the prior art that tolerance of Doppler is poor in the process of waveform design. The implementation steps of the method include the first step of fixing a set of waveform center frequency sequences of a radar and setting an initial phase to be zero, the second step of searching a frequency interval initial variable to find out a frequency interval initial value which satisfies the requirement that the 3dB bandwidth of an emitted energy function equals to the 3dB bandwidth of an ideal emitted energy function, the third step of setting a frequency interval according to the frequency interval initial value, the fourth step of inputting the frequency interval and the initial phase and outputting an emitted energy function difference value, emitting direction sidelobe gain and the sidelobe magnitude of signals after pulse synthesis to construct a cost function, the fifth step of calling the cost function and optimizing the cost function through the minimax method to obtain the frequency interval and the initial phase meeting the requirement, and the sixth step of obtaining a set of final signals according to the frequency interval and the initial phase and completing the design of the waveforms. The waveforms designed through the method have the advantages of being good in the tolerance of the Doppler and applicable to target detection.
Description
Technical field
The invention belongs to Radar Technology field, relate to a kind of method for designing of MIMO radar waveform, can be used for target detection.
Background technology
Multiinputoutput MIMO radar is a kind of emerging active detection technology.Its radiating antenna and receiving antenna can carry out flexible arrangement according to system requirements and each radiating element can be launched different signal waveforms.Transmitted waveform comprises orthogonal waveforms and simple crosscorrelation waveform, and common needs is simple crosscorrelation waveform.MIMO radar can utilize waveform diversity ability to design neatly simple crosscorrelation waveform, and electromagnetic energy is focused in interested target or on region, thereby effectively utilizes electromagnetic energy and suppress unnecessary noise jamming.In order to design better simple crosscorrelation waveform, current design cycle mainly: first according to required emitted energy function optimization, go out signal covariance matrix; According to this covariance matrix, utilize phase encoding comprehensively to go out required constant modulus signals waveform again.Its disclosed existing document has following two kinds:
1, the people such as Stoica Petre and Li Jian are at " On probing signal design for MIMO radar, IEEE Trans, on Signal Processing, 2007, Vol, 55 (8), 4151-4161 " in the covariance matrix that a kind of planning of positive semidefinite effectively (SDP) algorithm carrys out modelled signal is proposed, then at " Waveform Synthesis for Diversity-Based Transmit Beampattern Design, IEEE Trans, on Signal Processing, 2008, Vol, 56 (6), 2593-2598 " in round-robin algorithm has been proposed, by covariance matrix, designed the signal matrix of permanent mould.The peak sidelobe of the emitted energy function that this method is synthetic is higher, and in the situation that array element is more, computation complexity is high.
2, the people such as Hu Liangbing adopts positive semidefinite to plan that (SDP) optimizes the covariance matrix transmitting in " MIMO radar emission directional diagram coupling and waveform optimization method; Xian Electronics Science and Technology University's journal (natural science edition); 2009; 36 (6): 1021-1026 ", then based on this covariance matrix, the permanent mould that utilizes Chebyshev approximation to design continuous phase coding by optimizing waveform phase transmits.The peak sidelobe of this method emitted energy function is lower, but has equally the problem that computation complexity is higher.
Above-mentioned two kinds of methods are owing to all adopting phase encoding to cause its Doppler's tolerance poor, and particularly, when target velocity is larger, the synthetic impulse of target echo signal is processed seriously mismatch, can't detect target.
Summary of the invention
The object of the invention is to the deficiency for above-mentioned prior art, propose the method for designing of the good MIMO radar waveform of a kind of Doppler's tolerance, to reduce the synthetic impulse of target echo signal, process mismatch degree, realize the accurate detection to target.
The technical thought that realizes the object of the invention is: the good feature of Doppler's tolerance of utilizing linear FM signal, the frequency coding waveform of employing based on linear FM signal, by adjusting each waveform frequency interval and initial phase, meet the requirement that emitted energy function approaches the low distance side lobe of desirable emitted energy function, synthetic impulse figure and the low angle secondary lobe of transmit direction function, its technical scheme comprises the steps:
(1) the transmitting element number of array of establishing MIMO radar antenna is N, and carrier frequency is f
0, the pulse width of single waveform is Te, the total bandwidth of N waveform is B, and the bandwidth B s of each single waveform is identical, a fix N waveform centre frequency f
korder, i.e. f
1<f
2< ... <f
n-1<f
n, and by the first phase of N waveform
all be made as 0, k=1 wherein, 2,3 ..., N;
(2) define each waveform frequency and be spaced apart Δ f
m=f
m+1-f
m, establish each waveform frequency interval delta f
mbe equal to frequency interval initializaing variable Δ f
0, and between 0 to 1Te to frequency interval initializaing variable Δ f
0search for, find out the frequency interval initial value Δ f that three dB bandwidth that the three dB bandwidth that meets emitted energy function equals desirable emitted energy function requires
0', m=1 wherein, 2,3 ..., N-1;
(3) at frequency interval initial value Δ f
0near N-1 value of ' random generation, is assigned to respectively waveform frequency interval delta f
m, the centre frequency that obtains the 1st signal is f
1=f
0-sum (Δ f
mthe centre frequency f of i the signal in)/2, the
ifor f
i=f
i-1+ Δ f
m, i=2 wherein, 3 ..., N and m=i-1, sum () is summation operation;
(4) according to frequency interval Δ f
m, first phase
calculate bandwidth B s=B-sum (the Δ f of individual signals
m), then calculate chirp rate μ by μ=Bs/Te, thus baseband signal waveform s obtained
kfor:
Wherein, j is unit imaginary constant, and t is the L point sampling time in 0～Te, and exp () is for take the exponential function that natural constant e is the end;
(5) according to baseband waveform signal s
k, obtain signal covariance matrix R=ss
h, signal waveform matrix s=[s wherein
1; s
2; s
k; s
n], and then obtain emitted energy function and be:
P(θ)=a(θ)
TRa(θ)
*，
Wherein θ is the angle of-90～90 degree scopes, a (θ)=exp (j2 π ((1:N)
t-(N+1)/2) be * d.*G/c*sin (θ)) steering vector of angle θ, d is array element radius, and c is the light velocity, centre frequency vector G=[f
1, f
2..., f
k..., f
n]
t, ()
*for conjugate operation, ()
tfor transposition computing .* is point multiplication operation symbol;
(6) establishing transmit beam main lobe direction is β, according to signal covariance matrix R, obtains transmit direction function:
f(θ,β)=a(θ)
TRa(β)
*，
A (β)=exp (j2 π ((1:N) wherein
t-(N+1)/2) be * d.*G/c*sin (β)) steering vector of transmit beam main lobe direction;
(7) angle of establishing any one target in the angular range of desirable emitted energy function three dB bandwidth is φ, by signal waveform matrix s, obtains echoed signal sr=a (φ)
ts, carries out synthetic impulse to this echoed signal, and the signal st (φ) obtaining after synthetic impulse is processed is:
st(φ)=xcorr(a(φ)
Ts)，
A (φ)=exp (j2 π ((1:N) wherein
t-(N+1)/2) be * d.*G/c*sin (φ)) steering vector of target direction, xcorr () is auto-correlation computation;
(8) according to signal st (φ), transmit direction function f (θ after synthetic impulse processing, β), emitted energy function P (θ), calculate the maximum secondary lobe amplitude y2 of signal after synthetic impulse, the maximum difference p of the maximum side lobe gain y3 of transmit direction function and emitted energy function P (θ) and desirable emitted energy function Pt (θ), and then obtain row vector x:
x=[a*y2,b*y3,c*p]，
Wherein, a, b, c are greater than three different coefficients of 0, meet the different demands of waveform performance by adjusting its size;
(9) with frequency interval Δ f
m, first phase
as input, take row vector x as output, obtain cost function:
(10) with each waveform frequency interval delta f
m, first phase
for variable, call cost function
utilize Chebyshev approximation to be optimized, each waveform frequency interval delta f after being optimized
m' and first phase
make it meet low distance side lobe, low angle secondary lobe, and waveform emitted energy function approach the demand of desirable emitted energy function;
(11) with each waveform frequency interval delta f after optimizing
m' and first phase
try to achieve the final center frequency value f of each waveform
k', the final bandwidth value Bs ' of single waveform and final chirp rate value μ ', and then obtain final each signal s
k':
The present invention compared with prior art has the following advantages:
1, the present invention is due to employing linear FM signal and with frequency interval Δ f
m, first phase
for variable is optimized, thereby one group of designed linear FM signal, can not only meet low distance side lobe, low angle secondary lobe, and signal emitted energy function approaches the demand of desirable emitted energy function, and the synthetic impulse that has reduced target echo signal is processed mismatch degree, be that designed linear FM signal has good Doppler's tolerance, reduced in addition calculated amount.
Theoretical analysis and simulation result show, compared with prior art, Doppler's tolerance of waveform improves very large in the present invention, and secondary lobe and the emitted energy function of the distance side lobe of synthetic impulse figure, transmit direction function can meet the demands.
Accompanying drawing explanation
Fig. 1 is implementing procedure figure of the present invention;
Fig. 2 is Doppler's tolerance comparison diagram of signal of the present invention and the designed signal of prior art;
Fig. 3 is the designed signal emitted energy function comparison diagram of signal of the present invention and prior art;
Fig. 4 is the synthetic impulse functional arrangement of the designed signal of prior art;
Fig. 5 is the synthetic impulse functional arrangement of signal of the present invention;
Fig. 6 is the transmit direction function comparison diagram of signal of the present invention and the designed signal of prior art.
Embodiment
With reference to Fig. 1, concrete steps of the present invention are as follows:
The transmitting element number of array of supposing MIMO radar antenna is N, and carrier frequency is f
0, the pulse width of single waveform is Te, the total bandwidth of N waveform is B, and the bandwidth B s of each single waveform is identical; A fix N waveform centre frequency f
korder, i.e. f
1<f
2< ... <f
n-1<f
n, and by the first phase of N waveform
all be made as 0, k=1 wherein, 2,3 ..., N, N>=8.
Step 2, obtains the frequency interval initial value Δ f of centre frequency
0'.
2a) define each waveform frequency and be spaced apart Δ f
m=f
m+1-f
m, make each waveform frequency interval delta f
mbe equal to frequency interval initializaing variable Δ f
0, m=1 wherein, 2,3 ..., N-1;
2b) discrete value of getting uniformly Q point between 0 to 1/Te, and by its ascending [h that is arranged as
1', h
2' ..., h
n' ..., h
q'], n=1 wherein, 2 ..., Q, Q>=10;
2c) make frequency interval initializaing variable Δ f
0equal n discrete value h
n', according to frequency interval initializaing variable Δ f
0, the initial bandwidth of calculating individual signals: Bs
0=B-(N-1) Δ f
0; According to initial bandwidth Bs
0calculate the initial chirp rate of its correspondence: μ
0=Bs
0/ Te;
2d) according to frequency interval initializaing variable Δ f
0, calculate the initial center frequency η of each signal
k:
η
k=f
0-(N-1)Δf
0/2+(k-1)Δf
0。
2e) according to the initial center frequency η of each signal
k, initial chirp rate μ
0, carrier frequency f
0, the pulse width T e of individual signals and the first phase of each signal
obtaining initial baseband signal is gk:
Wherein, j is unit imaginary constant, and t is the L point sampling time in 0～Te, and exp () is for take the exponential function that natural constant e is the end;
2f) obtain initial transmissions energy function:
First, by each initial baseband signal g
kform initialize signal matrix g=[g
1; g
2; g
k; g
n];
Then, by initialize signal matrix computations, obtain initialize signal covariance matrix R
0=gg
h;
Finally, according to covariance matrix, obtain initial transmissions energy function: P
0(θ)=a ' is (θ)
tr
0a ' (θ)
*;
Wherein θ is the angle of-90～90 degree emitted energy coverages, and a ' is (θ) the initial steering vector of angle θ, a ' (θ)=exp (j2 π ((1:N)
t-(N+1)/2) * d.*G '/c*sin (θ)), d is array element radius, and c is the light velocity, initial center frequency vector G '=[η
1, η
2..., η
k..., η
n]
t, ()
*for conjugate operation, ()
tfor transposition computing .* is point multiplication operation symbol;
2g) find out initial transmissions energy function P
0(θ) maximal value Pm=max (abs (P
0(θ))), then obtain equation abs (P
0(θ) separate θ for two of)=0.5*Pm
1, θ
2, θ wherein
1< θ
2, abs () is for asking modular arithmetic;
2h) by above-mentioned two, separate θ
1, θ
2, calculate the three dB bandwidth θ that tries to achieve n initial transmissions energy function
w(n)=θ
2-θ
1;
2i) compare each frequency interval initializaing variable Δ f
0the three dB bandwidth θ of corresponding initial transmissions energy function
w(n) with the size of the three dB bandwidth of desirable emitted energy function, the corresponding frequency interval initializaing variable of the emitted energy function three dB bandwidth Δ f of desirable emitted energy function three dB bandwidth will be approached most
0value, as needed frequency interval initial value Δ f
0'.
Step 3, arranges each frequency interval Δ f
m, and calculate each centre frequency f
k.
3a) at frequency interval initial value Δ f
0' near, [max (0, Δ f
0'-0.3/Te), min (1/Te, Δ f
0'+0.3/Te)] interval interior random N-1 the value that produce, by its composition of vector [w
1, w
2..., w
m..., w
n-1], and make waveform frequency interval delta f
m=w
m;
3b) according to frequency interval Δ f
mwith carrier frequency f
0, the centre frequency that calculates the 1st signal is f
1=f
0-sum (Δ f
mthe centre frequency f of i the signal in)/2, the
ifor f
i=f
i-1+ Δ f
m, i=2 wherein, 3 ..., N and m=i-1, sum () is summation operation.
Step 4, procurement process baseband signal s
k.
4a) according to frequency interval Δ f
mwith total bandwidth B, calculate bandwidth B s=B-sum (the Δ f of individual signals
m), then by the bandwidth B s of individual signals and the pulse width T e of individual signals, calculate chirp rate: μ=Bs/Te;
4b) according to each centre frequency f
k, chirp rate μ and first phase
obtain the process baseband signal s in calculating process
kfor:
By process baseband waveform signal s
kform signal matrix s=[s
1; s
2; s
k; s
n]; According to signal matrix, calculate signal covariance matrix R=ss
h;
Utilize the steering vector a (θ) of signal covariance matrix and angle θ, obtain emitted energy function and be:
P(θ)=a(θ)
TRa(θ)
*，
Wherein a (θ) is the steering vector of angle θ, a (θ)=exp (j2 π ((1:N)
t-(N+1)/2) * d.*G/c*sin (θ)), centre frequency vector G=[f
1, f
2..., f
k..., f
n]
t.
Step 6, obtains transmit direction function f (θ, β).
If transmit beam main lobe direction is β, according to signal covariance matrix R, obtain transmit direction function:
f(θ,β)=a(θ)
TRa(β)
*，
A (β)=exp (j2 π ((1:N) wherein
t-(N+1)/2) be * d.*G/c*sin (β)) steering vector of transmit beam main lobe direction.
Step 7, obtains synthetic impulse and processes rear signal st (φ).
If the angle of any one target in the angular range of desirable emitted energy function three dB bandwidth is φ, by signal matrix s, obtain echoed signal sr=a (φ)
ts, carries out synthetic impulse to this echoed signal, and the signal st (φ) obtaining after synthetic impulse is processed is:
st(φ)=xcorr(a(φ)
Ts)，
A (φ)=exp (j2 π ((1:N) wherein
t-(N+1)/2) be * d.*G/c*sin (φ)) steering vector of target direction, xcorr () is auto-correlation computation.
Step 8, obtains row vector x.
8a) discretely in the angular range of desirable emitted energy function three dB bandwidth evenly get J angle on target value, by its ascending [φ that is arranged as
1', φ
2' ..., φ
z' ..., φ
j'], z=1 wherein, 2 ..., J, J>=10, and make angle on target φ=φ
z', the signal after synthetic impulse processing is st (φ
z'), and then obtain the signal matrix after synthetic impulse corresponding to all angle discrete values processed:
St=[st(φ
1′);st(φ
2′);…;st(φ
z′);…;st(φ
J′)]；
8b) the signal matrix St after synthetic impulse is processed presses row normalization, obtains normalized signal matrix Sg after synthetic impulse:
Sg=[abs(st(φ
1′))/max(abs(st(φ
1′)));abs(st(φ
2′))/max(abs(st(φ
2′)));…;abs(st(φ
J′))/max(abs(st(φ
J′)))]；
8c) main lobe of normalized signal matrix S g after synthetic impulse is partly taken as to 0, i.e. Sg (1:J, L+ (3:3))=0, and then obtain maximum secondary lobe amplitude corresponding to all discrete values of angle on target φ:
y2=max(max(Sg))；
8e) discretely in the angular range of desirable emitted energy function three dB bandwidth get M angle value, and by its ascending [β that is arranged as
1', β
2' ..., β
q' ..., β
m'], q=1 wherein, 2 ..., M, and make transmit beam main lobe direction β=β
q', transmit direction function f (θ, β
q')=a (θ)
tra (β
q')
*, obtain transmit direction Jacobian matrix corresponding to all discrete angle values of transmit beam main lobe direction β:
F(θ,β)=[abs(f(θ,β
1′)),abs(f(θ,β
2′)),…,abs(f(θ,β
M′))]；
8f) according to transmit direction Jacobian matrix F (θ, β), calculate normalized transmit direction Jacobian matrix and be:
Fg=[abs(f(θ,β
1'))/Fa(1),abs(f(θ,β
2'))/Fa(2),…,abs(f(θ,β
M'))/Fa(M)]，
Wherein, Fa=[max (abs (f (θ, β
1'))), max (abs (f (θ, β
2'))) ..., max (abs (f (θ, β
m')))] be the vector of the maximal value composition of each row of transmit direction Jacobian matrix F (θ, β);
8g) main lobe in transmit direction Jacobian matrix Fg is partly taken as to zero, is about to meet in matrix F g q row | θ-β
q| the element of <3 ° is taken as 0, and then obtains the maximum side lobe gain of the transmit direction function of all angle betas:
y3=max(max(Fg))；
8h) according to emitted energy function P (θ), calculate the maximal value of emitted energy function: Pa=max (abs (P (θ))); Calculate again the maximum difference p of emitted energy function P (θ) and desirable emitted energy function Pt (θ):
p=max(abs(abs(P(θ))/Pa-abs(Pt(θ))))；
8i), according to the maximum difference p of the maximum side lobe gain y3 of the maximum secondary lobe amplitude y2 of signal after synthetic impulse, transmit direction function and emitted energy function P (θ) and desirable emitted energy function Pt (θ), obtain row vector:
x=[a*y2,b*y3,c*p]，
Wherein, a, b, c are greater than three different coefficients of 0, meet the different demands of waveform performance by adjusting its size, and this example is got a=2, b=3, c=1.
Step 9, structure cost function.
With frequency interval Δ f
m, first phase
as input, take row vector x as output, obtain cost function:
With each waveform frequency interval delta f
m, first phase
for variable, call cost function
utilize Chebyshev approximation to be optimized, each waveform frequency interval delta f after being optimized
m' and first phase
make it meet low distance side lobe, low angle secondary lobe, and waveform emitted energy function approach the demand of desirable emitted energy function.
Described Chebyshev approximation, is under constraint condition, variable to be optimized, and the maximal value of a plurality of output valves is minimized, thereby makes all output valves minimum.This example is to be optimized variable by constrained optimization function f minimax in Calling MATLAB, detailed process can reference papers people's such as " Multipurpose Optimal Method Chebyshev approximation is inquired into " and Hu Liangbing " MIMO radar emission directional diagram mates and waveform optimization method; Xian Electronics Science and Technology University's journal (natural science edition); 2009,36 (6): 1021-1026 ".
Step 11, the baseband signal s of acquisition final design
k'.
11a) first according to each waveform frequency interval delta f after optimizing
m', calculate final bandwidth value Bs '=B-sum (Δ f of single waveform
m'); Then according to the pulse width T e of the final bandwidth value Bs ' of single waveform and signal, obtain final chirp rate value μ '=Bs '/Te;
11b) according to each waveform frequency interval delta f after optimizing
m', the final centre frequency that calculates the 1st signal is f
1'=f
0-sum (Δ f
m') the final centre frequency f of/2, a i signal
i' be f
i'=f
i-1'+Δ f
m', i=2 wherein, 3 ..., N and m=i-1;
11c) according to the final centre frequency f of each signal
k', final chirp rate μ ' and optimize after first phase
and then obtain one group of final signal s
k':
So far complete the design to MIMO radar waveform.
Effect of the present invention can further illustrate by following emulation experiment:
1. simulated conditions
The antenna that transmitting-receiving is put altogether adopts the uniformly-spaced even half-wavelength linear array of N=64 array element, carrier frequency f
0for 1.5GHz, interested target area is within the scope of 10～20 degree;
Emulation mode, adopts the signal of the present invention's design and the continuous phase coded signal designing by prior art.Wherein: wide Te=400us during signal pulse of the present invention, signal total bandwidth B=500kHz; The signal pulse width Te=100us of prior art design, signal total bandwidth B=2MHz; The time wide bandwidth of these two kinds of signals is long-pending identical.
2. emulation content
As seen from Figure 2, in the situation that speed is larger, it is very fast that Doppler's tolerance curve of the designed signal of prior art declines, and the decline of Doppler's tolerance curve of signal of the present invention is slower, therefore Doppler's tolerance of signal of the present invention is better.And wide length during the pulse of the wide signal more designed than prior art during due to signal pulse of the present invention, this more can prove the advantage of signal of the present invention on Doppler's tolerance.
Emulation 2, carries out emulation to the emitted energy function of the emitted energy function of signal of the present invention, the designed signal of prior art and desirable emitted energy function, and result is as Fig. 3.
As seen from Figure 3, the secondary lobe amplitude of the emitted energy function of the emitted energy function ratio signal of the present invention of the designed signal of prior art more approaches desirable emitted energy function, therefore signal of the present invention effect aspect emitted energy function is slightly poor, but still can meet the demands.
Emulation 3, the signal to the designed signal of prior art after synthetic impulse is processed carries out emulation, and result is as Fig. 4.
Emulation 4, the signal to signal of the present invention after synthetic impulse is processed carries out emulation, and result is as Fig. 5.
By Fig. 4, Fig. 5, can be found out, the distance side lobe amplitude of the signal of the distance side lobe Amplitude Ratio prior art design of the designed signal of the present invention after synthetic impulse is processed after synthetic impulse is processed is slightly high, but still can meet the demands.
As seen from Figure 6, the transmit direction secondary lobe amplitude of the transmit direction secondary lobe amplitude of signal of the present invention and the designed signal of prior art is close.
Claims (5)
1. a method for designing for MIMO radar waveform, comprises the steps:
(1) the transmitting element number of array of establishing MIMO radar antenna is N, and carrier frequency is f
0, the pulse width of single waveform is Te, the total bandwidth of N waveform is B, and the bandwidth B s of each single waveform is identical, a fix N waveform centre frequency f
korder, i.e. f
1<f
2< ... <f
n-1<f
n, and by the first phase of N waveform
all be made as 0, k=1 wherein, 2,3 ..., N, N>=8;
(2) define each waveform frequency and be spaced apart Δ f
m=f
m+1-f
m, establish each waveform frequency interval delta f
mbe equal to frequency interval initializaing variable Δ f
0, and between 0 to 1/Te to frequency interval initializaing variable Δ f
0search for, find out the frequency interval initial value Δ f that three dB bandwidth that the three dB bandwidth that meets emitted energy function equals desirable emitted energy function requires
0', m=1 wherein, 2,3 ..., N-1;
(3) at frequency interval initial value Δ f
0near N-1 value of ' random generation, is assigned to respectively waveform frequency interval delta f
m, the centre frequency that obtains the 1st signal is f
1=f
0-sum (Δ f
mthe centre frequency f of i the signal in)/2, the
ifor f
i=f
i-1+ Δ f
m, i=2 wherein, 3 ..., N and m=i-1, sum () is summation operation;
(4) according to frequency interval Δ f
m, first phase
calculate bandwidth B s=B-sum (the Δ f of individual signals
m), then calculate chirp rate μ by μ=Bs/Te, thus baseband signal waveform s obtained
kfor:
Wherein, j is unit imaginary constant, and t is the L point sampling time in 0～Te, and exp () is for take the exponential function that natural constant e is the end;
(5) according to baseband waveform signal s
k, obtain signal covariance matrix R=ss
h, signal waveform matrix s=[s wherein
1; s
2; s
k; s
n], and then obtain emitted energy function and be:
P(θ)=a(θ)
TRa(θ)
*，
Wherein θ is the angle of-90～90 degree scopes, a (θ)=exp (j2 π ((1:N)
t-(N+1)/2) be * d.*G/c*sin (θ)) steering vector of angle θ, d is array element radius, and c is the light velocity, centre frequency vector G=[f
1, f
2..., f
k..., f
n]
t, ()
*for conjugate operation, ()
tfor transposition computing .* is point multiplication operation symbol;
(6) establishing transmit beam main lobe direction is β, according to signal covariance matrix R, obtains transmit direction function:
f(θ,β)=a(θ)
TRa(β)
*，
A (β)=exp (j2 π ((1:N) wherein
t-(N+1)/2) be * d.*G/c*sin (β)) steering vector of transmit beam main lobe direction;
(7) angle of establishing any one target in the angular range of desirable emitted energy function three dB bandwidth is φ, by signal waveform matrix s, obtains echoed signal sr=a (φ)
ts, carries out synthetic impulse to this echoed signal, and the signal st (φ) obtaining after synthetic impulse is processed is:
st(φ)=xcorr(a(φ)
Ts)，
A (φ)=exp (j2 π ((1:N) wherein
t-(N+1)/2) be * d.*G/c*sin (φ)) steering vector of target direction, xcorr () is auto-correlation computation;
(8) according to signal st (φ), transmit direction function f (θ after synthetic impulse processing, β), emitted energy function P (θ), calculate the maximum secondary lobe amplitude y2 of signal after synthetic impulse, the maximum difference p of the maximum side lobe gain y3 of transmit direction function and emitted energy function P (θ) and desirable emitted energy function Pt (θ), and then obtain row vector x:
x=[a*y2,b*y3,c*p]，
Wherein, a, b, c are greater than three different coefficients of 0, meet the different demands of waveform performance by adjusting its size;
(9) with frequency interval Δ f
m, first phase
as input, take row vector x as output, obtain cost function:
(10) with each waveform frequency interval delta f
m, first phase
for variable, call cost function
utilize Chebyshev approximation to be optimized, each waveform frequency interval delta f after being optimized
m' and first phase
make it meet low distance side lobe, low angle secondary lobe, and waveform emitted energy function approach the demand of desirable emitted energy function;
(11) with each waveform frequency interval delta f after optimizing
m' and first phase
try to achieve the final center frequency value f of each waveform
k', the final bandwidth value Bs ' of single waveform and final chirp rate value μ ', and then obtain final each signal waveform s
k':
2. the method for designing of MIMO radar waveform according to claim 1, wherein step (2) described between 0 to 1Te to frequency interval initializaing variable Δ f
0search for, find out the frequency interval initial value Δ f that three dB bandwidth that the three dB bandwidth that meets emitted energy function equals desirable emitted energy function requires
0', carry out as follows:
2a) discrete value of getting uniformly Q point between 0 to 1/Te, and by its ascending [h that is arranged as
1', h
2' ..., h
n' ..., h
q'], n=1 wherein, 2 ..., Q;
2b) make frequency interval initializaing variable Δ f
0equal n discrete value h
n', according to frequency interval initializaing variable Δ f
0, by formula Bs
0=B-(N-1) Δ f
0calculate the initial bandwidth Bs of individual signals
0, and then obtain corresponding initial chirp rate: μ
0=Bs
0/ Te;
2c) according to frequency interval initializaing variable Δ f
0, calculate the initial center frequency η of each signal
k: η
k=f
0-(N-1) Δ f
0/ 2+ (k-1) Δ f
0thereby obtaining corresponding initial baseband signal waveform is g
k:
2d) according to each initial baseband signal waveform g
k, obtain corresponding initialize signal covariance matrix R
0=gg
h, initial matrix g=[g wherein
1; g
2; g
k; g
n], and then obtain initial transmissions energy function and be:
P
0(θ)=a′(θ)
TR
0a′(θ)
*，
Wherein a ' (θ)=exp (j2 π ((1:N)
t-(N+1)/2) be * d.*G '/c*sin (θ)) the initial steering vector of angle θ, initial center frequency vector G '=[η
1, η
2..., η
k..., η
n]
t;
2e) find out initial transmissions energy function P
0(θ) maximal value Pm=max (abs (P
0(θ))), then obtain equation abs (P
0(θ) separate θ for two of)=0.5*Pm
1, θ
2, θ wherein
1< θ
2, abs () is for asking modular arithmetic;
2f) by two, separate θ
1, θ
2, calculate the three dB bandwidth θ that tries to achieve n initial transmissions energy function
w(n)=θ
2-θ
1;
2g) compare each frequency interval initializaing variable Δ f
0the three dB bandwidth θ of corresponding initial transmissions energy function
w(n) with the size of the three dB bandwidth of desirable emitted energy function, the corresponding frequency interval initializaing variable of the emitted energy function three dB bandwidth Δ f of desirable emitted energy function three dB bandwidth will be approached most
0value, as needed frequency interval initial value Δ f
0'.
3. the method for designing of MIMO radar waveform according to claim 1, calculates the maximum secondary lobe amplitude y2 of signal after synthetic impulse in wherein said step (8), carry out as follows:
8a) discretely in the angular range of desirable emitted energy function three dB bandwidth evenly get J angle on target value, by its ascending [φ that is arranged as
1', φ
2' ..., φ
z' ..., φ
j'], z=1 wherein, 2 ..., J, and make angle on target φ=φ
z', the signal after synthetic impulse processing is st (φ
z'), and then obtain the signal matrix after synthetic impulse corresponding to all discrete values of angle on target processed:
St=[st(φ
1′);st(φ
2′);…;st(φ
z′);…;st(φ
J′)]；
8b) the signal matrix St after synthetic impulse is processed presses row normalization, obtains normalized signal matrix Sg after synthetic impulse:
Sg=[abs(st(φ
1′))/max(abs(st(φ
1′)));abs(st(φ
2′))/max(abs(st(φ
2′)));…;abs(st(φ
J′))/max(abs(st(φ
J′)))]，
Wherein max () is maximizing function, and abs () is for asking modular arithmetic;
8c) main lobe of normalized signal matrix S g after synthetic impulse is partly taken as to 0, i.e. Sg (1:J, L+ (3:3))=0, and then obtain maximum secondary lobe amplitude y2 corresponding to all discrete values of angle on target φ:
y2=max(max(Sg))。
4. the method for designing of MIMO radar waveform according to claim 1, calculates the maximum side lobe gain y3 of transmit direction function in wherein said step (8), carry out as follows:
8.1) discretely in the angular range of desirable emitted energy function three dB bandwidth get M angle value, and by its ascending [β that is arranged as
1', β
2' ..., β
q' ..., β
m'], q=1 wherein, 2 ..., M, and make transmit beam main lobe direction β=β
q', transmit direction function f (θ, β
q')=a (θ)
tra (β
q')
*, obtain transmit direction Jacobian matrix corresponding to all discrete angle values of transmit beam main lobe direction β:
F(θ,β)=[abs(f(θ,β
1′)),abs(f(θ,β
2′)),…,abs(f(θ,β
M′))]；
8.2), according to transmit direction Jacobian matrix F (θ, β), calculate normalized transmit direction Jacobian matrix and be:
Fg=[abs(f(θ,β
1'))/Fa(1),abs(f(θ,β
2'))/Fa(2),…,abs(f(θ,β
M'))/Fa(M)]，
Wherein, Fa=[max (abs (f (θ, β
1'))), max (abs (f (θ, β
2'))) ..., max (abs (f (θ, β
m')))] be the vector of the maximal value composition of each row of transmit direction Jacobian matrix F (θ, β);
8.3) main lobe in transmit direction Jacobian matrix Fg is partly taken as to zero, is about to meet in matrix F g q row | θ-β
q| the element of <3 ° is taken as 0, and then obtains the maximum side lobe gain of the transmit direction function of all angle betas:
y3=max(max(Fg))。
5. the method for designing of MIMO radar waveform according to claim 1, in wherein said step (8), calculate the maximum difference p of emitted energy function P (θ) and desirable emitted energy function Pt (θ), be first according to emitted energy function P (θ), calculate the maximal value of emitted energy function: Pa=max (abs (P (θ))); Calculate again the maximum difference p of emitted energy function P (θ) and desirable emitted energy function Pt (θ):
p=max(abs(abs(P(θ))/Pa-abs(Pt(θ)))),
Wherein, abs () is for asking modular arithmetic.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310499194.XA CN103592642B (en) | 2013-10-22 | 2013-10-22 | The method for designing of MIMO radar waveform |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310499194.XA CN103592642B (en) | 2013-10-22 | 2013-10-22 | The method for designing of MIMO radar waveform |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103592642A true CN103592642A (en) | 2014-02-19 |
| CN103592642B CN103592642B (en) | 2015-10-28 |
Family
ID=50082853
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310499194.XA Active CN103592642B (en) | 2013-10-22 | 2013-10-22 | The method for designing of MIMO radar waveform |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103592642B (en) |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104297732A (en) * | 2014-06-26 | 2015-01-21 | 电子科技大学 | Method and device for producing MIMO radar signal based on frequency-hopping technique |
| CN104569924A (en) * | 2014-12-12 | 2015-04-29 | 中国电子科技集团公司电子科学研究院 | Fuzzy-free region area expansion oriented incoherent MIMO radar waveform design method |
| CN105044683A (en) * | 2015-08-06 | 2015-11-11 | 西安电子科技大学 | Emission waveform design method of conformal MIMO radar system under multi-path condition |
| CN105068052A (en) * | 2015-08-19 | 2015-11-18 | 西安电子科技大学 | Design method for zero correlation zone like quadrature phase coding signals |
| CN105301572A (en) * | 2015-09-17 | 2016-02-03 | 西安电子科技大学 | Multi-peak low sidelobe phase encode signal and mismatched filter joint design method |
| CN105487068A (en) * | 2014-10-07 | 2016-04-13 | 罗伯特·博世有限公司 | Method and mimo radar device for determining a position angle of an object |
| CN106019237A (en) * | 2016-06-23 | 2016-10-12 | 哈尔滨工业大学（威海） | Radar LFM composite waveform design method |
| CN106597386A (en) * | 2016-08-01 | 2017-04-26 | 哈尔滨工业大学（威海） | Orthogonal coding waveform with discrete frequency FM gradient and design method thereof |
| CN108226878A (en) * | 2017-12-08 | 2018-06-29 | 西安电子科技大学 | MIMO radar transmitted waveform synthetic method based on coordinate descent algorithm |
| CN108490411A (en) * | 2018-03-13 | 2018-09-04 | 桂林电子科技大学 | A kind of method of color clutter in inhibition target detection |
| CN108885255A (en) * | 2016-04-15 | 2018-11-23 | 联发科技股份有限公司 | Radar chaff mitigates method and apparatus |
| CN109856603A (en) * | 2019-03-26 | 2019-06-07 | 森思泰克河北科技有限公司 | Radar anti-interference method, terminal device and storage medium |
| CN109946655A (en) * | 2019-03-18 | 2019-06-28 | 西安电子科技大学 | MIMO radar LFM signal zero setting waveform design method |
| CN110109065A (en) * | 2019-05-09 | 2019-08-09 | 西安电子科技大学 | MIMO radar part waveform correlation design method based on LFM signal |
| CN110398732A (en) * | 2019-06-21 | 2019-11-01 | 西北大学 | The target direction detection method of low calculation amount adaptive step iterative search |
| CN110412514A (en) * | 2019-08-02 | 2019-11-05 | 南京慧尔视智能科技有限公司 | A kind of linear frequency modulation continuous wave waveform optimization method under MIMO system |
| CN110632559A (en) * | 2019-09-17 | 2019-12-31 | 南京航空航天大学 | MIMO pulse radar waveform design and optimization method based on application characteristics |
| CN111478749A (en) * | 2020-02-16 | 2020-07-31 | 西安电子科技大学 | MIMO iterative detection method, system and application based on optimization initial value fast convergence |
| CN111693983A (en) * | 2020-06-12 | 2020-09-22 | 哈尔滨工业大学（威海） | Complementary waveform construction method and module, cognitive radar system and waveform emission method |
| CN111812612A (en) * | 2020-06-18 | 2020-10-23 | 西安电子科技大学 | Partial correlation waveform design method of MIMO radar based on subarray orthogonal LFM signal |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107329120B (en) * | 2017-06-29 | 2019-08-20 | 中国人民解放军信息工程大学 | The MIMO radar waveform design method differentiated towards approaching target |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6204797B1 (en) * | 1999-06-11 | 2001-03-20 | The United States Of America As Represented By The Secretary Of The Air Force | Waveform diversity for spatial-temporal denial of radar and communication system |
| CN101950014A (en) * | 2010-09-02 | 2011-01-19 | 西安电子科技大学 | MIMO (Multiple Input Multiple Output) radar wave shape online designing method |
| CN103018721A (en) * | 2012-12-11 | 2013-04-03 | 电子科技大学 | Method for generating Multiple Input Multiple Output (MIMO)-over the horizon (OTH) radar waveform |
-
2013
- 2013-10-22 CN CN201310499194.XA patent/CN103592642B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6204797B1 (en) * | 1999-06-11 | 2001-03-20 | The United States Of America As Represented By The Secretary Of The Air Force | Waveform diversity for spatial-temporal denial of radar and communication system |
| CN101950014A (en) * | 2010-09-02 | 2011-01-19 | 西安电子科技大学 | MIMO (Multiple Input Multiple Output) radar wave shape online designing method |
| CN103018721A (en) * | 2012-12-11 | 2013-04-03 | 电子科技大学 | Method for generating Multiple Input Multiple Output (MIMO)-over the horizon (OTH) radar waveform |
Non-Patent Citations (1)
| Title |
|---|
| 胡亮兵等: "MIMO雷达发射方向图匹配和波形优化方法", 《西安电子科技大学学报》, vol. 36, no. 06, 31 December 2009 (2009-12-31), pages 1021 - 1026 * |
Cited By (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104297732A (en) * | 2014-06-26 | 2015-01-21 | 电子科技大学 | Method and device for producing MIMO radar signal based on frequency-hopping technique |
| CN105487068A (en) * | 2014-10-07 | 2016-04-13 | 罗伯特·博世有限公司 | Method and mimo radar device for determining a position angle of an object |
| CN104569924A (en) * | 2014-12-12 | 2015-04-29 | 中国电子科技集团公司电子科学研究院 | Fuzzy-free region area expansion oriented incoherent MIMO radar waveform design method |
| CN104569924B (en) * | 2014-12-12 | 2017-01-25 | 中国电子科技集团公司电子科学研究院 | Fuzzy-free region area expansion oriented incoherent MIMO radar waveform design method |
| CN105044683A (en) * | 2015-08-06 | 2015-11-11 | 西安电子科技大学 | Emission waveform design method of conformal MIMO radar system under multi-path condition |
| CN105044683B (en) * | 2015-08-06 | 2017-05-24 | 西安电子科技大学 | Emission waveform design method of conformal MIMO radar system under multi-path condition |
| CN105068052A (en) * | 2015-08-19 | 2015-11-18 | 西安电子科技大学 | Design method for zero correlation zone like quadrature phase coding signals |
| CN105301572A (en) * | 2015-09-17 | 2016-02-03 | 西安电子科技大学 | Multi-peak low sidelobe phase encode signal and mismatched filter joint design method |
| CN105301572B (en) * | 2015-09-17 | 2017-11-10 | 西安电子科技大学 | Multi-peak low sidelobe phase-coded signal and mismatched filter co-design method |
| CN108885255A (en) * | 2016-04-15 | 2018-11-23 | 联发科技股份有限公司 | Radar chaff mitigates method and apparatus |
| CN106019237A (en) * | 2016-06-23 | 2016-10-12 | 哈尔滨工业大学（威海） | Radar LFM composite waveform design method |
| CN106019237B (en) * | 2016-06-23 | 2021-07-30 | 哈尔滨工业大学（威海） | Radar LFM composite waveform design method |
| CN106597386A (en) * | 2016-08-01 | 2017-04-26 | 哈尔滨工业大学（威海） | Orthogonal coding waveform with discrete frequency FM gradient and design method thereof |
| CN106597386B (en) * | 2016-08-01 | 2019-03-12 | 哈尔滨工业大学（威海） | A kind of the orthogonal coding waveform and its design method of discrete frequency chirp rate |
| CN108226878A (en) * | 2017-12-08 | 2018-06-29 | 西安电子科技大学 | MIMO radar transmitted waveform synthetic method based on coordinate descent algorithm |
| CN108226878B (en) * | 2017-12-08 | 2020-02-07 | 西安电子科技大学 | MIMO radar transmitting waveform synthesis method based on coordinate descent algorithm |
| CN108490411A (en) * | 2018-03-13 | 2018-09-04 | 桂林电子科技大学 | A kind of method of color clutter in inhibition target detection |
| CN108490411B (en) * | 2018-03-13 | 2021-08-03 | 桂林电子科技大学 | Method for inhibiting color clutter in target detection |
| CN109946655A (en) * | 2019-03-18 | 2019-06-28 | 西安电子科技大学 | MIMO radar LFM signal zero setting waveform design method |
| CN109946655B (en) * | 2019-03-18 | 2022-12-02 | 西安电子科技大学 | MIMO radar LFM signal zero setting waveform design method |
| CN109856603A (en) * | 2019-03-26 | 2019-06-07 | 森思泰克河北科技有限公司 | Radar anti-interference method, terminal device and storage medium |
| CN110109065B (en) * | 2019-05-09 | 2022-12-06 | 西安电子科技大学 | Partial correlation waveform design method of MIMO radar based on LFM signal |
| CN110109065A (en) * | 2019-05-09 | 2019-08-09 | 西安电子科技大学 | MIMO radar part waveform correlation design method based on LFM signal |
| CN110398732A (en) * | 2019-06-21 | 2019-11-01 | 西北大学 | The target direction detection method of low calculation amount adaptive step iterative search |
| CN110412514A (en) * | 2019-08-02 | 2019-11-05 | 南京慧尔视智能科技有限公司 | A kind of linear frequency modulation continuous wave waveform optimization method under MIMO system |
| CN110412514B (en) * | 2019-08-02 | 2023-05-23 | 南京慧尔视智能科技有限公司 | Linear frequency modulation continuous wave waveform optimization method under MIMO system |
| CN110632559A (en) * | 2019-09-17 | 2019-12-31 | 南京航空航天大学 | MIMO pulse radar waveform design and optimization method based on application characteristics |
| CN110632559B (en) * | 2019-09-17 | 2023-01-31 | 南京航空航天大学 | MIMO pulse radar waveform design and optimization method based on application characteristics |
| CN111478749A (en) * | 2020-02-16 | 2020-07-31 | 西安电子科技大学 | MIMO iterative detection method, system and application based on optimization initial value fast convergence |
| CN111693983A (en) * | 2020-06-12 | 2020-09-22 | 哈尔滨工业大学（威海） | Complementary waveform construction method and module, cognitive radar system and waveform emission method |
| CN111693983B (en) * | 2020-06-12 | 2023-04-28 | 哈尔滨工业大学（威海） | Complementary waveform construction method and module, cognitive radar system and waveform emission method |
| CN111812612A (en) * | 2020-06-18 | 2020-10-23 | 西安电子科技大学 | Partial correlation waveform design method of MIMO radar based on subarray orthogonal LFM signal |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103592642B (en) | 2015-10-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103592642B (en) | The method for designing of MIMO radar waveform | |
| CN103901417B (en) | Low-complexity space target two-dimensional angle estimation method of L-shaped array MIMO radar | |
| CN101369014B (en) | Bilateral constraint self-adapting beam forming method used for MIMO radar | |
| Huang et al. | Frequency-domain AoA estimation and beamforming with wideband hybrid arrays | |
| CN105912791B (en) | DOA estimation method based on local search in virtual relatively prime array | |
| CN105319545B (en) | Improve the MIMO-OFDM radar waveform design methods of STAP detection performances | |
| CN110927661A (en) | Single-basis expansion co-prime array MIMO radar DOA estimation method based on MUSIC algorithm | |
| CN109471064B (en) | Time modulation array direction-finding system based on pulse compression technology | |
| CN105467365A (en) | A low-sidelobe emission directional diagram design method improving DOA estimated performance of a MIMO radar | |
| CN107966677B (en) | Circular array modal domain orientation estimation method based on space sparse constraint | |
| CN103969633A (en) | Method for hierarchically designing transmission waveforms of MIMI (multiple input multiple output) radar for detecting targets in clutter | |
| CN103018732B (en) | MIMO (multi-input multi-output) radar waveform synthesis method based on space-time joint optimization | |
| CN103760527B (en) | Single base MIMO radar coherent source Wave arrival direction estimating method | |
| CN110109065B (en) | Partial correlation waveform design method of MIMO radar based on LFM signal | |
| CN108710112A (en) | Space and time code array space smoothing decorrelation LMS Power estimation method | |
| CN102841335A (en) | Iterative FFT-based quick MIMO radar waveform synthesis method | |
| CN109946655A (en) | MIMO radar LFM signal zero setting waveform design method | |
| CN114114188A (en) | FDA radar communication integrated waveform design method with low side lobe | |
| JP7150641B2 (en) | antenna system | |
| CN108896974A (en) | A kind of improved MIMO array high-resolution Estimation of Spatial Spectrum method | |
| CN104868946A (en) | Adaptive weighted interference suppression method of subarray level mixed MIMO-phased array system | |
| CN113281732A (en) | MIMO radar target positioning method and system based on space-time coding | |
| CN113238212A (en) | Space-time coding-based frequency diversity array radar range resolution enhancement method | |
| CN104459680B (en) | Method for rapidly estimating target direction through MIMO radar | |
| US10879609B1 (en) | Wave construction method for controlling, rotating, or shaping radio frequency or acoustic waves in free space or in a fluid |
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 | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |