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

用导电介质——UPML作为吸收边界写程序遇到的问题，迭代中为什么迭代不下去？ ei@3,{~5  
qstore=qx; iAK/d)bq  
    qx(:,2:je_tot,2:ke_tot)=C7ex(:,2:je_tot,2:ke_tot).*qx(:,2:je_tot,2:ke_tot)+... f3!n$lj  
                            C8ex(:,2:je_tot,2:ke_tot).*((hz(:,2:je_tot,2:ke_tot)-hz(:,1:je_tot-1,2:ke_tot))-... [a@B =E  
                            (hy(:,2:je_tot,2:ke_tot)-hy(:,2:je_tot,1:ke_tot-1)))./delta; +kM\ D~D1  
                         H.:9:I[n  
    px(:,2:je_tot,2:ke_tot)=C1ex(:,2:je_tot,2:ke_tot).*px(:,2:je_tot,2:ke_tot)+... kfXS_\@iW1  
                            C21ex(:,2:je_tot,2:ke_tot).*(qx(:,2:je_tot,2:ke_tot)-qstore(:,2:je_tot,2:ke_tot)); {gsW(T>)  
    pstore=px; >rKhlUD  
    ex(is:is+1,js,ks)=ex(is:is+1,js,ks)+J0.*(n-ndelay)*exp(-((n-ndelay)^2/tau^2));   J;~E<_"Hn  
   +$PFHXB  
    ex(:,2:je_tot,2:ke_tot)=C3ex(:,2:je_tot,2:ke_tot).*ex(:,2:je_tot,2:ke_tot)+... T8U[xu.>  
                            C41ex(:,2:je_tot,2:ke_tot).*(C5ex(:,2:je_tot,2:ke_tot).*px(:,2:je_tot,2:ke_tot)-... o5P&JBX<  
                     &n .. _ \l HI  
r&SO:#rOSM  

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