Нормальные векторы единиц триангуляции
F = faceNormal( возвращает единичные нормальные векторы ко всем треугольникам в 2-D триангуляции. TR)faceNormal функция поддерживает только 2-D триангуляции. F представляет собой матрицу из трех столбцов, где каждая строка содержит нормальные координаты единицы измерения, соответствующие треугольнику в TR.ConnectivityList.
Вычислите и постройте график единичных нормальных векторов к граням триангуляции на сферической поверхности.
Создайте набор точек на сферической поверхности.
rng default;
theta = rand([100,1])*2*pi;
phi = rand([100,1])*pi;
x = cos(theta).*sin(phi);
y = sin(theta).*sin(phi);
z = cos(phi);Триангулировать сферу с помощью delaunayTriangulation функция.
DT = delaunayTriangulation(x,y,z);
Найдите свободные граничные грани триангуляции и используйте их для создания 2-D триангуляции на поверхности.
[T,Xb] = freeBoundary(DT); TR = triangulation(T,Xb);
Вычислите нормали центров и граней каждой треугольной грани в TR.
P = incenter(TR); F = faceNormal(TR);
Постройте график триангуляции вместе с центрами и нормалями грани.
trisurf(T,Xb(:,1),Xb(:,2),Xb(:,3), ... 'FaceColor','cyan','FaceAlpha',0.8); axis equal hold on quiver3(P(:,1),P(:,2),P(:,3), ... F(:,1),F(:,2),F(:,3),0.5,'color','r');
