很多有关原子范数最小化的DOA估计的文章中提到“再对Toeplitz矩阵T(u)进行Vandermonde分解，即可以得到角度θ和功率P的估计结果”。我想请教一下实现“Vandermonde分解的编程思路。