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