天线原理大作业
学 生 张** 年 级 2010级 班 级 0210** 班 学 号 021012** 专 业 电子信息工程 学 院 电子工程学院
西安电子科技大学 2013年5月
1.1题目：
P63例2，两个半波振子齐平排列，利用matlab编程，画出E面、H面图和三维立体图 1.2程序代码
1. 绘制三维方向图的代码
%%%%%%% Project Name:The Pattern of Array Antenna %%%%%%% Author:WenJian Zhang %%%%%%% Date:5.15.2013 clc;clear all m=linspace(0,pi); n=linspace(0,2*pi); %生成一个等差数列 f=2*abs(cos(pi/2*sin(m).*cos(n))./sqrt(1-(sin(m).*cos(n)).^2)).*abs((cos(pi/4*(cos(m)-1)))); %二元阵方向函数 x1=(f.*sin(m))'*cos(n); y1=(f.*sin(m))'*sin(n); z1=(f.*cos(m))'*ones(size(n)); figure; surf(x1,y1,z1); axis equal title('二元阵三维方向图'); %极坐标转化为直角坐标 %天线阵三维方向图 %纵、横坐标采用等长刻度 2. 绘制E面，H面方向图的代码
%%%%%%% Project Name:The Pattern of Array Antenna %%%%%%% Author:WenJian Zhang %%%%%%% Date:5.15.2013 clc; clear all; lam=1; %波长 a=-pi/2; %二元阵振子间相位差 m=linspace(0,2*pi); %生成一个等差数列 n=linspace(0,2*pi); k=2*pi/lam; %自由空间相移常数 d=lam/4; %振子间隔 b=a+k*d.*cos(m); %振子在（m，n）方向上的相位差 fa_E=abs(cos((pi/4)*(cos(m)-1))); %半波振子在E面方向上的阵因子方向图 fa_H_1=abs(cos((pi/4)*(cos(pi/2)-1))); fa_H=fa_H_1*ones(size(n)); %半波振子在H面方向上的阵因子方向图 f1_E=2*abs(cos((pi/2).*sin(m))./cos(m)); %半波振子在E面方向上的方向函数 f1_H=2*abs(cos((pi/2).*cos(n))./sin(n)); %半波振子在H面方向上的方向函数 fE=fa_E.*f1_E; fH=fa_H.*f1_H; figure ; %二元阵E面方向函数 %二元阵H面方向函数 % 阵因子E面方向图 polar(m,fa_E,'r') title('阵因子E面方向图') figure polar(m,f1_E,'r') % 天线元在E面方向图 title('半波振子E面方向图') figure polar(m,fE,'r') figure ; % 天线阵在E面方向图 title('天线阵E面方向图') polar(n,fa_H,'r') title('阵因子H面方向图') figure % 天线元在H面方向图 % 阵因子H面方向图 polar(n,f1_H,'r') title('半波振子H面方向图') figure polar(n,fH,'r') % 天线阵在H面方向图 title('天线阵H面方向图') 1.3 运行结果