Open AccessArticle

The DOA Estimation Method for Low-Altitude Targets under the Background of Impulse Noise

Bin Lin

^1,2,

Guoping Hu

^1,*,

Hao Zhou

Guimei Zheng

¹ and

Yuwei Song

Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China

Graduate College, Air Force Engineering University, Xi’an 710051, China

Author to whom correspondence should be addressed.

Sensors 2022, 22(13), 4853; https://doi.org/10.3390/s22134853

Submission received: 7 June 2022 / Revised: 21 June 2022 / Accepted: 24 June 2022 / Published: 27 June 2022

(This article belongs to the Special Issue Advanced Signal Processing for Next Generation Wireless Communications)

Download

Browse Figures

Figure 1
Signal reception model in low-altitude multipath environment. "> Figure 2
Algorithm space spectrum under the background of GMM noise: (a) <math display="inline"><semantics> <mrow> <msub> <mi>c</mi> <mn>2</mn> </msub> </mrow> </semantics></math> = 0.1; (b) <math display="inline"><semantics> <mrow> <msub> <mi>c</mi> <mn>2</mn> </msub> </mrow> </semantics></math> = 0.3. "> Figure 3
Spatial spectrum under noise background with the algorithm: (a) <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.6</mn> </mrow> </semantics></math>; (b) <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.3</mn> </mrow> </semantics></math>. "> Figure 4
Analysis of mean square error and success rate under GMM noise: (a) Relationship between RMSE and SNR; (b) Relationship between success rate and SNR. "> Figure 5
Analysis of the mean square error and success rate of the algorithm under noise: (a) Relationship between RMSE and GSNR; (b) Relationship between success rate and GSNR. "> Figure 6
Influence of the number of snapshots on the performance of the algorithm: (a) GMM noise; (b) <math display="inline"><semantics> <mrow> <mi mathvariant="normal">S</mi> <mi>α</mi> <mi mathvariant="normal">S</mi> </mrow> </semantics></math> noise. "> Figure 7
Influence of characteristic exponent on algorithm performance: (a) GSNR = 10 dB; (b) GSNR = 20 dB. ">

Versions Notes

Abstract

Due to the discontinuity of ocean waves and mountains, there are often multipath propagation effects and obvious pulse characteristics in low-altitude detection. If the conventional direction of arrival (DOA) estimation method is directly used for direction finding, it will lead to a large error. In view of serious misalignment in the DOA estimation of multipath signals under the background of impulse noise, a DOA estimation method based on spatial difference and a modified projection subspace algorithm is proposed in this paper. Firstly, the covariance matrix of the received data vector is used for spatial difference to eliminate the multipath effects of low-altitude targets. Secondly, the modified projection matrix is constructed using the signal source estimated with the least squares criterion and then used for modifying the covariance matrix, thus eliminating the cross-covariance matrices that affect the estimation accuracy. Finally, the modified covariance matrix is used for the DOA estimation of targets. Simulations show that the proposed algorithm achieves a higher accuracy in the DOA estimation of low-altitude targets than conventional algorithms under two common impulse noise models, without requiring prior knowledge of impulse noise.

Keywords:

impulse noise; low-altitude targets; multipath effects; array signal processing; direction of arrival estimation

1. Introduction

As an important part of array signal processing, DOA estimation is a characteristic parameter commonly used in radar detection and localization. So far, a series of related studies have been carried out on DOA estimation. At the earliest, Schmidt [1] proposed the multiple signal classification (MUSIC) algorithm based on feature subspace. Later, scholars made improvements, most of which were based on the assumption of Gaussian noise. With the increasingly complex battlefield electromagnetic environment, there is more and more interference in radar detection, making it more difficult to detect targets [2,3,4]. Especially in the low-altitude environment, the diverse terrain environment will greatly weaken the energy of electromagnetic waves, always making the detection and tracking of low-altitude targets poor [5,6,7,8]. Recent studies have shown that due to the discontinuity of ocean waves and mountains, the output signal of the array exhibits obvious impulse characteristics, and the alpha-stable distribution with heavy smearing characteristics and the compound Gaussian model can effectively describe impulse noise [9,10]. Therefore, it is of great significance to develop a DOA estimation method for low-altitude targets under the background of impulse noise.

In order to solve the problem of low-altitude multipath effects, scholars have conducted in-depth research, and they have achieved many results. Li et al. [11] first proposed the use of spatial smoothing technology to deal with coherent signals in multipath effects, but in addition to aperture loss, it was only suitable for uniform arrays. In the literature [12,13,14], the spatial differencing algorithm was used for decoherence, where the accuracy was improved to some extent but the DOA estimation effect was not satisfactory under the condition of low signal-to-noise ratio (SNR). In the literature [15,16,17,18,19], a new method with a combination of time reversal technology and a generalized MUSIC was proposed for the DOA estimation of low-altitude targets, which improved the estimation accuracy of the algorithm to some extent.

Due to the absence of finite second-order moments in impulse noise, conventional higher-order cumulant algorithms cannot be used. For the purpose of direction finding under impulse noise, in the literature [20], a theory based on fractional lower-order statistics (FLOS) was first proposed. Based on this theory, Mendel et al. [21,22,23] proposed the robust covariation-based MUSIC (ROC-MUSIC) algorithm in combination with the idea of co-variation. In the literature [24,25,26], the FLOM-MUSIC algorithm based on the fractional lower-order moment (FLOM) was proposed. In the literature [27,28,29], a fractional lower-order cyclic MUSIC (FLOC-MUSIC) algorithm using a fractional lower-order cyclic covariance matrix was proposed. However, these FLOS-based algorithms require prior knowledge of the characteristic exponent of stable distribution. To this end, the infinite MUSIC (IN-MUSIC) algorithm for received data was proposed in the literature [30], which was robust for outliers, but it required large samples to obtain a satisfactory performance. In order to solve the problem of direction finding of coherent signals in the background of impulse noise, the literature [31] proposed an algorithm based on infinite norm normalization preprocessing and sparse representation. The algorithm does not need to predict the number of signal sources, and it can be well adapted to coherent sources. In the literature [32], suppression of impulse noise is achieved by using an infinite norm exponential kernel covariance matrix. In order to obtain the global optimal solution of the method, a quantum electromagnetic field optimization algorithm is designed, and the Monte Carlo experiment verifies the effectiveness of the method.

In order to improve the direction-finding accuracy of low-altitude targets under the background of impulse noise, a DOA estimation method based on a spatial difference algorithm and a modified projection subspace algorithm is proposed in this paper. This method restores the covariance matrix of the received signal to full rank through a spatial difference algorithm. The signal source vector is estimated according to the least squares criterion, whereby the modified projection matrix is constructed and the cross-covariance matrix of the signal and the noise is then calculated. The sample covariance matrix is modified to remove interference terms. Next, the best correction coefficient is obtained with the maximum likelihood criterion, followed by re-adjustment of the covariance matrix. Finally, DOA estimation is conducted using the subspace classification algorithm. Simulation results verify the effectiveness of the proposed method.

The advantages of the proposed method can be given as follows:

Compared with the algorithm based on fractional low-order moments, the algorithm in this paper does not require prior knowledge of the characteristic exponent of stable distribution, and it is more adaptable to the environment;
Compared with the algorithm in the literature [27,28,29,30,31], the algorithm in this paper does not require large sample sampling and loop iterations, so the computational complexity is low;
Compared with the conventional low-altitude target DOA estimation method that only deals with rank deficiencies, the proposed algorithm has better performance in low-altitude environment.

2. Signal Model in Low-Altitude Multipath Environment

The signal reception model of the array in the low-altitude multipath environment is shown in the following figure:

Suppose K far-field narrowband signals

s_{i} (t) (i = 1, 2, \dots, K)

are incident on a uniform linear array with an angle of

θ_{i}

. The spacing d of the array elements is one half of the signal wavelength. The array consists of M uniformly distributed array elements, where each element is omnidirectional. Then, the signal steering vector of the array can be expressed as:

A = [a_{d 1} + a_{r 1}, \dots, a_{d K} + a_{r K}]

(1)

where

a_{d i} (θ_{d i}) = {[1, e^{- j π \sin θ_{d i}}, \dots, e^{- j π (M - 1) \sin θ_{d i}}]}^{T}

and

a_{r i} (θ_{r i}) = δ_{i} a_{d i} (θ_{r i})

. As shown in Figure 1,

θ_{d i}

and

θ_{r i}

are the incident angles after direct incidence and reflection, respectively;

δ_{i}

is the reflection coefficient; (·)

^{T}

is the transposition of the matrix.

Due to the characteristics of low-altitude multipath effects,

θ_{d i} = - θ_{r i} = θ_{i}

. Therefore, the signal steering vector can be expressed as:

A = [a (θ_{1}) + δ_{1} a (- θ_{1}), \dots, a (θ_{K}) + δ_{K} a (- θ_{K})]

(2)

The received signals of array can be expressed as:

y (t) = A s (t) + n (t) = B x (t) + n (t)

(3)

where

s (t) = {[s_{1} (t), s_{2} (t), \dots, s_{K} (t)]}^{T}

;

B = {[a (θ_{1}), a (- θ_{1}), \dots, a (θ_{K}), a (- θ_{K})]}_{N \times 2 K}

;

x (t) = {[s_{1} (t), δ_{1} s_{1} (t), \dots, s_{K} (t), δ_{K} s_{K} (t)]}_{2 K \times 1}^{T}

; and

n (t)

is the impulse noise vector received by the array.

According to the above equations, it can be concluded that under ideal conditions, the covariance matrix of the received signal of array is:

R = E [y (t) y^{H} (t)] = B R_{x} B^{H} + E [n (t) n^{H} (t)] = B R_{x} B^{H} + R_{n} = B R_{x} B^{H} + R_{n}

(4)

where

R_{x}

is the covariance matrix of

x (t)

;

R_{n}

is the covariance matrix of noise

n (t)

; and (·)

^{H}

is the conjugate transpose of the matrix.

Equation (4) can be established only when the number of snapshots approaches infinity. In actual cases, after

y (t)

is expanded to

B x (t) + n (t)

, as shown in Equation (5) where the covariance of each term is calculated, it is found that the signal–noise cross-covariance calculated matrix

R_{xn}

and the noise–signal cross-covariance matrix

R_{nx}

are non-zero matrices:

\begin{array}{l} R = \frac{1}{N} \sum_{t = 1}^{N} y (t) y^{H} (t) = \frac{1}{N} \sum_{t = 1}^{N} (B x (t) + n (t)) {(B x (t) + n (t))}^{H} \\ = B [\frac{1}{N} \sum_{t = 1}^{N} x (t) x^{H} (t)] B^{H} + \frac{1}{N} \sum_{t = 1}^{N} n (t) n^{H} (t) + B [\frac{1}{N} \sum_{t = 1}^{N} x (t) n^{H} (t)] + [\frac{1}{N} \sum_{t = 1}^{N} x (t) n^{H} (t)] B^{H} \\ = R_{xx} + R_{nn} + R_{xn} + R_{nx} \end{array}

(5)

where N is the number of snapshots;

R_{xx}

is the signal covariance matrix; and

R_{nn}

is the noise covariance matrix.

3. DOA Estimation of Low-Altitude Targets under Impulse Noise

3.1. Decoherence with Spatial Difference Algorithm

In the low-altitude multipath environment, the direct signal and the reflected signal contained in received signals of array are correlated, as a result of which there will be angle loss in case of direct feature decomposition of

R

. Therefore, we use the spatial difference algorithm to divide the received signals

y (t)

of array into L sub-arrays

y_{l} (t)

that overlap each other, and each sub-array contains m array elements:

y_{l} (t) = B_{m} V^{l - 1} x (t) + n_{l} (t)

(6)

where

l = 1, 2, \dots, L ≜ M - m + 1

;

B_{m}

is the first m rows of B; and

n_{l} (t)

represents the noise received by the l-th subarray. Furthermore:

V = d i a g [v (θ_{1}), v (- θ_{1}), \dots, v (θ_{K}), v (- θ_{K})]

(7)

where

v (θ_{i}) = e^{- j π \sin θ_{i}}

;

i = 1, 2, \dots, K

The covariance matrix of submatrix

y_{l} (t)

is calculated. With Equation (5), the covariance submatrix is obtained:

\begin{array}{l} {\hat{R}}_{l} = \frac{1}{N} \sum_{t = 1}^{N} [y_{l} (t) y_{l}^{H} (t)] \\ = \frac{1}{N} \sum_{t = 1}^{N} [[B_{m} V^{l - 1} x (t) + n_{l} (t)] {[B_{m} V^{l - 1} x (t) + n_{l} (t)]}^{H}] \\ = R_{xx}^{l} + R_{nn}^{l} + R_{xn}^{l} + R_{nx}^{l} \end{array}

(8)

The spatial difference matrix consists of the difference between the first covariance submatrix and the backward covariance matrix [13]:

\hat{R} = \frac{1}{L} \underset{l = 1}{\sum^{L}} [{\hat{R}}_{1} - J {({\hat{R}}_{l})}^{*} J]

(9)

where

J

is the

m \times m

-dimensional anti-diagonal matrix; the elements on the anti-diagonal are all 1; the rest of the elements are 0; and (·)

^{*}

is the conjugate of the matrix.

Since it is difficult to estimate all target angles due to the influence of low-altitude multipath effects, the spatial difference algorithm is used to eliminate this influence. However, the impulse noise component in the covariance matrix still greatly affects the estimation accuracy.

3.2. Modified Projective Subspace Algorithm

As the cross-covariance component in the covariance matrix

\hat{R}

greatly reduces the accuracy of DOA estimation, we need to first estimate

x (t)

to remove the influence of these disturbance terms:

\hat{x} (t) = \arg \min_{y} {‖y (t) - \hat{B} (θ) x (t)‖}_{2}^{2}

(10)

where

\hat{x} (t)

is the estimated value of

x (t)

and

\hat{B} (θ)

is the estimated value of

B (θ)

. For the convenience of expression,

B (θ)

is denoted as

B

According to the least squares rule, formula (10) is minimized to obtain:

\hat{x} (t) = {({\hat{B}}^{H} \hat{B})}^{- 1} {\hat{B}}^{H} y (t)

(11)

Taking the difference between the observed value of the array signal and the estimated signal as an estimate of noise, we can get:

\hat{n} (t) = y (t) - \hat{B} \hat{x} (t)

(12)

Then, the estimated value of the third term in Equation (8) can be obtained:

\begin{array}{l} {\hat{R}}_{xn}^{l} = {\hat{B}}_{m} [\frac{1}{N} \sum_{t = 1}^{N} \hat{x} (t) {\hat{n}}^{H} (t)] \\ = {\hat{B}}_{m} [\frac{1}{N} \sum_{t = 1}^{N} {({\hat{B}}_{m}^{H} {\hat{B}}_{m})}^{- 1} {\hat{B}}_{m}^{H} y (t) (y^{H} (t) - y^{H} (t) {\hat{B}}_{m} {({\hat{B}}_{m}^{H} {\hat{B}}_{m})}^{- 1} {\hat{B}}_{m}^{H})] \\ = {\hat{B}}_{m} {({\hat{B}}_{m}^{H} {\hat{B}}_{m})}^{- 1} {\hat{B}}_{m}^{H} [\frac{1}{N} \sum_{t = 1}^{N} y (t) y^{H} (t) (I_{M N} - {\hat{B}}_{m} {({\hat{B}}_{m}^{H} {\hat{B}}_{m})}^{- 1} {\hat{B}}_{m}^{H})] \end{array}

(13)

where

{\hat{B}}_{m}

is the first m row of

\hat{B}

In order to facilitate subsequent calculation and representation, the following projection matrix is assumed:

{\hat{P}}_{B} = {\hat{B}}_{m} {({\hat{B}}_{m}^{H} {\hat{B}}_{m})}^{- 1} {\hat{B}}_{m}^{H}

(14)

{\hat{P}}_{B}^{⊥} = I_{M N} - {\hat{P}}_{B}

(15)

After simplification:

{\hat{R}}_{xn}^{l} = {\hat{P}}_{B} \hat{R} {\hat{P}}_{B}^{⊥}

(16)

Similarly, we can get:

{\hat{R}}_{nx}^{l} = {\hat{P}}_{B}^{⊥} \hat{R} {\hat{P}}_{B}

(17)

{\hat{R}}_{nn}^{l} = {\hat{P}}_{B}^{⊥} \hat{R} {\hat{P}}_{B}^{⊥}

(18)

Because the steering vector matrix

\hat{B}

is unpredictable,

{\hat{P}}_{B}

and

{\hat{P}}_{B}^{⊥}

cannot be calculated. Therefore, we perform eigen decomposition of the covariance matrix

\hat{R}

after spatial differentiation:

\hat{R} = {\hat{U}}_{x} {\hat{Σ}}_{x} {\hat{U}}_{x}^{H} + {\hat{U}}_{n} {\hat{Σ}}_{n} {\hat{U}}_{n}^{H}

(19)

where

{\hat{Σ}}_{x} = d i a g [λ_{1}, \dots, λ_{K}]

{\hat{Σ}}_{n} = d i a g [λ_{K + 1}, \dots, λ_{M}]

and

λ_{1} \geq \dots \geq λ_{K} > λ_{K} = λ_{K + 1} = \dots = λ_{M}

;

{\hat{U}}_{x} = [u_{1}, u_{2}, \dots u_{K}]

{\hat{U}}_{n} = [u_{K + 1}, u_{K + 2}, \dots u_{M}]

u_{i}

is the eigenvector corresponding to the eigenvalue

λ_{i}

Upon calculation, it is found that the subspace generated by the signal covariance matrices

{\hat{U}}_{x}

and

B_{m}

is the same, that is:

span \{u_{1}, u_{2}, \dots, u_{K}\} = span \{b_{1}, b_{2}, \dots, b_{K}\}

(20)

Hence, by assuming that there is a matrix

T

that makes

{\hat{U}}_{x} T = B_{m}

established, we get:

\begin{array}{l} {\hat{P}}_{B} = {\hat{B}}_{m} {({\hat{B}}_{m}^{H} {\hat{B}}_{m})}^{- 1} {\hat{B}}_{m}^{H} \\ = {\hat{U}}_{x} T {[{({\hat{U}}_{x} T)}^{H} {\hat{U}}_{x} T]}^{- 1} {({\hat{U}}_{x} T)}^{H} \\ = {\hat{U}}_{x} {({\hat{U}}_{x}^{H} {\hat{U}}_{x})}^{- 1} {\hat{U}}_{x}^{H} \end{array}

(21)

So,

{\hat{P}}_{B}^{⊥}

can be represented as:

{\hat{P}}_{B}^{⊥} = I_{M N} - {\hat{P}}_{B}

(22)

Formulas (21) and (22) are substituted into Formulas (17) and (18) to get:

{\hat{R}}_{xn}^{l} = {\hat{P}}_{B}^{⊥} \hat{R} {\hat{P}}_{B}

(23)

{\hat{R}}_{nx}^{l} = {\hat{P}}_{B}^{⊥} \hat{R} {\hat{P}}_{B}^{⊥}

(24)

{\hat{R}}_{l}

and

\hat{R}

after removal of disturbance terms can be expressed as:

{\hat{R}}_{l}^{'} = [{\hat{R}}_{l} - ε ({\hat{R}}_{x n}^{l} + {\hat{R}}_{n x}^{l})] {[{\hat{R}}_{l} - ε ({\hat{R}}_{x n}^{l} + {\hat{R}}_{n x}^{l})]}^{H}

(25)

{\hat{R}}^{'} = \frac{1}{L} \underset{l = 1}{\sum^{L}} [{\hat{R}}_{1}^{'} - J {({\hat{R}}_{l}^{'})}^{*} J]

(26)

in Equation (25),

ε \in [0, 1]

is used to correct the estimation error. If

{\hat{R}}_{x n}^{l}

{\hat{R}}_{n x}^{l}

are equal to

R_{x n}^{l}

R_{n x}^{l}

, respectively,

ε = 1

. In practice, since the error of simulation estimation cannot be accurately predicted, a specific uniform step size is substituted into the following equation for calculation:

f (ε) = \min_{ε} \ln \det ({\hat{P}}_{B}^{'} \hat{R} {\hat{P}}_{B}^{'} + \frac{Tr {{\hat{P}}_{B}^{⊥^{'}} \hat{R}}}{m - K} {\hat{P}}_{B}^{⊥^{'}})

(27)

according to the maximum likelihood criterion,

ε

is the best correction coefficient when

f (ε)

reaches the minimum value. In this equation,

{\hat{P}}_{B}^{'}

is the modified projection matrix of

{\hat{R}}^{'}

obtained by eigenvalue decomposition;

{\hat{P}}_{B}^{⊥}^{'} = I_{M N} - {\hat{P}}_{B}^{'}

and

Tr \{{\hat{P}}_{B}^{⊥}^{'} \hat{R}\}

is the trace of the matrix.

Finally, the DOA of the target is obtained through eigenvalue decomposition of

{\hat{R}}^{'}

using the MUSIC algorithm.

4. Basic Steps of the Algorithm and Complexity Analysis

4.1. The Basic Steps of the Algorithm

According to the above analysis, the algorithm in this paper is summarized as follows:

Calculate the data covariance matrix $R$ of the array element output vector $y (t)$ ;
As shown in Equations (6)–(9), spatial difference operation is performed to obtain $\hat{R}$ ;
Eigenvalue decomposition of $\hat{R}$ is performed to get ${\hat{U}}_{x}$ ;
Calculate the modified projection matrix ${\hat{P}}_{B}$ and ${\hat{P}}_{B}^{⊥}$ with Equations (21) and (22);
As shown in Equations (23)–(26), the cross-covariance matrices of ${\hat{R}}_{xn}^{l}$ and ${\hat{R}}_{nx}^{l}$ are constructed to correct the estimated value of ${\hat{R}}_{l}$ ;
As shown in Equation (27), the optimal correction coefficient $ε$ is obtained using the maximum likelihood criterion, and the estimated value of ${\hat{R}}_{l}$ is re-adjusted;
Finally, DOA estimation is performed for the adjusted ${\hat{R}}_{l}^{'}$ using MUSIC algorithm.

4.2. Algorithm Complexity Analysis

Let’s assume that the number of array elements is M; the number of signal sources is K; the number of sampling points is N; the number of angle searches is

N_{θ}

, and E is the number of

ε

divided by a uniform step size. The computational complexity of steps 1~5 of this algorithm is approximately

(E + 2 N + 1) {(M)}^{3} + (E N + K) {(M)}^{2} + K^{2} M + K^{3}

. The complexity of spatial spectrum calculation and search in steps 6~7 is approximately

O ((M - K) M N_{θ} (E + 1))

. Therefore, the total complexity of the algorithm herein is approximately

O (\begin{array}{l} (E + 2 N + 1) {(M)}^{3} + (E N + E N_{θ} + N_{θ} + K) {(M)}^{2} \\ + (K^{2} - K E N_{θ} - K N_{θ}) M + K^{3} \end{array})

5. Simulation Results and Analysis

5.1. Impulse Noise Model

Compared with Gaussian noise, the probability distribution of impulse noise is characterized by long tailing. Typical impulse noise models include a Gaussian mixture model (GMM), noise model, and

S α S

noise model.

The GMM noise model used in the simulation consists of two Gaussian components, and its probability density function (PDF) is:

p_{n} (x) = \underset{i = 1}{\sum^{2}} \frac{c_{i}}{π σ_{i}^{2}} \exp (- \frac{| x |^{2}}{σ_{i}^{2}})

(28)

where

i = 1, 2

0 \leq c_{i} \leq 1

is the probability of the i-th term and

c_{1} + c_{2} = 1

;

σ_{i}^{2}

is the variance of the i-th term; it is set in the simulation that

σ_{2}^{2} = 100 σ_{1}^{2}

. Meanwhile, the signal-to-noise ratio (SNR) is defined as:

S N R = \frac{σ_{s}^{2}}{σ_{1}^{2}}

(29)

the noise power should have been

c_{1} σ_{1}^{2} + c_{2} σ_{2}^{2}

, but since

σ_{1}^{2}

is set to a large noise background in this simulation, the presence of

σ_{2}^{2}

is considered an outlier.

As the most commonly used impulse noise distribution model, the

S α S

noise model PDF cannot usually be given in closed form, but its characteristic function is:

ϕ (x) = \exp (- γ^{α} | x |^{α})

(30)

where

α

is the characteristic exponent and

0 < α \leq 2

. The larger the

α

, the weaker the pulse characteristic. When

α = 2

, it degenerates into a Gaussian distribution. The dispersion coefficient is

γ

, and it is similar to the variance in a Gaussian distribution. Since the common SNR is nonsense for the

S α S

distribution with

α < 2

, the generalized SNR (GSNR) is used:

G S N R = \frac{σ_{s}^{2}}{γ^{α}}

(31)

A uniform linear array with 35 array elements is used in the simulation. Each subarray contains 30 array elements, that is, L = 6. The angles of two incoherent narrowband signals are 2° and 6°. The signal wavelength

λ = \frac{d}{2}

; the array element spacing d = 2; the source power is equal to

σ_{s}^{2}

; the reflection coefficient

δ_{i}

= 0.9; and the number of snapshots is 200.

5.2. Spatial Spectrum Estimation

Experiment one shows the spatial spectrum estimation effect with the simulation algorithm under the GMM noise model. Figure 2 shows the spatial spectrum of three independent experiments with the algorithm, when the probability of outliers is

c_{2}

= 0.1,

c_{2}

= 0.3, and SNR = 5 dB.

Figure 2 shows that under the background of GMM noise, when the probability of outliers is 0.1, the target direction can be accurately estimated with the algorithm. When the probability of outliers rises to 0.3, although specific spectral peaks are not steep, the algorithm still maintains a high resolution.

Experiment two shows the spatial spectrum estimation effect with the simulation algorithm under the

S α S

noise model. Figure 2 shows the spatial spectrum of three independent experiments with the algorithm, when the characteristic exponent

α = 1.6

and

α = 1.3

, where GSNR = 5 dB.

Figure 3 shows that under the background of

S α S

noise, when the characteristic exponent

α

decreases, the impact characteristic of the noise is enhanced and the position of the spectral peak deviates from the target direction to a certain extent. It can be learned from the spectral peak search that the algorithm can still maintain a high accuracy and a strong robustness.

5.3. Comparative Analysis of DOA Estimation Performance

In this section, simulation experiments are carried out to study the performance of this algorithm. Since the MUSIC algorithm fails to effectively identify coherent signal sources, the spatial smoothing method [11] is used in the following comparison algorithms. The comparison algorithms being compared include ROC-SSMUSIC [22], FLOM-SSMUSIC [24], FLOC-SSMUSIC [27], and IN-SSMUSIC [31].

The root mean square error (RMSE) of DOA estimation is defined as:

RMSE = \sqrt{1 / 2 N_{m o n t} \underset{i = 1}{\sum^{N_{m o n t}}} \underset{k = 1}{\sum^{K}} ({|θ_{i} - θ_{i k}|}^{2})}

(32)

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 actual value and the estimated value of the i-th signal azimuth in each experiment.

A successful estimation of each angle is defined as:

|θ_{i} - θ_{i k}| \leq 2 °, k = 1, 2, \dots K

(33)

namely, the difference between the estimated value of the angle and the actual value is less than 2°. The probability of success is the ratio of the number of successes to the total number of estimations.

In experiment three, the number of Monte Carlo experiments was 500, and other settings were the same as above. Different SNRs were set to compare the estimation performance under the background of GMM noise. Figure 4a shows the curve of the changes of RMSE with five different algorithms, when the SNR changes from −5 dB to 20 dB at a step of 5 dB. Figure 4b shows the curves of the success rate of the five algorithms with changes of SNR.

As can be seen from Figure 4, the theoretical method of edge asymptotic fitting of the covariant matrix adopted in the ROC-SSMUSIC algorithm is not very effective. Due to the difference between the impulse noise model and the selection of noise parameters, the estimation performance of the FLOM-SSMUSIC algorithm and the FLOC-SSMUSIC algorithm is significantly improved when the SNR is greater than −5 dB. With pre-treatment of infinite norm normalization to reduce the influence of impulse noise, the success rate of the IN-SSMUSIC algorithm is closest to the algorithm in this paper when the SNR is high. The proposed algorithm basically achieves the best estimation performance within the entire SNR simulation range, and especially when the SNR is greater than 0 dB, the estimation performance herein is significantly better than other algorithms.

In experiment four, while the simulation settings in this section remain unchanged, different GSNRs were set under the background of

S α S

noise to compare the algorithm’s performance, where

α

= 1.4. Figure 5a shows the curve of changes of RMSE with the five different algorithms when the GSNR changes from −5 dB to 20 dB at a step of 5 dB. Figure 5b shows the curves of the success rate of the five algorithms with changes of GSNR.

It can be seen from Figure 5 that with increased GSNR, the impact of impulse noise decreases and the success rate of DOA estimation is significantly improved. Low-order fractional moment algorithms such as FLOM-SSMUSIC and FLOC-SSMUSIC achieve better performance at low GSNR, but such algorithms depend too much on prior knowledge of noise. The performance of the proposed algorithm is better than the other algorithms being compared. However, in the case of a low GSNR, due to the strong impulse noise impact, its estimation accuracy decreases to a certain extent.

In Experiment five, the influence of the number of snapshots on the performance of the algorithm was analyzed. The performance of the algorithm was compared by setting a different number of snapshots under two different noise models. Figure 6a shows the curve of changes of RMSE with the five different algorithms when the number of snapshots changes from 50 to 300 at a step of 50 under the background of GMM noise with SNR = 5 dB. Figure 6b shows the change of RMSE with the five algorithms when the number of snapshots changes from 50 to 300 at a step of 50 under the background of

S α S

noise with GSNR = 5 dB and

α = 1.4

It can be seen from Figure 6 that with increased snapshots, the estimation performance of each algorithm is enhanced to a certain extent either under a GMM or an

S α S

noise model, and it gradually becomes stable. Regardless of the number of snapshots, the RMSE of the proposed algorithm is lower than other algorithms. Under the GMM noise, however, when the number of snapshots is less than 100, the performance of the proposed algorithm is poor. This is because when the number of snapshots is small, it is difficult for the sample-based estimation algorithm to correctly separate the signal and the noise subspaces, resulting in a larger error.

In Experiment six, the influence of the characteristic exponent

α

on the performance of the algorithm under the

S α S

noise background was analyzed. While the simulation in this section remains unchanged, different characteristic exponents are set. Figure 7 shows the curves of the success rate of five different algorithms when the characteristic exponent changes from 0.6 to 2 at a step of 0.2 when GSNR = 10 dB and GSNR = 20 dB.

Since the larger the characteristic exponent, the weaker the impulsiveness of the noise, the success rate of these algorithms is improved with increased characteristic exponent. It can be seen from the figure that the performance of ROC-SSMUSIC algorithm is the worst in an

S α S

noise environment, without a significant effect on the change of the characteristic exponent. The FLOM-SSMUSIC, FLOC-SSMUSIC, and IN-SSMUSIC algorithms are sensitive to the change of the impulse noise intensity to some extent, that is, the reduction of the noise intensity reduces RMSE. When the characteristic exponent of impulse noise contained in the output signal of the array changes from 1.6 to 2, the above three algorithms are insensitive to the change of impulse noise intensity, and they show strong robustness. However, the algorithm herein maintains a strong estimation performance in a 10 dB or a 20 dB

S α S

noise environment.

6. Conclusions

As to direction finding of low-altitude targets under the background of impulse noise, a DOA estimation method based on spatial difference and on a modified projection subspace algorithm is proposed. In this method, the covariance matrix is first constructed with the received signal vector, and then the spatial difference operation is performed to eliminate the influence of low-altitude multipath. Next, the array signal vector is estimated according to the least squares criterion, after which the modified projection matrix is constructed to calculate the cross-covariance matrix of the signal noise. Then, the covariance matrix of the samples is modified to remove the disturbance terms, thus completing the DOA estimation of the targets. The experimental results show that the proposed algorithm achieves satisfactory DOA estimation of low-altitude targets under different impulse noise backgrounds, with higher accuracy and more practical application compared with conventional algorithms. In the next step, on the basis of guaranteeing the performance of the algorithm, the extended study of the spatial difference method for coherent and incoherent mixed signals in the background of impulse noise will be carefully considered.

Author Contributions

This article was co-authored by B.L., G.H., H.Z., G.Z. and Y.S. 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. and Y.S.; 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 being obtained from the simulation of the signal models listed in the article and not from publicly available datasets.

Acknowledgments

We thank the college for providing us with an efficient simulation platform so that we could 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

Schmidt, R.; Schmidt, R.O. Multiple Emitter Location and Signal Parameter Estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef] [Green Version]
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] [PubMed] [Green Version]
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] [PubMed]
Zheng, G.; Song, Y. Signal Model and Method for Joint Angle and Range Estimation of Low-Elevation Target in Meter-Wave FDA-MIMO Radar. IEEE Commun. Lett. 2022, 26, 449–453. [Google Scholar] [CrossRef]
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] [Green Version]
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]
Zheng, G.; Song, Y.; Chen, C. Height Measurement with Meter Wave Polarimetric MIMO Radar: Signal Model And MUSIC-like Algorithm. Signal Process. 2021, 190, 1872–7557. [Google Scholar] [CrossRef]
Tan, J.; Nie, Z. Polarization Smoothing Generalized MUSIC Algorithm with Polarization Sensitive Array for Low Angle Estimation. Sensors 2018, 18, 1534. [Google Scholar] [CrossRef] [Green Version]
Liu, Q.; Gu, Y.; So, H. DOA Estimation in Impulsive Noise Via Low-Rank matrix Approximation and Weakly Convex OptimiZation. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 3603–3616. [Google Scholar] [CrossRef]
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]
Li, J. Improved Angular Resolution for Spatial Smoothing Techniques. Signal Process. 1992, 40, 3078–3081. [Google Scholar] [CrossRef]
Ma, X.; Dong, X.; Xie, Y. An Improved Spatial Differencing Method for DOA Estimation with The Coexistence of Uncorrelated and Coherent Signals. IEEE Sens. J. 2016, 16, 3719–3723. [Google Scholar] [CrossRef]
Liu, F.; Wang, J.; Sun, C. Spatial Differencing Method for DOA Estimation under The Coexistence of Both Uncorrelated and Coherent Signals. IEEE Trans. Antennas Propag. 2012, 60, 2052–2062. [Google Scholar] [CrossRef]
Al-Ardi, E.M.; Shubair, R.M.; Al-Mualla, M.E. Computationally Efficient High-Resolution DOA Estimation in Multipath Environment. Electron. Lett. 2014, 40, 908–910. [Google Scholar] [CrossRef]
Zeng, X.L.; Lei, M. Estimation of Direction of Arrival by Time Reversal for Low-Angle Targets. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 2675–2694. [Google Scholar] [CrossRef]
Xu, H.; Cui, W.; Du, Y. The Analysis of Using Spatial Smoothing for DOA Estimation of Coherent Signals in Sparse Arrays. Int. J. Antennas Propag. 2021, 2021, 4447562. [Google Scholar] [CrossRef]
Qi, B.; Liu, D. An Enhanced Spatial Smoothing Algorithm for Coherent Signals DOA Estimation. Eng. Comput. 2021, 39, 574–586. [Google Scholar] [CrossRef]
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]
Meng, F.X.; Yu, X.T.; Zhang, Z.C. Quantum Algorithm for Multiple Signal Classification. Phys. Rev. A 2021, 101, 23–34. [Google Scholar] [CrossRef]
Nikias, C.L.; Shao, M. Signal Processing with Alpha-Stable Distributions and Applications. Comput. Stat. Data Anal. 1996, 22, 20–55. [Google Scholar]
Li, L.; Qiu, T.S.; Song, D.R. Parameter Estimation Based on Fractional Power Spectrum under Alpha-Stable Distribution Noise Environment in Wideband Bistatic MIMO Radar System. AEU Int. J. Electron. Commun. 2013, 67, 947–954. [Google Scholar] [CrossRef]
Mendel, J.M. A Subspace-Based Direction Finding Algorithm Using Fractional Lower Order Statistics. IEEE Trans. Signal Process. 2001, 49, 1605–1613. [Google Scholar]
Tzagkarakisa, G.; Tsakalidesb, P.; Starcka, J.-L. Covariation-Based Subspace-Augmented MUSIC for Joint Sparse Support Recovery in Impulsive Environments. Signal Process. 2013, 93, 1365–1373. [Google Scholar] [CrossRef]
Hacioglu, I.; Harmanci, F.K.; Anarim, E. Time-of-Arrival Estimation with FLOM-MUSIC under Impulsive Noise. In Proceedings of the IEEE 14th Signal Processing & Communications Applications, Antalya, Türkiye, 16–19 April 2006. [Google Scholar]
Diao, M.; Yang, L.L.; Chen, C. DOA Estimation of Spatial-Temporal and Maximum-Likelihood Based on Fractional Lower-Order Statistics. Appl. Sci. Technol. 2010, 37, 15–18. [Google Scholar]
Li, S.; He, R.; Lin, B. DOA Estimation Based on Sparse Representation of The Fractional Lower Order Statistics in Impulsive Noise. IEEE J. Autom. Sin. 2018, 5, 860–868. [Google Scholar] [CrossRef]
Wang, H.B.; Long, J.B.; Zhou, Y.X. Frequency Spectrum Estimate Based on Decomposed FLOC Matrix under Stable Distribution Noise. Appl. Mech. Mater. 2014, 556–562, 4563–4567. [Google Scholar] [CrossRef]
Li, S.; Chen, X.; He, R. Robust Cyclic MUSIC Algorithm for Finding Directions in Impulsive Noise Environment. Int. J. Antennas Propag. 2017, 2017, 9038341. [Google Scholar] [CrossRef] [Green Version]
Liu, Y.; Wu, Q.; Zhang, Y. Cyclostationarity-Based DOA Estimation Algorithms for Coherent Signals in Impulsive Noise Environments. EURASIP J. Wirel. Commun. Netw. 2019, 2019, 81. [Google Scholar] [CrossRef]
Wang, J.T.; Sheng, J.; Jin, H.E. Subspace-Based DOA Estimation in Impulsive Noise Environments for MIMO Radars. J. Astronaut. 2009, 53, 89–98. [Google Scholar]
Gong, J.; Guo, Y. A Bistatic MIMO Radar Angle Estimation Method for Coherent Sources in Impulse Noise Background. Wirel. Pers. Commun. 2020, 116, 3567–3576. [Google Scholar] [CrossRef]
Yanan, D.U.; Gao, H.; Chen, M. Direction of Arrival Estimation Method Based on Quantum Electromagnetic Field Optimization in The Impulse Noise. Syst. Eng. Electron. Technol. Engl. Ed. 2021, 32, 11. [Google Scholar] [CrossRef]

Figure 1. Signal reception model in low-altitude multipath environment.

Figure 2. Algorithm space spectrum under the background of GMM noise: (a)

c_{2}

= 0.1; (b)

c_{2}

= 0.3.

Figure 2. Algorithm space spectrum under the background of GMM noise: (a)

c_{2}

= 0.1; (b)

c_{2}

= 0.3.

Figure 3. Spatial spectrum under noise background with the algorithm: (a)

α = 1.6

; (b)

α = 1.3

Figure 3. Spatial spectrum under noise background with the algorithm: (a)

α = 1.6

; (b)

α = 1.3

Figure 4. Analysis of mean square error and success rate under GMM noise: (a) Relationship between RMSE and SNR; (b) Relationship between success rate and SNR.

Figure 5. Analysis of the mean square error and success rate of the algorithm under noise: (a) Relationship between RMSE and GSNR; (b) Relationship between success rate and GSNR.

Figure 6. Influence of the number of snapshots on the performance of the algorithm: (a) GMM noise; (b)

S α S

noise.

Figure 6. Influence of the number of snapshots on the performance of the algorithm: (a) GMM noise; (b)

S α S

noise.

Figure 7. Influence of characteristic exponent on algorithm performance: (a) GSNR = 10 dB; (b) GSNR = 20 dB.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Lin, B.; Hu, G.; Zhou, H.; Zheng, G.; Song, Y. The DOA Estimation Method for Low-Altitude Targets under the Background of Impulse Noise. Sensors 2022, 22, 4853. https://doi.org/10.3390/s22134853

AMA Style

Lin B, Hu G, Zhou H, Zheng G, Song Y. The DOA Estimation Method for Low-Altitude Targets under the Background of Impulse Noise. Sensors. 2022; 22(13):4853. https://doi.org/10.3390/s22134853

Chicago/Turabian Style

Lin, Bin, Guoping Hu, Hao Zhou, Guimei Zheng, and Yuwei Song. 2022. "The DOA Estimation Method for Low-Altitude Targets under the Background of Impulse Noise" Sensors 22, no. 13: 4853. https://doi.org/10.3390/s22134853

APA Style

Lin, B., Hu, G., Zhou, H., Zheng, G., & Song, Y. (2022). The DOA Estimation Method for Low-Altitude Targets under the Background of Impulse Noise. Sensors, 22(13), 4853. https://doi.org/10.3390/s22134853

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

Article Menu

The DOA Estimation Method for Low-Altitude Targets under the Background of Impulse Noise

Abstract

1. Introduction

2. Signal Model in Low-Altitude Multipath Environment

3. DOA Estimation of Low-Altitude Targets under Impulse Noise

3.1. Decoherence with Spatial Difference Algorithm

3.2. Modified Projective Subspace Algorithm

4. Basic Steps of the Algorithm and Complexity Analysis

4.1. The Basic Steps of the Algorithm

4.2. Algorithm Complexity Analysis

5. Simulation Results and Analysis

5.1. Impulse Noise Model

5.2. Spatial Spectrum Estimation

5.3. Comparative Analysis of DOA Estimation Performance

6. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI