music_1.txt

来自「这里有三个music算法用于来波方向估计的算法」· 文本 代码 · 共 29 行

clear;
d=1;
M=8;
N=2;
lma=2;
sam=200;
q=[20,40];%信源方向
A=exp(-2*pi*j*d*[0:M-1]'*sin(q*pi/180)/lma);
n=1:sam;
v1=1000000000000;
v2=3000000000005;
SNR=10;
Am=10^(SNR/20);
signal=Am*[cos(v1*n);cos(v2*n)];
noise=randn(M,sam)+j*randn(M,sam);
u=A*signal+noise;
c=u*u'/sam;
[s,h]=eig(c);
Vn=s(:,[1:M-N]);
qlb=[-pi/2:pi/180:pi/2];
for n=1:length(qlb)
    qla(n)=qlb(n);
    Ala=exp(-2*pi*j*d*[0:M-1]'*sin(qla(n))/lma);
    Pmusic(n)=(Ala)'*Ala*(pinv((Ala)'*Vn*(Vn)'*Ala));
    T(n)=qla(n);
end
P1=10*log10(abs(Pmusic));
T1=T*180/pi;
plot(T1,P1);grid