二维弹性波数值模拟matlab

能够实现二维的弹性波数值模拟，是MATLAB程序。可以消除数值频散。

% 此程序用于模拟声波波场（未加边界条件和未消数值频散）
clear;clc;figure（2）;
% 雷克子波
tic;
dt=0.001;t=（0:40）*dt;f=30;
R=（1-2*（pi*f*t）.^2）.*exp（-（pi*f*t）.^2）;

% 模型参数设置
dx=2;dy=2;v=1000;%采样间隔和声波速度
x=-100:dx:100;y=-100:dy:100;%采样区间
x0=0;y0=0;t0=0.08;%炮点激发地点和时刻
m=length（x）;n=length（y）;

%数值计算与FCT消频散
data0=zeros（mn）;data1=data0;data2=data0;
m0=find（x==x0）;n0=find（y==y0）;l0=round（t0/dt）+1;
for k=1:l0
    if k<=length（t）;
       data1（m0n0）=R（k）;
    end
    for i=2:m-1
        for j=2:n-1
            data2（ij）=2*data1（ij）-data0（ij）+（v*dt）^2*...
               （（data1（i+1j）-2*data1（ij）+data1（i-1j））/dx^2+...
               （data1（ij+1）-2*data1（ij）+data1（ij-1））/dy^2）;
        end
    end
%     data2 = FCTforAW（data0data1data2）;
    data0=data1;data1=data2;
end

%数据成图
[XY]=meshgrid（xy）;
surf（XYdata2）;shading interp;
h=title（‘Numerical Stimulation without FCT for Acoustic Wave‘）;
set（h‘FontSize‘14‘Color‘‘b‘）;
xlabel（‘X（m）‘‘FontSize‘14‘Color‘‘b‘）;
ylabel（‘Y（m）‘‘FontSize‘14‘Color‘‘b‘）;
view（090）;%colormap（‘gray‘）;
axis square;

toc;

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

     文件       1148  2007-11-28 16:33  WaveFieldNumericalSimulation（StaggeredGrid）\AcousticWE_RGD_Simulation.m

     文件       1252  2007-11-14 19:36  WaveFieldNumericalSimulation（StaggeredGrid）\DCoef.m

     文件       3084  2007-11-14 19:36  WaveFieldNumericalSimulation（StaggeredGrid）\FCTforAW.m

     文件       2689  2007-11-14 19:36  WaveFieldNumericalSimulation（StaggeredGrid）\FCTforEW.m

     文件       5856  2008-01-24 09:27  WaveFieldNumericalSimulation（StaggeredGrid）\TotalWES.m

     目录          0  2008-01-25 10:49  WaveFieldNumericalSimulation（StaggeredGrid）

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

                14247                    7