music 算法

clear all;
N = 1000;
noise = （ randn（1  N） + 1i * randn（1  N）） / sqrt（2）;
%产生带噪声的信号样本u（n）
signal1 =  exp（ 1i * 0.5 * pi * （0 : N-1） + 1i * 2 * pi * rand）;
signal2 =  exp（ -1i * 0.3 * pi * （0 : N-1） + 1i * 2 * pi * rand）; 
un = signal1 + signal2 + noise;
%计算自相关矩阵
M = 8;
 for k=1:N-M
     xs（:k）=un（k+M-1:-1:k）.‘;
 end 
R=xs*xs‘/（N-M）;
%自相关矩阵的特征值分解
[UE] = svd（R）;
ev = diag（E）;

%根据AIC准则进行信号源个数的估计
for k = 1:M
    dec = prod（ev（k:M）.^（1/（M-k+1）））; %计算第一项中对数的自变量的分子
    nec = mean（ev（k:M））;%计算第一项中对数的自变量的分母
    lnv = （dec /nec）^（（M-k+1）*N）;%计算第一项中对数的自变量
    AIC（k） = -2 * log（lnv） + 2 * （k-1） * （2 * M - k + 1）;
end
[AminK] = min（AIC）;
N1 = K - 1;
%根据MDL准则进行信号源个数的估计
for k = 1:M