[求助]仿真结果和书上结果不一致

教材：时域有限差分方法及编程技巧 }C.M4{a\  
自己写了程序3.4的代码，用MATLAB绘图发现结果和图3.8 .. t QkEJ pj  
Z{RRhJ  

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

/* 1D FDTD simulation with an additive source. */ xMJF1O?3  
R< @o]p  
#include <stdio.h> L'=2Uk#.D  
#include <math.h> )iadu  
u38FY@U$  
#define SIZE 200 Dt?O_Bdv[  
2/S~l;x  
int main() 3?I^D /K^  
{ $+P9@Q$  
double ez[SIZE] = {0.}, hy[SIZE] = {0.}, imp0 = 377.0; K_j$iHqLF  
int qTime, maxTime = 200, mm; yyZH1A  
2Cn^<(F^4I  
char basename[80] = "sim", filename[100]; q+2yp&zF  
int frame = 0; oZ95)'L,  
FILE *snapshot; 7><ne|%  
B%!z7AT  
/* do time stepping */ K1#Y{k5D}  
for (qTime = 0; qTime < maxTime; qTime++) { :Z(?Ct&8  
yZ3/Ia>,  
/* update magnetic field */ Et~b^8$>  
for (mm = 0; mm < SIZE - 1; mm++) FrD.{(/~  
hy[mm] = hy[mm] + (ez[mm + 1] - ez[mm]) / imp0; mMmzi4HL  
2f\;#-  
/* update electric field */ f%9EZ+OP  
for (mm = 1; mm < SIZE; mm++) w#[Ul9=?6  
ez[mm] = ez[mm] + (hy[mm] - hy[mm - 1]) * imp0; Knsb`1"E^6  
HH+R47%*  
/* use additive source at node 50 */ rm2"pfs  
ez[50] += exp(-(qTime - 30.) * (qTime - 30.) / 100.); 6ae  
Jhu<^pjs  
/* write snapshot if time a multiple of 10 */ qQN&uBQ[  
if (qTime % 10 == 0) { j-|YE?AA  
sprintf(filename, "%s.%d", basename, frame++); dq~p]h~,H  
snapshot=fopen(filename, "w"); 'Y3
>+7bI  
for (mm = 0; mm < SIZE; mm++) 'TpW-r:  
fprintf(snapshot, "%g\n", ez[mm]); ;hR!j!3}  
fclose(snapshot); |`T3H5X>  
} "Q9S<O8)  
} /* end of time-stepping */ E 5}T_~-{  
WL-+;h@VQ  
return 0; Ge$cV}  
} mLX1w)=r  

图片:捕获.GIF

已经解决。

程序3.4中时间步数和图3.9并不对应，学习这本书的同学们注意下