%用零階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),'r');
title('圓陣的指向性：無相移');xlabel('信號入射方向'),ylabel('指向性函數值');
y2=4*pi*r/lamda*sin((x-pi/6)/2);
z2=bessel(0,y2);
subplot(212);plot(x,abs(z2),'r');
title('圓陣的指向性：有相移 pi/6');xlabel('信號入射方向'),ylabel('指向性函數值');