gongshishixian.m

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

%用公式法实现圆阵指向性图的绘制
close all
 clear all;
 N=60;
 c=1500;
 f=1500;
 alpha=2*pi/N;
 lamda=c/f;
r=lamda/0.4;
 x=0:pi/180:pi;
 theter=[0 pi/6];
 for j=1:2
     a=0;b=0;
 for i=1:60
    derta=r/c*[cos(x-(i-1)*alpha)-cos(theter(j)-(i-1)*alpha)];
     a=a+cos(2*pi*f.*derta);
     b=b+sin(2*pi*f.*derta);
 end
R=sqrt(a.^2+b.^2)/N;
 if j==1
     subplot(211);plot(x,abs(R));title('圆阵的指向性：无相移');xlabel('信号入射方向'),ylabel('指向性函数值');
 else 
subplot(212);plot(x,abs(R));title('圆阵的指向性：有相移');xlabel('信号入射方向'),ylabel('指向性函数值');
end
end