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