利用MATLAB给出仿真环境下100Hz声信号传播损失曲线，其中海深H=50m、声源深度z0=20m、接收深度z1=10m（z2=10m）、传播距离r=5000m。

在MATLAB中，我们可以使用声学传播模型如Rytov理论或者Rayleigh模型来模拟100Hz声信号在给定条件下的传播损失。首先，需要了解一些基本的声学参数，比如声速（通常取1500米每秒），密度（约为1000千克每立方米），以及声阻抗（声速乘以密度）。在这个场景中，我们有：

- 海深H = 50m
- 声源深度z0 = 20m
- 接收深度z1 = z2 = 10m (假设两个点位于同一水平线上)
- 传播距离r = 5000m

声信号的传播损失（也称衰减）可以用声压级变化来表示，通常是根据距离平方的比例来减少的。对于自由空间传播，可以使用贝塞尔定律（1/r^2），而对于海洋环境，要考虑吸收和散射效应，此时需要引入更多复杂的声学参数和修正因子。

在MATLAB中，你可以按照以下步骤操作：

1. 定义声学参数，并计算初始声压级。
2. 根据声速、深度和距离，估计声波在空气和水中传播的额外损耗。
3. 使用适当的函数或自定义函数，如`soundPressureLoss`或`rayleighRangeLoss`，计算总的声压损失。
4. 创建一个范围从0到r的距离数组，并对每个距离计算相应的声压损失。
5. 绘制声压损失随距离的变化曲线。

以下是部分伪代码示例：

% 初始化参数
c = 1500; % 声速（m/s）
rho = 1000; % 密度（kg/m^3）
Z0 = c * rho; % 声阻抗
frequency = 100; % 频率（Hz）

% 定义深度和距离
H = 50;
z0 = 20;
z1 = z2 = 10;
r = 5000;

% 计算声压初始值
p_initial = sqrt(2 * Z0 * frequency / (4 * pi));

% 自定义声压损失函数
loss_function(r) = ... % 这里填写实际的声压损失计算表达式，可能是贝塞尔定律加上海洋特定损耗项

% 范围数组
distance_range = linspace(0, r, 1000);

% 计算并绘制声压损失曲线
sound_pressure_loss = loss_function(distance_range);
plot(distance_range, sound_pressure_loss, 'LineWidth', 2);
xlabel('Distance (m)');
ylabel('Sound Pressure Loss (dB)');
title(['100Hz声信号在海深50m的传播损失曲线 (r=' num2str(r) 'm)']);