[求助]谁有色散介质的FDTD的仿真程序啊

色散介质中用H-P-D-E的推进，求D（n+1）中要用到n+1，n，n-1三个时刻的值，需要提前保存n，n-1两个时刻的值，怎么实现在推进循环中的即 .. wMa8HeBE\  
U-|]A\`)I  

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

这个很简单啊 wra0bS)4  
Fn-1=Fn |d0ZB_ci  
Fn=Fn+1

如果内存允许的话，楼主多定义个数组好了，去存放n-1时刻的数据， d?7?tL2  
等循环完了，在进行更新。 +vkmS  
preE[ ] = E[ ] l!*_[r  
preD[ ] = D[ ]

我用matlab写的， #,tT`{u1q  
  for k=upml+2:30 N,TV?Q5l7  
        for j=upml+2:upml+je ;TL.QN/l  
            for i=upml+2:upml+ie N|DfE{,  
         m?kiGC&m  
         pstore1=pstore; 529b. |  
         pstore=px; 'dwW~4|B  
         dstore1=dstore;
H9)@q3<  
         dstore=dx; hC2Fup1@  
         X%b1KG|#(  
    px(i,j,k)=C1ex(i,j,k).*  px(i,j,k)+... 7IxeSxXH  
                            C2ex(i,j,k).*((hz(i,j,k)-hz(i,j-1,k))-... sqKx?r72  
                                                        (hy(i,j,k)-hy(i,j,k-1)))./delta; k|-\[Yl.  
r,;ca6>5H  
    dx(i,j,k)=B1ex(i,j,k).*px(i,j,k)-B2ex.*pstore+B3ex.*pstore1+... CY o m  
                                                               B4ex.*dx(i,j,k)-B5ex.*dstore1;                                               AERJ]$\  
    ex(i,j,k)=C3ex(i,j,k).*ex(i,j,k)+... aDdxR:  
                            C4ex(i,j,k).*(C5ex(i,j,k).*dx(i,j,k)-... tLu&3<%  
                                                        C6ex(i,j,k).*dstore(i,j,k));             NifQsy)*%  
            end ]LB_ @#  
        end MhN;GMH  
    end h]<S0/  
想用pstore1 储存n-1时刻的值，pstore储存n时刻的值，从而实现更新，但是由于 pstore1=pstore;由于pstore开始没有定义，因此运行到这错误，决定用LS的建议试试，万分感谢