trin2d.m

细胞生长结构可视化工具箱-MATLAB Toolbox1999.zip

function [n12d,n22d,n32d,b1,b2] = trin2d(n1,n2,n3)

% Find two orthogonal basis vectors for the plane containing the triangle
b1 = n2-n1;

% find projection of n3-n1 onto b1
proj = (dot(n3-n1,b1)/dot(b1,b1))*b1;
% subtract this from n3-n1 to find second basis function orthogonal to b1
b2 = (n3-n1)-proj;

% can now express pts in traingle in form n1 +xb1 +yb2 => 2-D coords
n12d = [0,0];
n22d = [(n2-n1)/b1,(n2-n1)/b2];
n32d = [(n3-n1)/b1,(n3-n1)/b2];