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

我是新手，刚刚接触FDTD算法， oi m7=I0  
现在正从从sulivan书入手看， \=D+7'3  
看到ABC吸收边界处理的时候，很是迷惑啊 Fv7]1EO.  
[[HCP8Wk  
/*Absorbing Boundary Conditions*/ Zb=;\l*&  
             ex[0]      =       ex_low_m2; X!:J1'FE  
             ex_low_m2  =       ex_low_m1; :pM)I5MN[  
             ex_low_m1  =       ex[1]; \6
1H(,  
sj4\lpZ3h  
             ex[KE-1]   =       ex_high_m2; 6fo\z2  
             ex_high_m2 =       ex_high_m1; mKjTJzS  
             ex_high_m1 =       ex[KE-2]; ,ss"
s3  
1d~d1Rd  
low和high 分列三个公式，是为了表示空间步长的两个时间步， w[F})u]E  
low 和 high  ，我自己 .. >yr;Y4y7K  
-<g[P_#  

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

左右边界条件分开处理

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

左边界是用low表示 1gt 7My  
右边界是用high表示？ m#O; 1/P  
这个程序，并没有具体赋ex_low_m2  和 ex_low_m1值， kSCpr0c  
那我想程序使用ex_low_m2  和 ex_low_m1这两个变量，
T$<'ZC  
仅仅是为了表示取了两个网格，我这样理解对不对？

两临时变量，理解了好像不够深刻。 KMl3`+i  
9 *v14c%  
你把dt 改一下。 dg+"G|nr  
1）dt = dx / c W
>b\O">  
2) dt = dx / (2 * c) 5RY-.c4}  
3) dt = dx / (3 * c) {`>pigo  
p-!/p#  
在这些情况下，如何实现吸收边界。

我也不太懂到底ex[0]      =       ex_low_m2;              ex_low_m2  =       ex_low_m1; xd
4~[n\hm  
             ex_low_m1  =       ex[1]; 9!dG Xq  
是干啥用的，不懂，运行后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 f1d<xGx  
n(time step)     ex(0)     ex_low_m2      ex_low_m1    ex(1) CKDg3p';  
      1                      0               0                        1                    1 C{6m?6  
      2                      0               1                        2                    2 gX* &RsF  
      3                      1               2                        3                    3 W5&KmA  
      4                      2               3                        4                    4 V{rQ@7SE  
      5                      3               4                        5                    5 5)w;0{X!P  
      6                      4               5                        ..                    ..
?X7nM)  
      7                      5               ..                        ..                    ..  ~"h V-3U  
所以ex(1)的值在变量从传了两次才到ex(0)，进过了2*dt. m# ^).+  
zK*i:(>B  
oK&G  
就像gwzhao说的那样，把dt改成dt = dx / (3 * c) 改怎么实现,只需要再引入一个变量ex_low_m3就可以了，这样解释应该比较清楚了吧

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