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

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

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

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

不太清楚

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

%x 方向上的PML 参数设置 )anprhc  
for i = 1:npml+1 (jj=CLe  
xnum = npml -i + 1; %从npml 到0 zx:Qz  
xd = npml; #-T xhwYs  
xxn = xnum / xd; %辅助变量xxn，从1 到0 jz CA2N%  
xn = 0.33 * xxn^3; %成立方衰减 )fPN6x/e  
gi2(i) = 1 / ( 1 + xn ); SUv(MA&  
gi2(IE-i+1) = 1 / ( 1 + xn ); #0OW0:Q  
gi3(i) = (1 - xn) / (1 + xn); qg_M9xJ  
gi3(IE-i+1) = ( 1 - xn ) / ( 1 + xn ); R@ksYC3 F  
xxn = ( xnum - 0.5 ) / xd; 4:1URhE  
xn = 0.25 * xxn^3; ]7
0V  
fil(i) = xn;  ! @EZ  
fil(IE+1-i) = xn; %x(||cq  
fi2(i) = 1 / ( 1 + xn ); f]c{,LFvZ  
fi2(IE+1-i) = 1 / ( 1 + xn ); IRWVoCc9/\  
fi3(i) = ( 1-xn ) / ( 1 + xn ); [PXv8K%]p  
fi3(IE+1-i) = ( 1 - xn ) / ( 1 + xn ); !vd(WKq  
end; =3K}]3f  
for j = 1:npml+1 IIj :\?r  
xnum = npml - j + 1; >@Vap  
xd = npml; h=<x%sie  
xxn = xnum / xd; jtS-nQ|  
xn = 0.33 * xxn^3; MaMP7O|W  
gj2(j) = 1 / ( 1 + xn ); l1_hD,4  
gj2(JE+1-j) = 1 / ( 1 + xn ); be^+X[
 
gj3(j) = ( 1 - xn ) / ( 1 + xn ); bU$4"_eA B  
gj3(JE-j+1) = ( 1 - xn ) / ( 1 + xn ); V)|]w[(Y  
xxn = ( xnum - 0.5 ) / xd; k*xMe-  
xn = 0.25 * xxn^3; c, IAz  
fj1(j) = xn; {tE9m@[AF  
fj1(JE+1-j) = xn; 2c)Ez?  
fj2(j) = 1 / ( 1 + xn ); JDbRv'F:(  
fj2(JE+1-j) = 1 / ( 1 + xn ); d@d\9*mn  
fj3(j) = ( 1 - xn ) / ( 1 + xn );
-jXO9Q  
fj3(JE+1-j) = ( 1 - xn ) / ( 1 + xn ); BqoGHg4iq  
end; Mk-zeq<2z  
就是上面的部分……