使用matlab编程实现二维非均匀层状介质模型，模型网格数为200*200，第一层的纵波速度、横波速度和密度分别为1500m/s、500m/s和1000kg/m3；第二深度网格数为120，该层纵波速度、横波速度和密度分别为2500m/s、800m/s和1200kg/m3

在MATLAB中，你可以通过数值模拟库如PDE Toolbox或者Finite Element Toolbox来实现二维非均匀层状介质模型。这里是一个简单的概述步骤：

1. **初始化变量**：
- 定义网格尺寸（n = 200），深度网格（m = 120）以及各层的速度和密度数组。

n = 200;
m = 120;
v1z = 1500; % 第一层纵波速度 (m/s)
v1l = 500;  % 第一层横波速度 (m/s)
rho1 = 1000; % 第一层密度 (kg/m^3)
v2z = 2500; % 第二层纵波速度 (m/s)
v2l = 800;  % 第二层横波速度 (m/s)
rho2 = 1200; % 第二层密度 (kg/m^3)

2. **创建网格**：
- 创建一个二维矩阵表示网格，其中第一层的上半部分（1:100, :)是第一层，下半部分是第二层。

[x, y] = meshgrid(1:n, 1:m); % 网格坐标
depth = reshape(y, [m, n]);  % 深度向量
medium_type = depth <= (m/2 + 1); % 判断当前点属于哪一层

3. **设置属性**：
- 根据`medium_type`分配各个点的纵波速度、横波速度和密度。

[vx, vy] = deal(v1z, v1l);
vx(depth > m/2) = v2z; % 纵波速度切换
vy(depth > m/2) = v2l; % 横波速度切换
rho = ones(n*m, 1);     % 先假设所有点密度相同
rho(depth > m/2) = rho2; % 密度切换

4. **建立PDE模型**：
- 对于二维弹性波方程（例如波动方程），可以利用PDE Toolbox中的`pdepe`函数或其他适合的工具来描述。

5. **求解和可视化结果**：
- 解决PDE，得到位移场，然后可视化工件响应。

% 这里省略具体的求解步骤，因为实际的代码会更复杂，包括设定边界条件、时间步长等
[ux, uy] = pdepe(...); % 求解偏微分方程

% 可视化结果
surf(x, y, uy); % 或其他适合的二维图形展示纵波速度分布
xlabel('X');
ylabel('Y');
title('二维非均匀层状介质波速分布');