cisumyang.m

来自「关于MUSIC算法的源代码。很好用的」· M 代码 · 共 39 行

clc
clear
N=128;  %采样点%
n=1:N;
xn=sqrt(20)*sin(2*pi*0.2*n)+sqrt(2)*sin(2*pi*0.213*n)+randn(size(n));
Rx=xcorr(xn,'unbiased');
for i=1:128;
    for k=i:128;
        R1(i,k)=Rx(k-i+128);
    end
end
R2=R1';
for i=1:128;
    R1(i,i)=R1(i,i)-R1(i,i);
end
R=R1+R2;
p=4;m=N;
[V,D]=eig(R);
V=fliplr(V);
S=[V(:,1:p)];
G=[V(:,p+1:m)];
I=eye(128);
for r=1:1000
    w=2*pi*0.001*r;
    aw=1;
   for z=1:m-1
      aw=[aw;exp(-j*w*z)];
   end
        P=aw'*(I-S*S')*aw;
        Pw(r)=1/P;
end
w=linspace(0,1,1000);
plot(w,abs(Pw));
axis tight;
xlabel('w');
ylabel('Pw');
title('MUSIC算法');
grid on