求笛卡尔坐标系中两点的相对方位 matlab代码

在二维笛卡尔坐标系中，确定两点之间的相对方位通常涉及计算它们之间的角度差或方向。以下是一个简单的MATLAB代码示例，用于计算两个点(Ax, Ay)和(Bx, By)之间的方向角：

function direction = relative_direction(A, B)
    % A 和 B 是两个点的坐标，例如 A = [x1 y1], B = [x2 y2]
    
    Ax = A(1);
    Ay = A(2);
    Bx = B(1);
    By = B(2);

    % 计算向量AB的方向余弦值
    dx = Bx - Ax;
    dy = By - Ay;

    % 如果dx为0且dy为0，表示两个点重合，返回任意方向
    if abs(dx) < eps && abs(dy) < eps
        direction = 'Coincident';
    else
        angle = atan2(dy, dx); % 角度范围为(-π, π]
        
        % 根据需要调整角度至特定范围内，如0到2π
        if angle < 0
            direction = mod(angle + 2*pi, 2*pi);
        else
            direction = angle;
        end
        
        % 将角度转换为更直观的方向描述，如'North', 'Northeast', etc.
        directions = {'North', 'Northeast', 'East', 'Southeast', 'South', ' Southwest', 'West', 'Northwest'};
        direction = directions{mod(angle, 2*pi / length(directions)) + 1};
    end
end

% 示例用法
A = [0 0];
B = [1 1]; % 可以替换为任意两点坐标
direction = relative_direction(A, B);
disp(['Point B is ', direction, ' of Point A.']);