Direction-of-arrival (DOA) estimation performance of adaptive filtering-based methods, <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">e.g.</i> , least mean square (LMS), degrades significantly in adverse conditions ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">i.e.,</i> at low signal-to-noise ratio (SNR), small number of antenna array elements and snapshots, and in presence of two closely spaced signal sources). Moreover, these estimation methods require many regulating parameters to efficiently update the step size, which are very difficult to tune manually in practical scenarios. Although the subspace decomposition-based methods, <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">e.g</i> ., multiple signal classification (MUSIC), provide somewhat better performance in adverse conditions, they are based on the statistical properties of the signal such as spatial covariance matrix (SCM) and its eigenvalue decomposition (EVD), resulting in high computational complexity. In this paper, a robust DOA estimation method is proposed based on a variable step size LMS algorithm with low rank matrix approximation (LRMA), where the received signal de-noising problem is formulated first as a LRMA problem directly using the received signal observations instead of the SCM, and then the variable step size is calculated from the predicted instantaneous error and the estimated powers of reference and <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">de-noised</i> signals. The output spatial spectrum is determined by the reciprocal of the antenna array pattern, and the high peaks in the output spatial spectrum give the received signals estimated directions. The proposed method updates the step size efficiently without choosing any regulating parameter, and achieves improved performance even in adverse conditions due to the de-noising feature of LRMA with low computational complexity. Further, we also conduct the convergence and stability analyses of the proposed iterative estimation method with the <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">de-noised</i> signal. Numerical results demonstrate that the proposed method outperforms state-of-the-art methods especially those, which yielded based on adaptive filtering.