[讨论]公式中时间步问题（瞬态和等离子体中出现）

想问一下，如果一个公式中同时出现一个量，如Ex，它的n时间步和n+1时间步应怎么表达，我在瞬态外推中又一次应用了Ex的递推公式。但在PLjerc方法中，不能这么用，不知道还有什么别的 .. l^0 <a<P  
4dfe5\  

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

补充下，是n时间步和n+1时间步同时出现在一个公式里，如HHH=cout*（Ex（n+1）-Ex（n）），编程中该如何体现

在计算Ex(n+1)时，使用另外一个变量保存Ex(n)。 p6ZKyi  
例如： 6w@l#p  
Ex_temp(n)=Ex(n) 9h9Y:i*Gh5  
计算Ex(n+1) u|w[b9^r  
。。。。。。 12rr:(#%s  
。。。。。。 P|64wq{B8  
HHH=cout*(Ex(n+1)-Ex_temp(n))

我有点搞不清楚楼主问这个问题的目的。 lE'2\kxI?  
一个公式里面同时出现第n和n+1步的物理量，不就差不多是个对时间一次导么？而且从公式上看HHH=cout*（Ex（n+1）-Ex（n）），也没什么特殊的地方啊。 =g UOHH  
FDTD方法应用于电磁场数值计算中，就是通过同时处理对空间，时间的导数换算成离散差分的格式进行模拟的。 :tdN#m6&  
如对空间的一次导，有前向差分，后向差分，中间差分等等。 W2&o'(P\  
对空间的一次导，也是相应的保存上一个时刻的变量，然后差分就可以了啊。 Vo"Wr>F  
kZ>_m&g  
如公式里面HHH=cout*（Ex（n+1）-Ex（n）），要同时需要n和n+1时刻的电场才能求得磁场，所以简单的方法你把n时刻的电场用一个数组保存下来，然后用新求得的n+1时刻的电场差分一下不就可以了么？ B/0Xqyu  
whxTCI
V  
of659~EIW  

引用第1楼zl1369于2008-12-22 11:20发表的  : [0qe ?aI  
补充下，是n时间步和n+1时间步同时出现在一个公式里，如HHH=cout*（Ex（n+1）-Ex（n）），编程中该如何体现

引用第3楼gwzhao于2008-12-22 15:43发表的 : #
u}%r{T  
我有点搞不清楚楼主问这个问题的目的。 J
s
l2RdI  
一个公式里面同时出现第n和n+1步的物理量，不就差不多是个对时间一次导么？而且从公式上看HHH=cout*（Ex（n+1）-Ex（n）），也没什么特殊的地方啊。 y];-D>jk  
FDTD方法应用于电磁场数值计算中，就是通过同时处理对空间，时间的导数换算成离散差分的格式进行模拟的。 I`zd:o]  
如对空间的一次导，有前向差分，后向差分，中间差分等等。 W1 k]P.  
对空间的一次导，也是相应的保存上一个时刻的变量，然后差分就可以了啊。 0:f]&Ng  
.......

h2 KI  
我的整个程序的运行中，是基于时间步n的一个大循环。但是在E和H的存储中不涉及时间标量。 u\V^g  
=:aJZ[UU<2  
上面的公式是我随便举例，实际PLRC中是E的三个时间步间的关系。 v\GVy[Qyv  
|>yWkq  
如果按照回复所说，岂不是所有步骤的E都要保存，然后读取，太麻烦了吧

一次导要保留n时刻和n+1时刻的，二次导就要保留n-1,n和n+1时刻的。递推公式里面包含的东西总得有啊，除非你换种推导公式，把几个时间步长消掉一些。
V!W.P  

引用第4楼zl1369于2008-12-22 16:37发表的  :
我的整个程序的运行中，是基于时间步n的一个大循环。但是在E和H的存储中不涉及时间标量。
上面的公式是我随便举例，实际PLRC中是E的三个时间步间的关系。
.......

PLRC是指的什么？等离子体色散关系？

电流线性递归卷积方法 /:{%X(8  
球E（n+1）需要E（n）和fair（n） swKkY`g  
而fai（n）需要E（n）和E（n-1） 18X@0e  
所以问题变成了 在现有时间步下，如何表示前一时间步

多保存前几步的电场，肯定可以实现的，就要多些内存。 G=y~)B}  
我记得PLRC后来发展了很多改进的算法，不知道有没有节省内存的，提高速度的。 b}TvQ+W]2  
这个我倒没仔细看。   v4e4,Nt  

引用第7楼zl1369于2008-12-22 17:24发表的  : 3$hIc)  
电流线性递归卷积方法 cs6oD!h  
球E（n+1）需要E（n）和fair（n） B}[CU='P*  
而fai（n）需要E（n）和E（n-1） hh4R  
所以问题变成了 在现有时间步下，如何表示前一时间步

PLRC是片断线性递归卷积，是用于仿真色散介质的。 M(L6PyEa!Y  
-i?gYF!G  
如果在同一个循环中，可以只用一个字E_temp来保存E(n),而不是一个数组。 %16Lo<DPm  
_<RTes  
在计算电场E的循环中： S3M!"l  
c2NB@T9'v  
E_temp=E(n) bN-!&Td  
E(n+1)= .....H(...).... E(n)......fai(n)..... d<d3j9u(#  
fai(n)= ...E(n+1).......E_temp...... eP"B3Jw  
p-IJ':W  
OK，公式中系数和H分量被省略了，我想大家应该明白了吧！！！！ ^mn!;nu  
-TVwoK  
LZ应该把问题介绍仔细一些，才容易获得满意答案！！！！