作业帮我在用matlab编写一地震波数值模拟程序,运行时总显示有错,求高手指教.程序如下
来源：学生作业帮 编辑：拍题作业网作业帮 分类：综合作业 时间：2024/04/28 15:33:56
我在用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
V(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
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
V(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
问题出现V上面.
在第3行,程序明确V是一个数值而已
但是第19行到第24行,程序已经把V改变成为一个矩阵,行为685,列为750
这样直接导致第39行中(V*dt/dx).^2*(pp(x+1,z)+pp(x-1,z)-2*pp(x,z))该项的结果也是这样一个685*750的矩阵
这样的矩阵是不能直接赋值给第39行的左边的数值p(x,z)的.
在第3行,程序明确V是一个数值而已
但是第19行到第24行,程序已经把V改变成为一个矩阵,行为685,列为750
这样直接导致第39行中(V*dt/dx).^2*(pp(x+1,z)+pp(x-1,z)-2*pp(x,z))该项的结果也是这样一个685*750的矩阵
这样的矩阵是不能直接赋值给第39行的左边的数值p(x,z)的.
我在用matlab编写一地震波数值模拟程序,运行时总显示有错,求高手指教.程序如下
求高手帮我详细解释一下这个matlab程序,是关于地震波数值模拟的,事后另有重谢.
刚学matlab菜鸟一枚,为什么书上的程序运行不了啊,求大神指教.
编写一个程序,在同一行上显示数字1到4,相邻的数字用一个空格分开.按照如下要求编写该程序：a) 使用一
matlab中寻找回文数问题,程序如下,运行不了,求高手修改找错!
求一程序：用C语言编写,要求如下：
运行如下PHP程序 该程序运行结束时显示的结果
求编写一个matlab程序
matlab程序如下：为什么提示运行出错?
求软件高手帮我编写如下程序：求1到N的阶乘和
MATLAB的一个程序,运行出错,画不出图,求大牛!程序如下：
我的matlab运行的程序的收敛性曲线没有曲线,程序如下