已知三个点的经纬度，求第一个点分别和第二个第三个点连线的夹角的matlab代码

在MATLAB中，你可以使用地理坐标系计算两点之间的方位角。这里是一个简单的示例，假设你有三个点`lat1`, `lon1` (第一个点)、`lat2`, `lon2` (第二个点) 和 `lat3`, `lon3` (第三个点)，可以使用`atan2`函数来计算角度。注意，MATLAB的`atan2`返回的是从正X轴逆时针到向量的角度，范围是(-π, π]。

% 定义点的经纬度
lat1 = lat1; % 第一点纬度
lon1 = lon1; % 第一点经度
lat2 = lat2; % 第二点纬度
lon2 = lon2; % 第二点经度
lat3 = lat3; % 第三点纬度
lon3 = lon3; % 第三点经度

% 将经度转换为弧度
lon1_rad = deg2rad(lon1);
lon2_rad = deg2rad(lon2);
lon3_rad = deg2rad(lon3);

% 计算两点间的方位角
angle_to_second = deg2rad(atan2(lat2-lat1, lon2-lon1)); % 对于第二点
angle_to_third = deg2rad(atan2(lat3-lat1, lon3-lon1)); % 对于第三点

% 确保结果在0到360度之间（如果需要）
angle_to_second = wrap_to_degrees(angle_to_second); % 回绕角度处理
angle_to_third = wrap_to_degrees(angle_to_third); 

function angle = wrap_to_degrees(angle)
    while angle > pi
        angle = angle - 2*pi;
    end
    while angle < -pi
        angle = angle + 2*pi;
    end
end

这个代码片段会计算第一个点到第二个点和第三个点的相对方向角度。请注意，如果你需要其他特定的输出格式，例如四舍五入到特定小数位或数组形式，你可能需要添加额外的处理步骤。