%?此程序用于模擬聲波波場(chǎng)（未加邊界條件和未消數(shù)值頻散）
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）;

%?模型參數(shù)設(shè)置
dx=2;dy=2;v=1000;%采樣間隔和聲波速度
x=-100:dx:100;y=-100:dy:100;%采樣區(qū)間
x0=0;y0=0;t0=0.08;%炮點(diǎn)激發(fā)地點(diǎn)和時(shí)刻
m=length（x）;n=length（y）;

%數(shù)值計(jì)算與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

%數(shù)據(jù)成圖
[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;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????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