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

教材：时域有限差分方法及编程技巧 F.?01,J=1  
自己写了程序3.4的代码，用MATLAB绘图发现结果和图3.8 .. e[{mVhg4E  
'w.}2(  

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

/* 1D FDTD simulation with an additive source. */ 5 1dSFr<#  
Ri)uq\E/#  
#include <stdio.h> S3Y2O x  
#include <math.h> 4W;S=#1  
O{]9hm(tN  
#define SIZE 200 LFf`K)q  
qP1FJ89
H  
int main() tR)H~l7q  
{ dRaOGm)  
double ez[SIZE] = {0.}, hy[SIZE] = {0.}, imp0 = 377.0; L/Vx~r`P  
int qTime, maxTime = 200, mm; nidr\oFUIn  
X`JVR"=4  
char basename[80] = "sim", filename[100]; )9L/sKz  
int frame = 0; _#C()Ro*P  
FILE *snapshot; }6]0hWsN[  
F9flSeN  
/* do time stepping */ }]6f+  
for (qTime = 0; qTime < maxTime; qTime++) { j0mM>X HB  
p&Ed\aQ%z;  
/* update magnetic field */ p|[B =.c{  
for (mm = 0; mm < SIZE - 1; mm++) m3.sVI0I  
hy[mm] = hy[mm] + (ez[mm + 1] - ez[mm]) / imp0; Q5a)}6-5  
9P WY52!  
/* update electric field */ ig+4S[L~n  
for (mm = 1; mm < SIZE; mm++) *%gF2@=r8F  
ez[mm] = ez[mm] + (hy[mm] - hy[mm - 1]) * imp0; {whvTN1#dh  
|E3X  
/* use additive source at node 50 */ ~*.-  
ez[50] += exp(-(qTime - 30.) * (qTime - 30.) / 100.); M.[
rLJZ4  
:iUF7P1I  
/* write snapshot if time a multiple of 10 */  P_Hv%g  
if (qTime % 10 == 0) { 09h.1/  
sprintf(filename, "%s.%d", basename, frame++); xDLG=A%]z  
snapshot=fopen(filename, "w"); 5@!st  
for (mm = 0; mm < SIZE; mm++) ,wAz^cK|  
fprintf(snapshot, "%g\n", ez[mm]); Log|%P\  
fclose(snapshot); o{WyQ&2N  
} PXk?aJ  
} /* end of time-stepping */ .iwZ*b{  
bt/ =Kq#  
return 0; W+=o&V  
} <AoXEuD  

图片:捕获.GIF

已经解决。

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