DBF数字波束形成matlab代码实现

%DBF_scan

clear;close all;clc;

fc=10e9;%10GHz

c=3e8;%光速

lambda=c/fc;%载波波长

R=10e3;%发射机到天线阵列中心的距离

step=2;%发射机每次旋转的步进

%

d=lambda/2;%天线阵元间距,d<=半波长 符合空间采样定理

rev_num=12;%天线阵元数量

left_pos =[];

right_pos=[];

for i=1:rev_num/2

 left_pos_temp =-lambda/4-(i-1)*d;

 right_pos_temp= lambda/4+(i-1)*d;

 left_pos   =[left_pos left_pos_temp ];

 right_pos   =[right_pos right_pos_temp];

end

y_pos = zeros(1,rev_num);%全零，阵元的Y坐标

radar_position=[left_pos right_pos;y_pos];%N阵元雷达接收机(天线)坐标

for theta=30:150

 dis_radar2target=zeros(1,rev_num);%发射机到每个阵元的距离初始化

 rev_sig=zeros(1,rev_num);%每个阵元收到的信号_初始化

 target_position=[R*cosd(theta);R*sind(theta)];%发射机的位置

 for radar_idx=1:rev_num  

   dis_radar2target(radar_idx)=norm(radar_position(:,radar_idx)-target_position);%发射机到每个阵元的距离

   rev_sig(radar_idx)=exp(1i*2*pi*dis_radar2target(radar_idx)/lambda);%每个阵元收到的信号

 end

 w=conj(rev_sig);%导向矢量：角度为theta方向的天线阵列的导向矢量,每个阵元都有相应的矢量

 Rx=[];

 for i=0:step:360-1

   target_position=[R*cosd(i);R*sind(i)];

   for radar_idx=1:rev_num  

     dis_radar2target(radar_idx)=norm(radar_position(:,radar_idx)-target_position);

     rev_sig(radar_idx)=exp(1i*2*pi*dis_radar2target(radar_idx)/lambda)*w(radar_idx);%每个阵元收到的信号

   end

   Rev_temp=sum(rev_sig);

   Rx=[Rx Rev_temp];

%     %可视化

%     figure(1);

%     subplot(121)%直角坐标

%     plot(db(Rx));grid on;axis([0 180 -200 30]);

%     subplot(122)%极坐标

%     polarplot(abs(Rx));

 end  

 figure(1);

 subplot(121)%直角坐标

 plot(db(Rx));grid on;axis([0 180 -200 30]);

 subplot(122)%极坐标

 polarplot(abs(Rx));

end

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删