11.m

来自「谱估计 (1)修正协方差法 (2)多重信号分类（MUSIC）算法 (3)E」· M 代码 · 共 20 行

clear
m=sqrt(-1);
delta=0.101043;
a1=-0.850848;
sample=32;      %number of sample spot
p=10;           %number of sample spot in coef method;
f1=0.05; f2=0.40; f3=0.42;
fstep=0.01;
fstart=-0.5;
fend=0.5;  
f=fstart:fstep:fend;
nfft=(fend-fstart)/fstep+1;
%calculate xn
    for n=1:sample
        xn(n)=2*cos(2*pi*f1*(n-1))+2*cos(2*pi*f2*(n-1))+2*cos(2*pi*f3*(n-1));
    end
 z=pmusic(xn,p);
figure
semilogy(f,z);
title('人工计算协方差');