Модульные векторы нормали триангуляции
F = faceNormal(TR)F = faceNormal(TR,ID)F = faceNormal( возвращает модульные векторы нормали во все треугольники в 2D триангуляции. Функция TR)faceNormal поддерживает 2D триангуляции только. F является матрицей с тремя столбцами, где каждая строка содержит модульные нормальные координаты, соответствующие треугольнику в TR.ConnectivityList.
Вычислите и постройте модульные векторы нормали к фасетам триангуляции на сферической поверхности.
Создайте набор точек на сферической поверхности.
theta = gallery('uniformdata',[100,1],0)*2*pi; phi = gallery('uniformdata',[100,1],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');
