MatlabTools(1)mesh2tri

mutougeda

将结构的四边形网格转换为三角形网格。
Demo:
, E* Y) x. E5 j- ?1 s; N* Rclear all; close all; clc;

% 1. Create mesh grid
' c5 D( z4 j; `7 w  U[X,Y] = meshgrid(linspace(-10,10,25));
Z = sinc(sqrt((X/pi).^2+(Y/pi).^2));
0 c! m3 k* ~" l3 Qfigure('units','normalized','Position',[0 0 1 1],'Color','w'); colordef('white');
" T/ M: i' {1 E0 I  V% P
% 2. Plot the original rectangular mesh
  ~' M3 Z; W# \1 w9 ^* m3 esubplot(2,2,1);
0 Z6 o' A( p: V  |" z7 m* B# xsuRF(X,Y,Z); hold on;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
7 f2 |1 W0 }) w8 X% o9 X. C( N5 mtitle('Meshgrid','FontSize',20);
  \5 l  x7 g: x( I& O" ]& c" l+ q
, c) a/ u; E. v1 ^%3.1 forward slash division of quadrilateral

+ w1 X# l9 r7 x5 _) p6 S[F,V]=mesh2tri(X,Y,Z,'f');
C=V(:,3); C=mean(C(F),2);
subplot(2,2,2);
$ S/ h4 G5 q, ]! ?patch('Faces',F,'Vertices',V,'FaceColor','flat','CData',C); hold on;
4 f; n! F+ j5 U: y, E) uaxis tight; axis square; grid on; axis off; view(3); view(-30,70);
) S* l6 O, _4 |# l: rtitle('Forward slash','FontSize',20);
5 D. l7 F0 Z$ [: F+ Y  ?, H
%%3.2  back slash division of quadrilateral
[F,V]=mesh2tri(X,Y,Z,'b');
C=V(:,3); C=mean(C(F),2);
subplot(2,2,3);
%Example of using original meshgrid coordinates instead
trisurf(F,X,Y,Z);
axis tight; axis square; grid on; axis off; axis off; view(3); view(-30,70);
title('Back slash','FontSize',20);
' D3 f- `6 x, i3 N9 F9 C: v* o
%%3.3  Cross division of quadrilateral
[F,V]=mesh2tri(X,Y,Z,'x');
% Replace Z-coordinates of added points by interpolated values if desired
IND=(numel(X)+1):size(V,1);
3 H6 k3 Y0 e+ ]2 p# N' A2 yZI = interp2(X,Y,Z,V(IND,1),V(IND,2),'cubic');
! [7 v  d/ F% y+ C* E8 UV(IND,3)=ZI;
C=V(:,3); C=mean(C(F),2);
subplot(2,2,4);
- q; s! b) P8 v4 Zpatch('Faces',F,'Vertices',V,'FaceColor','flat','CData',C); hold on;
1 Z% a$ \' P# Saxis tight; axis square; grid on; axis off; view(3); view(-30,70);
title('Crossed','FontSize',20);

结果如下：

6 ^' y/ ]  u  T+ [, |  |# p


: Y- J0 L  l9 }8 w2 `  `& S

NingW

MatlabTools(1)mesh2tri