空间三维三角形的内切圆和外接圆绘制

最新推荐文章于 2022-03-28 17:03:42 发布

VIP文章 ThetaQing

最新推荐文章于 2022-03-28 17:03:42 发布

阅读量1.1k

点赞数

分类专栏：源码文章标签：几何学 matlab 算法

本文链接：https://blog.csdn.net/Thera_qing/article/details/123534349

版权

对于空间内一个三角形，绘制该三角形的外接圆和内切圆，如图所示。
在这里插入图片描述

三维空间的圆

函数plotCircle(c,r,n )

function [] = plotCircle(c,r,n )
% 画三维空间中圆
%n 法向量n
%r 圆的半径为1
%c 圆心的坐标
theta=(0:2*pi/100:2*pi)'; %theta角从0到2*pi
a=cross(n,[1 0 0]); %n与i叉乘，求取a向量
if ~any(a) %如果a为零向量，将n与j叉乘
    a=cross(n,[0 1 0]);
end
b=cross(n,a); %求取b向量
a=a/norm(a); %单位化a向量
b=b/norm(b); %单位化b向量

c1=c(1)*ones(size(theta,1),1);
c2=c(2)*ones(size(theta,1),1);
c3=c(3)*ones(size(theta,1),1);

x=c1+r*a(1)*cos(theta)+r*b(1)*sin(theta);%圆上各点的x坐标
y=c2+r*a(2)*cos(theta)+r*b(2)*sin(theta);%圆上各点的y坐标
z=c3+r*a(3)*cos(theta)+r*b(3)*sin(theta);%圆上各点的z坐标

plot3(x,y,z)

end