线性阵列自适应波束形成，干扰抑制

%%%%  这个程序是正确的事实上一旦考虑自适应，就根本不用考虑波束指向了
%%%%  波束指向只针对静态方向图有用
clear all;clc;
close all;
tic
N=24;            %%  the element number
c=pi;            %%  c=2*pi*d/lamda where d=0.5*lamda
thetaj=40;                      %%%%   干扰信号方向
theta0=10;                      %%%%   波束方向
j=sqrt（-1）;
index=0;
step=0.01
for w=-90:step:90
    index=index+1;
    aw=exp（-j*（0:N-1）*c*（sin（w*pi/180）-sin（theta0*pi/180）））;
    PSQ（index）=20*log10（abs（ones（1N）*aw‘））;            %% PSQ静态和方向图  P:pattern; S:sum;        Q:queisent;
end
w=-90:step:90;
figure;plot（wPSQ-max（PSQ））;xlabel（‘Angle （ \circ）‘）; ylabel（‘Power （dB）‘）; title（‘Queisent sum pattern‘）
temp1=sort（PSQ（find（diff（sign（diff（PSQ）））==-2）+1））;     %%%  静态和方向图所有极大值由小到大排列
mainlobeQS=temp1（length（temp1））;                        %%%  主瓣峰值
sidelobeQS=temp1（length（temp1）-1）;                      %%%  旁瓣峰值
if mainlobeQS==sidelobeQS
    sidelobeQS=t