面阵的波束形成

clear;
close all;
clc;

%均匀面阵方向图
%阵元间距波长一半
%16*16个阵元

freq =  1000;    
c = 1500;        
lambda = c/freq; 
d = lambda/2;    

N = 16;
M = 16;          %阵元个数



theta1 = linspace（0180181）; 
phi1 = linspace（09091）;


[XY] = meshgrid（phi1theta1）;

beta_x = 2*pi*d*cosd（Y）.*cosd（X）/lambda;
beta_y = 2*pi*d*sind（Y）.*cosd（X）/lambda;
beta_z = 2*pi*d*sind（X）/lambda;

G_row = 0;
for j = 1:1:N     
    G_row = G_row + exp（-i*（j-1）*beta_x）; 
end

G_col = 0;
for k = 1:1:M     
    G_col = G_col + exp