clc
clear
dx=10;??????????????????%差分網(wǎng)格10m*10m
dt=1/1000;??????????????%取樣間隔1ms
f0=10;??????????????????%頻率
c=2000;
u_0=zeros（400）;?????????%前一時(shí)刻的波場，相當(dāng)于0時(shí)刻波場
u_1=zeros（400）;?????????%當(dāng)前時(shí)刻的波場，相當(dāng)于1時(shí)刻波場，
u_x=zeros（400）;?????????%對當(dāng)前時(shí)刻的波場x方向的二階偏導(dǎo)矩陣
u_z=zeros（400）;?????????%對當(dāng)前時(shí)刻的波場z方向的二階偏導(dǎo)矩陣
u_2=zeros（400）;?????????%后一時(shí)刻的波場，待求波場
f=zeros（400）;???????????%震源矩陣
u_3=zeros（200）;?????????%無反射波場
u_4=zeros（1000400）;
u_5=[];?
z=1;
for?t=0:0.001:0.5?????????%時(shí)間共1s，每1ms取樣一次，共計(jì)1000個(gè)樣本?????
???a=pi*f0*（t-0.1）;?????%對震源進(jìn)行離散采樣?
???s=（1-2*a.^2）.*exp（-a.^2）*100;
???f（200200）=s;??????????%震源更新
???u_1（200200）=s;
???for?n=1:398??????
????????for?k=1:398
????????????u_x（nk+1）=（u_1（nk+2）-2*u_1（nk+1）+u_1（nk））/（dx^2）;
????????????%對x方向求偏導(dǎo)，行不變，列變?
????????????u_z（n+1k）=（u_1（n+2k）-2*u_1（n+1k）+u_1（nk））/（dx^2）;
????????????%對z方向求偏導(dǎo)，列不變，行變
????????end
???end
???
???for?n=1:400??????????????????
???????for?k=1:400
???????????u_2（nk）=（c^2）*（dt^2）*（u_x（nk）+u_z（nk）-f（nk））+2*u_1（nk）-u_0（nk）;
???????end
???end
???
???for?n=1:200?????????????????????
???????for?k=1:200
???????????u_3（nk）=u_1（n+100k+100）;
???????end
???end
%????for?i=1:400
%????????u_4（zi）=u_1（100i）;
%????end
%????z=z+1;
%?u_4（z）=u_1（100200）;
%?u_5（z）=u_1（150200）;
%?z=z+1;
???imagesc（u_3）
???set（gca‘XAxisLocation‘‘top‘）
???set（gca‘xticklabel‘{‘0‘;‘500‘;‘1000‘;‘1500‘;‘2000‘}‘xtick‘[0?50?100?150?200]）;
???set（gca‘yticklabel‘{‘0‘;‘500‘;‘1000‘;‘1500‘;‘2000‘}‘ytick‘[0?50?100?150?200]）;
???xlabel（‘位移/m‘）
???ylabel（‘深度/m‘）
???colorbar;
???drawnow
???u_0=u_1;????????????%更新時(shí)刻
???u_1=u_2;?
end
%??imagesc（u_4）
%???colorbar;
%???drawnow
%???title（‘距離/10m‘）
%???ylabel（‘時(shí)間/ms‘）
%?t=0:0.001:1;
%?figure（1）
%?plot（tu_4）;
%?title（‘勻速模型點(diǎn)（100200）處波形‘）;
%?xlabel（‘x軸‘）
%?ylabel（‘z軸‘）
%?figure（2）
%?plot（tu_5）;
%?title（‘勻速模型點(diǎn)（150200）處波形‘）;
%?xlabel（‘x軸‘）
%?ylabel（‘z軸‘）

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

?????文件???????2100??2018-06-08?09:57??modelling.m

?????文件????????239??2018-05-29?10:32??source.m

-----------?---------??----------?-----??----

?????????????????2339????????????????????2