[求助]UPML程序中不太理解 K 的设置，希望得到各位的点拨。

我的困惑出自以下这个三维UPML的程序，就是论坛里最常见的那个三维程序，其中由于程序比较长，我节选了我的困惑所在的那段程序。 D4$2'h  
%***********************************************************************   I*H($ a  
%     3-D FDTD code with UPML absorbing boundary conditions C KRb'kW  
%*********************************************************************** 7R$O~R3p  
~9n30j%]s  
%     Program author: Keely J. Willis, Graduate Student L   Tb}op XYK  
%                     UW Computational Electromagnetics Laboratory 1x^Vv;K  
%                           Director: Susan C. Hagness f-Zi!AGh>  
%                     Department of Electrical and Computer Engineering 5Kw?#  
%                     University of Wisconsin-Madison 40}7O<9*  
%                     1415 Engineering Drive m9 4O
Fv#$[  
%                     Madison, WI 53706-1691 OB 2ae"Sd!-2  
Ik-E_U2  
。。。。。。 ]D[\l$(  
。。。。。。             O2"5\@HfE  
%   x-varying material properties }/MmuPp  
delbc=upml*delta; d:=Z<Y?d/  
sigmam=-log(rmax)*(orderbc+1.0)/(2.0*eta*delbc); )Bu#ln"  
sigfactor=sigmam/(delta*(delbc^orderbc)*(orderbc+1.0)); vL8Rg} Jh4  
kmax=1; 8yk4#CZ  
kfactor=(kmax-1.0)/delta/(orderbc+1.0)/delbc^orderbc; for i=1:upml 3? F~H  
OxN[w|2\4  
% Coefficients for field components in the center of the grid cell 6o4Y]C2W{1  
%a']TX  
    x1=(upml-i+1)*delta; @;}vK=6L  
    x2=(upml-i)*delta;  Lr0:yo  
    sigma=sigfactor*(x1^(orderbc+1)-x2^(orderbc+1));   @W-0ybv  
    ki=1+kfactor*(x1^(orderbc+1)-x2^(orderbc+1)); D.gD4g_O/  
    facm=(2*epsr*epsz*ki-sigma*dt); i@mS8%|l  
    facp=(2*epsr*epsz*ki+sigma*dt); 2])e}&i  
    C5ex(i,:,:)=facp; b*1yvkX5  
    C5ex(ie_tot-i+1,:,:)=facp; &
pW2R}  
    C6ex(i,:,:)=facm; slAR<8  
    C6ex(ie_tot-i+1,:,:)=facm; *U +<Hv`C  
    D1hz(i,:,:)=facm/facp; c&C*'c-r  
    D1hz(ie_tot-i+1,:,:)=facm/facp; B[9y<FB+  
    D2hz(i,:,:)=2.0*epsr*epsz*dt/facp; wU|jw(  
    D2hz(ie_tot-i+1,:,:)=2.0*epsr*epsz*dt/facp; IYv.~IQO  
    D3hy(i,:,:)=facm/facp;   -^q;e]+J  
    D3hy(ie_tot-i+1,:,:)=facm/facp; Q hdG(`PY~  
    D4hy(i,:,:)=1.0/facp/mur/muz; _   0X5b32  
    D4hy(ie_tot-i+1,:,:)=1.0/facp/mur/muz; {@hJPK8  
iw%""q(`  
。。。 .. /[E2+g  
DMlr%)@{  

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

没有哪位朋友能给个建议或是解答么？ 顶一下，谢谢！！

如果你比较UPML和CPML。一个差别是K的值。好象K的值作用是伸展坐标（对应SIGMA的值是吸收波的能量。）K不为1的时候，可能回应想到对消逝波的吸收。但K直需要优化。 XUh&an$  
f/Y7@y  
这只是我的个人理解。一起学习哈。

谢谢你的回复，那UPML中，K一直设为1可以么，不按照talflove书中的公式，不影响结果？

taflove的那个公式？

就是第三版 talflove 书里第292页，7.6.2节中的 （7.60a,b)式的K迭代公式。如果Kmax设为1，那K不就一直为1了么？这样做可以么？

恩。你看他在书314页上给的UPML的例子，KMAX设为1。 ]E)gMf   
x_\e&"x  
对于CPML，KMAX ，SIGMA，ALPHA是需要优化的参数。 @cF aYI  
5X-cDY*|  
恩恩。觉得taflove的书在这个问题上的讨论，单薄了些。

嗯，谢谢你的解释。我就先当它为1吧，肯定有他的道理。