Модульные векторы нормали триангуляции
F = faceNormal( возвращает модульные векторы нормали во все треугольники в 2D триангуляции. TR)faceNormal функционируйте поддерживает 2D триангуляции только. 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);
Найдите свободные граничные фасеты триангуляции и используйте их, чтобы создать 2D триангуляцию на поверхности.
[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');
