求高手帮我详细解释一下这个matlab程序,是关于地震波数值模拟的,事后另有重谢.clear;%parameter form 参数表V=2563.5; %velocity dt=0.002; %steptime 时间步长time=0.5; %snaptime 快照时间f=30; %wav
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/16 23:59:27
求高手帮我详细解释一下这个matlab程序,是关于地震波数值模拟的,事后另有重谢.clear;%parameter form 参数表V=2563.5; %velocity dt=0.002; %steptime 时间步长time=0.5; %snaptime 快照时间f=30; %wav
求高手帮我详细解释一下这个matlab程序,是关于地震波数值模拟的,事后另有重谢.
clear;
%parameter form 参数表
V=2563.5; %velocity
dt=0.002; %steptime 时间步长
time=0.5; %snaptime 快照时间
f=30; %wavelet frequency 子波频率
wavelet=0.05; %wavelet length 子波长度
dx=20; %x方向步长
dz=20; %z方向步长
a=100;b=100; %source position
line=200;row=200;
p=zeros(line+1,row+1);
pp=zeros(line+1,row+1);
ppb=zeros(line+1,row+1);
fic=fopen('dizhen.dat','rb');
[v]=fread(fic,750*685,'float32');
k=1;
for i=1:685
for j=1:750
W(i,j)=v(k);
k=k+1;
end
end
%computate the wavelet
for t=0:dt:time
% t=0:0.001:0.05
if(t
答案可发至邮箱dgqb7@126@.com,谢谢
求高手帮我详细解释一下这个matlab程序,是关于地震波数值模拟的,事后另有重谢.clear;%parameter form 参数表V=2563.5; %velocity dt=0.002; %steptime 时间步长time=0.5; %snaptime 快照时间f=30; %wav
通过有限差分来解波动方程模拟波场,画的两个图一个是所用的子波,一个是最后时刻的波场快照.
k