[求助]新手关于sullivan书PML程序部分的疑问

Fd2d_3.2.c. 中设置gi，gj时，xn用的是0.33*pow（xxn，3.0），为什么计算fi，fj时用0.25呢？希望师兄们指点一下小弟，自己比较菜 .. !PA:#]J  
rRe5Q  

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

有书有code的，大家看一下。没书没code的我就闪了。

不太清楚

这个以前看过，好象是他们的实际位子是错了半个网格的， Z;QbqMj  
建议你在好好看以下网格模型。

%x 方向上的PML 参数设置 ^MGgFS]G  
for i = 1:npml+1 qqSf17sW  
xnum = npml -i + 1; %从npml 到0 gI qYIt  
xd = npml; :YU_ \EV  
xxn = xnum / xd; %辅助变量xxn，从1 到0 ?j{LE-(  
xn = 0.33 * xxn^3; %成立方衰减 jLf.qf8qm  
gi2(i) = 1 / ( 1 + xn ); spV/+jy{  
gi2(IE-i+1) = 1 / ( 1 + xn ); G:QaWqUb  
gi3(i) = (1 - xn) / (1 + xn); Wej8YF@  
gi3(IE-i+1) = ( 1 - xn ) / ( 1 + xn ); /,z4tf  
xxn = ( xnum - 0.5 ) / xd; hSq3LoHV  
xn = 0.25 * xxn^3; "3A.x1uQ  
fil(i) = xn; @O/-~,E68  
fil(IE+1-i) = xn; geL)v7t+#  
fi2(i) = 1 / ( 1 + xn ); \4`~J@5Y  
fi2(IE+1-i) = 1 / ( 1 + xn ); ([='LyH];z  
fi3(i) = ( 1-xn ) / ( 1 + xn ); V&E)4KBOs  
fi3(IE+1-i) = ( 1 - xn ) / ( 1 + xn ); YAqv:
  
end; k"V| f&  
for j = 1:npml+1 JZP>`c21y]  
xnum = npml - j + 1; {
mK=Vig  
xd = npml; ~1
Q$FgLk  
xxn = xnum / xd; ((TiBCF4  
xn = 0.33 * xxn^3; HgP9evz,0  
gj2(j) = 1 / ( 1 + xn ); c@%:aiEl  
gj2(JE+1-j) = 1 / ( 1 + xn ); $UAmUQg)}_  
gj3(j) = ( 1 - xn ) / ( 1 + xn ); ,25Qhz]  
gj3(JE-j+1) = ( 1 - xn ) / ( 1 + xn ); cbD&tsF  
xxn = ( xnum - 0.5 ) / xd; #=@( m.k:s  
xn = 0.25 * xxn^3; {i)k#`  
fj1(j) = xn; cWSiJr):r  
fj1(JE+1-j) = xn; `g&<7~\=A  
fj2(j) = 1 / ( 1 + xn ); J)iy6{0"  
fj2(JE+1-j) = 1 / ( 1 + xn ); ^9 gFW $]  
fj3(j) = ( 1 - xn ) / ( 1 + xn ); o })k@-oL  
fj3(JE+1-j) = ( 1 - xn ) / ( 1 + xn ); Lf|5miO  
end;
QXnL(z  
就是上面的部分……