?? capon.m
字號:
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);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -