有限差分法在matlab中绘制0-100km内地震一维波动传播动图，标明坐标物理量

以下是一个可能的 Matlab 代码示例：

% 定义模拟参数
dx = 0.1;  % 网格间距，单位为 km
dt = 0.01;  % 时间步长，单位为 s
v = 5;  % 地震波速，单位为 km/s
xmax = 100;  % 模拟区域长度，单位为 km
tmax = 20;  % 模拟时间长度，单位为 s

% 初始化模拟网格
nx = round(xmax/dx) + 1;
nt = round(tmax/dt) + 1;
u = zeros(nx, nt);

% 设置初始条件
u(:,1) = exp(-((0:dx:(xmax-dx))-50).^2/(5^2));

% 使用有限差分法进行模拟
for i = 2:nt
    for j = 2:nx-1
        u(j,i) = u(j,i-1) + (v*dt/dx)^2*(u(j+1,i-1) - 2*u(j,i-1) + u(j-1,i-1));
    end
end

% 绘制动图
figure;
for i = 1:nt
    plot(0:dx:(xmax-dx), u(:,i));
    xlabel('位置（km）');
    ylabel('位移（任意单位）');
    title(['时间 = ', num2str((i-1)*dt), ' s']);
    axis([0 xmax -1 1]);
    drawnow;
end

该代码使用了有限差分法模拟了一个长度为 100km、以速度 5 km/s 传播的地震波在 20 s 内的传播过程，并绘制了动态图。代码中的坐标物理量包括位置（km）和位移（任意单位），动态图中横轴表示位置，纵轴表示位移，每一帧的标题表示当前时间（s）。