相控阵天线面阵方向图（三种表示方法）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   均匀面阵方向图的三种不同的表达式
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clc;
close all;
clear all;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%             参数设置
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% theta是俯仰角 fine是方位角
c=3e8;       % 光速
f=3e10;     % 载波频率
lamda=c/f;  % 波长
d=lamda/2;   % 阵元间距
M=8;
N=8;     % 8*8的矩阵
% 波束指向角度(theta0,fine0)
theta0=30;
fine0=200;
bujing=0.5;
theta=0:bujing:90;   % 俯仰角取值范围
fine=0:bujing:360;   % 方位角取值范围%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%           第1种表达式
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%-----------权值
W =[];
for m=1:Mfor n=1:Nx = (m-1)*sind(theta0)*cosd(fine0)+(n-1)*sind(theta0)*sind(fine0);% 波程差W(m,n) = exp(1j*2*pi*f*d*x/c);    end
end
W1=reshape(W,[],1); % 将8*8的矩阵变成64*1的向量%---------方向矢量计算
E=[];
for p = 1:length(theta)for q = 1:length(fine)     for m=1:Mfor n=1:N%-------权值扫描x = (m-1)*d*sind(theta(p))*cosd(fine(q))+(n-1)*d*sind(theta(p))*sind(fine(q));V(m,n) = exp(1j*2*pi*f*x/c); %  方向矢量的求解endendV1=reshape(V,[],1); % 将8*8的矩阵变成64*1的向量E1(p,q)=W1'*V1;% 相控阵天线的方向图Beam_F1(p,q)=abs(E1(p,q));end
end
figure(4)
mesh(fine,theta,Beam_F1);
xlabel('方位角');ylabel('俯仰角');zlabel('方向图');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%           第二种表达式
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for p=1:length(theta)      % 俯仰角取值范围for q=1:length(fine)   % 方位角取值范围% 64个阵元功率叠加for m=1:M  for n=1:N  r(m,n)=exp(1i*2*pi/lamda*( (m-1)*d*(sind(theta(p))*cosd(fine(q))-sind(theta0)*cosd(fine0))...+(n-1)*d*(sind(theta(p))*sind(fine(q))-sind(theta0)*sind(fine0)) )  );endendE(p,q)=sum(sum(r)); % 将64个阵元的合成场强放到对应的扫描角度上end
end
Beam_F2=abs(E);      % 取实部
figure(2)
mesh(fine,theta,Beam_F2); hold on;
xlabel('方位角');ylabel('俯仰角');zlabel('方向图'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%           第三种表达式
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for p=1:length(theta)      % 俯仰角取值范围for n=1:length(fine)   % 方位角取值范围t1=2*pi/lamda*d*( sind(theta(p)) * cosd(fine(n)) - sind(theta0)*cosd(fine0) );beam1(p,n)=sin(M/2*t1)/sin(1/2*t1);t2=2*pi/lamda*d*( sind(theta(p)) * sind(fine(n)) - sind(theta0)*sind(fine0) );beam2(p,n)=sin(N/2*t2)/sin(1/2*t2);f3(p,n)= beam1(p,n)*beam2(p,n);end
end
Beam_F3=abs(f3);% 取实部
figure(3)
mesh(fine,theta,Beam_F3);
xlabel('方位角');ylabel('俯仰角');zlabel('方向图');max(max(Beam_F1-Beam_F2))
max(max(Beam_F1-Beam_F3))
max(max(Beam_F2-Beam_F3))

相控阵天线面阵方向图（三种表示方法）相关推荐

矩阵、坐标变换、相控阵天线拟合方向图代码
矩阵.坐标变换.相控阵天线拟合方向图代码文章目录矩阵.坐标变换.相控阵天线拟合方向图代码矩阵 Matrix.h Matrix.cpp 坐标变换 CConvert.h CConvert.cpp ...
tableau货架图制作_举个栗子！Tableau技巧（33）：实用漏斗图的三种制作方法
漏斗图适用于业务流程比较规范.周期长.环节多的流程分析,通过漏斗各环节业务数据的比较,能够直观地发现和说明问题所在. 零售电商行业的用户提出分析需求:访客数量→阅览数量→下单数量→付款数量→完成数量, ...
win 7计算机图标变了,win7系统图标变黑的三种解决方法（图文）
一些Win7旗舰版系统用户说桌面图标变黑了,也不知道怎么回事?严重影响桌面美观,这是由于下载一些有Bug的软件的缘故或者中了木马会导致图标缓存出错,具体表现为图标变黑的现象.那么有什么办法解决win7 ...
无限轮播图三种实现方式
无限轮播图三种实现方式:(源码地址) 1.首尾相接方式: https://github.com/jakajacky/DRScrollerView 2.imageView复用方式: https://gi ...
关于嘉立创拼板加工和三种拼板方法
关于嘉立创拼板加工和三种拼板方法一.嘉立创加工规则: 1. 拼板要求: 需嘉立创帮拼板的,如弧形和圆形板子不接受拼版,嘉立创统一是同一方向顺拼版,不处理倒扣旋转,或加邮票孔拼版,此类外形需麻烦客户提 ...
颗粒物检测仪常用的三种检测方法
颗粒物检测仪的测量原理:仪器由组装好的传感器和数据处理器组成.传感器是本装置数据采集的重要部件,原理是将激光束通过非球面镜组变为功率密度均匀分布的细小测量光束,在光束道的侧前方,前方焦点对准光束道,后 ...
视频剪辑完成，应该如何给视频配音？三种配音方法快来学
对短视频来说声音非常重要,因此为视频配音是每一位做短视频的作者都应当掌握的一个技能,我们很多时候看一段视频停不下来,就是因为视频中的bgm或者配音很好听,让我们忍不住把一个视频反复刷很多遍.所以要想视 ...
应用程序错误电脑黑屏_电脑运行程序出现APPCRASH错误的三种解决方法
好多小伙伴说电脑在运行程序时出现appcrash错误,导致软件无法正常运行,其实电脑是很复杂的程序代码设计的,遇到一些奇怪的问题也不奇怪,系统运行程序出现APPCRASH错误,很多新手用户都不知道怎么 ...
计算机网络学习19：IPV4的三种划分方法及应用规划
IPV4地址的概述后面对三种划分方法进行详细讨论. 由于32比特的IPV4地址不方便阅读记录等,采用点分十进制表示方法表示使用. 点分十进制表示方法 8位无符号二进制转10进制的方法分类编制的IP ...

相控阵天线面阵方向图（三种表示方法）

相控阵天线面阵方向图（三种表示方法）相关推荐

最新文章

热门文章