MUSIC算法的MATLAB实现2017年11月7日

clc
clear all
clf
M=64;  %天线阵列数
N=64;  %拍次       
x=zeros（MN）; 
%由于实验题目只给了128点，128/2=64故x最大规格为x（6464） 
for n=1:N
    for m=1:M
        nn=m+n-1;
        x（mn）=sqrt（20）*sin（2*pi*0.2*nn）+sqrt（2）*sin（2*pi*0.213*nn）+randn（11）;
    end
end

Rxx=x（:1）*x（:1）‘;
for n=2:N
    Rxx=Rxx+x（:n）*x（:n）‘;
end
Rxx=Rxx/N;
[vd]=eig（Rxx）;%求特征值，特征向量
p=0;
for n=1:M
    if d（nn）/d（MM）>0.05