MUSIC算法估计功率谱

clear all
N=256; 
n=0:1:N-1;                
f1=0.2;                              
f2=0.213;
M=10;
snr=3;
 
x1=5*exp（1i*2*pi*f1*n+pi/3）+10*exp（1i*2*pi*f2*n+pi/4）;  
xn=awgn（x1snr）;
 
    M=8;                                 %滤波器阶数
    for k=1:N-M                          
    xs（:k）=xn（k:1:k+M-1）.‘;
    end
    R=xs*xs‘/（N-M）;                      %估计样本相关矩阵
   
 
    [UE]=svd（R）;