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

三、计算的实例

相关推荐
2601_966949655 小时前
如何根据研究的起止时间截取分钟级行情?Python 实战与 QuantDash 数据处理指南
开发语言·python·pandas·量化交易·1024程序员节·股票数据·quantdash
晨曦中的暮雨7 小时前
限流算法——固定窗口、滑动窗口、漏桶、令牌桶
java·开发语言·算法
喵了几个咪8 小时前
GoWind CMS 拆解:在 Go 里给 TipTap 搭一套 Yjs 协同后端
开发语言·后端·golang
程序员清风8 小时前
数据处理脚本工程化:日志、配置、命令行与测试
服务器·开发语言·数据库
2601_9628857210 小时前
如何用 Python 把 K 线重采样成任意周期?
开发语言·python
Wang's Blog10 小时前
Java 项目实战: 交友项目 - 圈子功能实现:总结
java·开发语言·交友
下一站春风里11 小时前
位运算完整学习笔记:按位与(&)+ 异或(^)
开发语言·c++·算法
天外飞雨道沧桑12 小时前
我在浏览器里练了一个 792 参数的大模型,它学会了补全 JavaScript
开发语言·javascript·ai·ecmascript·ai编程
hhhjhl12 小时前
Yup 跨端三方库仓颉适配全流程
java·开发语言
运维行者_12 小时前
网络监控工具怎么选?从设备监控到智能运维的能力清单
运维·服务器·开发语言·网络·人工智能·支持向量机·php