写一个输入为入射角、双程传播时间、声速剖面，输出为水平距离与深度的声线跟踪matlab算法

以下是一个简单的声线跟踪Matlab算法，输入参数为入射角（单位：度）、双程传播时间（单位：秒）、声速剖面（包含深度和声速两列数据），输出为水平距离和深度两列数据。

function [x,z] = sound_trace(theta, t, sound_profile)
% theta: 入射角，单位：度
% t: 双程传播时间，单位：秒
% sound_profile: 包含深度和声速两列数据的矩阵
% x: 水平距离，单位：米
% z: 深度，单位：米

% 将角度转换为弧度
theta = theta * pi / 180;

% 获取声速剖面数据
z = sound_profile(:,1); % 深度
c = sound_profile(:,2); % 声速

% 计算声程
d = cumsum([0; diff(z)]); % 每个深度层的厚度
s = sin(theta); % 入射角的正弦值
t1 = d ./ (c * sqrt(1 - s^2)); % 到达每个深度层的时间
t2 = t - t1; % 返回每个深度层的时间
d2 = t2 .* (c * sqrt(1 - s^2)); % 返回每个深度层的距离

% 计算水平距离
x = cumsum(d2 .* cos(theta)); % 每个深度层的水平距离

注意，上述算法仅适用于水平层状介质中的声线跟踪，对于复杂的海底地形和不规则的声速剖面，可能需要使用更为复杂的算法。