[求助]关于葛老师书后的程序的问题

恩，我的意思就是这个，稳定性要求dt有个最大值，但没有规定最小值。 "x vizvR  
e#)NYcr6  
以前我做数值模型，包括等离子体FDTD的时候，经常出现发散，为了找到原因，我经常把dt 设置的比较小一些，：） 2M*i'K;;)P  
"BVp37m;?  
包括你要加源的问题，你先试想一下，如果要加个平面波，自己加的话，想怎么加才能保证传播出去的是平面波，然后再看一下书上是怎么加的，回味一下，然后就算在编程时出现问题，也能及时找到源头了。 h*R w^5,c  
1F_$[iIX]  
数值模拟解决的一些简单问题，比如平面波的传播，都是有数值解的，所以每一步，你都可以预先知道它大概是多少。 -p?&vQDo`  

引用第9楼wq_463于2009-04-07 19:39发表的  :
你跟我一样还没搞明白gwzhao师兄的问题,师兄的问题是为什么一维的时候也可以设置dt=dx/2c或者dt=dx/3c!! Y1ks'=c>  
我还没想明白,正在思考中

引用第10楼gwzhao于2009-04-07 21:42发表的 : v2E<~/|  
恩，我的意思就是这个，稳定性要求dt有个最大值，但没有规定最小值。
以前我做数值模型，包括等离子体FDTD的时候，经常出现发散，为了找到原因，我经常把dt 设置的比较小一些，：）
包括你要加源的问题，你先试想一下，如果要加个平面波，自己加的话，想怎么加才能保证传播出去的是平面波，然后再看一下书上是怎么加的，回味一下，然后就算在编程时出现问题，也能及时找到源头了。 mEbI\!}H0  
.......

^?Mp(o  
谢谢师兄指点，借这个帖子我还有些问题向师兄请教 Uns%6o  
上次师兄给的2维matlab平面波程序，里面有些地方我不是很明白 1C8xJ6F  
Isbias=0; B,cFvS  
Isbias=2*max(Itmax,Jtmax); bY2R/FNL=  
Isbias=Isource; 3i7E
F.  
T1=I*cp+Jtmin*sp+Isbias; w;gk=<_  
关于这个Isbias这个参数的意义和为什么要加在T1后面，请师兄给我解释一下，谢谢

这个程序问什么要延时两个Δt呢 还是不明白 请斑竹师兄稍微详细说明一下 兄弟愚钝 求您不吝赐教