Matlab之求直角坐标系下两直线的交点坐标

目的:在直角坐标系下,求两个直线的交点坐标

一、函数的参数说明

输入参数:

PointA:直线A上的点坐标;

AngleA:直线A的倾斜角,单位度;

PointB:直线B上的点坐标;

AngleB:直线B的倾斜角,单位度;

输出参数:

TargetPoint:交点坐标

二、函数的Matlab代码

cpp 复制代码
%% 双站测向交叉定位算法
%% 自拟算法,通过左除计算
function [TargetPoint]=DoubleStationCrossLocAlgorithm(PointA,AngleA,PointB,AngleB)
    A=[1,-tand(AngleA);1,-tand(AngleB)];
    b=[PointA(2)-tand(AngleA)*PointA(1);
       PointB(2)-tand(AngleB)*PointB(1)];
   
   TargetPoint=A\b;%左除运算,计算交叉点坐标
end

2.1 方向向量转换为真方位

direction_vector:方向向量

AOA:真方位

cpp 复制代码
% 归一化向量,方向向量
direction_vector = vector / norm(vector);% direction_vector 现在是表示P1到P2方向的单位向量

AOA=mod(atan2(direction_vector(1),direction_vector(2))*180/pi+360,360);

2.2 真方位转换为倾斜角

cpp 复制代码
%% 真方位转换为倾斜角
function [ElevationAngle]=Azimuth2ElevationAngle(Azimuth)

%% 第一步计算正弦和余弦值
sine_values = sind(Azimuth);
cosine_values = cosd(Azimuth);

% 构造斜对角线单位向量
% 斜对角线单位向量的x和y分量分别为正弦值和余弦值
vector_y = cosine_values;
vector_x = sine_values;

% 计算倾斜角对应的正切值
ElevationAngle =mod((180+atan2(vector_y, vector_x)*180/pi),180);

end

三、计算的实例

相关推荐
冻柠檬飞冰走茶2 小时前
字母移动游戏题解:贪心求最长可保留连续子串(附 C 实现)
c语言·开发语言·数据结构·算法·动态规划
FRing032 小时前
vue2干货总结7(自定义指令、插槽、单页应用程序、路由的使用(VueRouter))
开发语言·前端·vue.js·vue2
雪落漂泊2 小时前
C++算法 二分查找篇(上)
开发语言·c++·算法
波力海苔夹心脆6753 小时前
C# WPF 资源字典详解:ResourceDictionary、MergedDictionaries 与自定义控件案例
开发语言·c#·wpf·入门教程·resource·资源字典
不会算法的笨小孩4 小时前
异常的使用以及会出现的问题
开发语言·c++·c++11·异常
严肃小白兔6 小时前
MATLAB 数值计算中的精度与误差问题:浮点运算、误差传播与数值稳定性
人工智能·算法·matlab
Hhy_11076 小时前
《C++深度解构06》手撕 string —— 从底层内存管理到深浅拷贝彻底吃透
c语言·开发语言·c++·学习·visual studio
starzy19906 小时前
深入理解Scala函数定义:从方法到函数的类型跨越
开发语言·scala
starzy19906 小时前
深入理解Scala递归函数:从栈溢出到尾调用优化
开发语言·scala
三*一6 小时前
Mapbox‑GL JS 多边形折点编辑工具:实现原理、功能拆解与踩坑总结
开发语言·javascript·ecmascript