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

葛老师书后(第二版)的第二个程序中遇到一点问题请教大家 /c):}PJ^#7  
1.real EAC(4,2,0:1),括号里面的数字是什么意思? 4是不是代表有4个边连接边界? Zw wqSyuGf  
#97h6m?  
2,   real EBin(4) !入射波吸收边界条件临时变量 ioPUUUb)  
    Ein(IncidentStart)=EBin(4); 'l,V*5L  
    EBin(4)=EBin(3); r,_?F7  
    EBin(3)=Ein(IncidentStart+1);
b,8{ X<  
    Ein(IncidentEnd)=EBin(2); ] }f9JNf$  
    EBin(2)=EBin(1); _=|vgc  
    EBin(1)=Ein(IncidentEnd-1); 4Vq%N  
谁能给解释一下这段程序,另外括号里面的数字代表什么意思?? uYrfm:4S  
要是在三维里面的话,是不是入射波吸收边界也是这么设定(看了sullivan的程序发现他二维三维都是一样的程序设定入射波吸收边界)  d\|!Hg,  
3. T1=float(I)*cp+(float(Jtmin)-0.5)*sp (底部连接边界处HX .. !nAX$i~  
.jMq  

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

我自己顶一下,期待高人解答

hoping,waiting高人的解答

f/ ?_  
这个我觉得他程序里面都已经注明了啊。 $? m9")  
比如第二个，是存储临时变量的。 b*;Si7-  
括号里面的数字？ y6ECdVF  
A(2)就表示数组A里面的第二个数据，这个，这个。。。很难么，：（

引用第3楼gwzhao于2009-04-06 11:17发表的 : 6ezS{Q  
[表情] C4
hx@abA  
这个我觉得他程序里面都已经注明了啊。 zN{K5<7o  
比如第二个，是存储临时变量的。 )(tM/r4`c&  
括号里面的数字？ ynw5-aS3  
A(2)就表示数组A里面的第二个数据，这个，这个。。。很难么，：（

%H]lGN)  
Pb8^ b  
 ,这第二问题问的确实......也不知道自己但是脑子怎么想的.......我检讨! (tCUlX2  
第二段程序,问题我想明白了,一开始自己想偏了,因为一维入射波的dt=dx/c,而二维TM波的dt=dx/2c,所以传播同样的距离一维入射波需要2个时间步,传播同样的距离一维与二维相比有一个时间延迟,为了解决这个问题所以要这么设置,不知道这么理解对不对

如果你知道为什么一维吸收边界 lc_E!"1  
E(n) = E(n-1) | dt = dz/(c)， W5^.-B,(K  
E(n) = E(n-2) | dt = dx/(2*c) VdP`a(Yd;  
就更好了。 /t6X(*xoy  

引用第4楼wq_463于2009-04-06 12:09发表的  :
[表情] ,这第二问题问的确实......也不知道自己但是脑子怎么想的.......我检讨! {QbvR*gv  
第二段程序,问题我想明白了,一开始自己想偏了,因为一维入射波的dt=dx/c,而二维TM波的dt=dx/2c,所以传播同样的距离一维入射波需要2个时间步,传播同样的距离一维与二维相比有一个时间延迟,为了解决这个问题所以要这么设置,不知道这么理解对不对

引用第5楼gwzhao于2009-04-06 16:25发表的 : <!pQ  
如果你知道为什么一维吸收边界 cst}Ibfi  
E(n) = E(n-1) | dt = dz/(c)， 5j0{p$'9  
E(n) = E(n-2) | dt = dx/(2*c) z4@k$ L8  
就更好了。 cECi')  

=`Y.=RL+'n  
根据courant稳定性要求dT<=dx/sqrt(n)/c0,二维和三维我们取sqrt(n)=2,一维取sqrt＝1。 83cW=?UgA  
kkvtB<<Y  
E(n) = E(n-1) | dt = dz/(c)  这个我明白 rdnRBFt  
E(n) = E(n-2) | dt = dx/(2*c)  是不是可以理解为当一维入射波引入到二维的TM（TE）波时是这么设置呢？ w[_x(Ojq;  
请师兄指正

有点不对。 &h$|j  
就是一维时，我也可以设置dt = dx / (3 * c)或者 dt = dx / ( n * c)的啊。 MM{_Ur7Q  

引用第6楼wq_463于2009-04-06 20:54发表的  :
根据courant稳定性要求dT<=dx/sqrt(n)/c0,二维和三维我们取sqrt(n)=2,一维取sqrt＝1。
E(n) = E(n-1) | dt = dz/(c)  这个我明白 ]*%+H|l  
E(n) = E(n-2) | dt = dx/(2*c)  是不是可以理解为当一维入射波引入到二维的TM（TE）波时是这么设置呢？ h3u1K>R)  
.......

是这样的, 一维dt = dx/c0 W^g[L:s  
                 二维dt=dx/sqrt(2)/c0 OCyG_DLT$5  
                 三维dt=dx/sqrt(3)/c0 ,[48Mspp  
这是走过一个网格（无论几维）走对角线（最长距离）时的所需要的时间。因此，收敛条件就应该是 ^9"KTZc-*  
dt<=dx/sqrt(n)/c0，n是维数。那么即使当n=3时，取dt=dx/2/c0也是可以符合收敛条件的。可以看Talove的论文。 j!#OG  
我也不知道说得清楚不。

引用第8楼editwolf于2009-04-07 17:12发表的 : A[/I#Im7  
是这样的, 一维dt = dx/c0 A|x:UQlu  
                 二维dt=dx/sqrt(2)/c0 18!VO4u\I  
                 三维dt=dx/sqrt(3)/c0 9q4_j  
这是走过一个网格（无论几维）走对角线（最长距离）时的所需要的时间。因此，收敛条件就应该是 2wh#$zGy  
dt<=dx/sqrt(n)/c0，n是维数。那么即使当n=3时，取dt=dx/2/c0也是可以符合收敛条件的。可以看Talove的论文。 H30OUrD  
.......

-6EK#!+  

['Z{@9  
你跟我一样还没搞明白gwzhao师兄的问题,师兄的问题是为什么一维的时候也可以设置dt=dx/2c或者dt=dx/3c!! IS_Su;w>4  
我还没想明白,正在思考中