使用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)


相关推荐

  • java nextline next_java中的nextLine函数

    java nextline next_java中的nextLine函数今天在学习java异常处理的时候,下面这段程序中的nextLine()的用法怎么也看不明白。初学者看到这段代码会误以为程序中的input.nextLine()这句是多余的。其实,不使用这句的话,如果输入不是整数,程序会陷入死循环。当你不加input.nextLine()时,你输入小数,try块中给input.nextInt()就无法执行,因为小数无法被读取,我猜是小数一直留在键盘缓冲区。这时con…

  • 常用的网站后台编辑器

    常用的网站后台编辑器

  • 显示搜索dota2协调服务器,老司机教你处理搜索dota2游戏协调服务器中【操作流程】…[通俗易懂]

    显示搜索dota2协调服务器,老司机教你处理搜索dota2游戏协调服务器中【操作流程】…[通俗易懂]win7系统有很多人都喜欢使用,我们操作的过程中常常会碰到win7系统搜索dota2游戏协调服务器中的问题。如果遇到win7系统搜索dota2游戏协调服务器中的问题该怎么办呢?很多电脑水平薄弱的网友不知道win7系统搜索dota2游戏协调服务器中究竟该怎么解决?其实不难根据下面的操作步骤就可以解决问题1:DOTA2服务器蹦了之后,进入DOTA2,发现最顶端先是提示:“搜索DOTA2协调服务器中…”…

  • STM32中的NVIC详解[通俗易懂]

    STM32中的NVIC详解[通俗易懂]转载自https://blog.csdn.net/wuyuzun/article/details/72783152NVIC的全称是Nestedvectoredinterruptcontroller,即嵌套向量中断控制器。对于M3和M4内核的MCU,每个中断的优先级都是用寄存器中的8位来设置的。8位的话就可以设置2^8=256级中断,实际中用不了这么多,所以芯片厂商…

  • CSS控制文本超出指定宽度显示省略号和文本不换行

    一般的文字截断(适用于内联与块):.text-overflow{display:block;/*内联对象需加*/width:31em;/*何问起hovertree.com*/word

    2021年12月23日
  • Linux dos攻击服务器,Linux服务器如何防止DoS攻击「建议收藏」

    Linux dos攻击服务器,Linux服务器如何防止DoS攻击「建议收藏」对Linux系统所有的用户设置资源限制可以防止DoS类型攻击。如最大进程数和内存使用数量等。例如,可以在/etc/security/limits.conf中添加如下几行:*hardcore0*hardrss5000*hardnproc20然后必须编辑/etc/pam.d/login文件检查下面一行是否存在。sessionrequired/lib/security/pam_li…

发表回复

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

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