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

我的困惑出自以下这个三维UPML的程序，就是论坛里最常见的那个三维程序，其中由于程序比较长，我节选了我的困惑所在的那段程序。 3z: rUhA  
%***********************************************************************   !`,6E`Y#  
%     3-D FDTD code with UPML absorbing boundary conditions C ,{itnKJC  
%*********************************************************************** ?TE#4}p|  
9`.b   
%     Program author: Keely J. Willis, Graduate Student L   ix^:qw;  
%                     UW Computational Electromagnetics Laboratory (Tn*;Xjq  
%                           Director: Susan C. Hagness )rhKWg  
%                     Department of Electrical and Computer Engineering ?`\<t$M
 
%                     University of Wisconsin-Madison 4'+/R%jk"  
%                     1415 Engineering Drive m9 B/E1nBobC  
%                     Madison, WI 53706-1691 OB 3w{4G<I  
5#d"]7  
。。。。。。 rf;R"Uc  
。。。。。。             uVDB;6  
%   x-varying material properties (5-"5<-@R  
delbc=upml*delta; l/LRr.x  
sigmam=-log(rmax)*(orderbc+1.0)/(2.0*eta*delbc); 2K, 1wqf'  
sigfactor=sigmam/(delta*(delbc^orderbc)*(orderbc+1.0)); }
?*:uf  
kmax=1; FlVGi3  
kfactor=(kmax-1.0)/delta/(orderbc+1.0)/delbc^orderbc; for i=1:upml z^~uq:  
@[zPN[z.  
% Coefficients for field components in the center of the grid cell H
) q9.Jg  
ReKnvF~  
    x1=(upml-i+1)*delta; z2!NBOv  
    x2=(upml-i)*delta; 3KB)\nF#%  
    sigma=sigfactor*(x1^(orderbc+1)-x2^(orderbc+1));   R9S7p)B  
    ki=1+kfactor*(x1^(orderbc+1)-x2^(orderbc+1)); 9]|[z{v'>l  
    facm=(2*epsr*epsz*ki-sigma*dt); PKdM-R'Z  
    facp=(2*epsr*epsz*ki+sigma*dt); Js^ADUy  
    C5ex(i,:,:)=facp; OD>-^W t;%  
    C5ex(ie_tot-i+1,:,:)=facp; ]t0?,q.$7  
    C6ex(i,:,:)=facm; 1ErH \!  
    C6ex(ie_tot-i+1,:,:)=facm; F8b*Mt}p  
    D1hz(i,:,:)=facm/facp; "TtK!>!.  
    D1hz(ie_tot-i+1,:,:)=facm/facp; Ow//#:  
    D2hz(i,:,:)=2.0*epsr*epsz*dt/facp; }}grJh>tGg  
    D2hz(ie_tot-i+1,:,:)=2.0*epsr*epsz*dt/facp; AP3SOT3I  
    D3hy(i,:,:)=facm/facp;   3zWY%(8t4?  
    D3hy(ie_tot-i+1,:,:)=facm/facp; 6gakopZO  
    D4hy(i,:,:)=1.0/facp/mur/muz; _   BT d$n!'$n  
    D4hy(ie_tot-i+1,:,:)=1.0/facp/mur/muz; *PL+)2ob  
'r1&zw(  
。。。 .. Fhxg^  
5I T'u3V  

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

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

如果你比较UPML和CPML。一个差别是K的值。好象K的值作用是伸展坐标（对应SIGMA的值是吸收波的能量。）K不为1的时候，可能回应想到对消逝波的吸收。但K直需要优化。 6= 9  
E[_Z%zd^  
这只是我的个人理解。一起学习哈。

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

taflove的那个公式？

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

恩。你看他在书314页上给的UPML的例子，KMAX设为1。 e%Rg,dX  
\
qj4v^\  
对于CPML，KMAX ，SIGMA，ALPHA是需要优化的参数。 8b\XC%k  
Ro\ U T64  
恩恩。觉得taflove的书在这个问题上的讨论，单薄了些。

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