MATLAB 代码：计算坐标方位角（弧度制）

以下是一个用 MATLAB 代码计算坐标方位角的示例，该代码使用弧度制并可以根据输入的点坐标和参考点坐标计算两点之间的方位角。

% 输入点的坐标和参考点的坐标
x1 = input('请输入点的横坐标：');
y1 = input('请输入点的纵坐标：');
x2 = input('请输入参考点的横坐标：');
y2 = input('请输入参考点的纵坐标：');

% 计算两点之间的方位角
dx = x1 - x2;
dy = y1 - y2;
if dx > 0 && dy >= 0
theta = atan(dy/dx);
elseif dx > 0 && dy < 0
theta = 2*pi + atan(dy/dx);
elseif dx < 0
theta = pi + atan(dy/dx);
elseif dx == 0 && dy > 0
theta = pi/2;
elseif dx == 0 && dy < 0
theta = 3*pi/2;
elseif dx == 0 && dy == 0
theta = NaN;
end

% 输出结果
if ~isnan(theta)
    fprintf('点(%f,%f)相对于参考点(%f,%f)的方位角为%f弧度。\n', x1, y1, x2, y2, theta);
else
    fprintf('点(%f,%f)和参考点(%f,%f)重合，不存在方位角。\n', x1, y1, x2, y2);
end

该代码首先获取点的坐标和参考点的坐标。然后计算两点之间的水平和垂直距离（dx 和 dy），并根据这两个距离的符号来确定方位角的值。最后，代码输出方位角的结果，并在点和参考点重合的情况下输出一条特殊消息。