稳健Capon波束形成

稳健Capon波束形成
M = 10;%阵元数
j = sqrt(-1);
theta0 = 10;
A = zeros(10,8);
theta=[-75,-60,-45,-30,-10,25,35,50];
delta = 0.2;%角度增量
a0 = exp(-1i*(0:M-1)‘pisin(theta0));%原导向矢量
a_d = exp(-1i*(0:M-1)‘pisin(theta0 + delta));%信号方向增量的导向矢量
P_1=[];
P_2=[];
P_3=[];
for i = 1:8
A(:,i) = exp(-1i*(0:M-1)‘pisin(theta(i));
end
R0 = 100a0a0’+eye(10,10);
R1=100a_da_d’+eye(10,10);
sum = 0;
for i = 1:8
sum = sum + 10000*A(:,i)A(:,i)’;
R = R1 + sum;%接收信号有角度误差的自相关
C=3eye(10,10);%C矩阵
P1=1/(a_d’*R^(-1)a_d);%SOI功率
P1_dB = 10log10(P1);%SOI 标准功率
cvx_begin
variable a(M,1) complex;
minimize quad_form(a,R^(-1));
subject to
(a-a_d)'C^(-1)(a-a_d)<=1;
cvx_end
P2 = 1/(a’R^(-1)a);%SCB方法功率
P2_dB = 10log10(P2);
P3 = norm(a,2)^2 * P2/ M;%RCB功率（文中给出的功率）
P3_dB = 10log10(P3);
P_1(i) = P1_dB;
P_2(i) = P2_dB;
P_3(i) = P3_dB;
end
x = [1:8];
scatter(x,P_1,‘r’);
hold on
scatter(x,P_2,‘b’);
hold on
scatter(x,P_3,‘g’);
hold on
plot(x,P_1,‘r’);
hold on
plot(x,P_2,‘b’);
hold on
plot(x,P_3,‘g–’);
axis([1 8 19 29]);
xlabel(‘干扰数’);
ylabel(‘SOI能量/dB’);
title(‘delta = 0.2,e = 0.01’)