CN101902432B - Method for establishing and optimizing fully-complementary sequence in orthogonal MIMO radar system - Google Patents
Method for establishing and optimizing fully-complementary sequence in orthogonal MIMO radar system Download PDFInfo
- Publication number
- CN101902432B CN101902432B CN 201010237523 CN201010237523A CN101902432B CN 101902432 B CN101902432 B CN 101902432B CN 201010237523 CN201010237523 CN 201010237523 CN 201010237523 A CN201010237523 A CN 201010237523A CN 101902432 B CN101902432 B CN 101902432B
- Authority
- CN
- China
- Prior art keywords
- fully
- complementary sequence
- mimo radar
- complementary
- matrix
- 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.)
- Expired - Fee Related
Links
Images
Abstract
The invention discloses a method for establishing and optimizing a fully-complementary sequence in an orthogonal MIMO radar system, which is used in the orthogonal MIMO radar system. The method comprises the following steps of: firstly, setting an iteration number and an initial complementary sequence, and establishing the fully-complementary sequence by adopting an iteration method; secondly, selecting a suitable full-complementary sequence as a transmission signal according to antenna data in the orthogonal MIMO radar system; and finally, optimizing the selected full-complementary sequence according to clutter characteristics in the orthogonal MIMO radar system, and using the optimized full-complementary sequence as the final transmission signal. The method adopts the full-complementary sequence to construct orthogonal MIMO radar transmission signals, and finds the transmission signal capable of making mutual information amount be the largest through optimization so as to improve the mutual information amount.
Description
Technical field
The invention belongs to the Radar Technology field, the structure and the optimization method of fully-complementary sequence in particularly a kind of orthogonal MIMO radar.
Background technology
Multiple-input and multiple-output (MIMO; Multiple-Input Multiple-Output) technology is at first invented by the AT&T Labs the nineties in 20th century at wireless communication field, and mimo wireless communication is divided into the technology of utilizing Space Time Coding to improve signal to noise ratio and sends two kinds of main modes of technology that different information improve transmission rate simultaneously with many antennas.Along with the research of MIMO communication system, people have proposed the notion of MIMO radar again.
The notion of orthogonal MIMO radar is proposed by the Lincoln laboratory at first; The orthogonal MIMO radar is at the mutually orthogonal waveform of transmitting terminal emission; Compare with traditional phased array radar, have good advantages at aspects such as wide search wave beam formation, low probability of intercept (LPI), clutter inhibition.
The orthogonal MIMO radar has above-mentioned many advantages; Since coming out, just received the common concern of Chinese scholars; For this New System radar of orthogonal MIMO radar; Because the mutually orthogonal waveform collection of each transmission antennas transmit, receiving terminal recovers each component that transmits through matched filter processing, and the design that therefore transmits has directly influenced the systematic function of MIMO radar.In order to suppress to disturb and improve multiple target resolution, have good correlation function between requiring to transmit, preferably satisfy complete quadrature, promptly the non-periodic autocorrelation function secondary lobe is that zero-sum cross-correlation function aperiodic main lobe and secondary lobe are zero.The MIMO radar mainly adopts quadrature polyphase code and orthogonal frequency coding at present; Though the correlation function of above-mentioned two types of codings has lower side lobe performance; But still can not satisfy the complete quadrature between transmitting; Theoretical research shows that it is non-existent satisfying completely orthogonal sequence in traditional solid size field.Because fully-complementary sequence is made up of a plurality of subsequences; And between the subsequence relation of quadrature; Just in time coincide, thereby a new research direction has been opened up in the selection that appears as the orthogonal MIMO radar signal of fully-complementary sequence with the many signals and the orthogonality of the requirement of orthogonal MIMO radar.
Mutual information is important in an orthogonal MIMO theory system index, and the mutual information of system is big more, and the useful information amount of from system, extracting is also just many more.Because the complex characteristics of orthogonal MIMO Radar channel; Different clutter types can change the size of mutual information; In order to improve mutual information, we can improve the characteristic of channel from the angle of signal, promptly are optimized processing to transmitting to obtain optimum systematic function.
Summary of the invention
The invention provides fully-complementary sequence structure and optimization method in a kind of orthogonal MIMO radar system; Purpose is adopting fully-complementary sequence to construct aspect the orthogonal MIMO radar emission signal; And the fully-complementary sequence of being constructed has been realized the raising of channel capacity through optimizing.
The structure of fully-complementary sequence and optimization method in a kind of orthogonal MIMO radar system specifically may further comprise the steps:
Step 1, according to the demand of actual orthogonal MIMO radar, be provided for the iterations r that fully-complementary sequence makes up, require constructed complementary series number 2
R+1Be greater than or equal number of transmit antennas in the orthogonal MIMO radar system; Wherein, r is the integer greater than 0;
The number of phases of step 2, setting fully-complementary sequence is provided with initial complementary series { A
0, B
0, A in the set initial complementary series
0With B
0Complementary;
Step 3, employing iterative method make up fully-complementary sequence;
Described iterative method is the process of a recurrence, by { A
0, B
0Can construct two pairs of fully-complementary sequences: { A
0B
0A
0-B
0And
With
Again will be wherein each fully-complementary sequence proceed iteration and go down iteration r time;
Step 4, according to antenna data in the orthogonal MIMO radar system, select fully-complementary sequence; The fully-complementary sequence number of said selection is identical with number of transmit antennas;
According to following channel capacity formula fully-complementary sequence is optimized:
Wherein, I (h; Y) expression mutual information, h representes channel matrix,
y
1... y
NBe respectively the 1st to N the echo that receiver receives, N is the reception antenna number, and T representes the transposition computing; E is a unit matrix; Subscript H representes conjugate transpose; P is the distribution matrix of noise signal v; L representes the code element number of sequence, equals the sub-sequence length of complementary series;
The variance yields of expression channel; M is the number of transmit antennas of orthogonal MIMO radar;
The matrix S that is made up of fully-complementary sequence is:
When noise signal v is relevant, according to (ψ S)
HP
-1(ψ S)=S
H(ψ
HP
-1ψ) S=Λ, the ψ S that obtains are the fully-complementary sequence matrixes after optimizing;
That wherein, hints obliquely at matrix ψ asks method following:
Because matrix P is a positive definite matrix, can decompose as follows clutter matrix P:
Then:
That is:
When noise signal v is white Gaussian noise,
arranged according to
then S be exactly the fully-complementary sequence matrix after optimizing; Wherein
representes the variance yields of noise signal, and Λ is a diagonal matrix; With the fully-complementary sequence matrix S after optimizing as transmitting in the orthogonal MIMO radar.
Advantage of the present invention and good effect are:
(1) to utilize alternative manner to construct number abundant in the present invention, and the sufficiently long fully-complementary sequence of code length can satisfy the application of actual orthogonal MIMO radar system fully;
(2) different distributions of pair of orthogonal MIMO radar clutter of the present invention has been carried out optimum processing to transmitting, and the Optimal Signals that obtains has enlarged mutual information, has improved systematic function.
Description of drawings
Fig. 1 is the flow chart of steps of fully-complementary sequence structure of the present invention and optimization method;
Fig. 2 is the sketch map that fully-complementary sequence makes up in the step 2 of the present invention;
Fig. 3 is that the present invention optimizes the back signal and the former mutual information that transmits compares sketch map.
Embodiment
To combine accompanying drawing and embodiment that the present invention is done further detailed description below.
Fully-complementary sequence structure and optimization method are realized through following steps in the orthogonal MIMO radar system provided by the invention, and be as shown in Figure 1:
Step 1, iterations r is set.Choosing iterations in the embodiment of the invention is 3.According to the demand of actual orthogonal MIMO radar, sequence length can not be oversize can not be too short, oversize complexity height a bit, having lacked very much systematic function can be a little bit poorer, general iterations is between 3~5.Require constructed complementary series number 2
R+1Be greater than or equal number of transmit antennas in the orthogonal MIMO radar system.
The number of phases of step 2, setting fully-complementary sequence is provided with initial complementary series { A
0, B
0, A in the selected initial complementary series
0With B
0Complementary.The number of phases of choosing fully-complementary sequence in the embodiment of the invention is the structure that 5 initial complementary series is explained fully-complementary sequence: A
0=(1i-i-1i), B
0=(111i-i).In practical application, the number of phases of fully-complementary sequence generally is set at 2 or 4.The number of phases of fully-complementary sequence is the same with the initial complementary series number of phases.
Initial complementary series according to selected can obtain: the length L of initial complementary series
0=5; Initial complementary series number of phases P_n=4; The length L that needs the fully-complementary sequence of structure
r=L
0* 2
r=40, this length is meant the sub-sequence length of each complementary series.
Step 3, employing iterative method make up fully-complementary sequence.
Described iterative method is the process of a recurrence, by { A
0, B
0Can construct two pairs of fully-complementary sequences: { A
0B
0A
0-B
0And
With
Wherein, { A
0B
0A
0-B
0With
It is complementary fully,
With
Complementary fully.In like manner, again will be wherein each fully-complementary sequence proceed iteration and go down, iteration r time, the number of the fully-complementary sequence that constructs are 2
r, the number of complementary series is 2
R+1, the length L of fully-complementary sequence
r=L
0* 2
rIf initial complementary series number of phases P_n>2 then need conjugation to handle in building process, building process is as shown in Figure 2, wherein
Expression B
0The opposite sequence of sequence,
Expression B
0Conjugation.
As shown in Figure 2, iteration once, by { A
0, B
0Construct { A
1, B
1, { A
1', B '
1, { (AA)
1, (BB)
1, { (AA)
1', (AA) '
1, wherein, { A
1, B
1Be { A
0B
0A
0-B
0, { A
1', B '
1Do
{ (AA)
1, (BB)
1Do
{ (AA)
1', (BB) '
1Do
Iteration for the second time, with { A
1, B
1Be example, construct { A
2, B
2, { A
2', B '
2, { (AA)
2, (BB)
2, { (AA)
2', (BB) '
2, wherein, { A
2, B
2Be { A
1B
1A
1-B
1, { A
2', B '
2Do
{ (AA)
2, (BB)
2Do
{ (AA)
2', (BB) '
2Do
Iteration for the third time, with { A
2, B
2Be example, construct { A
3, B
3, { A
3', B '
3, { (AA)
3, (BB)
3, { (AA)
3', (BB) '
3, wherein, { A
3, B
3Be { A
2B
2A
2-B
2, { A
3', B '
3Do
{ (AA)
3, (BB)
3Do
{ (AA)
3', (BB) '
3Do
In like manner, can carry out r time iteration.
The constructed fully-complementary sequence that goes out is a quadrature.Specify as follows:
The a pair of fully-complementary sequence that constructs through the r rank is { A
r, B
rAnd A '
r, B '
r, as satisfying
With
Wherein * representes related operation, proves that then sequence is a fully-complementary sequence, is quadrature, forms one type of fully-complementary sequence collection by these sequences of same length.
Iteration is 3 times in the embodiment of the invention, and the fully-complementary sequence collection that constructs comprises 8 pairs of fully-complementary sequences, i.e. 16 complementary seriess, and each complementary series comprises two sub-sequence, and each sub-sequence length is 40, and a pair of fully-complementary sequence wherein is:
Wherein,
Can verify that thus constructed sequence is a fully-complementary sequence.
Step 4, according to antenna data in the orthogonal MIMO radar system, select suitable fully-complementary sequence.According to the number of transmitting antenna, require to construct complementary series number 2
R+1Be greater than or equal number of transmit antennas, from 2
R+1Choose complementary series number in the middle of the individual complementary series with the transmitting antenna similar number.
Number of transmit antennas is 2 in the embodiment of the invention, 2 of reception antenna data bit.Choose following fully-complementary sequence as transmitting:
Number according to orthogonal MIMO radar system transmitting antenna and reception antenna; Choose the fully-complementary sequence that makes up in the step 3 suitably; According to the type of clutter in the real system, be foundation with the maximize channel capacity, the fully-complementary sequence of launching is optimized processing.
The clutter distribution matrix P of orthogonal MIMO radar system in the embodiment of the invention
K, q=0.4
| k-q|, k, q be the capable ordinal sum row ordinal number of representing matrix P respectively; Embodiment P is 2L rank matrixes, and then k, the equal value of q are 1 here ..., 2L.I.e.
Because the number of transmit antennas of orthogonal MIMO radar is M, the fully-complementary sequence that then makes up from step 3 is concentrated and is chosen M to complementary series { A
m, B
m(m=1 ..., M) as transmitting S, that is:
Wherein,
{ S
Am, S
BmBe the complementary series of choosing in the step 4,1≤m≤M; T representes the transposition computing, sets reception antenna number N simultaneously, orthogonal MIMO Radar channel parameter h
n=[h
N1..., h
NM]
T, v
nExpression clutter item, then n the echo y that receiver receives
nFor:
y
n=Sh
n+v
n n=1,…,N (2)
The signal that each antenna is received makes up,
simultaneously
X represent the signal combination of all antennas; Channel matrix
then
y=Xh+v (3)
V representes the clutter noise matrix;
formula (3) can get the orthogonal MIMO radar system capacity C according to this formula to be for the orthogonal MIMO radar system emission receives equation:
C=max{I(h;y)}=max{H(y)-H(y|h)}(4)
Wherein, I (h; Y) be mutual information, H (y) is the entropy of matrix y, and H (y|h) is the conditional entropy of y under the given h, because the independence of h and y, then formula (4) can be write as:
I(h;y)=H(y)-H(v|h)=H(y)-H(v)(5)
H (is the entropy of noise matrix v v), supposes that channel matrix h obeys average E (h)=0, variance
Distribution,
The variance yields of expression channel, E
NMExpression E is NM rank matrixes, the variance of definition v
L representes the code element number of sequence, the just sub-sequence length of complementary series; Shannon's theorems according to communication theory can get:
H(v)=lg[det(∑)](7)
The objective of the invention is to find the X that transmits that can make mutual information maximum; Because so S (i.e. the fully-complementary sequence of emission) is determining X, formula (8) can be write as
:
Be S
HP
-1S is a diagonal matrix, at this moment I (h; Y) obtain maximum, this moment, fully-complementary sequence was optimal sequence, and when noise signal v was relevant, we need transform ψ S to the S that transmits, and make (ψ S)
HP
-1(ψ S) is diagonal matrix, improves channel capacity:
(ψS)
HP
-1(ψS)=S
H(ψ
HP
-1ψ)S=Λ(11)
At this moment, ψ S is the sequence after the optimization.
That hints obliquely at matrix ψ asks method following:
Because matrix P is a positive definite matrix, can decompose as follows clutter matrix P:
Then further can obtain formula (13):
Further obtain formula (14) again:
Calculate ψ from formula (14) according to concrete clutter distribution matrix.
Fully-complementary sequence after optimizing is applied in the orthogonal MIMO radar.
Fully-complementary sequence constructed among the present invention is when noise signal v is relevant; Its different-effect that before and after optimizing, is applied in the orthogonal MIMO radar is as shown in Figure 3; Abscissa SCR representes signal to noise ratio; Ordinate is represented mutual information, can see that the fully-complementary sequence that adopts after optimizing than not adopting the fully-complementary sequence of optimizing in application, can improve the mutual information of MIMO radar.
Claims (3)
1. the structure and the optimization method of fully-complementary sequence in the orthogonal MIMO radar is characterized in that this method may further comprise the steps:
Step 1, according to the demand of actual orthogonal MIMO radar, be provided for the iterations r that fully-complementary sequence makes up, require constructed complementary series number 2
R+1Be greater than or equal number of transmit antennas in the orthogonal MIMO radar system; Wherein, r is the integer greater than 0;
The number of phases of step 2, setting fully-complementary sequence is provided with initial complementary series { A
0, B
0, A in the set initial complementary series
0With B
0Complementary;
Step 3, employing iterative method make up fully-complementary sequence;
Described iterative method is the process of a recurrence, by { A
0, B
0Construct two pairs of fully-complementary sequences: { A
0B
0; A
0-B
0And
With
Again will be wherein each fully-complementary sequence proceed iteration and go down iteration r time; Wherein
Expression B
0The opposite sequence of sequence,
Expression B
0Conjugation;
Step 4, according to antenna data in the orthogonal MIMO radar system, select fully-complementary sequence; The fully-complementary sequence number of said selection is identical with number of transmit antennas;
Step 5, the fully-complementary sequence of step 4 being chosen according to the noise performance of noise signal v in the orthogonal MIMO radar system are optimized;
According to following mutual information formula fully-complementary sequence is optimized:
Wherein, I (h; Y) expression mutual information, h representes channel matrix,
y
1... Y
NBe respectively the 1st to N the echo that receiver receives, T representes the transposition computing; E is a unit matrix; Subscript H representes conjugate transpose; P is the distribution matrix of noise signal v; L representes the sub-sequence length of complementary series;
The variance yields of expression channel; M is the number of transmit antennas of orthogonal MIMO radar;
The matrix of S for constituting by fully-complementary sequence:
Wherein
{ S
Am, S
BmBe the complementary series of choosing in the step 4,1≤m≤M;
When noise signal v is relevant, according to (ψ S)
HP
-1(ψ S)=S
H(ψ
HP
-1ψ) S=Λ asks for and hints obliquely at matrix ψ, and the ψ S that further obtains again is exactly the fully-complementary sequence matrix after optimizing, with ψ S as transmitting in the orthogonal MIMO radar;
When noise signal v is white Gaussian noise, according to
Have
Then S is exactly the fully-complementary sequence matrix after optimizing; Wherein
The variance yields of expression noise signal, Λ is a diagonal matrix; With the fully-complementary sequence matrix S after optimizing as transmitting in the orthogonal MIMO radar.
2. the structure and the optimization method of fully-complementary sequence is characterized in that in a kind of orthogonal MIMO radar according to claim 1, and the iterations described in the step 1 is 3~5.
3. the structure and the optimization method of fully-complementary sequence is characterized in that the number of phases of the fully-complementary sequence described in the step 2 is set at 2 or 4 in a kind of orthogonal MIMO radar according to claim 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010237523 CN101902432B (en) | 2010-07-27 | 2010-07-27 | Method for establishing and optimizing fully-complementary sequence in orthogonal MIMO radar system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010237523 CN101902432B (en) | 2010-07-27 | 2010-07-27 | Method for establishing and optimizing fully-complementary sequence in orthogonal MIMO radar system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101902432A CN101902432A (en) | 2010-12-01 |
CN101902432B true CN101902432B (en) | 2012-12-05 |
Family
ID=43227641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010237523 Expired - Fee Related CN101902432B (en) | 2010-07-27 | 2010-07-27 | Method for establishing and optimizing fully-complementary sequence in orthogonal MIMO radar system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101902432B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104639473B (en) * | 2015-02-01 | 2018-01-19 | 中国传媒大学 | A kind of mimo channel method of estimation based on fully-complementary sequence and compressed sensing |
JP6566396B2 (en) * | 2015-08-06 | 2019-08-28 | パナソニック株式会社 | Radar equipment |
CN106209704B (en) * | 2016-07-11 | 2019-06-11 | 中国传媒大学 | Time domain mimo channel estimation method based on fully-complementary sequence |
CN107677999B (en) * | 2017-09-25 | 2020-09-08 | 西北工业大学 | Sequence set design method for accurately controlling correlation side lobe |
CN110568409B (en) * | 2019-08-09 | 2023-02-03 | 南京航空航天大学 | Subcarrier allocation and waveform joint optimization design method for radar communication integrated system |
CN111817758B (en) | 2020-07-21 | 2022-03-01 | 上海交通大学 | Discrete modulation signal MIMO transmission method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101241180A (en) * | 2008-01-29 | 2008-08-13 | 电子科技大学 | Orthonormal discrete frequency coding design method possessing relative low self correlation performance |
CN101251597A (en) * | 2008-04-08 | 2008-08-27 | 西安电子科技大学 | Method for self-correction of array error of multi-input multi-output radar system |
-
2010
- 2010-07-27 CN CN 201010237523 patent/CN101902432B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101241180A (en) * | 2008-01-29 | 2008-08-13 | 电子科技大学 | Orthonormal discrete frequency coding design method possessing relative low self correlation performance |
CN101251597A (en) * | 2008-04-08 | 2008-08-27 | 西安电子科技大学 | Method for self-correction of array error of multi-input multi-output radar system |
Non-Patent Citations (2)
Title |
---|
"完全互补序列在MIMO雷达中的应用";李树锋 等;《北京航空航天大学学报》;20100531;第36卷(第5期);全文 * |
李树锋 等."完全互补序列在MIMO雷达中的应用".《北京航空航天大学学报》.2010,第36卷(第5期),全文. |
Also Published As
Publication number | Publication date |
---|---|
CN101902432A (en) | 2010-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101902432B (en) | Method for establishing and optimizing fully-complementary sequence in orthogonal MIMO radar system | |
CN101369014B (en) | Bilateral constraint self-adapting beam forming method used for MIMO radar | |
CN102540187B (en) | Orthogonal waveform designing method for formation flying satellites SAR (synthetic aperture radar) | |
CN105891771B (en) | It is a kind of improve estimated accuracy based on continuously distributed angle estimating method and equipment | |
CN102175989B (en) | Method for measuring incoherently distributed signal two-dimensional DOA (direction of arrival) | |
CN107728118B (en) | Low sidelobe transmission beam pattern design method without fitting covariance matrix | |
CN101799535A (en) | Method for estimating target direction by multiple input multiple output (MIMO) radar | |
CN104199029B (en) | Measurement matrix design method for improving target imaging performance of compressed sensing radar | |
CN113179231B (en) | Beam space channel estimation method in millimeter wave large-scale MIMO system | |
CN106772224A (en) | A kind of L-type array estimating two-dimensional direction-of-arrival algorithm of use time frequency analysis | |
CN104977558A (en) | Distributed source center direction-of-arrival estimation method based on Bayesian compressed perception | |
CN106646387A (en) | MIMO radar method capable of resisting active interference based on emission wave beam domain | |
CN111044979A (en) | Blind source separation-based main lobe interference cancellation and target angle estimation method | |
CN110113088B (en) | Intelligent estimation method for wave arrival angle of separated digital-analog hybrid antenna system | |
CN103983952A (en) | Low-complexity receiving and transmitting angle joint estimation method for non-circular signal double-base MIMO radar | |
CN103684700A (en) | 3D (three-dimensional) MU-MIMO (multiple user-multiple input multiple output) precoding method based on orthogonal joint codebook set | |
CN104615854A (en) | Beam broadening and sidelobe suppression method based on sparse constraint | |
CN101241180A (en) | Orthonormal discrete frequency coding design method possessing relative low self correlation performance | |
CN111193679B (en) | Channel estimation method and system based on co-prime array system | |
CN103323810B (en) | L-array azimuthal angle and pitch angle paired signal processing method | |
CN101483280A (en) | Weight solving method for stable wave beam synthesizer | |
CN112965034B (en) | Method for improving Doppler tolerance of slow time phase coded signal of sky-wave radar | |
CN109490846B (en) | Multi-input multi-output radar waveform design method based on space-time joint optimization | |
CN104868946A (en) | Adaptive weighted interference suppression method of subarray level mixed MIMO-phased array system | |
CN105182292A (en) | Multi-waveform phase coding method based on mode search algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121205 Termination date: 20190727 |
|
CF01 | Termination of patent right due to non-payment of annual fee |