music.m

array processing code for communication engineering

clear all;
close all;
A=1/10;
N =10;
n = (-(N-1)/2:(N-1)/2)';
Vm = ones(N,1);
u=1:100;
SNR=10.^(20/10);
ui=0.15;
us=0;
Vi = exp(j*n*pi*ui);
Vs=exp(j*n*pi*us);
Sx=zeros(10,10);
for i=1:2
X=sqrt(SNR)*exp(j*2*pi*randn)*(randn+j*randn)*Vi/sqrt(2)+exp(j*2*pi*randn)*(randn+j*randn)'*sqrt(SNR)*Vs/sqrt(2)+(randn+j*randn)/sqrt(2);

Sx=(X*X')+Sx;
end
Sx=Sx/2;
[U,S,V]=svd(Sx);
Us=zeros(10,2);
Us(:,1)=U(:,1);
Us(:,2)=U(:,2);
i=1;
for u=-1:0.001:1
  V=exp(j*n*pi*u);
  P(i)=1/(V'*(eye(10,10)-Us*Us')*V);
  i=i+1;
end
u=-1:0.001:1;
P= 20*log10(abs(P));
plot(u,P);
  M1=P(1100);
for i=1100:1:1200
    if P(i)>M1
    M1=P(i)
    j1=i
    end
end
M2=P(950);
for i=950:1:1050
    if P(i)>M2
    M2=P(i)
    j2=i
    end
end
U1=-1+(j1-1)*0.001;
U2=-1+(j2-1)*0.001;