capon.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);
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:100
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/100;
i=1
for u=-1:0.01:1
  V=exp(j*n*pi*u)
  P(i)=1/(V'*inv(Sx)*V)
  PP(i)=(V'*Sx*V)/(V'*V)
  i=i+1
end
u=-1:0.01:1;

subplot(2,1,1);
plot(u,P);
subplot(2,1,2);
plot(u,PP);