linearbeamdirect1.m

DBF 数字波束直接形成 MATLAB 仿真程序

%**********M阵元-线阵波束*********矩阵法
%********邓益群写于2008年10月31日******

clear all;

M=16;                      % 阵元数
c=3e8;                     % 光速
fc=6e9;                    % 工作频率（hz）
l=c/fc;                    % 波长
d=l/2;                     % 阵元间距
r=-90:0.1:90           % 入射角(rad)
r0=20                      % 波束方向
W=zeros(M,1);
A=zeros(M,1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
k=(2*pi)/l                 % 波数
u=0.5*k*d*(sin(r*pi/180)-sin(r0*pi/180))     % 空间角频

%计算方向图
for m=1:M
    A(m,1)=exp(j*(m-1)*2*pi*d/l*sin(r0*pi/180));  % 加权矢量 
end

for i=1:length(r);
    for m=1:M
        W(m,1)=exp(j*(m-1)*2*pi*d/l*sin(r(i)*pi/180)); % 控制矢量
    end
    y(i)=W'*A;%方向图函数
end
s=abs(y);
S=s/max(s);
Gain=10*log10(S);
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure 
subplot(2,2,1);plot(r,y);xlabel('theta-度');ylabel('y');title('阵列输出');axis([-90 90 -5 20]);
subplot(2,2,2);plot(r,s);xlabel('theta-度');ylabel('s=abs(y)/dB');title('波束图');axis([-90 90 0 20]);
subplot(2,2,3);plot(r,S);xlabel('theta-度');ylabel('S/dB');title('归一化波束图');axis([-90 90 0 1]);
subplot(2,2,4);plot(r,Gain);xlabel('theta-度');ylabel('Gain/dB');title('增益波束图');axis([-90 90 -60 0]);
hold off