矩量法分析八木天线.docx

矩量法分析八木天线
如何用基于hallen方程的MoM分析八木天线，
clear
lamda=1;
a=;
me=-12;
mu=4*pi*(1e-7); arg=2*pi*(3e8)/lamda; k=2n(theta)*sin(pha(m))+A(n).z*cos (theta)))*dL*(-sin(thet
a));
end
F21二abs(Fll);%方向图函数
F21二F21/max(F21);%归一化方向图函数
figure(3)
polar(pha,F21,'mo');
Rin=real(Zin);
G=(l/(480*pi"2)) *arg"2 *m「2 *max(F厂2/(Rin *abs(I(fedp)厂2);
T=(Zin-50)/(Zin+50); SWR=(1+abs(T))/(1-abs(T)); 谢谢哦~~我已经知道 了。，会产生奇点的。
矩量法分析八木天线
clear
lamda=1;% 波长
a=;%振子的半径
me=-12;% 介电常数
mu=4 *pi *(le-7);% 磁导率
arg=2 *pi *( 3e8)/lamda;% 角频率
k=2 *pi/lamda;% 波数
N2=51;
L1=*lamda;
L2=*lamda;
L3=*lamda;
d=*lamda;
dL=L2/(N2+1);
N1=fix(L1/dL);
dL1=rem(L1,dL);
N3=fix(L3/dL); dL3=rem(L3,dL);
l1=L1/2-dL1/2; l2=L2/2-dL/2;
l3=L3/2-dL3/2; lzintal1=-l1:dL:l1; lzm1=lzintal1(1:N1)+dL/2;
lzintal2=-l2:dL:l2;
lzm2=lzintal2(1:N2)+dL/2; lzintal3=-l3:dL:l3;
lzm3=lzintal3(1:N3)+dL/2; M=N1+N2+N3;
for km=1:M if km<=N1 A(km).y=-d;
A(km).z=lzm1(km); elseif km<=N1+N2 A(km).y=0;
A(km).z=lzm2(km-N1); else
A(km).y=d; A(km).z=lzm3(km-N1-N2); end
for m=1:M
for n=1:M
if n==m
Fmnmm=(1/(2*pi*dL))*log(dL/a)-j*k/(4*pi);
Fmnee=(1/(2*pi*dL))*log(dL/a)-j*k/(4*pi);
Fmnss=(1/(2*pi*dL))*log(dL/a)-j*k/(4*pi);
Fmnse=exp(-j*k*dL)/(4*pi*dL);
Fmnes=exp(-j*k*dL)/(4*pi*dL);
elseif abs(n-m)==1
Rmm=sqr t( (A(m).y—A(n).y厂2+(A(m).z—A(n).z厂2);
Fmnmm=exp(-j*k*Rmm)/(4*pi*Rmm);
Fmnee=exp(—j*k*Rmm)/(4*pi*Rmm);
Fmnss=exp(—j*k*Rmm)/(4*pi*Rmm);
if m<n
Rmnse二sqr t((A(m).y-A(n).y厂2+((A(m).z-dL/2)-(A(n).z+dL/2)厂2);
Fmnse=exp(—j*k*Rmnse)/(4*pi*Rmnse);
if (m==N1)||(m==N1+N2)
Rmnes=sqr t((A(m).y-A(n).y厂2+((A(m).z+dL/2)-(A(n).z-dL/2)厂2);
Fmnes=exp(—j*k*Rmnes)/(4*pi*Rmnes);
else
Fmnes=(1/(2*pi*dL))*log(dL/a)—j*k/(4*pi);
end
elseif m>n
Rmnes=sqr t((A(m).y-A(n).y厂2+((A(m).z+dL/2)-(A(n).z-dL/2)厂2);
Fmnes=exp(-j*k*Rmnes)/(4*pi*Rmnes);
if (n==N1)||(n==N1+N2)
Rmnse二sqr t((A(m).y-A(n).y厂2+((A(m).z-dL/2)-(A(n).z+dL/2)厂2);
Fmnse=exp(-j*k*Rmnse)/(4*pi*Rmnse);
else
Fmnse=(1/(