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 结论
方位维的方向图关于波束指向对称;
俯仰维的方向图不具备对称性。