%交錯網格時間二階空間二階各向同性均勻介質%
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