同轴馈电的微带天线的FDTD分析

我参照“E型宽带微带天线的时域有限差分法分析”一文中提到的
思路，用矩形网格模拟同轴线。
1）在计算同轴线中的电压和电流时（采用传输线模型），需要用到同轴线
的特性阻抗Za，请问Za应该怎么计算？
2）计算Zin时，取什么位置，我取在激励源处，妥不妥？
3）程序附
% 同轴线口径处电场
ex(30,50:52,1,1)=U(M+1)/dx;
ex(33,50:52,1,1)=-U(M+1)/dx;
ey(30:32,50,1,1)=U(M+1)/dy;
ey(30:32,53,1,1)=-U(M+1)/dy;
%****************************(同轴线内电压和电流的计算）
U1=U; %保留上一次计算的U
I(M+2)=-(((hy(33,50,2,1)+hy(33,51,2,1)+hy(33,52,2,1)-hy(30,50,2,1)-
         hy(30,51,2,1)-hy(30,52,2,1))*dy+(hx(30,50,2,1)+hx(30,51,2,1)+
         hx(30,52,2,1)-hx(30,53,2,1)-hx(31,53,2,1)-hx(32,53,2,1))*dx));
U(ks)=exp(-((ct*dt-3.0*T)^2/(T^2))); %激励源
I( 2:M+1)=I(2:M+1)-1/(2*Za)*(U(2:M+1)-U(1:M));
U(2:M+1)=U(2:M+1)-Za/2*(I(3:M+2)-I(2:M+1));
U(1)=U1(2)+(C*dt-dz)/(C*dt+dz)*(U(2)-U1(1));  %吸收边界在k'=1

申明：网友回复良莠不齐，仅供参考。如需专业解答，请学习本站推出的微波射频专业培训课程。