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

#V,~d&_k  
w,#W&>+&  
麻烦各位大神帮我看看下面的程序，是我摘取的一个时间步内关于同轴馈电部分的更新程序，我的参数更新步骤对吗？为什么我的程序运行的结果和书上有差别？书上计算的金属平板上圆柱天线问题的反射电压结果是这样的 j#>![km Mu  

图片:IMG_20140611_160908.jpg

xr3PO?:  
]fS~N9B  
GZ[h`FJg/  
Py0i%pZ  
下面的程序仿真结果是这样的，不知道哪里出问题了呢。。。。 }fMFQA)  
*>*/|  
b0=AQ/:  

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

$ 5-2cL  
lGBdQc]IL  
[)6E)E`_e  
G`lhvpifG  
%更新激励源 tsC|R~wW  
uin = exp(-((current_time - t_0)/(sqrt(2)*tau)).^2);   Q pc^qP^-  
w8#>xV^~  
%更新同轴线口径处的磁场 ~zC fan/  
z>g& ?vo2  
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); kQ'xs%Fw  
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); kJpHhAn4  
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); 9@*>$6  
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); n.UM+2G  
{Qmb!`F  
%更新同轴线电流 19{?w6G<k  
[a*>@IR  
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); Wv/5#
_  
I(nk1+1)=dy*(Hy(is,js,ks)-Hy(is-1,js,ks))-dx*(Hx(is,js,ks)-Hx(is,js-1,ks)); #_sVB~sn@  
I(nks)=I(nks)-1/(2*Za)*(U(nks+1)-U(nks))+1/Za*c*dt/dz*uin; CIvT5^}  
<).qe Z  
%再次更新激励源 (gC^5&11  
uin = exp(-((current_time - t_0)/(sqrt(2)*tau)).^2); fTtSx_}3H  
aNU%OeQA  
%更新同轴线电压
x(N}^Hu  
X.Y)'qSf  

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

为什么都没人啊

还是没有人 Bm^vKzp  
求解答

还是没有人 1x V~EX  
求解答