使用Matlab计算两条线的交点及三角形垂心

使用Matlab计算两条线的交点及三角形垂心1、我用了12个点,每两个点一条线,能组成6条线;2、每两条线相交于一点,一共有3个点;3、3个点构成了三角形的三个顶点,然后再根据三个顶点计算三角形的垂心。

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

Matlab版本:R2016A

操作系统:Win-8

为了让整个流程比较完整,

1、我用了12个点,每两个点一条线,能组成6条线;

2、每两条线相交于一点,一共有3个点;

3、3个点构成了三角形的三个顶点,然后再根据三个顶点计算三角形的垂心。

1、现在先给出12个点的坐标(坐标可以随意设置,只要构成的线不是平行没有交点即可)

(294.94, 2959.12) (668.61, 1464.42)

(1678.70, 2848.19) (1731.25, 1914.00)
(540.16, 1978.22) (1731.25, 1919.83)

(335.80, 2778.12) (1696.22, 2678.86)
(1106.51, 892.23) (1404.28, 1832.25)

(1351.74, 880.56) (1591.12, 1820.58)

上述点是我从一幅图片里面点选出来的,将每个点的橫坐标放在一个数组里,纵坐标放在一个数组里,然后进行绘制:

X = [294.94, 668.61, 1678.70, 1731.25, 540.16, 1731.25, 335.80, 1696.22, 1106.51, 1404.28, 1351.74, 1591.12];  
Y = [2959.12, 1464.42, 2848.19, 1914.00, 1978.22, 1919.83, 2778.12, 2678.86, 892.23, 1832.25, 880.56, 1820.58];
plot(X, Y, '.');

Jetbrains全家桶1年46,售后保障稳定

绘制出来的结果如图所示:

使用Matlab计算两条线的交点及三角形垂心

将12个点组成6条线:

for i = 1: 6
    plot([X(2*i-1), X(2*i)], [Y(2*i-1), Y(2*i)]);
    hold on;
end

依次绘制出的线如图所示:

使用Matlab计算两条线的交点及三角形垂心

2、下面计算每两条线之间的交点(即上图所示的1,2一组,3,4一组,5,6一组计算其交点)
两条直线相交,必然需要求出两条直线的表达式,每条直线的表达式可以用 y = ax + b来表示,下面用matlab求解a和b两个变量,a即为直线的斜率,b为直线的截距。

a = zeros(6, 1);
b = zeros(6, 1);
for i = 1: 6
    a(i, 1) = (Y(2*i) - Y(2*i-1)) / (X(2*i) - X(2*i-1));
    b(i, 1) = Y(2*i) - a(i, 1) * X(2*i);
end

上述公式即求出了全部直线的方程式,下面计算它们两两之间的交点

points = zeros(3, 2);
for i = 1: 3
    % 两个直线方程之差的函数
    fun = @(x) (a(2*i-1, 1) * x + b(2*i-1, 1)) - (a(2*i, 1) * x + b(2*i, 1));
    
    % 两个直线方程之差函数值为0时的x值,由得到的x值计算y值
    points(i, 1) = fzero(fun, 2);
    points(i, 2) = a(2*i) * points(i, 1) + b(2*i, 1);
end

得到直线交点如图所示:

使用Matlab计算两条线的交点及三角形垂心

3、将三个点连线绘制,可以得到三角形

hold on;
plot([points(:, 1); points(1, 1)], [points(:, 2); points(1, 2)])

绘制得到结果如图所示:

使用Matlab计算两条线的交点及三角形垂心

计算三角形每条边的垂线,需要先计算每条线的斜率:

% 计算各个顶点连线的方程式
slope1_2 = (points(2, 2) - points(1, 2)) / (points(2, 1) - points(1, 1));
slope2_3 = (points(3, 2) - points(2, 2)) / (points(3, 1) - points(2, 1));
slope1_3 = (points(3, 2) - points(1, 2)) / (points(3, 1) - points(1, 1));
b1_2 = points(2, 2) - slope1_2 * points(2, 1);
b2_3 = points(3, 2) - slope2_3 * points(3, 1);
b1_3 = points(1, 2) - slope1_3 * points(1, 1);
slopevertical1_2 = -1 / slope1_2;
slopevertical2_3 = -1 / slope2_3;
slopevertical1_3 = -1 / slope1_3;

verticalb1_2 = points(3, 2) - slopevertical1_2 * points(3, 1);
verticalb2_3 = points(1, 2) - slopevertical2_3 * points(1, 1);
verticalb1_3 = points(2, 2) - slopevertical1_3 * points(2, 1);

这样就求出了三个边的垂线,下面对这三个垂线进行绘制,按照上面的方程计算出两个点,两个点的横坐标一个选择很小,一个选择很大,这样可以保证画出来的直线足够长,方便观察(这个方法似乎有些简单粗暴~~~~~~~~,不过还是能够绘制出来的)

% 设置预绘制直线的两个端点的X值
verticalLinePoints = zeros(6, 2)
for i = 1: 3
    verticalLinePoints(2*i-1, 1) = 0;
    verticalLinePoints(2*i, 1) = 10000;
    if i == 1
         slopevertical = slopevertical1_2;
         verticalb = verticalb1_2;
    elseif i == 2
         slopevertical = slopevertical2_3;
         verticalb = verticalb2_3;
    else
	 slopevertical = slopevertical1_3;
	 verticalb = verticalb1_3;
    end
    
   % 计算两个端点的y值
    verticalLinePoints(2*i-1, 2) = slopevertical * verticalLinePoints(2*i-1, 1) + verticalb;
    verticalLinePoints(2*i, 2) = slopevertical * verticalLinePoints(2*i, 1) + verticalb;
end

然后绘制所得结果(为了使图更加简洁,开始的点线都去掉了,仅保留了三角形)

hold on;	% 为了保留前面绘制的图像,使用hold on
plot(verticalLinePoints(1:2, 1), verticalLinePoints(1:2, 2));
plot(verticalLinePoints(3:4, 1), verticalLinePoints(3:4, 2));
plot(verticalLinePoints(5:6, 1), verticalLinePoints(5:6, 2));
axis equal;	% 使横纵轴单位统一

绘制得到结果如图所示:

使用Matlab计算两条线的交点及三角形垂心

由图或者基础几何可知,三条垂线交于一点,称之为垂心(中心),求垂心的方法与前面求两条线的交点的方法一样

fun = @(x) (slopevertical1_2 * x + verticalb1_2) - (slopevertical1_3 * x + verticalb1_3);
finalX = fzero(fun, 2);
finalY = slopevertical1_2 * finalX + verticalb1_2;
hold on;
plot(finalX, finalY, 'o');

至此,计算垂心的工作已经完成,对垂心进行绘制,绘制的结果如图所示:

使用Matlab计算两条线的交点及三角形垂心

代码下载地址:http://download.csdn.net/detail/zhangwellyear/9796085

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/204267.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • SqlSessionTemplate探究

    SqlSessionTemplate探究问题就是:无论是多个dao使用一个SqlSessionTemplate,还是一个dao使用一个SqlSessionTemplate,SqlSessionTemplate都是对应一个sqlSession,当多个web线程调用同一个dao时,它们使用的是同一个SqlSessionTemplate,也就是同一个SqlSession,如何保证线程安全,关键就在于代理:(1)首先,通过如下代码创建代理类

  • linux安装配置SVN并设置钩子

    linux安装配置SVN并设置钩子

    2021年10月29日
  • hdu 4455 Substrings (DP 预处理思路)「建议收藏」

    hdu 4455 Substrings (DP 预处理思路)

  • Win10搭建ftp(含设置用户名和密码)测试成功!!!

    Win10搭建ftp(含设置用户名和密码)测试成功!!!Win10搭建ftp(含设置用户名和密码)测试成功1、打开ftp服务方法:win+R输入control打开控制面板点击程序与功能→启动或关闭Windows功能,选择一下选项,打开ftp服务2、打开ISS管理器方法:win+R输入inetmgr打开iss管理器1)、网站—>添加FTP站点…—>站点信息2)、输入IP和端口号,IP就是自己电脑的ip,端口号最好改一下,21…

  • Servlet-JSP入门

    Servlet-JSP入门JSP作用:主要作用是代替Servlet程序回传html页面的数据。因为Servlet程序回传html页面数据是一件非常繁琐的事情,开发成本和维护成本都很高。一般的Servlet回传页面如下protected void doGet(HttpServletRequest req, HttpServletResponse resp) throws ServletException, IOException { resp.setContentType(“text/html;charset=U

  • 循环队列 基本概念「建议收藏」

    循环队列 基本概念「建议收藏」循环队列是队列的一种特殊形式。首先介绍队列,然后引申出循环队列。队列又称为“先进先出”(FIFO)线性表限定插入操作只能在队尾进行,而删除操作只能在队首进行队列也可以采用顺序存储结构或链表结构来实现,分别称为顺序队列和链队列队列的顺序表示—顺序队列用一组连续的存储单元依次存放从队首到队尾的元素,附设两个指针head和tail分别指向队首元素和队尾元素的位置,(有的地方用f…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号