Staggered-grid-finite-difference 很简单的交错网格有限差分正演程序

很简单的交错网格有限差分正演程序，程序通俗易懂，非常适合初学者学习

%交错网格时间二阶空间二阶各向同性均匀介质%
clear;clc;
dt=0.002;
dh=6;
p=0.002;
l=4000;
m=2000;
ml=8000;
a1=dt/p/dh;
a2=dt*ml/dh;
a3=l*dt/dh;
a4=dt*m/dh;
u2=zeros（102102）;
v2=u2;
r1=u2;t1=u2;h1=u2;
u0=u2;v0=u2;
r3=u2;t3=u2;h3=u2;

for k=1:120
    for i=2:100
        for j=2:100
            if i==50&j==50
                u2（ij）=u0（ij）+a1*（r1（ij）-r1（i-1j）+h1（ij+1）-h1（ij））+（1-2*（pi*30*（k-30）/1000）^2）*exp（-（pi*30*（k-30）/1000）^2）;
                v2（ij）=v0（ij）+a1*（h1（i+1j）-h1（ij）+t1（ij）-t1（ij-1））+（1-2*（pi*30*（k-30）/1000）^2）*exp（-（pi*30*（k-30）/1000）^2）;
              
            else
                u2（ij）=u0（ij）+a1*（r1（ij）-r1（i-1j）+h1（ij+1）-h1（ij））;
                v2（ij）=v0（ij）+a1*（h1（i+1j）-h1（ij）+t1（ij）-t1（ij-1））;
            end
        end
    end
    for

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1496  2014-01-08 21:21  交错\jiaocuo1.m

     文件         98  2014-01-08 21:21  交错\zibo.m

     文件      30784  2014-01-08 21:21  交错\交错100ms.png

     文件       3163  2014-01-08 21:21  交错\子波.png

     目录          0  2014-03-16 22:05  交错

----------- ---------  ---------- -----  ----

                35541                    5