三角形兩角的角平分線就能確定內切圓。
結果如下:
matlab程式碼如下:
clear all;close all;clc; p=rand(3,2); %(x,y) v12=(p(2,:)-p(1,:))/norm(p(2,:)-p(1,:)); %一個頂點兩邊的歸一化向量 v13=(p(3,:)-p(1,:))/norm(p(3,:)-p(1,:)); v21=(p(1,:)-p(2,:))/norm(p(1,:)-p(2,:)); %另一個頂點兩邊的歸一化向量 v23=(p(3,:)-p(2,:))/norm(p(3,:)-p(2,:)); v1=v12+v13; %過頂點一的角平分線 k1=v1(2)/v1(1); b1=p(1,2)-k1*p(1,1); v2=v21+v23; %過頂點二的角平分線 k2=v2(2)/v2(1); b2=p(2,2)-k2*p(2,1); x0=-(b1-b2)/(k1-k2); %求兩直線交點,圓心 y0=-(-b2*k1+b1*k2)/(k1-k2); k=(p(1,2)-p(2,2))/(p(1,1)-p(2,1)); b=p(1,2)-k*p(1,1); r=(k*x0-y0+b)/sqrt(k^2+1); %求圓心到一條邊的距離,半徑 hold on; plot(p(:,1),p(:,2)); p=circshift(p,1); plot(p(:,1),p(:,2)); theta=0:0.01:2*pi; x=x0+r*cos(theta); y=y0+r*sin(theta); plot(x,y,'-',x0,y0,'.'); axis equal