Matlab绘制射线声学中的声线的代码

请参考以下Matlab代码实现射线声学中的声线绘制：

%定义声速和密度
c = 1500; % m/s，声速
rho = 1000; % kg/m^3，密度

%定义发射角度和接收角度
theta_i = 30; % 初始入射角度，单位为度
theta_r = -30; % 接收角度，单位为度

%定义声线的参数
t = 0 : 0.0001 : 0.15; % 定义时间范围
x = (c * t) / 2; % 计算声线的路径长度
y = (2 * rho * x) ./ cosd(theta_i); % 计算声线的深度

%绘制声线
figure;
plot(x, y, '-r', 'LineWidth', 2); % 坐标轴为声线的长度和深度，绘制红色线条
xlabel('Ray Length (m)');
ylabel('Depth (m)');
title('Acoustic Ray Path');
grid on;

该代码的基本思路是使用浅水声学模型计算声线在给定入射和接收角度时的路径，然后将其绘制出来。需要注意的是，在实际应用中需要根据具体情况设置声速和密度等参数，并且还需要考虑海底地形和介质等因素对声线的影响。