A Robust DOA Tracking Method Using a Nested Array in Impulse Noise

Wang, Yufeng; Zhao, Dasheng; Zhang, Yidu; Wang, Tiebin

doi:https://doi.org/10.1155/2022/6227794

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest References Copyright Related Articles

Special Issue

Potential Problems and Recent Advances in Direction of Arrival Estimation

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 6227794 | https://doi.org/10.1155/2022/6227794

A Robust DOA Tracking Method Using a Nested Array in Impulse Noise

Yufeng Wang,¹Dasheng Zhao,²Yidu Zhang,²and Tiebin Wang³

Academic Editor: Yanan Du

Received06 Jun 2022

Accepted14 Jul 2022

Published09 Aug 2022

Abstract

To achieve the underdetermined direction of arrival (DOA) tracking in impulse noise, we propose a robust DOA tracking method in this paper. Firstly, an infinite norm difference covariance (INDC) matrix is introduced to suppress the impulse noise. Secondly, using a nested array, we present a maximum likelihood DOA tracking equation based on the INDC matrix. Furthermore, a quantum-inspired multiverse algorithm is proposed to maximize efficiently the proposed tracking equation. The simulation results show that our method has better robustness and superiority compared to other DOA tracking methods, which can achieve underdetermined DOA tracking in impulse noise as well as in the multipath environment.

1. Introduction

In the past few decades, many high-resolution direction of arrival (DOA) estimation algorithms [1–3] have been proposed and applied in many fields such as radar and wireless communication, most of which assume that the targets are stationary, whereas the targets are usually time-varying in fact. For time-varying DOA tracking problem, many researchers have made great efforts for its development. In [4], the authors proposed a projection approximation subspace tracking (PAST) algorithm to track the subspace. In [5], a particle filter (PF) DOA tracker was proposed, which utilized a partitioned state-vector method to achieve multiple targets tracking. In [6], a fast approximated power iteration subspace tracking method was proposed, which offered a faster tracking response. In [7], on the basis of low-rank and sparse recovery, the authors proposed a novel DOA tracking method. In [8], a multisource DOA tracking approach using a superposition model was proposed, which utilized the PF to achieve the DOA tracking.

The above methods are derived in Gaussian noise, whereas in real world, there exists some impulse noises, including the lightning and the low-frequency atmospheric noise, which can be described well as symmetric α-stable (SαS) distribution [9], and these forms of impulse noise exist a long tail, which deteriorates the performance of the existing methods. Many methods, including fractional lower-order moment [10], correntropy technique [11, 12], infinite norm [13], and kernel method [14, 15], have been proposed to achieve accurate estimates in impulse noise. DOA tracking in impulse noise has also made great progress in recent years. In [16], the authors proposed a robust PAST (RPAST) algorithm for tracking subspace in impulse noise. In [17], a PF method for DOA tracking in impulse noise was proposed, which offered robustness in impulse noise but cannot obtain excellent performance in fast time-varying scenario. In [18], a correntropy method was proposed to achieve subspace tracking in impulse noise.

Furthermore, the above algorithms are based on uniform linear array (ULA), which will be invalid when the number of targets exceeds the number of antennas. Many nonuniform arrays have been proposed to achieve underdetermined DOA estimation [19, 20]. Some researchers also applied them in DOA tracking problems. In [21], a PF method using coprime array was proposed for DOA tracking, which was based on the propagator method. In [22], a spatial smoothing projection approximation subspace tracking (SSPAST) method was proposed, which utilized the difference coarray of nested array (NA) and coprime array to achieve underdetermined DOA tracking. In [23], a DOA tracking method using NA was proposed, which was based on offset compensation to obtain accurate tracking results.

The noted methods are only an improvement on one aspect of the DOA tracking problem; in order to address the above problems simultaneously, a robust DOA tracking method using nested array in impulse noise is proposed in our work. First, we propose an infinite norm difference covariance (INDC) matrix to obtain robustness in impulse noise, and on this basis, a maximum likelihood (ML) DOA tracking equation based on NA is utilized, which involves a multidimension joint optimization problem requiring enormous computational complexity. Fortunately, intelligent optimization algorithms, such as the sea lion optimization algorithm [24] and the Archimedes optimization algorithm [25], can allow for this cost function. However, the above algorithms are easy to trap in local optimum. Therefore, we propose a quantum-inspired multiverse algorithm (QMVA) to avoid this problem, which is inspired by quantum computation [26] and theory of cosmology [27]. The resulting method is termed as QMVA-INDC-ML-NA. The simulation results demonstrate that our method offers better robustness and effectiveness compared to other approaches.

The main contributions are as follows:(1)An INDC matrix is proposed to suppress the strong impulsive noise.(2)A ML tracking equation based on INDC, using NA, is proposed to achieve underdetermined DOA tracking in impulse noise.(3)A quantum-inspired multiverse algorithm is proposed to solve efficiently the considered ML tracking equation.

The following is the rest of this paper. The DOA tracking model in impulse noise is introduced in Section 2. The DOA tracking method based on the QMVA is presented in Section 3. The simulation results are shown in Section 4. Finally, the conclusions of our work are given in Section 5.

2. DOA Tracking Model in Impulse Noise

Assume that a nested array consists of N ULAs and M isotropic antennas and denotes the distance between the mth antenna and the first antenna, where , . The minimum spacing of the antenna elements is , and then the coordinates of the antenna elements arewhere are integers, and and denote the coordinates of the nth ULA and the first ULA, respectively, where and . A set is a continuous or nearly continuous set of natural numbers.

Consider that P narrowband signals with wavelength impinging on a nested array with M isotropic antenna elements thus the receiving vector is described aswhere is the array manifold, the pth steering vector is , , , is the signal vector, and is the complex impulse noise vector.

The characteristic function of zero-location SαS distribution is given bywhere and denote the characteristic exponent and the scale. In impulse noise, one usually use generalized signal-to-noise ratio (GSNR)where denote the expectation.

As conventional second-moment-based methods will deteriorate or even be invalid in impulse noise, we employ the infinite norm normalization preprocessing method, and on this basis, we propose an infinite norm difference covariance (INDC) matrix for the kth snapshot, and the INDC matrix is given bywhere

The ith row and jth column element of is , where ; denotes the difference constant, ; and denotes the number of snapshots.

Next, virtualize the INDC of the nested array into an extended INDC of a virtual ULA with more antenna elements, and virtualize the array manifold into a virtual array manifold. The maximum correlation delay calculated by the nested array is , the number of antenna elements of virtual ULA is , the pth virtual steering vector is , and the virtual INDC is given bywhere , , , , , and .

For the (k + 1)th snapshot, the updated INDC is obtained aswhere is the updated INDC of the kth snapshot, is the virtual INDC of the (k+1)th snapshot, is the update constant, and for the first snapshot, .

The maximum likelihood (ML) tracking equation can be represented bywhere denotes the projection matrix of , , and denotes the trace of the matrix.

3. DOA Tracking Method Based on the QMVA

3.1. Quantum-Inspired Multiverse Algorithm

The quantum-inspired multiverse algorithm (QMVA) is inspired by quantum computation [26] and theory of cosmology [27]. Assume that denotes the number of quantum universes and G denotes the maximum number of iterations. The quantum state of the quantum universe at the th iteration is , where with B denoting the number of variables of the optimization problems, and the actual state of the qth quantum universe can be obtained bywhere , and denote the bth dimensional upper bound and lower bound, respectively, and denotes the fitness of .

In the exploration stage, quantum universes are sorted according to their fitness at each iteration, and on this basis, a quantum universe is chosen by the roulette wheel, and the specific formula is given bywhere denotes the normalized fitness of the qth quantum universe at the th iteration, where denotes the set of the fitness at the th iteration, and denotes the Euclidean norm of the vector, denotes a uniformly random number in [0, 1], and denotes the quantum state of the quantum universe at the th iteration selected by the roulette wheel. This roulette wheel mechanism can guarantee the exploration capacity of the QMVA.

In the exploitation stage, generate a wormhole existence probability at first, which is given bywhere and denote the user defined maximum and minimum of , respectively. Then, a travelling distance rate is given bywhere denotes the exploitation constant.

For the qth quantum universe at the th iteration, generate a uniformly random number in [0, 1]. If , the bth quantum rotational angle of the qth quantum universe at th iteration is given bywhere and are random numbers that obeys a uniform distribution in the range [0,1]. Then, the quantum state of the qth quantum universe is updated bywhere denotes the bth quantum state of the quantum universe with best fitness until the th iteration. The wormhole existence probability increases adaptively during entire iterations, which emphasizes the exploitation capacity of the QMVA. The travelling distance rate increases the accuracy of local search during entire iterations. Therefore, the balance between exploration and exploitation of the proposed QMVA is guaranteed by the effective combination of two stages.

3.2. DOA Tracking Based on the QMVA

To reduce the computational cost of DOA tracking, we propose the dynamic upper and lower bounds of the search space, which will continue to decrease as the number of snapshots increases, and the dynamic upper and lower bounds are defined as , where and denote the bth dimension upper and lower bounds of the search space for the kth snapshot, respectively.

At the first snapshot, the upper and lower bounds are the definition domain of the search space. For the kth snapshot, the upper and lower bound of the search space are updated bywhere denotes the convergence constant, denotes the centre value of the bth dimension search space for the (k − 1)th snapshot, , denotes the genetic factor, and denotes the bth dimension estimated value for the (k − 1)th snapshot.

Moreover, the maximum number of iterations is , where is a positive integer, and denotes the round down operation.

In the QMVA, the quantum states of the initial quantum universes are randomly generated in [0,1], and for the proposed DOA tracking method, the fitness function is defined aswhere is corresponding to the estimation of DOAs of the P moving targets; thus, B = P. Therefore, the proposed DOA tracking method can be summarized as follows:

Step 1. Obtain the first snapshot data, initialize

Step 2. Initialize the upper and lower bounds of the search space, and initialize the parameters of the QMVA: the number of quantum universes, the minimum and maximum of wormhole existence probability, the exploitation constant, and the maximum number of iterations .

Step 3. Initialize randomly the quantum states of the initial quantum universes, calculate the fitness of all quantum universes and store the quantum state of the quantum universe with best fitness

Step 4. In the exploration stage, quantum universes are sorted according to their fitness at each iteration, and on this basis, a quantum universe is chosen by the roulette wheel, for details, see equation (11). In the exploitation stage, quantum universes are evolved through equations (14) and (15)

Step 5. Calculate the fitness of all quantum universes, update the quantum state of the quantum universe with best fitness

Step 6. Examine whether is reached, if not, let , and then go back to step 4; otherwise, stop the loop iteration, output the actual state of the quantum state with best fitness for the kth snapshot and go to the next step

Step 7. Examine whether K_P is reached, if not, obtain the next snapshot data, update the INDC through (8), let , and go back to step 2; otherwise, output the DOA tracking results according to the actual state of the quantum universe with best fitness of all snapshots

4. Simulation Results

The root mean square error (RMSE) is defined aswhere denotes the ith true DOA at kth snapshot, denotes the ith estimated DOA at kth snapshot in the run, and N_e denotes the number of Monte-Carlo runs.

The probability of convergence (PROC) is defined aswhere

We consider that NA consists of two ULAs, , , difference constant , update constant , , , , and . For the QMVA, , , , and .

We execute 300 Monte-Carlo runs in the numerical simulations, and the proposed method is compared with three alternative methods: PAST [4], RPAST [16], and SSPAST [22].

4.1. The Independent Sources Tracking Scenario

Consider that two independent time-varying sources withwhere , and and are in degrees. Figures 1 and 2 show the RMSE curves and the PROC curves via α = 1.2, respectively. Figure 3 shows the PROC curves via GSNR = 10 dB. Figure 4 plots the DOA tracking results of four methods via GSNR = 10 dB and α = 1.2. From Figure 1, the proposed method has highest tracking accuracy in the low GSNR scenario compared with alternative algorithms. From Figures 2 and 3, the PROC of the proposed method exceeds 90% in terms of GSNR and characteristic exponent. From Figure 4, we can conclude that the tracking curve of the proposed method is smoother; in other words, the proposed method can achieve accurate tracking in the considered scenario compared with other methods.

(a)

(b)

(c)

(d)

4.2. The Coherent Sources Tracking Scenario

Consider that two coherent time-varying sources with the identical time-varying DOAs to the first simulation. Figures 5 and 6 show the RMSE curves and the PROC curves via α = 1.2, respectively. Figure 7 shows the PROC curves via GSNR = 10 dB. Figure 8 plots the DOA tracking results of four methods via GSNR = 10 dB and α = 1.2. From the figures, the PAST and the RPAST fail to work in the coherent sources tracking scenario. The SSPAST can achieve DOA tracking in the considered scenario, whose tracking performance is relatively weak. Furthermore, it is obvious that the proposed method yields robust tracking performance in the coherent sources tracking scenario.

(a)

(b)

(c)

(d)

4.3. The Underdetermined DOA Tracking Scenario

Consider that seven independent time-varying sources withwhere , and the corresponding , , and are in degrees. In the underdetermined DOA tracking scenario, the PAST and the RPAST will be invalid; thus, we plot the tracking results of the proposed method and the SSPAST method. Figures 9 and 10 plot the DOA tracking results of the proposed method and the SSPAST method via GSNR = 10 dB and α = 1.2. From the tracking results, although both the SSPAST and the proposed method can achieve underdetermined DOA tracking, it is obvious that our tracking method has better DOA tracking performance in the underdetermined DOA tracking scenario.

5. Conclusions

In this paper, we propose a robust DOA tracking method using nested array for achieving the underdetermined DOA tracking in the impulse noise. Simulation results demonstrate that our method offers better robustness and effectiveness both in independent and coherent sources tracking scenarios and obtains better DOA tracking results in the underdetermined DOA tracking scenario. In the future, we will try to generalize it to other complex DOA tracking problems.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

J. Shen, J. Yi, X. Wan, D. Xie, and F. Cheng, “Robust DOA estimation for passive radar with target signals mixed in the reference channel,” IEEE Geoscience and Remote Sensing Letters, vol. 18, no. 3, pp. 456–460, 2021.
View at: Publisher Site | Google Scholar
S. Li, Y. Liu, L. You, W. Wang, H. Duan, and X. Li, “Covariance matrix reconstruction for DOA estimation in hybrid Massive MIMO systems,” IEEE Wireless Communications Letters, vol. 9, no. 8, pp. 1196–1200, 2020.
View at: Publisher Site | Google Scholar
X. Han, M. Liu, S. Zhang, and J. Lan, “A passive DOA estimation algorithm of underwater multipath signals via spatial time-frequency distributions,” IEEE Transactions On Vehicular Technology, vol. 70, no. 4, pp. 3439–3455, 2021.
View at: Publisher Site | Google Scholar
B. Yang, “Projection approximation subspace tracking,” IEEE Transactions on Signal Processing, vol. 43, no. 1, pp. 95–107, 1995.
View at: Publisher Site | Google Scholar
V. Cevher, R. McClellan, and J. H. McClellan, “Acoustic multitarget tracking using direction-of-arrival batches,” IEEE Transactions on Signal Processing, vol. 55, no. 6, pp. 2810–2825, 2007.
View at: Publisher Site | Google Scholar
R. Badeau, B. Richard, and G. Richard, “Fast approximated power iteration subspace tracking,” IEEE Transactions on Signal Processing, vol. 53, no. 8, pp. 2931–2941, 2005.
View at: Publisher Site | Google Scholar
B. Lin, J. Liu, and M. ZhuJ. Zhu, “Direction-of-arrival tracking via low-rank plus sparse matrix decomposition,” IEEE Antennas and Wireless Propagation Letters, vol. 14, pp. 1302–1305, 2015.
View at: Publisher Site | Google Scholar
A. Masnadi-Shirazi and B. D. Rao, “A covariance-based Superpositional CPHD filter for multisource DOA tracking,” IEEE Transactions on Signal Processing, vol. 66, no. 2, pp. 309–323, 2018.
View at: Publisher Site | Google Scholar
Q. Liu, Y. So, and H. C. So, “DOA estimation in impulsive noise via low-rank matrix approximation and weakly convex optimization,” IEEE Transactions on Aerospace and Electronic Systems, vol. 55, no. 6, pp. 3603–3616, 2019.
View at: Publisher Site | Google Scholar
G. Jian, H. Wang, and Y. Guo, “An angle estimation method for monostatic MIMO radar based on RCC-FLOM algorithm,” Wireless Communications and Mobile Computing, vol. 2019, Article ID 9710656, 2019.
View at: Google Scholar
W. Peng, T. Qiu, F. Ren, and A. Song, “A robust DOA estimator based on the correntropy in alpha-stable noise environments,” Digital Signal Processing, vol. 60, pp. 242–251, 2017.
View at: Google Scholar
F. Ma, J. He, and X. Zhang, “Median-difference correntropy for DOA under the impulsive noise environment,” Wireless Communications and Mobile Computing, vol. 2019, Article ID 8107176, 2019.
View at: Publisher Site | Google Scholar
J. He, Z. Liu, and K. T. Wong, “Snapshot-instantaneous ||. ||_infin normalization against heavy-tail noise,” IEEE Transactions on Aerospace and Electronic Systems, vol. 44, no. 3, pp. 1221–1227, 2008.
View at: Google Scholar
Y. Du, H. Gao, and M. Chen, “Direction of arrival estimation method based on quantum electromagnetic field optimization in the impulse noise,” Journal of Systems Engineering and Electronics, vol. 32, no. 3, pp. 527–537, 2021.
View at: Google Scholar
H. Gao, M. Chen, Y. Du, and A. Jakobsson, “Monostatic MIMO radar direction finding in impulse noise,” Digital Signal Processing, vol. 117, Article ID 103198, 2021.
View at: Google Scholar
S.-C. Chan, Yu Wen, and Ka-L. Ho, “A robust past algorithm for subspace tracking in impulsive noise,” IEEE Transactions on Signal Processing, vol. 54, no. 1, pp. 105–116, 2006.
View at: Publisher Site | Google Scholar
X. Zhong, A. B. Madhukumar, and A. S. Madhukumar, “Particle filtering for acoustic source tracking in impulsive noise with alpha-stable process,” IEEE Sensors Journal, vol. 13, no. 2, pp. 589–600, 2013.
View at: Publisher Site | Google Scholar
J.-feng Zhang and T.-shuang Qiu, “A robust correntropy based subspace tracking algorithm in impulsive noise environments,” Digital Signal Processing, vol. 62, pp. 168–175, 2017.
View at: Publisher Site | Google Scholar
S. Ren, W. Dong, X. Li, W. Li, and X. Li, “Extended nested arrays for consecutive virtual aperture enhancement,” IEEE Signal Processing Letters, vol. 27, pp. 575–579, 2020.
View at: Publisher Site | Google Scholar
X. Wang and X. Wang, “Hole identification and filling in $k$-Times extended Co-prime arrays for highly efficient DOA estimation,” IEEE Transactions on Signal Processing, vol. 67, no. 10, pp. 2693–2706, 2019.
View at: Publisher Site | Google Scholar
F. Dong, L. Li, and X. Li, “Particle filter algorithm for DOA tracking using co-prime array,” IEEE Communications Letters, vol. 24, no. 11, pp. 2493–2497, 2020.
View at: Publisher Site | Google Scholar
Z. Zheng, Y. Huang, and H. C. So, “Spatial smoothing PAST algorithm for DOA tracking using difference coarray,” IEEE Signal Processing Letters, vol. 26, no. 11, pp. 1623–1627, 2019.
View at: Publisher Site | Google Scholar
Q. Zhang, J. Li, and P. Li, “A DOA tracking method based on offset compensation using nested array,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 69, no. 3, pp. 1917–1921, 2022.
View at: Publisher Site | Google Scholar
Raja Masadeh and A. Basel, “Mahafzah and ahmad sharieh, “sea lion optimization algorithm”,” International Journal of Advanced Computer Science and Applications, vol. 10, no. 5, pp. 388–395, 2019.
View at: Publisher Site | Google Scholar
F. A. Hashim, K. Hussain, E. H. Houssein, and W. Al Atabany, “Archimedes optimization algorithm: a new metaheuristic algorithm for solving optimization problems,” Applied Intelligence, vol. 51, no. 3, pp. 1531–1551, 2021.
View at: Publisher Site | Google Scholar
Z. Zhang, H. Gao, J. Ma, W. Shihao, and S. Helin, “Blind source separation based on quantum slime mould algorithm in impulse noise,” Mathematical Problems in Engineering, vol. 2021, Article ID 1496156, 145 pages, 2021.
View at: Publisher Site | Google Scholar
S. Mirjalili, S. M. Mirjalili, and A. Hatamlou, “Multi-Verse Optimizer: a nature-inspired algorithm for global optimization,” Neural Computing & Applications, vol. 27, no. 2, pp. 495–513, 2016.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Yufeng Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

194

Downloads

459

Citations