[求助]FDTD 计算出来的电场是发散的

FDTD程序计算出来的电场是发散的，dt=dx/2c,可能是什么原因呢？会不会是激励源的缘故。我把激励源改小的话可不可以消除色散呢？ X;i~<Tq  
激励源设置如下： gu0j.XS^  
E=sin(omega0*(n-1)*dt)*exp(-((n-1981)/660)**2). \9cG36  
我的意思是我把1981改成更大的值，660换成更小的值能不能消除色散呢。色散 .. eM$
sv9?  
JBXrFC;  

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

发散和色散是两个不同的概念，发散可能是你的算法或程序有问题，或是你的计算时间步还不够，还没达到收敛。而色散是FDTD 方法本身所固有的特性，可以通过减小空间步来减小数值色散，色散不能消除，空间步必须满足足色散条件，时间步必须满足稳定性条件。

先看看FDTD方面的书，将基本概念搞清楚。

如果你计算出来的电场是发散的，那说明代码有问题，肯定哪里出错了。 j[`j9mM8  
)c8rz
[i  
至于色散，一般在你研究色散介质的时候才需要详细了解，而一般的数值色散，在精度足够的前提下，你可以先不用考虑这个问题。 :r{<zd>;  
D{GfLib"U  

激励源：E=sin(omega0*(n-1)*dt)*exp(-((n-1981)/660)^2). [$mHv,~  
式中，1981为高斯型脉冲的峰值所对应的位置n，660与高斯脉冲的束腰宽度相关[sqrt(2)联系] ?4+9fE<Q  
调整这两个数值并不能消除色散

dt=dx/2c,?? yg`j-9[8  
dt=dx/(2c)

以我的经验也是代码或者你公式推错了，特别是那些带符号的系数最容易错了，错一个就会发散。