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

想问一下，如果一个公式中同时出现一个量，如Ex，它的n时间步和n+1时间步应怎么表达，我在瞬态外推中又一次应用了Ex的递推公式。但在PLjerc方法中，不能这么用，不知道还有什么别的 .. QqM[W/&R  
A
mZXUb  

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

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

在计算Ex(n+1)时，使用另外一个变量保存Ex(n)。 :T-DxP/  
例如： eiK_JPFA-  
Ex_temp(n)=Ex(n) Odwe1q&  
计算Ex(n+1) Z6I|Y5#H  
。。。。。。 oSCaP,P  
。。。。。。 )oIh?-WL  
HHH=cout*(Ex(n+1)-Ex_temp(n))

我有点搞不清楚楼主问这个问题的目的。 oxxuw Dcl  
一个公式里面同时出现第n和n+1步的物理量，不就差不多是个对时间一次导么？而且从公式上看HHH=cout*（Ex（n+1）-Ex（n）），也没什么特殊的地方啊。 .~i|kc]Ue  
FDTD方法应用于电磁场数值计算中，就是通过同时处理对空间，时间的导数换算成离散差分的格式进行模拟的。 b6-N2F1Fs  
如对空间的一次导，有前向差分，后向差分，中间差分等等。 z$Jm1l  
对空间的一次导，也是相应的保存上一个时刻的变量，然后差分就可以了啊。 6q^$}eOt  
t'U=K>7  
如公式里面HHH=cout*（Ex（n+1）-Ex（n）），要同时需要n和n+1时刻的电场才能求得磁场，所以简单的方法你把n时刻的电场用一个数组保存下来，然后用新求得的n+1时刻的电场差分一下不就可以了么？ 6ld4'oM  
lE
k@I"  
` @QZK0Ox  

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

引用第3楼gwzhao于2008-12-22 15:43发表的 : h`6 (Oo|  
我有点搞不清楚楼主问这个问题的目的。 E&8Nh J  
一个公式里面同时出现第n和n+1步的物理量，不就差不多是个对时间一次导么？而且从公式上看HHH=cout*（Ex（n+1）-Ex（n）），也没什么特殊的地方啊。 <q7o"NI6FZ  
FDTD方法应用于电磁场数值计算中，就是通过同时处理对空间，时间的导数换算成离散差分的格式进行模拟的。 it}h8:^<  
如对空间的一次导，有前向差分，后向差分，中间差分等等。 V#Wy` ce  
对空间的一次导，也是相应的保存上一个时刻的变量，然后差分就可以了啊。 <D|&)/#  
.......

'maX  
我的整个程序的运行中，是基于时间步n的一个大循环。但是在E和H的存储中不涉及时间标量。 &GLDoLk6[  

-_ 9k+AV  
上面的公式是我随便举例，实际PLRC中是E的三个时间步间的关系。 e* [wF}))  
`jeATxWv  
如果按照回复所说，岂不是所有步骤的E都要保存，然后读取，太麻烦了吧

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

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

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

电流线性递归卷积方法 |H8C4^1Rq  
球E（n+1）需要E（n）和fair（n） yL<u>S0  
而fai（n）需要E（n）和E（n-1） 9pi{)PDJ  
所以问题变成了 在现有时间步下，如何表示前一时间步

多保存前几步的电场，肯定可以实现的，就要多些内存。 d|$-l:(J  
我记得PLRC后来发展了很多改进的算法，不知道有没有节省内存的，提高速度的。 nZkMyRk  
这个我倒没仔细看。   +D#.u^  

引用第7楼zl1369于2008-12-22 17:24发表的  : <``krPi  
电流线性递归卷积方法 yTU'voE.|  
球E（n+1）需要E（n）和fair（n） (FNX>2Mv  
而fai（n）需要E（n）和E（n-1） BlfW~l'mx  
所以问题变成了 在现有时间步下，如何表示前一时间步

PLRC是片断线性递归卷积，是用于仿真色散介质的。 $^"_Fox]A\  
I~I$/j]e`  
如果在同一个循环中，可以只用一个字E_temp来保存E(n),而不是一个数组。 vE+OL8V  
\%:]o-+"I  
在计算电场E的循环中： "t)|N dZm  
{V9}W<  
E_temp=E(n) X~Vr}  
E(n+1)= .....H(...).... E(n)......fai(n)..... ^5![tTJ  
fai(n)= ...E(n+1).......E_temp...... f:K`MW  
HxH=~B1"P  
OK，公式中系数和H分量被省略了，我想大家应该明白了吧！！！！ ;Cqjg.wkB  
vxC,8Z  
LZ应该把问题介绍仔细一些，才容易获得满意答案！！！！