[求助]关于FDTD的系数

最近在做FDTD有关方面的科研项目。在编程时发现如果直接按照书上的递推公式所示的系数的计算方法来运算的话发现得到的值，比如说，ex的值，加入的入射波时的绝对值为一，但经过几个时间不的递推运算后会变到几千或者更大。虽然是稳定的值但是便得太大。不知道是怎么回事呢？ 0SQrz$y  
特此求助。 1'N<ITb  
另外我在看别人给的程序（葛德彪，闫玉波）的时发现其递推系数直接使用了0.5这个值，不知道这有什么原理呢。希望得到高人的 .. v*OV\h.  
EJ(36h  

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

编写的问题，是计算发散了 %1:caa@_p  
检查公式、代码、下标对应

谢谢指教 TYb$+uY  
我回去好好查查看

su1lv#
 
最近在做FDTD有关方面的科研项目。在编程时发现如果直接按照书上的递推公式所示的系数的计算方法来运算的话发现得到的值，比如说，ex的值，加入的入射波时的绝对值为一，但经过几个时间不的递推运算后会变到几千或者更大。虽然是稳定的值但是便得太大。不知道是怎么回事呢？

4SJ aAeIZ  
检查代码 .Q,"gsY  

{,Q)D$i  
另外我在看别人给的程序（葛德彪，闫玉波）的时发现其递推系数直接使用了0.5这个值，不知道这有什么原理呢。希望得到高人的帮助 P3&s<mh  

FN-/~Su~J  
没发现葛老的书上直接使用了0.5，用0.5做系数的书有一本，是Sullivan的，他把电场做了下转换，然后去dt=delta/2c，系数恰好是0.5

好像在数值稳定性里有讨论，至于那份文献就知道，应该一般详细讨论FDTD的数值稳定性的文献都会讲到~大概是：在二维情况下，要求时间步长小于delta/sqrt(2)/c，三维情况下时间步长小于delta/sqrt(3)/c，其中delta为空间步长~

谢谢，今天刚仔细看了SULIIAN的书，才明白那个0.5的意思。