Transmit/Receive Spatial Smoothing with Improved Effective Array Aperture for Angle and Mutual Coupling Estimation in Bistatic MIMO Radar

. We proposed a transmit/receive spatial smoothing with improved effective aperture approach for angle and mutual coupling estimation in bistatic MIMO radar. Firstly, the noise in each channel is restrained, by exploiting its independency, in both the spatial domain and temporal domain. Then the augmented transmit and receive spatial smoothing matrices with improved effective aperture are obtained, by exploiting the Vandermonde structure of steering vector with uniform linear array. The DOD and DOA can be estimated by utilizing the unitary ESPRIT algorithm. Finally, the mutual coupling coefficients of both the transmitter and the receivercanbefiguredoutwiththeestimatedanglesofDODandDOA.Numericalexamplesarepresentedtoverifytheeffectiveness oftheproposedmethod.


Introduction
A novel array radar named as multiple-input multiple-output (MIMO) radar has improved the progress of array signal processing [1,2].Because the MIMO radar is more flexible than the phased array radar [3], it has been researched in many fields [4].In general, the MIMO radar can be divided into two classes: statistical MIMO radar and collocated MIMO radar.It is statistical MIMO radar where the transmit and receive antennas are widely distributed [5], which owns the spatial diversity.The transmit and receive antennas are closely located in collocated MIMO radar [6].Because the antennas transmit totally or partially noncoherent waveforms in the transmitter, a virtual array with larger aperture in the receiver can be formed and higher resolution can be obtained for the waveform diversity.In this paper, we focus on the collocated MIMO radar.
The direction of departure (DOD) and direction of arrival (DOA) estimation are one of the most important aspects in bistatic MIMO radar with collocated antennas.And a lot of algorithms have been presented for this issue.In [7], the Capon based algorithm for DOD and DOA estimation in bistatic MIMO radar is presented.An estimation of signal parameters via rotational invariance technique (ESPRIT) method [8] is proposed by exploiting the invariance property of the transmit array and the receive array.The unitary ESPRIT algorithm for joint angles estimation is presented in [9], which has comparable angle estimation performance as ESPRIT with lower computational complexity.The multiple signal classification (MUSIC) is utilized for estimation of DOD and DOA in [10].As the tensor is used widely in multidimensional signal processing, a three-dimensional tensor decomposition method, parallel factor (PARAFAC) analysis, is used to estimate DOD and DOA [11], which has a better performance than the other methods.There are a lot of methods for joint DOD and DOA estimation in bistatic MIMO radar, and we have just presented representative ones.
In recent years, although many institutes and researchers have been studying this novel radar, only a few of institutes have built up physical systems (e.g., the ONERA in France).In [12], it points out that the mutual coupling is a major cause of descending the performance of radar by real data experiment.In order to eliminate the effect of mutual coupling, a MUSIC-Like method is proposed in [13], and the mutual coupling coefficients can be estimated.But this method loses partial effective aperture and it needs lots of snapshots for 2 International Journal of Antennas and Propagation performance guarantee.In [14], an ESPRIT-Like method performs better than MUSIC-Like method.However, in the case of small snapshots, the accuracy of angle estimation using both of the above approaches will degrade remarkably.By exploiting the multidimensional structure of the received data, a three-order tensor is constructed [15], which are DOD, DOA, and temporal dimensions, respectively.And a realvalued subspace approach is proposed; it computes the subspace utilizing the higher order singular value decomposition (HOSVD).Due to use of the forward-backward averaging technique, this approach is suitable for coherent targets and small snapshots.However, the tensor-based real-valued subspace approach just employs partial aperture; it loses a lot of information.When there exists coherent targets, the performance of that approach will degrade, and it cannot deal with more than two coherent targets.Spatial smoothing technique is an effective approach to deal with the situation of small snapshots and coherent targets.In [16], a spatial smoothing with improved aperture (SSIA) method is proposed to estimate DOA for the phased array radar; it improves the effective array aperture twice larger than the conventional spatial smoothing approaches.This technique also can be used to estimate DOD and DOA for bistatic MIMO radar.
In this paper, we proposed a transmit/receive spatial smoothing with improved effective aperture (TRSSIA) method for joint DOD and DOA estimation in bistatic MIMO radar with unknown mutual coupling.Firstly, the white Gaussian noise is restrained using its dependency in both the spatial and temporal domains.Then the TRSSIA is used to construct the transmit spatial smoothing matrix and the receive spatial smoothing matrix.Due to the Vandermonde structure of steering vector with uniform linear arrays (ULA), the transmit and receive augmented spatial smoothing matrices are constructed, and these two matrices can improve the effective aperture two times larger than conventional ones.Thirdly, by using the centro-Hermitian structure of augmented matrices, the real-valued subspace methods (e.g., unitary ESPRIT) can be used to estimate DOD and DOA.Finally, an additional DOD and DOA pairing technique is proposed and the mutual coupling coefficients are estimated.The proposed approach restrains white Gaussian noise and takes full advantage of the received data, so it provides better angle estimation performance.And it can deal with more than two coherent targets.
The remainder of the paper is organized as follows.In Section 2, the bistatic MIMO radar signal model is introduced.The coupling calibration approach is demonstrated in Section 3. In Section 4, simulations are employed to verify the analytical derivations.Finally, Section 5 gives the conclusion.
Notation.(⋅) H , (⋅) T , (⋅) −1 , and (⋅) * denote conjugate transpose, transpose, inverse, and conjugate, respectively.⊗ and ⊙ denote the Kronecker product operation and Khatri-Rao product operation.Toeplitz(c) denotes the symmetric Toeplitz matrix constructed by the vector c. diag(⋅) denotes the diagonalization operation.mat(⋅) denotes the matrixing operation.Re{⋅} and Im{⋅} denote abstracting the real part and imaginary part of complex number, respectively.e min (⋅) is an operator of getting the minimum eigenvector.mean(⋅) is used to compute the mean of numbers.I  denotes the ×M identity matrix, and 0 × is the × zero matrix.A[ : , :] denotes the th to the th rows of A.

Bistatic MIMO Radar Signal Model
Consider a bistatic MIMO radar system equipped with  transmit antennas and  receive antennas, both of which are ULA with half-wavelength spacing antennas.The transmit antennas emit  orthogonal waveforms S = [s 1 , s2 , . . ., s ] T , and SS H / = I  , where  is the number of samples per pulse period.All the targets are modeled as point-scatters in the farfield, and we assume that there exists  targets in the same range bin.Consider the effect of mutual coupling in both the transmitter and receiver; the received data of the th pulse is shown as where X() ∈ C × is the received data during the th pulse period and  is the number of pulses.C  and C  are the mutual coupling matrices of the transmit and receive arrays, respectively, which can be expressed as banded symmetric Toeplitz matrices [17]:  [18,19].After matched filtering, the output data at the receiver can be expressed as where SH /).It assumes that the noise is i.i.d.complex white Gaussian noise, and we can obtain the equation written as I  .Furthermore, according to the properties of Kronecker product [20,21], we can obtain the following equation:

Proposed Algorithm
3.1.Restraining Noise.According to (4), the noise of all  channels are independent.Meanwhile, the noise of each channel is also independent in the temporal domain.So we can restrain the noise both in the spatial and temporal domain.For the ULA, the received data of each channel is expressed as where  , () means the ( − 1) + th row of x(),  , () is the ( − 1) + th row of w(),and     and     are the th element and the th element of a   and a   , respectively.When ), the correlation coefficient can be expressed as where is the correlation coefficient between   1 () and   2 ( + Δ).For the white Gaussian noise, the correlation coefficient of two adjacent noises and the correlation coefficient of signal and noise are both zeros in the temporal domain.And in the spatial domain, we can get the correlation efficient between any channel and the first channel which can be shown as As we known from (4), the noise is independent in spatial domain.So the correlation coefficient of noise between any channel and the first channel equals zero.
Factually, we can only get limited snapshots, so we obtain the asymptotic correlation coefficients: where ε2 (1,1) and ε2 (,) are the correlation coefficients of noise, which nearly equal zeros.
According to (8a)-(8b), we obtain a new vector, which can be written as where (1,1) , ε2 (1,2) , . . ., ε2 (,) ] T .Consider the effect of mutual coupling in both the transmitter and receiver; we construct two selection matrices: Then we use these two selection matrices on the received data; it can get the selected data: where and s Meanwhile, we can get another selected data vector, which can be expressed as where where

Transmit/Receive Spatial
Smoothing with Improved Aperture.After restraining the noise, we obtain the new received data vector r (Δ), r (Δ).In the following, we omit Δ and write the received data as r , r .In [16], a spatial smoothing with improved aperture (SSIA) method with single snapshot is proposed.It is suitable for coherent sources and improves the effective aperture.This approach performs well for DOA estimation in the phased array radar.In this paper, it will prove that this technique can work well in the MIMO radar.
Firstly, we define a  ×  exchange matrix Π  with ones on its antidiagonal and zeros elsewhere.Then, the left-Π-real matrix Q can be expressed as [22] Equation ( 14) is a left-Π-real matrix of odd order.The 2 order one is obtained from Q 2+1 by dropping the center row and center column.
For r and 1 =  sub +   − 1, we define And ( 15) can be expressed as where Λ = diag(ρ  ), Ȧ and Ä are the first  sub rows and the first   rows of A  , respectively.
Then we obtain an augmented matrix: where ).The steering matrix A taug is twice larger than Ȧ ⊙ Ã .Although the noise matrix V is not white Gaussian, its effect on Xtaug is weak, as every element nearly equals zero.See Appendix for a brief illustration of the rotational invariance property of Xtaug .The subspace approaches, such as MUSIC and ESPRIT, can be used to estimate DOD.We note that Xtaug is a centro-Hermitian matrix, it satisfies the following identity: So the real-valued space methods, for example, unitary MUSIC and unitary ESPRIT, are suitable for Xtaug to estimate angles.In this paper, we estimate angles by utilizing the unitary ESPRIT algorithm.Then the complex matrix can be transformed into real-valued matrix as follows: The real-valued signal subspace can be obtained by making SVD on ( Xtaug ).The property of rotational invariance in real-valued subspace is shown as K (1)  taug E taug Υ taug ≈ K (2)  taug E taug , where E taug contains the  dominant left singular vectors of ( Xtaug ).Meanwhile, K (1)  taug and K (2)  taug are the transformed selection matrices; they are both obtained from J (2)  taug in the following way: where J (2)  taug = I 2 ⊗ J 2 , J 2 = [0 ( sub −1)× , I ( sub −1) ].A class of least squares (LS) approaches [20], for example, structured least squares (SLS) and total least squares (TLS), can be used to solve (20) for Υ taug .Ultimately, the angles of DOD can be figured out by {  }  =1 = asin(2atan(eig(Υ taug ))/ ).

Pairing DOD and DOA.
Because the proposed TRSSIA approach tries to make full use of the received data, the angles of DOD and DOA are figured out from transmit augmented matrix and receive augmented matrix, respectively.By exploiting the relationship between the steering vectors a  , a  and the signal subspace E taug , E raug , then {  }  =1 can be paired with {  }  =1 correctly.Note that their relationship can be expressed as where A  , A  are the first 1 and 1 rows of A  and A  and T 1 and T 2 are both  ×  nonsingular matrices.The selection matrices abstract parts of signal subspace are shown as Then the pairing matrix can be written as follows: where P is a  ×  matrix.ȧ  (  ), a  (  ) are the first  sub , 1 rows of a  (  ), respectively.ȧ  (  ), a  (  ) are the first  sub , 1 rows of a  (  ), respectively.If P(, ) is the largest in the th row, then we pair   with   .

Mutual Coupling Estimation.
In order to calibrate the antennas with unknown mutual coupling, the mutual coupling coefficients need to be estimated.
Lemma 1.For any  × 1 complex vector x and any  ×  banded complex symmetric Toeplitz matrix A, we have (Lemma 3 in [17]) where the  × 1 vector a is given by and  is the highest superdiagonal that is different from zero.
The  ×  matrix Q(x) is given by the sum of the two  ×  following matrices: By Lemma 1, we can get The steering matrix and the signal subspace own the relationship as follows: where T 3 is a  ×  nonsingular matrix.In theory, and it also can be written as , where i  is a  sub ×1 vector, whose all elements are ones.U ⊥ 3 is the orthogonal complement space of U 3 .Then the mutual coupling coefficients in transmitter can be obtained as where . By exploiting the structure of c taug , we get a matrix C = mat(c taug ) = c  i T  .Then the estimated mutual coupling coefficient can be obtained with the operations: The mutual coupling coefficients in receiver also can be figured out in the same way.By exploiting the relationship between steering matrix and signal subspace ( Ȧ ⊙ [C  A  ])T 4 = J (4)  taug (Q 2 sub  E taug ) = U 4 , where J (4)  taug = [I  sub  , 0  sub  ].The vector c raug can be obtained and the estimated mutual coupling coefficients can be obtained by {ĉ  }   =1 = mean( C ( + 1, :)).
In the first simulation, we investigate the angle estimation performances of MUSIC-Like, ESPRIT-Like, and tensorbased real-valued subspace (in this paper, we call it HOSVD) methods and our proposed TRSSIA method.There are two cases: (1)   =   = 1, c  = [1, 0.2 + 0.0061]T , and c  = [1, 0.15 + 0.0251]T : the number of pulses is  = 64; and (2)   =   = 2, c  = [1, 0.7 + 0.002, 0.2 + 0.0061]T , and c  = [1, 0.6+0.0121,0.15+0.0251]T : the number of pulses is  = 128.In Figure 1(a), it shows that TRSSIA estimates angles with more accuracy than the other approaches.The MUSIC-Like and ESPRIT-Like methods perform well in the high SNR region for case (1), as shown in Figure 1(b), and MUSIC-Like fails to work in case (2) owing to angle ambiguity, which is explained in [13].At both cases, the TRSSIA method is better than the other methods.Although the HOSVD methods make use of the multidimensional structure of the received data, it just selects  − 2  elements of transmit array and  − 2  elements of receive array for angles estimation.The TRSSIA algorithm restrains the noise of received data, improves the effective aperture, and uses all the elements of transmitter and receiver.
The second simulation is carried out to show the RMSE of angle estimation versus number of pulses with SNR = 5 dB. Figure 2(a) shows that the performances of MUSIC-Like and ESPRIT-Like approaches highly depend on the number of  pulses, and the TRSSIA method can give out more precise angle estimation even with small snapshots.In Figure 2(b), it analyses the ratio of the power of signal to the power of noise. 2 diag ,  2 off , and  2  represent the sum of diagonal elements of autocorrelation matrix of signal, the sum of off-diagonal elements of autocorrelation matrix of signal, and the sum of all elements of autocorrelation matrix of noise, respectively.In our proposed algorithm, the newly autocorrelation matrices of signal are Λ and Λ .For the other methods, we call them conventional methods; it is ∑  =1 s()s() H /. According to (1), the  2 off of the TRSSIA method is zero.Because the TRSSIA method restrains the noise, the  2 diag / 2  is tremendously higher than conventional methods.Though the  2 diag / 2  of conventional methods is constant, the  2 off / 2  gets lower as the number of pulses increases, which improves the accuracy of signal subspace estimation and angle estimation.
In the third simulation, it compares the angle estimation performance of the TRSSIA method with the performance of the HOSVD methods in the scenario of coherent targets.There are three cases: (1) f  = [100, 110, 200] T Hz; (2) f  = [100, 100, 150] T Hz; and (3) f  = [100, 100, 100] T Hz.As Figure 3 shows, the angle estimation performance of the HOSVD methods degrade, when all the targets move in a narrow speed zone.This method even cannot work when more than two coherent targets exist.Equation (1) proves that the rank is always  whether the targets are coherent or not, so the angle estimation performance of the TRSSIA method does not degrade even if there are more than two coherent targets.
As we know, the forward-backward (FB) smoothing is a preferable technique to cope with coherent signals [23].It should be notable that our proposed approach can perform as well as the joint transmit and receive diversity smoothing  (TRDS) algorithm with FB smoothing in [24], as it is shown in Figure 4(a).But the TRSSIA needs to compute covariance matrices only once.We present an evaluation of computational complexity using TIC and TOC instruction in MAT-LAB.In Figure 4  to measure the performance.In [13,14], a technique of estimating mutual coupling is proposed; we call it conventional method in this paper.Figure 5 shows that our proposed approach can estimate mutual coupling more accurately.There are two reasons.On one hand, the angle estimation of our proposed method is more accurate.On the other hand, our proposed method computes the mean of every mutual coupling coefficient.

Conclusion
This paper has proposed an algorithm for angle estimation with unknown mutual coupling in both the transmitter and receiver.The preliminary work is to restrain the white Gaussian noise of each channel by computing the correlation coefficients, because the noise is independent in both the spatial and temporal domains.In order to use more information of the received data, we do spatial smoothing in both the transmit array and the receive array and construct the augmented steering matrices with improved aperture.The TRSSIA algorithm adopts more elements from the transmitter and the receiver to estimate angles.So more information improves the angle estimation.For restraining noise, improving aperture and the spatial smoothing technique, the TRSSIA method proves better angle estimation than MUSIC-Like, ESPRIT-Like, and tensor-based real-valued approaches at small number of pulses and low SNR cases, and its angle estimation performance does not descend even for more than two coherent targets.Based on the more accurately estimated angles and computing the mean of every mutual coupling efficient, the mutual coupling estimation is more accurate than the other methods.The simulation results verify the advantage of the proposed method.
sub  rows of U  , respectively.Then the above equation can be written as Then EVD can be employed to solve Φ  ; furthermore, DOD can be got.The above deducing procedure proves the rotational invariance property of X taug .