Нормальные векторы модуля триангуляции
F = faceNormal( возвращает единичные нормальные векторы ко всем треугольникам в 2-D триангуляции. The 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');
