《矩量法分析八木天线.doc》由会员分享，可在线阅读，更多相关《矩量法分析八木天线.doc（7页珍藏版）》请在金锄头文库上搜索。
1、如何用基于hallen方程的MoM分析八木天线？clear W$ 2C47i lamda=1; (oy N7NT a=0.0025; opVpbB me=8.85e-12; :ak mu=4*pi*(1e-7); mh75( arg=2*pi*(3e8)/lamda; M |UvO k=2*pi/lamda; Q_X.rUL0w L1=0.499*lamda;%为什么取0.499而不是0.5 ouFYvtFg L2=0.460*lamda; y)!5R3b L3=0.419*lamda; |NC*7/ d=0.20*lamda; 0dh aAqk N1=59; YLv43PL N2=59; F
2、!Q u N3=59; rpaD2& dL1=L1/(N1+1); Fc;_GG dL2=L2/(N2+1); ( _)jkI dL3=L3/(N3+1); kDN:ep/ lz1=-L1/2:dL1:L1/2; $?Z!*x lzm1=lz1(2:N1+1); DSk/q-u lz2=-L2/2:dL2:L2/2; S ;?z lzm2=lz2(2:N2+1); 4?7W+/& lz3=-L3/2:dL3:L3/2; QQ - lzm3=lz3(2:N3+1); 5bjRj6 M=N1+N2+N3; 4K(xd&Q for km=1:M muxI za if km=N1 xmN U A(km
3、).y=-d; (wH Pzq A(km).z=lzm1(km); $mE9a6 elseif (kmN1) xJZuTN A(km).y=0; 3GHj A(km).z=lzm2(km-N1); yA(K=?sq else ,JN8fag A(km).y=d; u%_ 1 A(km).z=lzm3(km-N1-N2); KiOcu=F end 0&SI5E end !MI& for m=1:M g7#w,o if m=N1 H0SQ? flagm=1; sx51Xd Vim=0; ZMLg;-T.&4 hm.y=-d; *8IWg hm.z=L1/2; njnDWSnb elseif m=N
4、1+N2 TdGnf flagm=2; ?zJOh Vim=1; zC:wNzzK hm.y=0; fH ?ha hm.z=L2/2; VL.h else U(;&(WM flagm=3; l-$uHHyu* Vim=0; W$YcE ; hm.y=d; % 9#gB hm.z=L3/2; |AwN end hwa:kD for n=1:M *cg( ?yg if n=N1 JwX_?ln flagn=1; /,CbKi,+ dL=dL1; 6k/U3&R elseif nr #% Fhmn=exp(-j*k*Rhmn)/(4*pi*Rhmn)*dL; r+i=SQax else Gr(|Ra
5、 . if flagm=flagn ELh3 Rhmn=sqrt(hm.z-A(n).z)2+a2); nD,3B# Rmn=sqrt(A(m).z-A(n).z)2+a2); 9|Sub else g3?U#7i Rhmn=sqrt(hm.y-A(n).y)2+(hm.z-A(n).z)2); C/v#cLD Rmn=sqrt(A(m).y-A(n).y)2+(A(m).z-A(n).z)2); hWcTIv end .2A9_s Fhmn=exp(-j*k*Rhmn)/(4*pi*Rhmn)*dL; &G0hCT! Fmn=exp(-j*k*Rmn)/(4*pi*Rmn)*dL; inT_
6、g end /Ij7nDl% Z(m,n)=cos(k*hm.z)*Fmn-cos(k*A(m).z)*Fhmn; u+% tPe end DA LQF?JB Zin=1/I(fedp) /wC6z ! grid on D9ywg/Q91 +xGWkltz theta=0:pi/100:2*pi;%画E面方向图 I./L)0= pha=pi/2; 0Uaem for m=1:length(theta) Q2$/e+ F(m)=0; g*r/u; for n=1:M Q=XAR if nN1&n)- polar(theta,F1,mo) kzWz4f FBR=20*log10(F(51)/F(1
7、51); M3)v- Q1K% theta=pi/2;%画H面方向图 gjO *h3 pha=0:pi/100:2*pi; FuV Yl for m=1:length(pha) p.qrf7N$ F11(m)=0; p63 for n=1:M YmdsI+DbIu if nN1&n=N1+N2 Ky =(urAd dL=dL2; 1uKD&k%q else Ubn dL=dL3; w04S figure(3) oD=XUQw( polar(pha,F21,mo); z. Ve# zSN#v$ Rin=real(Zin); r(8A:#d G=(1/(480*pi2)*arg2*mu2*max(
8、F)2/(Rin*abs(I(fedp)2); i1.5z T=(Zin-50)/(Zin+50); AH ?MJKYZ SWR=(1+abs(T)/(1-abs(T); 谢谢哦我已经知道了。那个如果反射线取0.500，会产生奇点的。矩量法分析八木天线clear rTy*W L lamda=1;%波长 emCM|NQg& a=0.0025;%振子的半径 0g0i4IV me=8.85e-12;%介电常数 WZ.UN, mu=4*pi*(1e-7);%磁导率 RhncBKm*M arg=2*pi*(3e8)/lamda;%角频率 zgxYD*v k=2*pi/lamda;%波数 : Zlwp6 N
