[求助]求助：ElectroMagnetic Simulation Using The FDTD Method 书中程序3.2为什么没有出书上的图形

如题：哪位前辈能帮我看看下面的程序有什么错误，为什么出不了书上的图形，下面是我自己稍微改的程序。万分感谢啊！！！ z6>ZV6(d2^  
附件也是c程序 ePcI^}{  
0g*r!aa  
#include<stdio.h> (Q p]0  
#include<math.h> uD*s^  
#include<stdlib.h> v%PWr5]  
#define IE 150 W;AWO0+  
#define JE 150 Gv2./<{#  
main() ~[/c'3+4qn  
{ t-gNG!B  
  float ga[IE][JE],dz[IE][JE],ez[IE][JE],hx[IE][JE],hy[IE][JE]; :@807OYzy  
  int l,n,i,j,ic,jc,nsteps,npml; UFY~D"%/  
  float ddx,dt,T,epsz,pi,epsilon,sigma,eaf; Appz1q  
  float xn,xxn,xnum,xd,curl_e; {*r$m>HpM  
  float t0,spread,pulse; $6x:aG*F  
  float gi2[IE],gi3[IE],gj2[JE],gj3[JE]; ^HN  
  float fi1[IE],fi2[IE],fi3[IE]; D"XQ!1B%  
  float fj1[JE],fj2[JE],fj3[JE]; |>fS"u  
  float ihx[IE][JE],ihy[IE][JE]; 1?#p !;&  
  FILE *fp; `h{mj|~  
  ic=IE/2; 4a50w:Jy]  
  jc=JE/2; F p=Q$J|  
  ddx=0.1; ^3@a0J=F  
  dt=6e-10; $izpH  
  epsz=8.8e-12; .AI'L|FQ%c  
  pi=3.14159;  #=~1hk  
  GBOmVQ $Hb  
  for(j=0;j<JE;j++) ?)QBJ9F  
  { -3%)nV  
   printf("%2d",j); $m0x8<7nu  
    for(i=0;i<IE;i++) 6XCX#4'i%  
 { Bst>9V&R  
   dz[j]=0.0; Tq_X8X#p  
      hx[j]=0.0; }Z$G=;3#  
      hy[j]=0.0; L<>NL$CrN  
      ihx[j]=0.0; }eRG$)'  
      ihy[j]=0.0; adPU)k_j:  
      ga[j]=1.0; SIVLYi  
      printf("%5.2f",ga[j]); (G
Orfr  
   } 92ww[+RQ@  
printf("\n"); iwx0V  
} Hh=::Bi  
   for(i=0;i<IE;i++) ~W2&z]xD  
  { \
#c+vfq  
    gi2=1.0; yacGJz^f=  
    gi3=1.0; MxA'T(Ay  
    fi1=0.0; v?=VZ~`O(  
    fi2=1.0; P\0%nyOG(%  
    fi3=1.0; i1\ /\^  
  } bc}OmPE  
for(j=0;j<JE;j++) T@.m^|~  
 { t>u9NZt G  
  gj2[j]=1.0; V>j`  
  gj3[j]=1.0; f]65iE?x  
  fj1[j]=0.0; jY6=+9Jz5  
  fj2[j]=1.0; rd~W.b_b  
  fj3[j]=1.0; g;Q^
_4@  
  } -#M~NbI,  
printf("Number of PML cells-->"); l'8TA~  
scanf("%d",&npml); $ttr_4=  
for(i=0;i<npml;i++) (G!J=
=  
{ q x }fn/:  
  xnum=npml-i; -"fq34v  
  xd=npml; j?)`VLZ  
  xxn=xnum/xd; 
aq| [g  
  xn=0.33*pow(xxn,3.0); BCB/cBE  
  printf("%d %7.4f %7.4f \n",i,xxn,xn); 6Tmb@<I_  
   gi2=1.0/(1.0+xn); 3t6'5{  
   gi2[IE-1-i]=1.0/(1.0+xn); =C2C~Xd  
   gi3=(1.0-xn)/(1.0+xn); r>@/XYK&\  
   gi3[IE-i-1]=(1.0-xn)/(1.0+xn); T>o# *{qn  
 xxn=(xnum-0.5)/xd; )F2tV ]k\  
 xn=0.25*pow(xxn,3.0); =9MH
  
  fi1=xn; ^K[tO54  
  fi1[IE-2-i]=xn; H pfI  
  fi2=1.0/(1.0+xn); 1OeDWEcB  
  fi2[IE-2-i]=1.0/(1.0+xn); {oeQK  
  fi3=(1.0-xn)/(1.0+xn); zXPj7K*  
  fi3[IE-2-i]=(1.0-xn)/(1.0+xn); Q'M Ez  
 } OB@t(KNx*P  
 for(j=0;j<npml;j++) .HJHJ.Js8X  
 { lY?d*qED  
   xnum=npml-j; ]}c=U@D,9  
   xd=npml; 2xj`cFT  
   xxn=xnum/xd; \{mJO>x  
   xn=0.33*pow(xxn,3.0); /YR$#&N2  
   printf("%d %7.4f %7.4f\n",i,xxn,xn); $^ 3 f}IzA  
    gj2[j]=1.0/(1.0+xn); )q-!5^ak  
    gj2[JE-1-j]=1.0/(1.0+xn); IsL/p3|  
    gj3[j]=(1.0-xn)/(1.0+xn); GQNiBsV  
    gj3[JE-j-1]=(1.0-xn)/(1.0+xn); ?OdA`!wE  
   xxn=(xnum-0.5)/xd; bcs(#  
   xn=0.25*pow(xxn,3.0); ^yZEpQN_  
     fj1[j]=xn; I2Rp=L:z5  
     fj1[JE-2-j]=xn; U<gw<[>f  
     fj2[j]=1.0/(1.0+xn); ]gk1h=Y~h  
     fj2[JE-2-j]=1.0/(1.0+xn); rnaDo\5  
     fj3[j]=(1.0-xn)/(1.0+xn); }PuO$ L  
     fj3[JE-2-j]=(1.0-xn)/(1.0+xn); @SB+u+mOS  
  } >GR
uS\B  
printf("gi+fi\n"); E/ )+hK&  
for(i=0;i<IE;i++) )j*qGsOg  
  { 8r,%!70  
   printf("%2d %5.2f %5.2f\n",i,gi2,gi3); }H.vH  
   printf("%5.2f %5.2f %5.2f\n",fj1[j],fj2[j],fj3[j]); EHjhez  
} j(2T,WM  
  t0=40.0; je%12DM  
  spread=12.0; e:R[  
  T=0;  {`  
  nsteps=1; o`0H(\en  
  31Du@h8YX  
while(nsteps>0) (k45k/PAP  
{ ?J-KB3Uv3  
  printf("nsteps-->"); x{5*%}lX8  
  scanf("%d",&nsteps); -l
Eh}r  
  printf("%d\n",nsteps); gH.^NO5\'  
  SQx):L)P6  
 for(n=1;n<=nsteps;n++) [i_x 1  
 { ?K#$81;[  
  T=T+1; J n/=v\K@  
for(j=0;j<JE-1;j++) y9<Fv|Ric  
{ 9$i`B>C
~  
  for(i=0;i<IE;i++) I78pul8!  
   { N*W.V,6yH  
    curl_e=ez[j]-ez[j+1]; gQwmYe  
    ihx[j]=ihx[j]+curl_e; $Yc9><i  
    hx[j]=fj3[j]*hx[j]+fj2[j]*0.5*curl_e+fi1*ihx[j]; T]`" Xl8  
    } C3gz)!3  
} _=#mmZkq  
for(j=0;j<=JE-1;j++) ( PlNaasV  
{  ,Fiiw  
  for(i=0;i<IE-1;i++) +(0eOO'\M  
  { 2Mp;/b!  
    curl_e=ez[i+1][j]-ez[j]; X=~QE}x  
    ihy[j]=ihy[j]+curl_e; IgJC>;]u  
    hy[j]=fi3*hy[j]+fi2*0.5*curl_e+fj1[j]*ihy[j]; [8#l~ |U  
    } 6 [E"  
  }  `SrVMb(  
  for(j=1;j<IE;j++) @RW%EXKt  
 { IuD<lMeJJ  
   for(i=1;i<JE;i++) SF7 Scd  
    { HS9U.G>  
      dz[j]=gi3*gj3[j]*dz[j]+gi2*gj2[j]*0.5*(hy[j]-hy[i-1][j]-hx[j]+hy[j-1]); RPwSo.c4
 
   } #u^d3 $Nj  
} *^()
el,d  
pulse=exp(-0.5*(pow((t0-T)/spread,2.0))); ]ghPbS@  
dz[ic][jc]+=pulse; jRc#>;dN  
for(j=0;j<JE;j++) -`&;3 7  
{ i9tM]/SP  
 for(i=0;i<IE;i++) 4LLCb7/5lP  
  { (z{xd  
   ez[j]=ga[j]*dz[j]; f5Oh#  
   } )-gyDA
 
} jT',+  
printf("%3f %6.2f\n",T,ez[ic][jc]); 6WZp&pO  
for(j=0;j<JE-1;j++) ;{RQ+ZX'[  
{ rD gl@B
3  
  ez[0][j]=0.0; >{Q2S  
  ez[IE-1][j]=0.0; Ns(F%zkm  
 } j/&7L@Y  
for(i=0;i<IE-1;i++) <LY+" Y  
 { `PfC:L  
   ez[0]=0.0; YR2/`9s\QJ  
   ez[JE-1]=0.0; ML-g"wv  
} ?*=Jq  
hrK^oa_[W  
} *6DKUCA/  
for(j=0;j<JE;j++) ,@xZuq+K<  
{ X,~C&#  
  printf("%2d",j); u#<]>EtbB  
  for(i=0;i<IE;i++) !WXSrICX[  
   { A:(uK>5{Kk  
    printf(" % .. 2z:9^a/]Na  
F
44")fY  

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

有没有人给看看啊，求助啊

具体说下那个图啊，问题具体在什么地方，不然没事喜欢看枯燥的code

就是electromagnetic simulation using fdtd method (Sullivan Dennis)书中程序3.2出的图，我传到附件里了结果自己做出的却是第二张图，完全不一样啊

大概一个月前看的这本书，发现代码有些问题的。我改了一下就好了

我当时是把calculate the Hy field那段并入到Hx那段里一起去算然后就没事了。。。 Z ZMz0^V  
也不知道是什么原因，可能是我编译器有问题吧。 GY~Q) Z  
lz你先运行一两步看看结果吧，慢慢分析应该是可以找出原因的。 Hy*_4r  
我还要去上课，实在是懒得看代码了~当时在准备期末被它弄得晕死了