DOA Estimation for Sources with Large Power Differences

Algorithm 1

Pseudocode of the proposed algorithm.

(i)	Step 1. Obtain the output of the uniform linear antenna array from equation (1) and compute the sample covariance matrix as in equation (3).
(ii)	Step 2. Perform the EVD of to obtain the noise subspace eigenvalue matrix .
(iii)	Step 3. Apply diagonal loading on to obtain through equation (14). is equal to , which is the largest eigenvalue of the signal subspace. The noise subspace eigenvalues of are .
(iv)	Step 4. Suppose that the virtual source impinges on the antenna array from direction ; then, obtain the matrix as in equation (18), where , denotes the trace of and M is the number of array elements.
(v)	Step 5. Perform the EVD of to obtain the eigenvalues in descending order as .
(vi)	Step 6. For each , compute the spatial zero spectrum through equation (23). The DOAs of the sources are P sets of that enable to reach P maximum values, as in equation (24).