dongmain.m

MUSIC 算法的原程序,希望对大家有帮助

clear;
N=10000;
M=20;
P=5;
Am=[ones(1,2)];
Sp=sum(Am.*Am);
Sn=Sp/(50*2);
m=(1:M)';
c=1500;
freq=1000;
% d=0.45;
theta=pi/2;
phi=[90,90.5];
% t=(1:N)/200000;
noise=randn(M,N)+j*randn(M,N);
X=signal(Am,N);
% S=dm(d,phi, m);
Tor=timedelay(phi, M);
S=exp(-j*2*pi*freq*Tor);
% S(1,:)=exp(j*35/180)*S(1,:);
% S(2,:)=exp(j*35/180)*S(2,:);
% S(3,:)=exp(j*35/180)*S(3,:);

% S(1,:)=1/2*S(1,:);
% S(2,:)=1/2*S(2,:);
% S(3,:)=1/2*S(3,:);

Y=S*X+noise*(sqrt(Sn));
R=Y*Y'/N;
[U,V] = eig(R);
[lamda,index] = sort(diag(V));
G=U(:,index(1:(M-P)));
phi2=79:0.001:82;
NN=length(phi2);
for k=1:NN
    Tor2=timedelay(phi2(k),M);
    S2=exp(-j*2*pi*freq*Tor2);
% AA=dm(d,phi2(k),m);
pw(k)=1/abs((S2'*G*G'*S2));
end
F=10*log10(pw);
plot(phi2,F);
xlabel('DOA  (Degree)');
ylabel('F (dB)');
% axis([-2 2 max(F)-10 max(F)]);
grid;