besselshixian.m

来自「均匀圆阵的指向性仿真和模拟」· M 代码 · 共 19 行

%用零阶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('指向性函数值');