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

Fd2d_3.2.c. 中设置gi，gj时，xn用的是0.33*pow（xxn，3.0），为什么计算fi，fj时用0.25呢？希望师兄们指点一下小弟，自己比较菜 .. c`<2&ke  
)9H5'Wh#  

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

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

不太清楚

这个以前看过，好象是他们的实际位子是错了半个网格的， @Kr)$F  
建议你在好好看以下网格模型。

%x 方向上的PML 参数设置 U.F65KaKF  
for i = 1:npml+1 bN&da [K  
xnum = npml -i + 1; %从npml 到0 VT7NWTJ,  
xd = npml; ksaC[G;}:  
xxn = xnum / xd; %辅助变量xxn，从1 到0 iP%=Wo.  
xn = 0.33 * xxn^3; %成立方衰减 RGLA}|  
gi2(i) = 1 / ( 1 + xn ); !~l%6Z5  
gi2(IE-i+1) = 1 / ( 1 + xn ); 7dJaWD:&   
gi3(i) = (1 - xn) / (1 + xn); @#tSx  
gi3(IE-i+1) = ( 1 - xn ) / ( 1 + xn ); Xt84Evo  
xxn = ( xnum - 0.5 ) / xd; `@D4?8_  
xn = 0.25 * xxn^3; ^e,RM_.  
fil(i) = xn; %nh'F6bNgv  
fil(IE+1-i) = xn; 5w1[KO#K|  
fi2(i) = 1 / ( 1 + xn ); U
G_0Y8$  
fi2(IE+1-i) = 1 / ( 1 + xn ); 9p4U\hx  
fi3(i) = ( 1-xn ) / ( 1 + xn ); lpUtNy  
fi3(IE+1-i) = ( 1 - xn ) / ( 1 + xn ); fNJ;{&#  
end; 8!SiTOzR?  
for j = 1:npml+1 K-u/q6ufK  
xnum = npml - j + 1; -JENY|6  
xd = npml; k#) .E X  
xxn = xnum / xd; (FYJ^o  
xn = 0.33 * xxn^3; @
GtZK  
gj2(j) = 1 / ( 1 + xn ); WCu%@hh=h  
gj2(JE+1-j) = 1 / ( 1 + xn ); *~uuCLv_  
gj3(j) = ( 1 - xn ) / ( 1 + xn ); rfi`Bp  
gj3(JE-j+1) = ( 1 - xn ) / ( 1 + xn ); A%2}?Ds  
xxn = ( xnum - 0.5 ) / xd; |]cDz  
xn = 0.25 * xxn^3; b]#d04]  
fj1(j) = xn; ~Wm}M  
fj1(JE+1-j) = xn; \AzcW;03g[  
fj2(j) = 1 / ( 1 + xn );
rtx]dc1m  
fj2(JE+1-j) = 1 / ( 1 + xn );  vB*oI~<  
fj3(j) = ( 1 - xn ) / ( 1 + xn ); c7IR06E  
fj3(JE+1-j) = ( 1 - xn ) / ( 1 + xn ); " G6jUTt  
end; I}IW!K  
就是上面的部分……