[求助]FDTD计算中遇到的问题

用导电介质——UPML作为吸收边界写程序遇到的问题，迭代中为什么迭代不下去？
-^+fZBU;  
qstore=qx; ml@;ngmp.  
    qx(:,2:je_tot,2:ke_tot)=C7ex(:,2:je_tot,2:ke_tot).*qx(:,2:je_tot,2:ke_tot)+... LI*=T   
                            C8ex(:,2:je_tot,2:ke_tot).*((hz(:,2:je_tot,2:ke_tot)-hz(:,1:je_tot-1,2:ke_tot))-... {8>g?4Q#  
                            (hy(:,2:je_tot,2:ke_tot)-hy(:,2:je_tot,1:ke_tot-1)))./delta; $Q"D>Qf{G  
                         D6-R>"}  
    px(:,2:je_tot,2:ke_tot)=C1ex(:,2:je_tot,2:ke_tot).*px(:,2:je_tot,2:ke_tot)+... ;.'?(iEB  
                            C21ex(:,2:je_tot,2:ke_tot).*(qx(:,2:je_tot,2:ke_tot)-qstore(:,2:je_tot,2:ke_tot)); ulE5lG0c  
    pstore=px; `*6|2  
    ex(is:is+1,js,ks)=ex(is:is+1,js,ks)+J0.*(n-ndelay)*exp(-((n-ndelay)^2/tau^2));   #N"m[$;QR  
   t W+"/<U  
    ex(:,2:je_tot,2:ke_tot)=C3ex(:,2:je_tot,2:ke_tot).*ex(:,2:je_tot,2:ke_tot)+... IdqCk0lVD  
                            C41ex(:,2:je_tot,2:ke_tot).*(C5ex(:,2:je_tot,2:ke_tot).*px(:,2:je_tot,2:ke_tot)-... ]"}BqS0  
                     &n .. /k8I6  
,p{naT%R  

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