[求助]葛德彪书上同轴馈电程序结果有差别是怎么回事

!u0|{6U  
( q8uB  
麻烦各位大神帮我看看下面的程序，是我摘取的一个时间步内关于同轴馈电部分的更新程序，我的参数更新步骤对吗？为什么我的程序运行的结果和书上有差别？书上计算的金属平板上圆柱天线问题的反射电压结果是这样的 bIs@CDB  

图片:IMG_20140611_160908.jpg

LS(J%\hMDm  
JIJ79HB  
1UyQ``v/  
Pgdv)i3  
下面的程序仿真结果是这样的，不知道哪里出问题了呢。。。。 #}y(D{zc  
w]-,X`  
`-,y
J  

图片:FC[FDV2XTUVPE6KO`17}T{N.jpg

)$Tcip`  
ItLP&S=  
)wFr%wNe  
T1Gy_
G/  
%更新激励源 :Ca]
/]]  
uin = exp(-((current_time - t_0)/(sqrt(2)*tau)).^2);   -hn~-Sy+  
{Kdr-aC  
%更新同轴线口径处的磁场 U`25bb1Wj  
I{rW+<)QGC  
Hy(is,js,ks)=Hy(is,js,ks)-CQ0*(Ex(is,js,ks+1)-2/log(dx/a)/dx*U(nk1+1))+CQ0*2/log(dx/a)*Ez(is+1,js,ks); 4|yZA*Q^  
Hx(is,js,ks)=Hx(is,js,ks)+CQ0*(Ey(is,js,ks+1)-2/log(dy/a)/dy*U(nk1+1))-CQ0*2/log(dy/a)*Ez(is,js+1,ks); fQ.{sQ$@h  
Hy(is-1,js,ks)=Hy(is-1,js,ks)-CQ0*(Ex(is-1,js,ks+1)+2/log(dx/a)/dx*U(nk1+1))-CQ0*2/log(dx/a)*Ez(is-1,js,ks); RCCI}ovU  
Hx(is,js-1,ks)=Hx(is,js-1,ks)+CQ0*(Ey(is,js-1,ks+1)+2/log(dy/a)/dy*U(nk1+1))+CQ0*2/log(dy/a)*Ez(is,js-1,ks); .PAR  
7E4Xvg
+c  
%更新同轴线电流 k2axGq  
Qx{[#[Da  
I(1:nk1)=I(1:nk1)-1/(2*Za)*(U(2:nk1+1)-U(1:nk1)); %+1/Za*c*dt/dz*uin(1:2:2*nk1-1); L~Hgf/%5  
I(nk1+1)=dy*(Hy(is,js,ks)-Hy(is-1,js,ks))-dx*(Hx(is,js,ks)-Hx(is,js-1,ks)); ]Sl]G6#Iwv  
I(nks)=I(nks)-1/(2*Za)*(U(nks+1)-U(nks))+1/Za*c*dt/dz*uin; )'pc1I  
e 9p+  
%再次更新激励源 XwerQwO=  
uin = exp(-((current_time - t_0)/(sqrt(2)*tau)).^2); X!'nfN  
'OERW|BO  
%更新同轴线电压
6e%|.}U  
PR{y84$  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

为什么都没人啊

还是没有人 lt4UNJ3w  
求解答

还是没有人 %aBJ+V F  
求解答