[求助]mpl的设置问题

附件中的程序中，mpl是这样设置的    $\\lx_)  
}NmNanW^  
 for ( j=0; j < n_pml; j++ ) { V[hK2rVH.  
        xxn = (npml-j)/npml; (GU9p>2  
        xn  = .33*pow(xxn,3.); ;e s^R?z  
        fj0[j] = 1.; &a)vdlZSE=  
        fj0[JE-j-1] = 1.; WS0JS'  
        fj1[j] = xn; 1~5trsB+5  
        fj1[JE-j-1] = xn; g`,AaWlF  
        gj2[j] = 1./(1.+xn); }GZbo kWg.  
        gj2[JE-j-1] = 1./(1.+xn); 'fW#7W  
        gj3[j] = (1.-xn)/(1.+xn); \2!!L=&4G  
        gj3[JE-j-1] = (1.-xn)/(1.+xn); \7 a4uc  
       xxn = (npml-j-.5)/npml; Vb4;-?s_  
        xn  = .33*pow(xxn,3.); nkY@_N
 
        gj0[j] = 1.; r6.N4eW.L  
        gj0[JE-j-2] = 1.; Je7RrCz  
        gj1[j] = xn; JdNF-64ky  
        gj1[JE-j-2] = xn; M'5'O;kn  
        fj2[j] = 1./(1.+xn); j$@tK0P  
        fj2[JE-j-2] = 1./(1.+xn); \L]|-f(4  
        fj3[j] = (1.-xn)/(1.+xn); Dgi~rr1`'s  
        fj3[JE-j-2] = (1.-xn)/(1.+xn); hU {-a`  
      } 9y]J/1#  
H7&xLYQ2  
可以看出，他完全没有处理总场的结构参数，默 .. (.X]F_*sc  
I3y9:4  

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

就相当于通用的PML设置，可以根据介质、损耗进行公式修正，介质参数根据PML的无反射的等效条件

随便找本参考书，把公式重新推导一次就可以了。