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

如题：哪位前辈能帮我看看下面的程序有什么错误，为什么出不了书上的图形，下面是我自己稍微改的程序。万分感谢啊！！！ ::g"dRS<v  
附件也是c程序 X1]&j2WR  
FAjO-T4(  
#include<stdio.h> x1Q}B   
#include<math.h> <^UB@'lCm  
#include<stdlib.h> k;3Bv
6  
#define IE 150 D
IB Az s  
#define JE 150 =   
main() Cfyas'  
{ f-y4V}  
  float ga[IE][JE],dz[IE][JE],ez[IE][JE],hx[IE][JE],hy[IE][JE]; ?1xBhKq  
  int l,n,i,j,ic,jc,nsteps,npml; DifRpj I-0  
  float ddx,dt,T,epsz,pi,epsilon,sigma,eaf; nZfTK>)A0  
  float xn,xxn,xnum,xd,curl_e; o$'Fz[U  
  float t0,spread,pulse; #^gn,^QQ  
  float gi2[IE],gi3[IE],gj2[JE],gj3[JE]; (zWzF_v  
  float fi1[IE],fi2[IE],fi3[IE]; 95_?F7}9  
  float fj1[JE],fj2[JE],fj3[JE]; Cnd*%CPZ  
  float ihx[IE][JE],ihy[IE][JE]; *qAF
#  
  FILE *fp; 82&JYx  
  ic=IE/2; ?g}n$%*5y!  
  jc=JE/2; jZ7/p^c5R  
  ddx=0.1;  @zEEX9U  
  dt=6e-10; qBpv[m  
  epsz=8.8e-12; %/,PY>:|  
  pi=3.14159; 9xB^dKM3  
  'd]9u9u  
  for(j=0;j<JE;j++)
e!l!T@ pf  
  { kWF4k  
   printf("%2d",j); GIWgfE?  
    for(i=0;i<IE;i++) q#pBlJ.LK  
 { 79~,KFct  
   dz[j]=0.0; q =b.!AZy  
      hx[j]=0.0; ZFOYYht  
      hy[j]=0.0; N:)`+}  
      ihx[j]=0.0; qDqIy+WR  
      ihy[j]=0.0; ~=Z&l  
      ga[j]=1.0; &jl'1mZ  
      printf("%5.2f",ga[j]); RecA?-0  
   } -3/:Dk`3  
printf("\n"); G&0&*mp  
} G-G!c2o  
   for(i=0;i<IE;i++)  V-}d-Y  
  { %;!@\5$  
    gi2=1.0; Ig'Y]%Z0  
    gi3=1.0; owCQ71Q  
    fi1=0.0; aj20, w  
    fi2=1.0; yrdJX  
    fi3=1.0; t^Aios~F  
  } r[eZV"  
for(j=0;j<JE;j++) wHBkaPO!  
 { =9\=5_V  
  gj2[j]=1.0; oZ& ns!#  
  gj3[j]=1.0; W:3u$LTf*f  
  fj1[j]=0.0; >L$y|8O  
  fj2[j]=1.0; H #J"'  
  fj3[j]=1.0; DvG.G+mo#  
  } m1gJ"k6 `j  
printf("Number of PML cells-->"); O*z x{a6  
scanf("%d",&npml); YZ7rs]A  
for(i=0;i<npml;i++) gu/eC  
{ [NFg9y;{h  
  xnum=npml-i; t+qLQY}=  
  xd=npml; H'<9;bD -  
  xxn=xnum/xd; %\}|&z6  
  xn=0.33*pow(xxn,3.0); 3LG}x/l  
  printf("%d %7.4f %7.4f \n",i,xxn,xn); lT@5=ou[  
   gi2=1.0/(1.0+xn); nqBZp N^  
   gi2[IE-1-i]=1.0/(1.0+xn); +IuV8XT2(  
   gi3=(1.0-xn)/(1.0+xn); TIW6v4  
   gi3[IE-i-1]=(1.0-xn)/(1.0+xn); s Ce7ni  
 xxn=(xnum-0.5)/xd; 2K.. ;A$  
 xn=0.25*pow(xxn,3.0); 8doT`rI1  
  fi1=xn; OZbwquF@  
  fi1[IE-2-i]=xn; HoT5 5v!o  
  fi2=1.0/(1.0+xn); V{HZ/p_Y  
  fi2[IE-2-i]=1.0/(1.0+xn); x1/Usupi  
  fi3=(1.0-xn)/(1.0+xn); ~1S7\e7{  
  fi3[IE-2-i]=(1.0-xn)/(1.0+xn); Z"lL=0rY/  
 } ?p. dc~tZ  
 for(j=0;j<npml;j++) j5R= K*y  
 { q+~z#
jFX  
   xnum=npml-j; ynw^nmM  
   xd=npml; $sU?VA'h  
   xxn=xnum/xd; Jb6rEV>  
   xn=0.33*pow(xxn,3.0); ~r]ZD)  
   printf("%d %7.4f %7.4f\n",i,xxn,xn); {uM0J$P:  
    gj2[j]=1.0/(1.0+xn); L-Xd3RCD  
    gj2[JE-1-j]=1.0/(1.0+xn); Umv_{n`  
    gj3[j]=(1.0-xn)/(1.0+xn); b]g}h
 
    gj3[JE-j-1]=(1.0-xn)/(1.0+xn); <eO 7b6_  
   xxn=(xnum-0.5)/xd; -v&Q'a  
   xn=0.25*pow(xxn,3.0); D,mFme  
     fj1[j]=xn; +=$\7z>s  
     fj1[JE-2-j]=xn; j*I0]!-  
     fj2[j]=1.0/(1.0+xn); wC{=o`v  
     fj2[JE-2-j]=1.0/(1.0+xn); ({m["d  
     fj3[j]=(1.0-xn)/(1.0+xn); (*/P~$xIj  
     fj3[JE-2-j]=(1.0-xn)/(1.0+xn); 6"|PJ_@P  
  } Sj+gf~~  
printf("gi+fi\n"); yZb@  
for(i=0;i<IE;i++) ;p BXAl  
  { '$ s:cS`=  
   printf("%2d %5.2f %5.2f\n",i,gi2,gi3); /Or76kE  
   printf("%5.2f %5.2f %5.2f\n",fj1[j],fj2[j],fj3[j]); {AO3o<-h  
} <9"s&G@  
  t0=40.0; KFAB  
  spread=12.0; c[sC 2  
  T=0; ,T|iA/c
 
  nsteps=1; %p/Qz|W  
  ;1NZY.pyc  
while(nsteps>0) x4*8q/G=D  
{ }C/}8<  
  printf("nsteps-->"); 56>Zqtp*  
  scanf("%d",&nsteps); 8ji!FZf  
  printf("%d\n",nsteps); _L_SNjA_  
  L^RyJ;^c  
 for(n=1;n<=nsteps;n++) JGn@)!$+/  
 { &'WgBjP  
  T=T+1; FDbx"%A  
for(j=0;j<JE-1;j++) 9`{Mq9J  
{ ~SV;"e2N.  
  for(i=0;i<IE;i++) =WyAOgy}  
   { O_-.@uo./(  
    curl_e=ez[j]-ez[j+1]; >Eik>dQ a  
    ihx[j]=ihx[j]+curl_e; xO/44D  
    hx[j]=fj3[j]*hx[j]+fj2[j]*0.5*curl_e+fi1*ihx[j]; klT6?'S  
    } VEpIAC4  
} 5f/[HO)  
for(j=0;j<=JE-1;j++) a+A/
l  
{  ;r49H<z  
  for(i=0;i<IE-1;i++) bkmX@+Pe  
  { _]j=[|q 9  
    curl_e=ez[i+1][j]-ez[j]; b7aAP*$  
    ihy[j]=ihy[j]+curl_e; y|.wL=;  
    hy[j]=fi3*hy[j]+fi2*0.5*curl_e+fj1[j]*ihy[j]; /%=#*/E7  
    } AW+4Vm_!l  
  } W;~^3Hz6  
  for(j=1;j<IE;j++) mrFMdpaHl%  
 { 7; TS  
   for(i=1;i<JE;i++) +ywWQ|V  
    { 8Ral%
I:gr  
      dz[j]=gi3*gj3[j]*dz[j]+gi2*gj2[j]*0.5*(hy[j]-hy[i-1][j]-hx[j]+hy[j-1]); 8}xU]N#EV  
   } M[<O]p6  
} {4}Sl^kn*  
pulse=exp(-0.5*(pow((t0-T)/spread,2.0))); 6@H&S  
dz[ic][jc]+=pulse; |g8 ]WFc  
for(j=0;j<JE;j++) ~i))Zc3,g\  
{ '!yyg#  
 for(i=0;i<IE;i++) GS8,mQ8l*l  
  { -|
J?-  
   ez[j]=ga[j]*dz[j]; 0Yfk/}5  
   } tDy1Gh/c  
} d`({
z]W;  
printf("%3f %6.2f\n",T,ez[ic][jc]); f,t[`0 va  
for(j=0;j<JE-1;j++) xS,):R  
{ +8Peh9"  
  ez[0][j]=0.0; \Q^\z   
  ez[IE-1][j]=0.0; h=`rZC  
 } A javV

 
for(i=0;i<IE-1;i++) uBl&|yvxB  
 { ;wW6x  
   ez[0]=0.0; qpp:h_E  
   ez[JE-1]=0.0; *$`r)pV%AK  
} p7.@ez ;  
kZUuRB~om  
} %{zM> le9  
for(j=0;j<JE;j++) #n2GW^x  
{ u:g(x+u4:  
  printf("%2d",j); A<>W^ow  
  for(i=0;i<IE;i++) 7T4rx53  
   { kwWO1=ikz@  
    printf(" % .. .?AtW:<*I  
KC}B\~ +  

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

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

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

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

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

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