有限差分法在matlab中绘制0-100km内地震传播动图，一维波动传播模拟

以下是一个简单的 MATLAB 代码示例，使用有限差分法模拟地震波传播过程，并绘制传播动画。

% 定义模拟参数
nx = 1000;       % 空间网格数
nt = 500;        % 时间步数
dx = 100;        % 空间步长，单位：m
dt = 0.01;       % 时间步长，单位：s
vp = 3000;       % 纵波速度，单位：m/s

% 初始化波场
u = zeros(nx, nt);

% 初始化边界条件
u(1,:) = 0;
u(nx,:) = 0;

% 初始化震源
src = zeros(1, nt);
src(1) = 1;

% 循环迭代计算波场
for i = 2:nx-1
    for j = 2:nt
        u(i,j) = (vp*dt/dx)^2 * (u(i+1,j-1) - 2*u(i,j-1) + u(i-1,j-1)) + 2*u(i,j-1) - u(i,j-2);
    end
    % 添加震源
    u(i,1) = u(i,1) + src(1);
end

% 绘制传播动画
figure;
for i = 1:nt
    plot(u(:,i))
    axis([0 nx -1 1])
    pause(0.01)
end

在此代码中，我们使用了二阶有限差分法来模拟地震波传播过程，其中 `vp` 为纵波速度，`src` 为震源，`u` 为波场。我们通过循环迭代计算每个时刻的波场，最后使用 `plot` 函数绘制传播动画。

请注意，在实际地震模拟中，需要考虑更加复杂的物理现象，比如横波、衰减等，需要使用更为复杂的模拟方法和模型。