MATLAB 代码计算坐标方位角

以下是一个计算坐标方位角的基本 MATLAB 代码：

function [azimuth] = calculate_azimuth(x1,y1,x2,y2)
% 计算两个坐标点之间的方位角
% 输入参数：x1、y1、x2、y2 分别为两个坐标点的横纵坐标
% 输出参数：azimuth 为从第一个坐标点指向第二个坐标点的方位角，单位为度

if x1 == x2
    if y1 > y2
        azimuth = 180;
    else
        azimuth = 0;
    end
elseif y1 == y2
    if x1 > x2
        azimuth = 270;
    else
        azimuth = 90;
    end
else
    if x1 > x2
        azimuth = atan((y2-y1)/(x2-x1))*180/pi + 270;
    else
        azimuth = atan((y2-y1)/(x2-x1))*180/pi + 90;
    end
end

end

使用方法：

假设有两个坐标点，分别为 (1,2) 和 (3,4)，则可以使用以下命令计算它们之间的方位角：

azimuth = calculate_azimuth(1,2,3,4)

输出结果为：

azimuth = 135

这表示从第一个坐标点 (1,2) 指向第二个坐标点 (3,4) 的方位角为 135 度。