%用零階bessel函數近似描述圓陣的指向性（必須滿足條件：N>=4*pi*r/lamda+2） 
clear all;
close all;
N=60;
alpha=2*pi/N;
c=1500;
f=1500;
lamda=c/f;
r=lamda/0.4;
x=0:pi/180:pi; 
y1=4*pi*r/lamda*sin(x/2);
z1=bessel(0,y1);
subplot(211)
plot(x,abs(z1));title('圓陣的指向性：')
y2=4*pi*r/lamda*sin((x-pi/6)/2);
z2=bessel(0,y2);
subplot(212)
plot(x,abs(z2));