设为首页收藏本站

EPS数据狗论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

查看: 1625|回复: 0

matlab练习程序(三角形内切圆)

[复制链接]

30

主题

333

金钱

471

积分

入门用户

发表于 2019-7-12 16:04:20 | 显示全部楼层 |阅读模式

三角形两角的角平分线就能确定内切圆。
结果如下:
1.jpg

matlab代码如下:
  1. clear all;close all;clc;

  2. p=rand(3,2);    %(x,y)

  3. v12=(p(2,:)-p(1,:))/norm(p(2,:)-p(1,:));   
  4. v13=(p(3,:)-p(1,:))/norm(p(3,:)-p(1,:));
  5. %一个顶点两边的归一化向量

  6. v21=(p(1,:)-p(2,:))/norm(p(1,:)-p(2,:));  
  7. v23=(p(3,:)-p(2,:))/norm(p(3,:)-p(2,:));
  8. %另一个顶点两边的归一化向量

  9. v1=v12+v13;            
  10. k1=v1(2)/v1(1);
  11. b1=p(1,2)-k1*p(1,1);
  12. %过顶点一的角平分线

  13. v2=v21+v23;            
  14. k2=v2(2)/v2(1);
  15. b2=p(2,2)-k2*p(2,1);
  16. %过顶点二的角平分线

  17. x0=-(b1-b2)/(k1-k2);            
  18. y0=-(-b2*k1+b1*k2)/(k1-k2);
  19. %求两直线交点,圆心

  20. k=(p(1,2)-p(2,2))/(p(1,1)-p(2,1));  
  21. b=p(1,2)-k*p(1,1);
  22. r=(k*x0-y0+b)/sqrt(k^2+1);     
  23. %求圆心到一条边的距离,半径

  24. hold on;
  25. plot(p(:,1),p(:,2));
  26. p=circshift(p,1);
  27. plot(p(:,1),p(:,2));

  28. theta=0:0.01:2*pi;
  29. x=x0+r*cos(theta);
  30. y=y0+r*sin(theta);
  31. plot(x,y,'-',x0,y0,'.');
  32. axis equal
复制代码

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

关闭

站长推荐上一条 /1 下一条

客服中心
关闭
在线时间:
周一~周五
8:30-17:30
QQ群:
653541906
联系电话:
010-85786021-8017
在线咨询
客服中心

意见反馈|网站地图|手机版|小黑屋|EPS数据狗论坛 ( 京ICP备09019565号-3 )   

Powered by BFIT! X3.4

© 2008-2028 BFIT Inc.

快速回复 返回顶部 返回列表