直线阵和圆阵数字波束形成MatlAB程序

clear all; close all; %程序用于均匀圆阵天线的数字波束形成 %设定参数：波长（lambda），主瓣方向（theta0）， %圆阵直径（d），天线数目（N） %圆心到1号天线的射线为0度参照方向，顺时针递增 %以圆心为相位基准点 lambda=10; %波长，（m） d=lambda*2; %天线间距，（m） theta0=pi; %主瓣方向，（rad） N=16; %天线数目 k=2*pi/lambda; %波数，（/m） r=d/2; %半径 alpha=2*pi/N; %两相邻天线到圆心的夹角 theta=(0:pi/180:2*pi)'; %间隔为1度 L=length(theta); echo=zeros(L,N); %回波 for n=1:N echo(:,n)=exp(1i*k*r*cos(theta-(n-1)*alpha)); end n=(1:N)'; weigh=exp(-1i*k*r*cos(theta0-(n-1)*alpha)); %加权值 BF=echo*weigh; %波束形成 plot(theta,abs(BF)/N); xlabel('Direction of arrival,rad'); ylabel('Normalized amplitude'); xlim([0,2*pi]); ylim([0,1]);