clear?all;
close?all;
clc;
f=1000;???????????%陣元發(fā)射信號(hào)頻率
c=1500;???????????%聲速
lamta=c/f;????????????%波長(zhǎng)
d=lamta/2;????????????%陣元間距為半波長(zhǎng)
%%%%%%%%%%%%%%%%%均勻線(xiàn)陣%%%%%%%%%%%%%%%
theta=-pi:pi/360:pi;
phai=-pi:pi/360:pi;?%方位角范圍
thetad=30*pi/180;??????????????%陣列指向
betad=2*pi*d*sin（thetad）/lamta;
m1=8;??????????????????????????%陣元個(gè)數(shù)
beta?=?2*pi*d*sin（theta）/lamta;
G11=sin（m1*（beta-betad）/2）;
G12=m1*sin（（beta-betad）/2）;
G1=abs（G11./G12）;
GL=20*log（G1/max（G1））;
figure
plot（theta*180/piGL）;grid;
axis（[-180?180?-100?0]）;
title（‘陣元M=8thetad=30均勻線(xiàn)陣方向圖‘）;
xlabel（‘方位角/度‘）;
ylabel（‘G（\theta）/dB‘）;
figure
polar（thetaG1）;title（‘陣元M=8thetad=30均勻線(xiàn)陣方向圖‘）;
xx=inline（‘a(chǎn)bs（（sin（8*（pi*sin（theta））/2））./（8*sin（（pi*sin（theta））/2）））*sin（theta）*cos（phai）‘）;
yy=inline（‘a(chǎn)bs（（sin（8*（pi*sin（theta））/2））./（8*sin（（pi*sin（theta））/2）））*sin（theta）*sin（phai）‘）;
yz=inline（‘a(chǎn)bs（（sin（8*（pi*sin（theta））/2））./（8*sin（（pi*sin（theta））/2）））*cos（theta）‘）;
figure
ezmesh（xxyyyz[0?2*pi?0?pi]100）;
%%%%%%%%%%%%%%%%%平面陣%%%%%%%%%%%%%%%%%
x=16;?
y=16;???????????????????%16*16面陣
z=16;
theta1=0:pi/360:pi;????%方位角
fai1=0:pi/360:pi;??????%俯仰角
[fai1theta1]=meshgrid（fai1theta1）;
beta1=2*pi*d*cos（theta1）.*cos（fai1）/lamta;
beta2=2*pi*d*sin（theta1）.*cos（fai1）/lamta;
beta3=2*