[求助]mpl的设置问题

附件中的程序中，mpl是这样设置的    khQ@DwO*\=  
Rc2JgV  
 for ( j=0; j < n_pml; j++ ) { *o}7&Hw#9f  
        xxn = (npml-j)/npml; +aZcA#%  
        xn  = .33*pow(xxn,3.); (b#4Z  
        fj0[j] = 1.; W!o|0u!D  
        fj0[JE-j-1] = 1.; 1%$Z%?  
        fj1[j] = xn; G/N1[)  
        fj1[JE-j-1] = xn; )N3XbbV  
        gj2[j] = 1./(1.+xn); =OamN7V=  
        gj2[JE-j-1] = 1./(1.+xn); OVUJiBp  
        gj3[j] = (1.-xn)/(1.+xn); r0/aw  
        gj3[JE-j-1] = (1.-xn)/(1.+xn); :ZsAWe{%,J  
       xxn = (npml-j-.5)/npml; P<C=9@`!  
        xn  = .33*pow(xxn,3.); ;2}wrX  
        gj0[j] = 1.; qtlcY8!  
        gj0[JE-j-2] = 1.; .X\9vVJ  
        gj1[j] = xn; rr^?9M*{V  
        gj1[JE-j-2] = xn; M97MIku~9  
        fj2[j] = 1./(1.+xn); C0gO^A.d  
        fj2[JE-j-2] = 1./(1.+xn); wBI>H 7A  
        fj3[j] = (1.-xn)/(1.+xn); XP^[,)E  
        fj3[JE-j-2] = (1.-xn)/(1.+xn); V{C{y5  
      } ddzMwucjp  
W4p4[&c|  
可以看出，他完全没有处理总场的结构参数，默 .. K};~A?ET,h  
ngOGo =
 

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

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

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