Mathematical Problems in Engineering
Volume 2017 (2017), Article ID 7904356, 7 pages
https://doi.org/10.1155/2017/7904356

Research Article

MVDR Algorithm Based on Estimated Diagonal Loading for Beamforming

Yuteng Xiao,¹ Jihang Yin,¹ Honggang Qi,² Hongsheng Yin,¹ and Gang Hua¹

¹China University of Mining and Technology, Xuzhou 221116, China
²University of Chinese Academy of Sciences, Beijing 101408, China

Correspondence should be addressed to Hongsheng Yin; xuzhouyhs@sina.com

Received 1 June 2017; Revised 24 August 2017; Accepted 11 September 2017; Published 12 October 2017

Academic Editor: Thomas Schuster

Copyright © 2017 Yuteng Xiao 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.

Abstract.

Beamforming algorithm is widely used in many signal processing fields. At present, the typical beamforming algorithm is MVDR (Minimum Variance Distortionless Response). However, the performance of MVDR algorithm relies on the accurate covariance matrix. The MVDR algorithm declines dramatically with the inaccurate covariance matrix. To solve the problem, studying the beamforming array signal model and beamforming MVDR algorithm, we improve MVDR algorithm based on estimated diagonal loading for beamforming. MVDR optimization model based on diagonal loading compensation is established and the interval of the diagonal loading compensation value is deduced on the basis of the matrix theory. The optimal diagonal loading value in the interval is also determined through the experimental method. The experimental results show that the algorithm compared with existing algorithms is practical and effective.

1. Introduction

The beamformer is used for the purpose of detecting the desired signals and suppressing the interference signals. It filters the interference signals and outputs the desired signals whose mainlobe points to desired direction. Beamformer, which is the key technology in many fields, is widely used in radar, communications, sonar, GPS, and other fields. For example, weather radar system needs the algorithm to get the desired signals [1] from the meteorological information. In MIMO communication system, the beamformer is used to achieve large array gains [2].

MVDR (Minimum Variance Distortionless Response) is a typical beamforming (beamformer) algorithm, which was proposed by Capon in 1967 [3]. It makes the outputting power with minimum interference and noise in the desired direction through adjusting a weight factor. However, the performance of the beamforming algorithm will decline dramatically with the inaccuracy in the steering vector and the covariance matrix [4]. To reduce the dispersive degree of the eigenvalues of the covariance matrix, Carlson proposed diagonal loading by pulsing a constant diagonal loading compensation value to the diagonal elements of the covariance matrix [5]. On this basis of constant diagonal loading compensation, the literature [6] selects the loading value which is higher than the background noise. However, when the snapshots are small, there is a large inaccuracy in the sample covariance matrix, and the constant diagonal loading compensation value cannot correct the inaccuracy. Literatures [7–9] reconstruct the covariance matrix which is also one of the diagonal loadings. In literature [10], an optimal algorithm employing variable diagonal loading of the sample covariance matrix eigenvalues is proposed. In literature [11], an algorithm is used to compute the diagonal loading level automatically from the given data with a shrinkage method of enhancing covariance matrix. Literature [12] uses beam-to-reference ratio which is estimated as a weighting factor for variable diagonal loading. These algorithms have solved the problem that the diagonal loading value is not easy to determine and cannot be adjusted automatically with the sample changing. However, these algorithms are more complex and less efficient in the calculation.

In this paper, we will adopt the diagonal loading method to improve MVDR algorithm. It uses the adaptive diagonal loading compensation method, adjusting the diagonal elements in sample covariance matrix. The sharp decline, which is caused by the inaccuracy in the covariance matrix, is improved in performance.

The rest of this paper is organized as follows. In Section 2, the beamforming array signal model is introduced. In Section 3, the model with beamforming MVDR algorithm is proposed and the model is solved by Lagrange multiplier. Then, the model with MVDR optimization based on diagonal loading in theory is discussed in Section 4. In Section 5, the method for estimating the diagonal loading value is discussed and analyzed. The experimental results are shown in Section 6. Finally, we conclude this paper in Section 7.

2. Beamforming Array Signal Model

Figure 1 presents uniform linear array signal model. There are narrow band multisignal sources (containing both desired signals and interference signals) radiating at the same time, whose wavelengths are all . The received antenna is an linear array with distance and . Taking the weather radar system as an example, the desired signals are usually the meteorological information with different desired directions and the interference signals are the various clutter information reflected by the ground. According to the Shannon communication system model, the noise is concentrated on the beamformer. In practice, the desired signals are much smaller than the interference signals and the noise, whose Signal-Noise Ratio and Signal/Interference-Noise Ratio are often up to −20 dB. Thus, the beamformer is used to receive narrow band multisignal source. It makes the desired signals without attenuation through adjusting weight factor . And at the same time, the interference signals and noise are greatly suppressed.

Figure 1: Uniform linear array signal model.

Referring to Figure 1, single signal is first analyzed. is radiated to the array parallel, whose desired direction is . The array filters snapshot data . The snapshot data received have a time delay due to the distance of the array, which can be expressed as

where is speed of light. The signal filtered by the th array is expressed as

After the single signal and the noise are got through the beamformer together, the output is

where , is column vector with noise, and , .

When the multisignal source radiates at the same time, the output matrix of array is

where ; ; it is dimension matrix; could be further divided into two parts: output matrix of array with desired signals and output matrix of array with interference signals .

After the multisignal source and the noise are got through the beamformer together, the output is

3. Beamforming MVDR Algorithm

According to the literature [3], MVDR makes the outputting power with minimum interference and noise in the desired direction through adjusting a weight factor, while ensuring the output desired signal without minimum distortion.

From (5), the output power of interference signal adding noise can be expressed as

For convenience, (6) , , where is the interference signal adding noise. Mathematically, is equal to the covariance matrix of .

Similarly, from (5), in order to guarantee the nondistortion output of the desired signals, we can deduce the necessary and sufficient condition:

where is directional vector with normalizing factor .

Therefore, beamforming MVDR algorithm can be expressed as the following minimization problem:

Using Lagrange multiplier, (8) can be changed into the unconstraint equation

The weighted vector is obtained as

4. MVDR Algorithm Based on Diagonal Loading Compensation

According to literature [5], beamforming MVDR algorithm requires an accurate covariance matrix . In practice, the performance of beamforming MVDR algorithm declines dramatically due to the inaccuracy with the covariance matrix based on sampling and the theoretical covariance matrix . Therefore, literature [5] corrects the inaccuracy of by adding a constant value to the diagonal elements of . As a result, the performance is improved and the robustness is enhanced.

In practice, can be estimated using snapshots data with finite order. The method is as below:

where is the inaccuracy term of covariance matrix of sampling data. It is a random matrix with zero mean and unit variance, all the elements of which are far less than 1. is the sampling data of interference signals adding noise in a moment. is the sampling snapshots.

from (8) can be expressed as covariance matrix with diagonal loading:

where is covariance matrix with diagonal loading, is the unit matrix, and is the factor of diagonal loading, which is used to control the diagonal loading value.

Therefore, to get the MVDR optimization problem based on diagonal loading, we amend (8) as

Similarly, using Lagrange multiplier in (13), the weighted vector is obtained as

5. The Estimated Method with Diagonal Loading Compensation

In practice, from (11) is determined manually [6]. It needs strong experience. The covariance matrix on sampling has irreconcilable inaccuracy especially in the small sampling snapshots. The constant value of diagonal loading regulated by experience often fails to work. It also cannot correct the inaccuracy. This paper presents an adaptive diagonal loading method which can change the diagonal loading value with the difference of input signals. It decreases the interval of diagonal loading value depending on the relation between the inaccuracy of covariance matrix on sampling and the diagonal loading value. Therefore, this paper provides theoretical guidance for researchers to select diagonal loading value effectively.

From (12), due to , the inverse of covariance matrix on diagonal loading can be expressed as

The is unfolded as

Because all the elements of from (16) are far less than the elements of matrix , . Then

Equation (15) can be approximately expressed as

from (18) can be expressed as

where is signal power corresponding to the th interference signal, is the noise power, is the number of interference signals, is the directional matrix with interference signal, and is diagonal matrix corresponding to the interference signal power.

Using the matrix inversion lemma [13] and substituting (19) into (18), there is

where

The value of first bracket from (20) is close to the covariance matrix ; the diagonal loading value should be far less than the diagonal elements of covariance matrix ; that is,

It is the reason that performance of MVDR beamforming is degraded by the in the curly bracket from (19). And all the elements of are far less than 1 and close to 0. To get the optimal performance of MVDR, it should meet the condition . And it can be explained as

From (22) and (23), there is

where is the average of diagonal elements in the covariance matrix of sampling data, . Here is the number of arrays; is the trace of matrix. is the inaccuracy of covariance matrix of sampling data, which can be expressed by obtaining the standard deviation of , . Here is the diagonal elements of matrix; is the standard deviation of matrix. Moreover, (24) shows that , will change with the difference of input signals, and the diagonal loading value will change relatively. Thus, the algorithm is adaptive.

6. Discussion and Simulation

6.1. Evaluation Methodology of Experimental Result

6.1.1. The Output Signal/Interference-Noise Ratio (SINR)

The ratio with output power of the desired signals to the interference adding noise signals is defined as the output SINR. When the unit of SINR is dB, it can be rewritten as

where is the output power of desired signals.

6.1.2. Beamforming Array Pattern

According to Figure 1, the relation between the absolute value of the output array and the direction angle is called beampattern, which is used to evaluate the performance of beamformer.

From (3), the ratio of output to input arrays is defined as ; that is,

From (26), to get the gain of beampattern in unit dB, it can be rewritten as

where is the maximum beampattern. Ideally, when the direction corresponding to its mainlobe is from desired signal, gets the max; the rest of directions are sidelobes. The smaller the sidelobes are, the stronger the suppression for interference adding noise is.

6.2. Simulation and Analysis

It can be observed from (24) that the diagonal loading value has been limited to an interval. In this paper, with the value increasing gradually in the interval, we try to observe the performance of MVDR on diagonal loading beamformer with the output SINR as the index and get the better diagonal loading value. Then we use (27) to evaluate MVDR algorithms on both constant diagonal loading and estimated diagonal loading , where is optimal value from the first experiment. In this paper, we design two experiments.

Experiment 1 (optimization and estimation of diagonal loading ). In this experiment, we employ a uniform 16-element linear array. The distance of elements is half of the wavelength, the direction of the desired signal is 0 degrees, and the direction of the interference signal is 40 degrees. We select three representative input scenarios: the first scenario is that the power of desired signal is far less than the noise and the power of interference signal is greater than the noise. At this point, the input SNR is −5 dB and the input INR is 10 dB. Therefore, input SINR is −15.414 dB. The second scenario is that the power of desired signal is far less than the noise and the power of interference signal is equal to the noise. At this point the input SNR is −5 dB and the input INR is 0 dB. Therefore, input SINR is −8.01 dB. The third scenario is that the power of desired signal is greater than the noise and the power of interference is equal to the noise. The input SNR is 5 dB and the input INR is 0 dB. Therefore, input SINR is 1.990 dB. From (24), 10 points are selected in the interval . And we further optimize the diagonal loading by using the output SINR. Considering the relations between diagonal loading value and the output SINR in the three scenarios, the snapshots are set as 5, 25, and 50, and the diagonal loading value is set as the unit on . Figures 2, 3, and 4 show that no matter whether snapshots are 5, 25, or 50 in various scenarios, when , the output SINR gets the maximum value, and the performance gets the best. The results were in line with general cognition: diagonal loading factor is a compensation factor. When the compensation factor is equal to the error of the covariance matrix, the disturbance problem on the performance can be counteracted preferably, and the factor is independent of the input SINR.

Figure 2: The relations between diagonal loading value and output SINR in different snapshots (input: SNR = −5 dB, INR = 10 dB, and SINR = −15.414 dB).

Figure 3: The relations between diagonal loading value and output SINR in different snapshots (input: SNR = −5 dB, INR = 0 dB, and SINR = −8.01 dB).

Figure 4: The relations between diagonal loading value and output SINR in different snapshots (input: SNR = 5 dB, INR = 0 dB, and SINR = 1.990 dB).

Experiment 2 (comparison of algorithm performance). The experiment selects the input parameters in the most unfavorable situation (the first scenario) in Experiment 1. The input SNR is −5 dB, INR is 10 dB, and input SINR = −15.414 dB. Under the results of Experiment 1, , analyzing the MVDR algorithm, MVDR on constant diagonal loading [6] and the MVDR on estimated diagonal loading with the output SINR and the beamforming array pattern in different snapshots (5, 15, 25, and 50). The results are shown in Table 1 and Figures 5 and 6–9.
Table 1 and Figure 5 show that, under the same input conditions, with the increase of the snapshots, the suppression on the interference adding noise of three methods is getting better. When the snapshots are small, the traditional MVDR algorithm has less suppression on the interference adding noise. And when the snapshots are big, the suppression is tolerable. But it shows poor stability. The algorithm with the paper () and the algorithm with literature [6] () are obviously better than MVDR algorithm in the stability and the suppression on the interference adding noise. However, the algorithm with this paper is superior to the algorithm with literature [6] in the stability and the suppression on the interference adding noise.
Figures 6 and 7 show that when the snapshots are big (), the amplitudes of mainlobe in the algorithm with literature [6] and the algorithm with this paper are equal; the amplitude of sidelobe in the algorithm with this paper is below −15 dB, which is slightly superior to the algorithm with literature [6]; the above two methods are both superior to MVDR algorithm.
Figures 8 and 9 show that when snapshots are small (), the algorithm with this paper has obvious advantages of the same direction in pattern. It mainly shows that the amplitudes of sidelobe are superior to the algorithm with literature [6], which means that the suppression on the interference and noise is the best.

Table 1: Output SINR comparison of three algorithms in different snapshots (input: SNR = −5 dB, INR = 10 dB, and SINR = −15.414 dB).

Snapshots
Output SINR/dB

MVDR

5
−10.5200
1.1889
1.6289

15
−8.6626
2.3498
2.7657

25
−0.4494
2.8155
3.1409

50
3.0710
3.5075
3.6187

Figure 5: Output SINR comparison of three algorithms in different snapshots (input: SNR = −5 dB, INR = 10 dB, and SINR = −15.414 dB).

Figure 6: Comparison of three algorithms for beamforming array pattern ().

Figure 7: Comparison of three algorithms for beamforming array pattern ().

Figure 8: Comparison of three algorithms for beamforming array pattern ().

Figure 9: Comparison of three algorithms for beamforming array pattern .

7. Conclusion

To enhance the robustness of the existing beamforming MVDR algorithms, we propose the MVDR algorithm based on estimated diagonal loading for beamforming. It is on the basis of learning the uniform linear array signal model and beamforming MVDR algorithm. The paper establishes the MVDR optimization model on diagonal loading compensation solved by Lagrange multiplier. The interval of the diagonal loading compensation is deduced on the basis of the matrix theory. And the optimal diagonal loading is also determined through the experimental method. The experimental results show that when the snapshots are small, the traditional MVDR algorithm has less suppression on the interference adding noise. And when the snapshots are big, the suppression is tolerable. It shows poor stability. The algorithm with the paper and the algorithm with literature [6] are obviously better than MVDR algorithm in the stability and the suppression on the interference adding noise. However, the algorithm with this paper is superior to the algorithm with literature [6] in the stability and the suppression on the interference adding noise.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is partly supported by the National Natural Science Foundation of China (61379100, 61472388, and 51574232).

References

C. Fulton, J. L. Salazar, Y. Zhang et al., “Cylindrical Polarimetric Phased Array Radar: Beamforming and Calibration for Weather Applications,” IEEE Transactions on Geoscience and Remote Sensing, vol. 55, no. 5, pp. 2827–2841, 2017.

N. Song, T. Yang, and H. Sun, “Overlapped subarray base hybrid beamfroming for millimeter wave multiuser massive MIMO,” IEEE Signal processing Letters, vol. 24, no. 5, pp. 550–554, 2017.

J. Capon, “High-resolution frequency-wavenumber spectrum analysis,” Proceedings of the IEEE, vol. 57, no. 8, pp. 1408–1418, 1969.

F. Shen, F. Chen, and J. Song, “Robust adaptive beamforming based on steering vector estimation and covariance matrix reconstruction,” IEEE Communications Letters, vol. 19, no. 9, pp. 1636–1639, 2015.

B. D. Carlson, “Covariance matrix estimation errors and diagonal loading in adaptive arrays,” IEEE Transactions on Aerospace and Electronic Systems, vol. 24, no. 4, pp. 397–401, 1988.

H. Cox, R. Zeskind, and M. Owen, “Robust adaptive beamforming,” IEEE Transactions on Acoustics, Speech and Signal Processing, vol. 35, no. 10, pp. 1365–1376, 1987.

Y. Gu and A. Leshem, “Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation,” IEEE Transactions on Signal Processing, vol. 60, no. 7, pp. 3881–3885, 2012.

X. Yang, L. Yan, Y. Sun, and T. Zeng, “Improved orthogonal projection approach utilising interference covariance matrix reconstruction for adaptive beamforming,” IEEE Electronics Letters, vol. 50, no. 20, pp. 1446-1447, 2014.

Z. Lu, Y. Li, M. Gao, and Y. Zhang, “Interference covariance matrix reconstruction via steering vectors estimation for robust adaptive beamforming,” IEEE Electronics Letters, vol. 49, no. 22, pp. 1373-1374, 2013.

J. Gu and P. J. Wolfe, “Robust adaptive beamforming using variable loading,” in Proceedings of the 4th IEEE Sensor Array and Multichannel Signal Processing Workshop Proceedings, SAM 2006, pp. 1–5, Waltham, MA, USA, July 2006.

L. Du, J. Li, and P. Stoica, “Fully automatic computation of diagonal loading levels for robust adaptive beamforming,” IEEE Transactions on Aerospace and Electronic Systems, vol. 46, no. 1, pp. 449–458, 2010.

C. Wu, Y. Guo, Y. Na, X. Wang, Q. Fu, and Y. Yan, “Robust beamforming using beam-to-reference weighting diagonal loading and Bayesian framework,” IEEE Electronics Letters, vol. 51, no. 22, pp. 1772–1774, 2015.

D. J. Tylavsky and G. R. L. Sohie, “Generalization of the Matrix Inversion Lemma,” Proceedings of the IEEE, vol. 74, no. 7, pp. 1050–1052, 1986.