子阵均匀波束形成

clc;
clear;
clc;
%这个程序给的角度和最后形成的方向图主瓣零点角度刚好相反
%%  一、阵元基本参数
NN = 1000; %观测数据样本个数
N = 40;   %阵元总个数
L = 4;    %子阵内阵元个数
dn = 1/2; %阵元间距为半波长
M = N/L ; %子阵个数e
noise = （randn（NNN）+1i*randn（NNN））/sqrt（2）; %噪声功率为1
theta_sig = 10;% 期望信号存在时期望信号入射角度
theta_jam =50;%干扰信号入射角度
SNR_sig = 0; %期望信号存在时信噪比
SNR_jam = 20; %干扰信号信噪比
sigma= 10;    %波束形成角度

%% 二、产生期望信号和干扰信号的导向矢量，构造降维矩阵T
A_sig_ele = exp（-1i*2*pi*dn*sin（sigma*pi/180）*（0:N-1）‘）; %期望信号的导向矢量
A_jam_ele = exp（-1i*2*pi*dn*sin（theta_jam*pi/180）*（0:N-1）‘）; %干扰信号的导向矢量
C = exp（1i*（0:L-1）*2*pi*dn*sin（sigma*pi/180））;
B = eye（M）;
T0 = kron（BC）;
T=T0‘;
A_sig_sub = T0*A_sig_ele;
A_jam_sub = T0*A_jam_ele;


%% 三、产生期望信号和干扰信号
As_sig = 10^（SNR_sig/10）;  %期望信号的幅度
%As_sig= 0 ;           %无期望信号时
As_jam = 10^（SNR_jam/10）;   %干扰信号的幅度
sig = As_sig*exp（1i*0.15*2*pi*（0:NN-1）+1i*2*pi*rand）; %产生期望信号
jam = As_jam*exp（1i*0.1*2*pi*（0:NN-1）+1i*2*pi*rand）;  %产生干扰信号