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

色散介质中用H-P-D-E的推进，求D（n+1）中要用到n+1，n，n-1三个时刻的值，需要提前保存n，n-1两个时刻的值，怎么实现在推进循环中的即 .. ,b b/ $   
As5*)o"&  

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

这个很简单啊 1^R:[L4R`  
Fn-1=Fn /2jw]ekQ'  
Fn=Fn+1

如果内存允许的话，楼主多定义个数组好了，去存放n-1时刻的数据， X
A%a7Xtni  
等循环完了，在进行更新。 ".=LzjE<gv  
preE[ ] = E[ ] Z
O!  
preD[ ] = D[ ]

我用matlab写的， qXQ/M]  
  for k=upml+2:30 I )LO@  
        for j=upml+2:upml+je mm5y'=#  
            for i=upml+2:upml+ie /KWR08ftp  
         L`!M3c@u  
         pstore1=pstore; +z2+z  
         pstore=px; +.RC{o,  
         dstore1=dstore; l6C^,xU~IX  
         dstore=dx; vFL\O  
         vj23j[!|  
    px(i,j,k)=C1ex(i,j,k).*  px(i,j,k)+... |4F3Gu  
                            C2ex(i,j,k).*((hz(i,j,k)-hz(i,j-1,k))-... kK]^q|vb6  
                                                        (hy(i,j,k)-hy(i,j,k-1)))./delta; {D(_"  
vGT#BS%  
    dx(i,j,k)=B1ex(i,j,k).*px(i,j,k)-B2ex.*pstore+B3ex.*pstore1+... bBFdr  
                                                               B4ex.*dx(i,j,k)-B5ex.*dstore1;                                               08!pLE  
    ex(i,j,k)=C3ex(i,j,k).*ex(i,j,k)+... j+_pF<$f:  
                            C4ex(i,j,k).*(C5ex(i,j,k).*dx(i,j,k)-... 8] BOq:
 
                                                        C6ex(i,j,k).*dstore(i,j,k));             $YyN-C  
            end J}035  
        end t`XYY  
    end L,XWX8  
想用pstore1 储存n-1时刻的值，pstore储存n时刻的值，从而实现更新，但是由于 pstore1=pstore;由于pstore开始没有定义，因此运行到这错误，决定用LS的建议试试，万分感谢