[求助]新手求助sulivan书中一维ABC吸收边界问题

我是新手，刚刚接触FDTD算法， ~XyW&@  
现在正从从sulivan书入手看， pl
5P2&k  
看到ABC吸收边界处理的时候，很是迷惑啊 R)M_|ca  
-+7uy.@cS  
/*Absorbing Boundary Conditions*/ d"5oD@JG:  
             ex[0]      =       ex_low_m2; ]W Zq^'q.  
             ex_low_m2  =       ex_low_m1; e1>aTu@  
             ex_low_m1  =       ex[1]; Z.W66\8~}^  
7_.11$E=H  
             ex[KE-1]   =       ex_high_m2; sf OHl  
             ex_high_m2 =       ex_high_m1; Aub]IO~  
             ex_high_m1 =       ex[KE-2]; i^_#%L  
Xg=x7\V  
low和high 分列三个公式，是为了表示空间步长的两个时间步， r'dr9"-{  
low 和 high  ，我自己 .. M t*6}Cl  
`((Y
c]:7  

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

左右边界条件分开处理

用户被禁言,该主题自动屏蔽!

左边界是用low表示 b-@\R\T  
右边界是用high表示？ P20|RvE  
这个程序，并没有具体赋ex_low_m2  和 ex_low_m1值， ?@R")$  
那我想程序使用ex_low_m2  和 ex_low_m1这两个变量， DlyMJ#a  
仅仅是为了表示取了两个网格，我这样理解对不对？

两临时变量，理解了好像不够深刻。 J(h3]J/Yw  
uhn%lV]  
你把dt 改一下。 cfoYnM  
1）dt = dx / c FMu!z  
2) dt = dx / (2 * c) o;J;*~g  
3) dt = dx / (3 * c) PDw{R]V+  
K5 vNhA  
在这些情况下，如何实现吸收边界。

我也不太懂到底ex[0]      =       ex_low_m2;              ex_low_m2  =       ex_low_m1; URw!7bTz  
             ex_low_m1  =       ex[1]; \|{*arS  
是干啥用的，不懂，运行后ex(2)延迟三步有赋给ex(1)为什么

有没有人能够详细讲解一下呢，还是终究无法理解！关于两个临时变量的作用，怎么实现相差两个时间步的？

必须满dx=1/2*c*dt

这个问题只需要把每一次循环后ex(0),ex_low_m1,ex_low_m2,ex(1)的值列出来就清楚了，书上的程序里面没有设定初始值，需要将ex_low_m1, ex_low_m2都设定成0. 然后开始循环，下面是每一次循环后变量的值. 假设ex(1)在每一步的值是 1 2 3 4 5 kX ~-g  
n(time step)     ex(0)     ex_low_m2      ex_low_m1    ex(1) [HC8-N^.}  
      1                      0               0                        1                    1 6Tm Rc  
      2                      0               1                        2                    2 o2[vM$]  
      3                      1               2                        3                    3 &CB.*\0  
      4                      2               3                        4                    4 [dSDg2]  
      5                      3               4                        5                    5 icKg7-$N  
      6                      4               5                        ..                    .. PfwI@%2  
      7                      5               ..                        ..                    .. x>1iIpBv^  
所以ex(1)的值在变量从传了两次才到ex(0)，进过了2*dt. h|(ZXCH  
-13P 2<i+  
oTp
lxF1  
就像gwzhao说的那样，把dt改成dt = dx / (3 * c) 改怎么实现,只需要再引入一个变量ex_low_m3就可以了，这样解释应该比较清楚了吧

我想问一下··ABC代表了什么意思啊···