matlab練習程式(三角形內切圓)

Dsp Tian發表於2014-09-27

三角形兩角的角平分線就能確定內切圓。

結果如下:

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

 

相關文章