t=1:n*N; figure(1);
plot(t,abs(si)),ylabel('I'),title('电流分布') in=U'*(Z\\V); i=V'*si; Zin=1/i
theta=(-pi:pi/100:pi)+eps; for m=1:length(theta)
E1=-j*w*u*exp(-j*k).*exp(j*k.*sqrt(mid(:,1).^2+mid(:,2).^2).*cos(abs(atan(mid(:,1)./(mid(:,2)+eps))-theta(m)))).*mid(:,3).*sin(theta(m))/(4*pi); Etheta(m)=E1'*si; end
Etheta=Etheta./max(Etheta); figure(2);
polar(theta,abs(Etheta)/max(abs(Etheta))),title('E平面方向图 (\\Phi = 0)'); Lo=find((abs(Etheta-1/sqrt(2))<0.05)==1);
G=abs(4*pi.*Etheta.*conj(Etheta)/(y*real(Zin).*si((N+1)/2).*conj(si((N+1)/2)))); Gmax=max(G)
phi=(0:pi/100:2*pi)+eps; for m=1:length(theta)
E2=-j*w*u*exp(-j*k).*exp(j*k.*sqrt(mid(:,1).^2+mid(:,2).^2).*cos(abs(atan(mid(:,1)./(mid(:,2)+eps))-theta(m)))).*mid(:,3).*sin(theta(m))/(4*pi); Ephi(m)=E2'*abs(si); end figure(3);
polar(theta-pi/2,(abs(Etheta)/120/pi)/max(abs(Etheta)/120/pi)),title('H平面方向图 (\\theta = \\pi / 2)');
16
5.2 MATLAB仿真结果
图(5) ?=0时，E面方向图
图（5） ???/2时，H面方向图
17
图（6） 电流分布
该天线输入阻抗为：
Zin = 1.4271e+003 -2.5439e+002i 该天线最大增益为：
Gmax =115.2401
六、设计总结
从仿真结果可以看出，该5元八木天线基本满符合天线的基本功能。
本次课程设计是一个艰难的过程，因为这是一个综合知识的运用，并且要使用MATLAB编程、仿真，这就要求我们要有坚实的课本知识，也要对MATLAB软件有比较深入的了解。刚开始动手的时候，不知从何下手，但通过自己查阅资料，以及在老师的知道下，才有了思路。
通过这次课程设计，使我对天线的原理有了更深的认识，更重要的是在设计的过程中收获了许多感悟。例如，动手做课程设计前一定要先仔细审题，提出切实可行的方案，如果这一点都没做好，那么发再多时间做设计都是徒劳。生活中做任何事情都应该三思而后行。还要，要多虚心请教老师，在老师或者是能力比你强的人的指导下，会让你少走很多弯路。
18
七、参考资料
[1] 邓华. MATLAB通信仿真及应用实例详解. 人民邮电出版社. 2003 [2] 宋铮、张建华、黄冶. 天线与电波传播. 西安电子科技大学出版社.2011
19