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

教材：时域有限差分方法及编程技巧 _$_CR\$  
自己写了程序3.4的代码，用MATLAB绘图发现结果和图3.8 .. SK}g(X7IWH  
%c2i.E/G  

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

/* 1D FDTD simulation with an additive source. */ iiB$<b.((I  
8[IifF1M=&  
#include <stdio.h> &"n9,$  
#include <math.h> q8]k]:r  
@6`@.iZ  
#define SIZE 200 P*sCrGO%  
a8c]B/  
int main() f#b[KB^Z,2  
{ d_B5@9e#  
double ez[SIZE] = {0.}, hy[SIZE] = {0.}, imp0 = 377.0; K#m o+n5-;  
int qTime, maxTime = 200, mm; KphEw[4/  
El}z^e  
char basename[80] = "sim", filename[100]; `d#_66TLr  
int frame = 0; th|Q NG  
FILE *snapshot; `"4EE}eQc
 
DO{4n1-U  
/* do time stepping */ .?:~s8kB  
for (qTime = 0; qTime < maxTime; qTime++) { QV1%Zou  
[wYQP6Cyy  
/* update magnetic field */ >@
iV!!  
for (mm = 0; mm < SIZE - 1; mm++) o 4wKu  
hy[mm] = hy[mm] + (ez[mm + 1] - ez[mm]) / imp0; 6nhB1Aei  
:gv#_[k  
/* update electric field */ >X)G`N@!  
for (mm = 1; mm < SIZE; mm++) wyM3|%RZ  
ez[mm] = ez[mm] + (hy[mm] - hy[mm - 1]) * imp0; \zyGJyy.  
b=EZtk6>  
/* use additive source at node 50 */ `rQl{$9IC  
ez[50] += exp(-(qTime - 30.) * (qTime - 30.) / 100.); D~>P/b)v{j  
FgaBwd^W  
/* write snapshot if time a multiple of 10 */ mJT m/C  
if (qTime % 10 == 0) { 1t'\!  
sprintf(filename, "%s.%d", basename, frame++); ~=*_I4,+r  
snapshot=fopen(filename, "w"); YdaJ&  
for (mm = 0; mm < SIZE; mm++) =v9;HPiO  
fprintf(snapshot, "%g\n", ez[mm]); Y [S^&pF  
fclose(snapshot); 5Px_vtqP  
} Pw{+7b$  
} /* end of time-stepping */ H;E{Fnarv  
'5\1uB PKW  
return 0; {eV_+@dT  
} VvKH]>*  

图片:捕获.GIF

已经解决。

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