DBF波束形成

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（