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

色散介质中用H-P-D-E的推进，求D（n+1）中要用到n+1，n，n-1三个时刻的值，需要提前保存n，n-1两个时刻的值，怎么实现在推进循环中的即 .. N{SQ(%V  
s'aip5P  

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

这个很简单啊 #,{+3Y&5-+  
Fn-1=Fn /ywD{*  
Fn=Fn+1

如果内存允许的话，楼主多定义个数组好了，去存放n-1时刻的数据， H1!iP$1#V  
等循环完了，在进行更新。 >-E<n8  
preE[ ] = E[ ] ,_!6U  
preD[ ] = D[ ]

我用matlab写的， IB.yU,v  
  for k=upml+2:30 +@X5!S6  
        for j=upml+2:upml+je v;{{ y-  
            for i=upml+2:upml+ie I[o*RKT'"  
         x0a.!  
         pstore1=pstore; T~sTBGcv  
         pstore=px; wLuv6\E  
         dstore1=dstore; 6VD1cb\lF  
         dstore=dx; 
]L2Oz  
         }~Q"s2  
    px(i,j,k)=C1ex(i,j,k).*  px(i,j,k)+... C@o%J.9"#  
                            C2ex(i,j,k).*((hz(i,j,k)-hz(i,j-1,k))-... iq?#rb P#I  
                                                        (hy(i,j,k)-hy(i,j,k-1)))./delta; Oc9#e+_&  
l? #xAZx&_  
    dx(i,j,k)=B1ex(i,j,k).*px(i,j,k)-B2ex.*pstore+B3ex.*pstore1+... *cWmS\h|  
                                                               B4ex.*dx(i,j,k)-B5ex.*dstore1;                                               dsK&U\ej}  
    ex(i,j,k)=C3ex(i,j,k).*ex(i,j,k)+... 0ZAj=u@O  
                            C4ex(i,j,k).*(C5ex(i,j,k).*dx(i,j,k)-... {>0V[c[~  
                                                        C6ex(i,j,k).*dstore(i,j,k));             cIXwiC8t  
            end PXP`ZLF  
        end 'NZGQebK  
    end 1|,Pq9  
想用pstore1 储存n-1时刻的值，pstore储存n时刻的值，从而实现更新，但是由于 pstore1=pstore;由于pstore开始没有定义，因此运行到这错误，决定用LS的建议试试，万分感谢