wide_band.m

基于CSM自适应波束形成——酉聚焦 ％ 8阵列； ％信噪比为10dB 干扰噪声比为20dB 干扰信号与需要信号不相关. ％干扰和信号频率所在范围为80—120Hz

% 基于CSM自适应波束形成——酉聚焦
% 8阵列；
%信噪比为10dB;干扰噪声比为20dB;干扰信号与需要信号不相关.
%干扰和信号频率所在范围为80—120Hz，相对带宽为40％ ，入射角度分别为seta_S_reg、seta_J_reg;
close all;
clear all;
M=8;
fl=80;
fh=120;
seta_S_reg=10*pi/180;%信号DOA为10度方向
seta_J_reg=20*pi/180;%干扰DOA
seta_reg=[seta_S_reg seta_J_reg];
Power_s=10;
Power_j=100;
d_l=1/2;
N=eye(M,M);
fc=100;%
fr=fc;
num_BW=40;
%生成聚焦频点的阵列流型
A_f_fc=exp(j*pi*(0:M-1).'*fr/fc*sin(seta_reg));
T_f_sub=zeros(M,M,num_BW);%（阵元数，阵元数，频点数）
R=zeros(M,M);
%生成每个频点的阵列流型
for i=1:num_BW
   f_sub(i)=fl+(fh-fl)*i/40;
   a_s_angle=pi*f_sub(i)/fc*sin(seta_S_reg);
   a_j_angle=pi*f_sub(i)/fc*sin(seta_J_reg);
   a_s=exp(j*a_s_angle*(0:M-1)).';
   a_j=exp(j*a_j_angle*(0:M-1)).';
   Rx_sub=Power_s*a_s*a_s'+Power_j*a_j*a_j'+N;
   A_f_sub=exp(j*pi*(0:M-1).'*f_sub(i)/fc*sin(seta_reg));
   E=A_f_sub*A_f_fc';
   [U,segma,V]=svd(E);
   %[segma,index_segma]=sort(diag(segma),'descend');
   T_f_sub(:,:,i)=V*U';
   R_sub=T_f_sub(:,:,i)*Rx_sub*T_f_sub(:,:,i)';
   R=R+R_sub;
end
a_s_angle_fc=pi*fr/fc*sin(seta_S_reg);
a_s_fc=exp(j*a_s_angle_fc*(0:M-1)).';
%a_s_fc=exp(j*pi*sin(seta_S_reg)*(0:M-1)).';
R=R/num_BW;
W_opt=inv(R)*a_s_fc./(a_s_fc'*inv(R)*a_s_fc);
D_angle=linspace(-pi/2,pi/2,2000);
angle=pi*sin(D_angle);
f=abs(W_opt'*(exp(j).^((0:M-1).'*angle)));
F=20*log10(f/max(f));
figure
plot(D_angle/pi*180,F);
grid