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）;??%產生干擾信號