圆环阵方向图

1 单层圆环阵

close all;
clear;
clc;

%%   均匀圆阵的方向图
c=3e8;
fC=1e9;         % 信号频率
lamda=c/fC;     % 信号波长

M=10;           % 阵列规模
m=(0:1:M-1)';   % 第m个阵元  列矢量

theta0 = 20;    %  波束指向
fine0=180;       %  
bujing=1;
theta=0:bujing:90;
fine=0:bujing:360;


%-----同心圆环阵的阵元个数，层数等信息
R=2*lamda;          % 圆环半径
gama=(360*m/M);     % 第m个阵元和圆心O之间的连线与x轴的夹角,单位为度


%% 权值
W =exp(1j*(2*pi/lamda)*R*sind(theta0)*cosd(fine0 - gama));   % 列矢量

            
%% 方向图            
Beam_F=zeros(length(theta),length(fine));
V=[];
for p=1:length(theta)
    for q=1:length(fine)
        V =exp(1j*(2*pi/lamda)*R*sind(theta(p))*cosd( fine(q)-gama) ); % 列矢量  
        Beam_F(p,q)=W'*V;  
    end      
end

Beam_F=abs(Beam_F);
Beam_F_max=max(max(Beam_F));
Beam_F_dB=20*log10(Beam_F/Beam_F_max);
figure(1)
mesh(fine,theta,Beam_F_dB);zlim([-60 0]);
xlabel('方位角/度');ylabel('俯仰角/度');zlabel('功率/dB');


Beam_F_dB_fine=Beam_F_dB(theta0/bujing+1,:);
Beam_F_dB_theta=Beam_F_dB(:,fine0/bujing+1);
figure(2)
plot(theta,Beam_F_dB_theta);grid on;
xlabel('俯仰角/度');ylabel('功率/dB');

figure(3)
plot(fine,Beam_F_dB_fine);grid on;
xlabel('方位角/度');ylabel('功率/dB');

2 同心圆环阵

close all;
clear;
clc;


%%   同心圆环阵的方向图
c=3e8;
fc=1e9;          % 信号频率
lamda=c/fc;      % 信号波长
d=lamda/2;
M=5;             % 圆环层数
m=(0:1:M-1)';

theta0 = 45;    %  波束指向
fine0 = 80;        
bujing=1;
theta=0:bujing:90;
fine=0:bujing:360;


%-----同心圆环阵的阵元个数，层数等信息
rou_0=0;       	% 圆心
delta_rou=d;    % 相邻圆环之间的间距
rou_m=m*delta_rou+rou_0; % 第m层圆环的半径
rou_m(1)=0;
N_m=round(2*pi*rou_m/d); % 第m层圆环的阵元个数
N_m(1)=1;
Num_array=sum(N_m);      % 均匀圆环阵的阵元总个数

for m=1:M           % 第m层
   for n=1:N_m(m)   % 第n个阵元 
        angle(m,n)=360/N_m(m)*(n-1);% 第m层圆环的第n个阵元的角度,单位为度
    end
end


%%       权值
alfar=[sind(theta0)*cosd(fine0);...
       sind(theta0)*sind(fine0);cosd(theta0)];% 方向矢量

for m=1:M           % 第m层
   for n=1:N_m(m)   % 第n个阵元  
        position=[rou_m(m)*cosd(angle(m,n));rou_m(m)*sind(angle(m,n));0];
        W(m,n) =exp(1j*2*pi*fc*(alfar.'*position)/c);   % 矩阵
   end
end
W1=reshape(W,[],1); % 将8*8的矩阵变成64*1的向量

% fff
 
%% 方向图            
%---------方向矢量计算
E=[];
for p = 1:length(theta)
    for q = 1:length(fine)   
        alfar=[sind(theta(p))*cosd(fine(q)); ...
               sind(theta(p))*sind(fine(q));cosd(theta(p))];
        for m=1:M           % 第m层
           for n=1:N_m(m)   % 第m层的第n个阵元  
                position=[rou_m(m)*cosd(angle(m,n));rou_m(m)*sind(angle(m,n));0];
                V(m,n) =exp(1j*2*pi*fc*(alfar.'*position)/c);   % 矩阵
           end
        end
        V1=reshape(V,[],1); % 将8*8的矩阵变成64*1的向量
        
        Beam_F(p,q)=abs(W1'*V1);% 相控阵天线的方向图
    end
end

Beam_F_dB=20*log10(Beam_F/max(max(Beam_F)));
figure(1)
mesh(fine,theta,Beam_F_dB);zlim([-60 0]);
xlabel('方位角/度');ylabel('俯仰角/度');zlabel('功率/dB');


Beam_F_dB_fine=Beam_F_dB(theta0/bujing+1,:);
Beam_F_dB_theta=Beam_F_dB(:,fine0/bujing+1);
figure(2)
plot(theta,Beam_F_dB_theta);grid on;
xlabel('俯仰角/度');ylabel('功率/dB');

figure(3)
plot(fine,Beam_F_dB_fine);grid on;
xlabel('方位角/度');ylabel('功率/dB');



%% 极坐标图
% Beam_F_fine=Beam_F(theta0/bujing+1,:);
% Beam_F_theta=Beam_F(:,fine0/bujing+1);
% yuzhi=-70;
% for m=1:length(Beam_F_dB_fine)
%     if(Beam_F_dB_fine(m)<yuzhi)
%         Beam_F_dB_fine(m)=yuzhi;
%     else
%         Beam_F_dB_fine(m)=Beam_F_dB_fine(m);
%     end
% end
% for m=1:length(Beam_F_dB_theta)
%     if(Beam_F_dB_theta(m)<yuzhi)
%         Beam_F_dB_theta(m)=yuzhi;
%     else
%         Beam_F_dB_theta(m)=Beam_F_dB_theta(m);
%     end
% end
% 
% figure(4)
% polarplot(theta/180*pi,Beam_F_dB_theta);grid on;
% rlim([yuzhi 0]);
% thetalim([0,90]); % 修改此处，决定扇形区域
% 
% figure(5)
% polarplot(fine/180*pi,Beam_F_dB_fine);grid on;
% rlim([yuzhi 0]);
% % thetalim([-90,90]); % 修改此处，决定扇形区域

3 结论

方位维的方向图关于波束指向对称；
俯仰维的方向图不具备对称性。