clc;
clear?all;
close?all;
imag=sqrt（-1）;
element_num=48;%陣元數為8
d_lamda=1/2;%陣元間距d與波長lamda的關系
d1_lamda=1/4;
theta=linspace（-pi/3pi/348）;
theta0=0;%來波方向


w1=exp（imag*2*pi*d_lamda*sin（theta0）*[0:element_num-1]）;
m1=exp（imag*2*pi*d1_lamda*cos（30/180*pi））;
m2=exp（imag*2*pi*d1_lamda*2*cos（30/180*pi））;
m3=exp（imag*2*pi*d1_lamda*3*cos（30/180*pi））;

w2=w1*m1;
w3=w1*m2;
w4=w1*m3;
w=[w1;w2;w3;w4];
%%%%%%%%DBF處理%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%對列處理%%%%%%%%%%%%%%%%
w=[0;m1;m2;m3]‘*w;

for??j=1:length（theta）
????a=exp（imag*2*pi*d_lamda*sin（theta（j））*[0:element_num-1]）;
???p（