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

(:|rCZC  
n8.Tag(#  
麻烦各位大神帮我看看下面的程序，是我摘取的一个时间步内关于同轴馈电部分的更新程序，我的参数更新步骤对吗？为什么我的程序运行的结果和书上有差别？书上计算的金属平板上圆柱天线问题的反射电压结果是这样的 \c\z 6;j  

图片:IMG_20140611_160908.jpg

(7*((  
th^&wp  
~(x"Y\PEu  
G_E \p%L>]  
下面的程序仿真结果是这样的，不知道哪里出问题了呢。。。。 0% zy 6{  
QFg{.F?3q>  
3+WmM4|  

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

W @]t  
2EC<8}CG  
&Cm$%3  
>mW*K _~  
%更新激励源 Koi-b  
uin = exp(-((current_time - t_0)/(sqrt(2)*tau)).^2);   nrM_ay  
Vpe\Okt:  
%更新同轴线口径处的磁场 2uzy]faM  
?'^dYQ4  
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); :>*0./hG  
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); QB<~+dW  
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); 3Hi[Y[O`%P  
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); oU/{<gs  
v/GZByco>  
%更新同轴线电流 <i ";5+  
1~
ZFkcV_C  
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); N Wf IRL  
I(nk1+1)=dy*(Hy(is,js,ks)-Hy(is-1,js,ks))-dx*(Hx(is,js,ks)-Hx(is,js-1,ks)); w!r
w%  
I(nks)=I(nks)-1/(2*Za)*(U(nks+1)-U(nks))+1/Za*c*dt/dz*uin; >|rU*+I`  
Fql|0Fq  
%再次更新激励源 L#`Vr$  
uin = exp(-((current_time - t_0)/(sqrt(2)*tau)).^2); 3y}8|ML  
FT Ytf4t  
%更新同轴线电压
G0lg5iA<fC  
r E&}B5PN=  

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

为什么都没人啊

还是没有人 IV|})[n*  
求解答

还是没有人 xZp`Ke
!  
求解答