Next Article in Journal
Fragmentation Attacks and Countermeasures on 6LoWPAN Internet of Things Networks: Survey and Simulation
Previous Article in Journal
Fringe Detection and Displacement Sensing for Variable Optical Feedback-Based Self-Mixing Interferometry by Using Deep Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coherent Signal DOA Estimation for MIMO Radar under Composite Background of Strong Interference and Non-Uniform Noise

1
Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China
2
Graduate College, Air Force Engineering University, Xi’an 710051, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(24), 9833; https://doi.org/10.3390/s22249833
Submission received: 31 October 2022 / Revised: 8 December 2022 / Accepted: 13 December 2022 / Published: 14 December 2022
(This article belongs to the Special Issue Signal Processing for Bi- and Multistatic Radar Systems)

Abstract

:
To address the problems of low accuracy and low robustness of the conventional algorithm in estimating the direction of arrival (DOA) of coherent signals against a composite background of strong interference and non-uniform noise, in this paper, a coherent signal DOA estimation algorithm based on fixed projection blocking is proposed in conjunction with a multi-input multi-output (MIMO) radar. The covariance matrix of the received signal is first decomposed by eigenvalues, and a fixed projection matrix orthogonal to the interference guidance vector is constructed as the interference blocking matrix. Then, the received array signal is pre-processed to re-form the covariance matrix, and this matrix is rendered decoherent through a Toeplitz reconstruction. Finally, the reconstructed covariance matrix is estimated by DOA using the propagation operator algorithm to reduce the complexity. The simulation verifies that the proposed algorithm has a better robustness and higher accuracy than conventional algorithms for the DOA estimation of coherent signals in composite backgrounds.

1. Introduction

The MIMO radar is a new radar system with many advantages over traditional phased array radars, including a low interception probability, strong anti-jamming capability and virtual arrays with large apertures. However, with the increasing complexity of the battlefield electromagnetic environment, MIMO radar detection also faces an increasing amount of interference, including coherent signals from multipath effects that interfere with each other, non-uniform noise caused by differences in antenna arrays and received signal channels, and strong interference signals from enemy jammers that are close to the radar or have a high effective interference power [1,2,3,4]. The effects of these disturbances will lead to problems of reduced accuracy, reduced stability and even failure of traditional direction-finding algorithms [5,6,7,8,9,10], which need to be addressed by effective technical means. If the sensor nodes can be connected to each other, the DV-hop algorithm proposed in the literature [11] can also help with localization.
In addition, to overcome the problem of a directional measurement under non-uniform noise, the authors of [12] pre-whitened received signal data to achieve a higher accuracy estimation, but the problem of aperture loss led to a limited number of measurement sources. The authors of [13] used the maximum likelihood estimation method to gradually iterate to obtain the signal subspace and noise covariance matrix, which improved the accuracy, because of the need for iterative operations; thus, the complexity was high. The authors of [14] used the centrosymmetric property of uniform line array to transform complex matrices into real-valued matrices and proposed the unitary low-rank matrix decomposition method to improve DOA estimation performance and reduce complexity. The authors of [15] constructed a spatial difference matrix to suppress noise and used a rotation-invariant subspace-like algorithm to solve the model sparsely, which has a strong robustness. To address the problem of coherent signal orientation under non-uniform noise, a generalized covariance matrix was constructed in [16] using a forward-backward spatial smoothing algorithm and the difference between the complex conjugate of the data matrix, which fully eliminated spatial non-uniform noise and achieved the estimation of coherent signals, but there was a problem of aperture loss. The authors of [17] proposed a DOA estimation based on the tail minimization method, which exploits the sparsity of the signal subspace and has the advantage of expanding the virtual array aperture. The authors of [18] used a spatial smoothing algorithm to reconstruct a full-rank signal covariance matrix, which is beneficial for the DOA estimation of uncorrelated or coherent signals, but with high complexity.
For the current problem of direction-finding in the context of strong interference, research ideas can be broadly divided into two types: one involves sieving out the strong interference signal and the weak signal after estimating them at the same time, and the other involves estimating the weak signal after removing the strong interference signal. For the former idea, the RELAX algorithm [19] method is mainly used; such an algorithm separates the array output part of all signals into multiple data blocks, and then, according to the characteristics of the interference signal, finds out the data blocks containing an interference signal and reject, so as to achieve the purpose of target identification; however, such an algorithm is problematic due to high complexity. The second idea is the subarray interference suppression algorithm proposed in the literature [20,21], which shows a good performance by improving the beam formation method, so that it can suppress interference before weighting the subarray for the DOA estimation of weak signals. The authors of [22] use the angular information of the signal matrix to construct an interference blocking matrix to suppress strong interference signals, but this method requires precisely known interference directions, and different blocking matrices need to be constructed for different array structures. Gong et al. [23] analyzed in detail the direction-preserving orthogonality in array signal data under the background of strong interference and performed simulations for a variety of scenarios, verifying that the algorithm was suitable for situations with low signal-to-noise ratios (SNR) and small angular spacing. The authors of [24] propose a broadband signal localization method based on sparse spectral fitting and a zero-value matrix filter with a high localization accuracy in a strong interference environment. For the DOA estimation of coherent sources in a strong interference background, Li et al. [25] proposed an improved interference blocking algorithm to render the coherent signal decoherent by Toeplitz matrix reconstruction, and then applied the interference blocking matrix method to suppress the interfering signal, but the algorithm was not robust.
The above algorithms suffer from a low estimation accuracy, poor robustness and high complexity in direction-finding against a composite background of strong interference and non-uniform noise. This paper presents a coherent weak-target DOA estimation method based on a fixed projection blocking algorithm and Toeplitz matrix reconstruction. The method obtains the interference subspace by decomposing the eigenvalues of the received signal covariance matrix, then constructs a fixed projection matrix orthogonal to the interference-oriented vector as the interference blocking matrix, multiplies it with the received array signal for pre-processing and re-forming the covariance matrix, performs a Toeplitz reconstruction to suppress non-uniform noise, and finally uses the propagator method (PM) for DOA estimation, avoiding the spectral peak search.

2. MIMO Radar Signal Model

As shown in Figure 1, it is assumed that the number of array elements of the transmitting array is M, the number of array elements of the receiving array is N, the array element spacing is d, the measured target is a far-field target, the signal reflected from the target is a narrow-band signal, and the wave front arriving at each array element is a plane wave.
K are far-field narrowband signals, and J are strong interference signals s i ( t ) ( i = 1 , 2 , , J , , K + J ) incident to the uniform line array; the first J signals are strong interference signals whose power satisfies σ 1 2 > σ 2 2 > > σ J 2 >> σ J + 1 2 > > σ K + J 2 ; the received K are far-field narrowband signals containing Q independent signals; L are the group of coherent signals in the presence of P coherent signals, where the number of coherent signals in the l-th group is P l ; the radar transmit and receive guide vectors are:
a t ( θ ) = [ 1 , e j 2 π d 0 sin θ / λ , , e j 2 π ( M 1 ) d 0 sin θ / λ ] T
a r ( θ ) = [ 1 , e j 2 π d 0 sin θ / λ , , e j 2 π ( N 1 ) d 0 sin θ / λ ] T
Based on the orthogonality of the transmit waveform, the received signal at the m-th array element of the array at moment t after matching and filtering the signal received at the receiver end of the MIMO radar can be expressed as:
x m ( t ) = i = 1 J ( a r ( θ i ) a t ( θ i ) ) s i ( t ) + k = 1 K ( a r ( θ k ) a t ( θ k ) ) s k ( t ) + n m ( t )
where n m ( t ) is the non-uniform noise of the m-th array element at moment t.
Rewriting the received signal into vector form, one obtains:
x ( t ) = j = 1 J ( a r ( θ j ) a t ( θ j ) ) s j ( t ) + i = 1 Q ( a r ( θ i ) a t ( θ i ) ) s i ( t ) + l = 1 L k = 1 P l ( a r ( θ l k ) a t ( θ l k ) ) b l k s c ( t ) + n ( t ) = A p s p ( t ) + A u s u ( t ) + A c s c ( t ) + n ( t ) = A s ( t ) + n ( t )
where A p = [ a r ( θ 1 ) a t ( θ 1 ) , a r ( θ 2 ) a t ( θ 2 ) , , a r ( θ J ) a t ( θ J ) ] denotes the array flow pattern of the strong interference signal, s p ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s J ( t ) ] T denotes the incident signal vector of J strong interfering signals, A u = [ a r ( θ 1 ) a t ( θ 1 ) , a r ( θ 2 ) a t ( θ 2 ) , , a r ( θ Q ) a t ( θ Q ) ] indicates the array flow pattern of the independent signal, s u ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s Q ( t ) ] T denotes the incident signal vector of Q independent signals, A c = [ A c 1 b 1 , , A c L b L ] , b l = [ b l 1 , b l 2 , , b l P l ] T , b l k denotes the coherence coefficient of the k-th signal in the l-th set of coherent signals, s c ( t ) = [ s J + Q + 1 ( t ) , s J + Q + 2 ( t ) , , s J + Q + L ( t ) ] T denotes the incident signal vector of the L-group coherent signal, and A c l = [ a r ( θ l 1 ) a t ( θ l 1 ) , a r ( θ l 2 ) a t ( θ l 2 ) , , a r ( θ l P l ) a t ( θ l P l ) ] , A = [ A p , A u , A c ] , s ( t ) = [ s p ( t ) , s u ( t ) , s c ( t ) ] T , n ( t ) is the non-uniform noise at the receiver.
The covariance matrix R x of the signal is expressed as:
R x = E [ x ( t ) x H ( t ) ] = A R s A H + R n
where R s is the covariance matrix of the signal vector s ( t ) , and R n = diag { σ 1 2 , , σ M N 2 } is the covariance matrix of the non-uniform noise n ( t ) , where σ 1 2 , , σ M N 2 is the unequal noise power.

3. Coherent Signal DOA Estimation in A Composite Background

3.1. Fixed Projection Blocking Algorithm Rejects Interference

R x is the positive definite Hermitain matrix, for which the eigenvalue decomposition is:
R x = U Σ U H = i = 1 M λ i e i e i H
where U = [ e 1 , e 2 , e M N ] H , Σ = diag { λ 1 , λ 2 , , λ M N } , e i is the eigenvector of matrix R x , and λ i is the eigenvalue of matrix R x , sorted as follows:
λ 1 > λ 2 > > λ J > > λ J + 1 > > λ J + K > > λ J + K + 1 > > λ MN
The eigenvalues are ordered to obtain the corresponding eigenvectors as e 1 , e 2 , , e M N , where the eigenvectors corresponding to the J large eigenvalues of the matrix R x form the strong interference subspace:
E J = [ e 1 , e 2 , , e J ]
The eigenvectors corresponding to the K larger eigenvalues form the signal subspace:
E K = [ e J + 1 , e J + 2 , , e J + K ]
Other eigenvectors corresponding to M N K J small eigenvalues form the noise subspace:
E N = [ e K + J + 1 , e K + J + 2 , , e M N ]
Construct the fixed projection blocking matrix G :
G = I E J ( E J H E J ) 1 E J
which is orthogonal to the strong interference subspace, where I is the unit matrix of M N × M N dimensions.
The subspace E = [ E J , E K , E N ] is pre-transformed by the matrix G :
D = G E = G [ E J , E K , E N ] = [ 0 J , D J + K + 1 , D J + K + 2 , D M N ]
From the pre-transformation, it follows that the transformed subspace D is only relevant for the target signal and the noise, that is the fixed projection blocking matrix G is able to suppress strong interference signals.
Similarly, the data vector x ( t ) of the signal is subjected to interference rejection:
y ( t ) = G x ( t )
Subsequently, the new covariance matrix is obtained from the pre-processed signal data vector:
R y = G A R s A H G H + G R n G H
For the new covariance matrix, although the influence of interfering signals on the target signal is eliminated, the presence of coherent signals as well as non-uniform noise information in the matrix still affects the estimation accuracy of the signal.

3.2. Toeplitz Matrix Reconstruction

The matrix R y can be minimized by Toeplitzisation of the difference with the noiseless covariance matrix, that is:
min R T S T R T R
where R T is the Toeplitz reconstruction of R y , R is the noiseless covariance matrix, and S T is the set of Toeplitz matrices. Essentially, the covariance matrix is averaged over the diagonal elements of the covariance matrix.
Define each element of the matrix as:
R y = [ r y ( 1 , 1 ) r y ( 1 , 2 ) r y ( 1 , 3 ) r y ( 1 , M N ) r y ( 2 , 1 ) r y ( 2 , 2 ) r y ( 2 , 3 ) r y ( 2 , M N ) r y ( M N , 1 ) r y ( M N , 2 ) r y ( M N , M N ) ]
Construct a new matrix by extracting the elements of the first row of the matrix R y :
R ˜ y = [ r y ( 1 , 1 ) r y ( 1 , 2 ) r y ( 1 , 3 ) r y ( 1 , M N ) 0 r y ( 1 , 1 ) r y ( 1 , 2 ) r y ( 1 , M N 1 ) 0 0 r y ( 1 , 1 ) r y ( 1 , M N 2 ) 0 0 r y ( 1 , 1 ) ]
The correlation vector of Equation (17) contains all the angular information of the incident source, and a Toeplitz matrix can be constructed from this vector with the following structure:
R ˜ ˜ y = [ r y ( 1 , 1 ) r y ( 1 , 2 ) r y ( 1 , 3 ) r y ( 1 , M N ) r y * ( 1 , 2 ) r y ( 1 , 1 ) r y ( 1 , 2 ) r y ( 1 , M N 1 ) r y * ( 1 , 3 ) r y * ( 1 , 2 ) r y ( 1 , 1 ) r y ( 1 , M N 2 ) r y * ( 1 , M N ) r y * ( 1 , M N 1 ) r y ( 1 , 1 ) ]
where R ˜ ˜ y is the M N × M N -dimensional Hermitian matrix. As the effect of non-stationary noise is concentrated in the main diagonal of the covariance matrix and the Toeplitz matrix is reconstructed with the main diagonal elements of the covariance matrix averaged, the effect of non-stationary noise is reduced and decoherence is achieved, but the matrix also suffers from excessive complexity compared to conventional algorithms.

3.3. PM Algorithm DOA Estimation

According to the literature [26], the matrix R ˜ ˜ y is divided into two submatrices:
R ˜ ˜ y = [ R K 1 T , R K 2 T ] T
where R K 1 consists of the first K rows of the matrix R ˜ ˜ y (K is the number of target signals), and R K 2 consists of the next M N K rows. Since R K 1 , R K 2 are full rank vandermond matrices, there exists an K × ( M N K ) -dimensional linear operator P, so that:
P = ( R K 1 R K 1 H ) 1 R K 1 R K 2 H
For ease of calculation, define Ω = [ P T , I M N K ] T ( I M N K is the M N K -dimensional unit matrix). Similarly, define a ( θ ) = a r ( θ ) a t ( θ ) , and the estimate θ of θ ¯ is the value taken when minimizing the cost function f ( θ ) :
f ( θ ) = a H ( θ ) Π a ( θ )
where Π = Ω ( Ω H Ω ) - 1 Ω H .

4. Basic Steps of the Algorithm

Based on the above analysis, the operational steps of the Toeplitz reconstruction algorithm based on fixed projection blocking can be summarized as follows:
  • Calculate the data covariance matrix R of the array element output vector x ( t ) ;
  • The covariance matrix R x is decomposed to obtain the interference subspace E J , the signal subspace E K , and the noise subspace E N ;
  • To eliminate interference, construct a fixed projection blocking matrix G using Equation (11);
  • Pre-processing the received signal by Equation (13) gives y ( t ) ;
  • The covariance matrix R y is reconstructed according to Equation (18) and replaced to obtain the matrix R ˜ ˜ y ;
  • After partitioning the matrix R ˜ ˜ y , construct the linear operator P according to Equation (20);
  • The DOA is obtained by minimizing the cost function x using Equation (21).

5. Algorithm Performance Analysis

5.1. Source Overload Capacity Analysis

For transceiver arrays consisting of M and N arrays, the proposed algorithm is able to identify M N 1 coherent signals without loss of aperture using the Toeplitz matrix reconstruction algorithm. If the forward spatial smoothing algorithm is used to identify M N / 2 coherent signals at most, the forward and backward spatial smoothing algorithm can identify 2 M N / 3 coherent signals at most. Compared with the conventional algorithm, the proposed algorithm has an advantage in the number of coherent sources that are identified.

5.2. Algorithmic Complexity Analysis

Assuming that the number of sources measured is K, the number of samples is L, and the number of angular searches is N θ , the complexity of the fixed projection blocking-Toeplitz (FPB-To) algorithm proposed in this section is mainly derived from the covariance matrix formation, the eigenvalue decomposition and the PM algorithm, and the computational effort is approximated by O ( 2 L ( M N ) 2 + ( M N ) 3 + 10 K 3 + 4 K M N ) . In the case of N θ > L > M N > K , the comparison is slightly lower than for the jamming jam method (JJM) algorithm, the extended noise subspace (ENS) algorithm, the Toeplitz-jamming jam method (To-JJM) algorithm [25] with complexity O ( L ( M N ) 2 + ( M N ) 3 + N θ ( M N ) 2 ) , and the fixed projection blocking (FPB) algorithm with complexity O ( 2 L ( M N ) 2 + 2 ( M N ) 3 + N θ ( M N ) 2 ) .

6. Simulation Results and Analysis

The simulation uses a uniform line array with the number of array elements M and N. R n = diag { σ 1 2 , σ M N 2 } , σ m 2 ( m = 1 , , M N ) is the different noise power of each array element; then, the signal-to-noise ratio is defined as SNR = 10 log 10 ( 1 M N m = 1 M N σ k 2 σ m 2 ) , and the interference signal ratio (ISR) is defined as 10 log 10 ( σ j 2 / σ k 2 ) , where σ k 2 is the average power of the target signal and σ j 2 is the average power of the strong interference signal.
The root mean square error (RMSE) of the DOA estimate is:
RMSE = 1 / 2 N m o n t i = 1 N m o n t k = 1 K ( | θ i θ i k | 2 )
where N m o n t is the number of Monte–Carlo experiments, K is the number of signals, and θ i and θ i k are the true and estimated values of the i-th signal angle in each trial, respectively.

6.1. Spatial Spectrum Estimation

The simulation sets M = N = 12; there is one set of coherent signals with a coherence coefficient of 0.9, incidence angles of −60° and –40°, three non-coherent signals with incidence angles of –20°, 0° and 20°, and two strong interfering signals with incidence angles of 40° and 60°; the SNR is 0 dB, the ISR is 30 dB, and the number of snapshots is 350. Figure 2 gives the spatial spectral curves of the JJM algorithm, ENS algorithm, FPB algorithm, To-JJM algorithm and the spatial spectrum curve of the FPB-To algorithm proposed in this paper.
Figure 2 shows that the JJM algorithm, ENS algorithm, and FPB algorithm all have specific interference cancellation mechanisms, which can suppress strong interference signals and estimate the DOA of non-coherent signals better, but they are less effective in estimating coherent signals. The To-JJM algorithm can better estimate the DOA of coherent targets in a composite background but requires an accurate prediction of the direction of strong interference signals, which is not conducive to the engineering application of the algorithm. The FPB-To algorithm is the most effective and can accurately detect the target DOA without predicting the direction of interference.

6.2. Comparative Error Analysis of Different Algorithms

The number of Monte–Carlo experiments performed is 500, and the other settings are the same as in Section 6.1. Figure 3a provides the curves of RMSE versus snapshots for the JJM algorithm, ENS algorithm, FPB algorithm, To-JJM algorithm, and FPB-To algorithm by varying the snapshots from 50 to 500 in steps of 50. Figure 3b shows the curve of RMSE with SNR for the five algorithms by varying the SNR from –10 dB to 10 dB in steps of 2 dB.
As can be seen from Figure 3, the RMSE of these methods decreases as the SNR and the number of snapshots increase. As the number of snapshots and the SNR increase to a certain level, the RMSE of several algorithms level off. The To-JJM algorithm outperforms several other algorithms at a higher SNR due to the decoherence of the target and non-uniform noise suppression. The algorithm in this paper is pre-transformed in the subspace, which reduces the sensitivity to noise to a certain extent and still provides a good directional accuracy at a low SNR.

6.3. Comparative Error Analysis with Different Numbers of Interference Sources

We change the number of interference sources on the basis of the settings in Section 6.1. For a strong interference number of 1, the angle of incidence is –40°; for a strong interference number of 2, the angles of incidence are –40° and –60°; for a strong interference number of 3, the angles of incidence are –40°, –60°, and –80°. Figure 4a shows the curve of the To-JJM algorithm SNR varying from –10 dB to 10 dB in steps of 2 dB, with RMSE varying with SNR for three numbers of interference sources. Figure 4b shows the curve of RMSE versus SNR for the FPB-To algorithm SNR varying from –10 dB to 10 dB in steps of 2 dB for three sources of interference.
As can be seen from Figure 4, both algorithms still do a good job of suppressing interference and accurately estimating weak signals against the composite background. The To-JJM algorithm performs better when suppressing one interfering signal than when suppressing three, the reason being that for every interfering signal suppressed by the algorithm, a certain amount of aperture loss is incurred, which affects the direction-finding accuracy. However, the proposed FPB-To algorithm does not make a significant difference as the number of interference sources increases, and still achieves a high accuracy in direction finding.

6.4. Algorithm Robustness Analysis

To measure the difference in noise intensity received in different arrays, the worst noise power ratio (WNPR) is defined as:
WNPR = σ m max 2 / σ m min 2
This is the ratio of the power of the strongest noise to the weakest noise received by each physical array element.
The simulation setup is the same as in Section 6.1. Figure 5a provides the curves of the algorithm WNPR varying from 10 to 100 in steps of 10, with the RMSE varying with WNPR for different algorithms.
In addition, the simulation is set up with three different noise covariance matrices:
N 1 = diag { 1 , 1 , 1 , 2 , 2 , 2 , 2 , 3 , 3 , 3 } N 2 = diag { 5 , 9 , 15 , 2 , 3 , 4 , 2 , 4 , 4 , 1 } N 3 = diag { 1 , 2 , 3 , 4 , 5 , 1 , 2 , 3 , 4 , 5 }
Other simulation conditions are the same as in Section 6.1. Figure 5b shows the spatial spectral estimation performance of the proposed algorithm for different noise covariance matrices.
As can be seen from Figure 5, the WNPR characterizes the inhomogeneity of the non-uniform noise, and the larger the WNPR, the stronger the inhomogeneity; additionally, the algorithm performances are all degraded to different degrees. Compared with other algorithms, the algorithm in this section is able to maintain an estimation accuracy under different WNPRs and noise covariance matrices, and has a strong robustness.

7. Conclusions

In answer to the problem of coherent source vectoring in the composite background of strong interference and non-uniform noise, this paper proposes a fixed projection blocking algorithm based on Toeplitz reconstruction. The algorithm first constructs the covariance matrix from the received signal vector, while an eigenvalue decomposition is then performed, the resulting strong interference vector is used to construct the orthogonal projection matrix of the interference subspace, and the covariance matrix is obtained again after pre-processing the received array signal. This is followed by a Toeplitz matrix reconstruction for decoherence and finally a PM algorithm to obtain the DOA information of the faint signal. The experimental results show that the proposed algorithm is able to perform a DOA estimation of coherent targets in strong interference and non-uniform noise environments with a higher accuracy and robustness than conventional algorithms. It is conducive to a better use in future complex and changing electromagnetic environments [27,28,29,30,31].

Author Contributions

This article was co-authored by B.L., G.H., H.Z. and G.Z. The major individual contributions are as follows: conceptualization, B.L. and G.H.; methodology, B.L., G.H. and H.Z.; software, B.L.; validation, H.Z.; formal analysis, G.Z.; investigation, B.L.; resources, G.H.; data curation, H.Z.; writing—original draft preparation, G.Z.; writing—review and editing, B.L.; supervision, G.H.; project administration, G.H.; funding acquisition, G.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of China, grant number 62071476.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available, due to the data in this paper not coming from publicly available datasets but being obtained from the simulation of the signal models listed in the article.

Acknowledgments

We thank the college for providing us with an efficient simulation platform so that we can complete the experimental simulation as scheduled. Funding from the National Science Foundation of China (No. 62071476) is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Salari, S.; Chan, F.; Chan, T. DOA Estimation Using Compressive Sampling-Based Sensors in the Presence of Interference. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 4395–4405. [Google Scholar] [CrossRef]
  2. Chen, C.; Zhou, J.; Tang, M. Direction of Arrival Estimation in Elliptical Models via Sparse Penalized Likelihood Approach. Sensors 2019, 19, 2356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, B.; Liu, S.; Zhao, D. DOA Estimation of Unknown Emitter Signal Based on Time Reversal and Coprime Array. Sensors 2019, 19, 1398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zheng, G.; Song, Y.; Chen, C. Height Measurement with Meter Wave Polarimetric MIMO Radar: Signal Model and MUSIC-like Algorithm. IEEE Signal Process. 2021, 190, 1872–7557. [Google Scholar] [CrossRef]
  5. Sun, K.; Liu, Y.; Meng, H.; Wang, X. Adaptive Sparse Representation for Source Localization with Gain/Phase Errors. Sensors 2011, 11, 4780–4793. [Google Scholar] [CrossRef] [Green Version]
  6. Fu, W.; Hu, Z.; Li, D. A Sorting Algorithm for Multiple Frequency-Hopping Signals in Complex Electromagnetic Environments. Circuits Syst. Signal Process. 2020, 39, 245–267. [Google Scholar] [CrossRef]
  7. Nicola, M.; Falco, G.; Morales Ferre, R.; Lohan, E.-S.; de la Fuente, A.; Falletti, E. Collaborative Solutions for Interference Management in GNSS-Based Aircraft Navigation. Sensors 2020, 20, 4085. [Google Scholar] [CrossRef]
  8. Muhammad, M.; Li, M.; Abbasi, Q.; Goh, C.; Imran, M.A. A Covariance Matrix Reconstruction Approach for Single Snapshot Direction of Arrival Estimation. Sensors 2022, 22, 3096. [Google Scholar] [CrossRef]
  9. Zhang, C.; Han, Y.; Zhang, P. Research on Modern Radar Emitter Modelling Technique under Complex Electromagnetic Environment. J. Eng. 2019, 20, 7134–7138. [Google Scholar] [CrossRef]
  10. Yao, Y.; Lei, H.; He, W. A-CRNN-Based Method for Coherent DOA Estimation with Unknown Source Number. Sensors 2020, 20, 2296. [Google Scholar] [CrossRef]
  11. Gui, L.; Xiao, F.; Zhou, Y.; Shu, F.; Val, T. Connectivity Based DV-Hop Localization for Internet of Things. IEEE Trans. Veh. Technol. 2020, 69, 8949–8958. [Google Scholar] [CrossRef]
  12. Wu, Y.; Hou, C.; Liao, G. Direction-of-Arrival Estimation in the Presence of Unknown Nonuniform Noise Fields. IEEE J. Ocean. Eng. 2006, 31, 504–510. [Google Scholar] [CrossRef] [Green Version]
  13. Liao, B.; Chan, S.; Huang, L. Iterative Methods for Subspace and DOA Estimation in Nonuniform Noise. IEEE Trans. Signal Process. 2016, 64, 3008–3020. [Google Scholar] [CrossRef]
  14. Zhu, Y.; Wang, X.; Wan, L. Unitary Low-Rank Matrix Decomposition for DOA Estimation in Nonuniform Noise. In Proceedings of the 2018 IEEE 23rd International Conference on Digital Signal Processing, Shanghai, China, 19–21 November 2018. [Google Scholar]
  15. Wei, J.; Cao, K.; Ning, F. DOA Estimation Algorithm for Sparse Representation under Nonstationary Noise. J. Beijing Univ. Posts Telecommun. 2020, 43, 116–121. [Google Scholar]
  16. Qi, C.; Chen, Z.; Wang, Y. DOA Estimation for Coherent Sources in Unknown Nonuniform Noise Fields. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1195–1204. [Google Scholar]
  17. Zheng, C.; Fang, Y.; Li, J. Sparse Representation Based Direction-of-Arrival Estimation in Nonuniform Noise via Tail Minimisation. IET Radar Sonar Navig. 2021, 15, 788–795. [Google Scholar] [CrossRef]
  18. Wan, J.; Liao, B.; Guo, C. Spatial Smoothing Based Methods for Direction-of-Arrival Estimation of Coherent Signals in Nonuniform Noise. Digit. Signal Process. 2017, 67, 116–122. [Google Scholar] [CrossRef]
  19. Tian, Y.; Liang, G.; Fu, J. Direction Estimation of Coherent Signals Based on the Symmetry of Uniform Linear Array. Proc. Meet. Acoust. 2019, 33, 1–9. [Google Scholar]
  20. Chai, L.; Luo, J. A Combined Algorithm Both for Restraining Intensive Interferences and for Weak Signals DOA Estimation. In Proceedings of the ICMMT 4th International Conference on, Proceedings Microwave and Millimeter Wave Technology, Beijing, China, 18–21 August 2004; pp. 243–246. [Google Scholar]
  21. Chai, L.; Luo, J. A New Method for DOA Estimation of Weak Signals under Strong Interference Conditions. J. Electron. Inf. 2005, 27, 1517–1520. [Google Scholar]
  22. Su, C.; Luo, J. A Uniform Circular Array Subarray Interference Suppression DOA Estimation Algorithm. Signal Process. 2010, 26, 1355–1360. [Google Scholar]
  23. Gong, J.; Lou, S. DOA Estimation Method of Weak Sources for An Array Antenna under Strong Interference Conditions. Int. J. Electron. 2018, 105, 1931–1944. [Google Scholar] [CrossRef]
  24. Yang, Y.; Zhang, Y.; Long, Y. Wideband Sparse Spatial Spectrum Estimation Using Matrix Filter with Nulling in A Strong Interference Environment. J. Acoust. Soc. Am. 2018, 143, 3891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Li, H.; Guo, Y.; Gong, J. DOA Estimation Method for Coherent Weak Sources under Strong Interference Background. Mod. Radar 2012, 34, 45–49. [Google Scholar]
  26. Al-Ardi, E.; Shubair, R.; Al-Mualla, M. Computationally Efficient High-Resolution DOA Estimation in Multipath Environment. Electron. Lett. 2004, 40, 908–910. [Google Scholar] [CrossRef]
  27. Liu, C.; Wu, R.; He, Z. Modeling and Analyzing Interference Signal in A Complex Electromagnetic Environment. EURASIP J. Wirel. Commun. Netw. 2016, 1, 1–9. [Google Scholar] [CrossRef] [Green Version]
  28. Wan, F.; Duvail, F.; Savatier, X. Electromagnetic Interference Detection Method to Increase The Immunity of A Icrocontroller-Based System in A Complex Electromagnetic Environment. Sci. Meas. Technol. 2012, 6, 254–260. [Google Scholar] [CrossRef]
  29. Chen, Z.; Ding, Y.; Ren, S.; Chen, Z. A Novel Noncircular MUSIC Algorithm Based on the Concept of the Difference and Sum Coarray. Sensors 2018, 18, 344. [Google Scholar] [CrossRef] [Green Version]
  30. Wu, S.Y.; Zhao, J.; Dong, X.D.; Xue, Q.T.; Cai, R.H. DOA Tracking Based on Unscented Transform Multi-Bernoulli Filter in Impulse Noise Environment. Sensors 2019, 19, 4031. [Google Scholar] [CrossRef] [Green Version]
  31. Qi, B.; Liu, D. An Enhanced Spatial Smoothing Algorithm for Coherent Signals DOA Estimation. Eng. Comput. 2021, 39, 574–586. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the MIMO radar architecture.
Figure 1. Schematic diagram of the MIMO radar architecture.
Sensors 22 09833 g001
Figure 2. Comparison of spatial spectral curves.
Figure 2. Comparison of spatial spectral curves.
Sensors 22 09833 g002
Figure 3. RMSE analysis: (a) RMSE in relation to snapshots; (b) RMSE in relation to SNR.
Figure 3. RMSE analysis: (a) RMSE in relation to snapshots; (b) RMSE in relation to SNR.
Sensors 22 09833 g003
Figure 4. Variation curve of algorithm performance with SNR under different numbers of interference sources: (a) To-JJM algorithm; (b) FPB-To algorithm.
Figure 4. Variation curve of algorithm performance with SNR under different numbers of interference sources: (a) To-JJM algorithm; (b) FPB-To algorithm.
Sensors 22 09833 g004
Figure 5. Algorithm robustness analysis: (a) Comparison of errors at different WNPRs; (b) Spatial spectral curves for different noise parameters.
Figure 5. Algorithm robustness analysis: (a) Comparison of errors at different WNPRs; (b) Spatial spectral curves for different noise parameters.
Sensors 22 09833 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, B.; Hu, G.; Zhou, H.; Zheng, G. Coherent Signal DOA Estimation for MIMO Radar under Composite Background of Strong Interference and Non-Uniform Noise. Sensors 2022, 22, 9833. https://doi.org/10.3390/s22249833

AMA Style

Lin B, Hu G, Zhou H, Zheng G. Coherent Signal DOA Estimation for MIMO Radar under Composite Background of Strong Interference and Non-Uniform Noise. Sensors. 2022; 22(24):9833. https://doi.org/10.3390/s22249833

Chicago/Turabian Style

Lin, Bin, Guoping Hu, Hao Zhou, and Guimei Zheng. 2022. "Coherent Signal DOA Estimation for MIMO Radar under Composite Background of Strong Interference and Non-Uniform Noise" Sensors 22, no. 24: 9833. https://doi.org/10.3390/s22249833

APA Style

Lin, B., Hu, G., Zhou, H., & Zheng, G. (2022). Coherent Signal DOA Estimation for MIMO Radar under Composite Background of Strong Interference and Non-Uniform Noise. Sensors, 22(24), 9833. https://doi.org/10.3390/s22249833

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop