MatlabTools(1)mesh2tri

mutougeda

将结构的四边形网格转换为三角形网格。
Demo:
7 c# _% ]/ ]3 E% @clear all; close all; clc;
$ d6 E# t) {( ~/ E' g
0 N4 Y* b( M8 |, V% 1. Create mesh grid
[X,Y] = meshgrid(linspace(-10,10,25));
Z = sinc(sqrt((X/pi).^2+(Y/pi).^2));
figure('units','normalized','Position',[0 0 1 1],'Color','w'); colordef('white');
: _% o6 v" m. N
6 p; @' A$ P7 `0 ?% X% g/ u, [2 i8 K% 2. Plot the original rectangular mesh
. `. x& M% k% N0 p, Ksubplot(2,2,1);
suRF(X,Y,Z); hold on;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
  p3 s) G9 M3 Q7 U) C+ Ftitle('Meshgrid','FontSize',20);

%3.1 forward slash division of quadrilateral
5 S3 b: g! \- P4 d
[F,V]=mesh2tri(X,Y,Z,'f');
C=V(:,3); C=mean(C(F),2);
subplot(2,2,2);
" d0 S3 |( v' h: h, l& opatch('Faces',F,'Vertices',V,'FaceColor','flat','CData',C); hold on;
9 a, m. `" S: B0 M- vaxis tight; axis square; grid on; axis off; view(3); view(-30,70);
8 r8 Y% T/ B; q0 B3 Q1 u% Ptitle('Forward slash','FontSize',20);

%%3.2  back slash division of quadrilateral
[F,V]=mesh2tri(X,Y,Z,'b');
2 D  x$ j! |( e0 H: [, A2 xC=V(:,3); C=mean(C(F),2);
, E8 h; Y% ~0 X  Dsubplot(2,2,3);
5 W* I4 f3 x- \) t6 C( i- {%Example of using original meshgrid coordinates instead
& E+ v/ v9 s5 E9 m# N. k, j; itrisurf(F,X,Y,Z);
. W* S* H* u9 }; T3 V/ Naxis tight; axis square; grid on; axis off; axis off; view(3); view(-30,70);
title('Back slash','FontSize',20);

%%3.3  Cross division of quadrilateral
( m6 `% f& J1 L" H( `1 X. S# L# X[F,V]=mesh2tri(X,Y,Z,'x');
$ \$ H+ |% n6 B6 W6 m$ V% Replace Z-coordinates of added points by interpolated values if desired
IND=(numel(X)+1):size(V,1);
ZI = interp2(X,Y,Z,V(IND,1),V(IND,2),'cubic');
V(IND,3)=ZI;
C=V(:,3); C=mean(C(F),2);
subplot(2,2,4);
patch('Faces',F,'Vertices',V,'FaceColor','flat','CData',C); hold on;
) S# ?+ A* ^6 e# J8 J8 Caxis tight; axis square; grid on; axis off; view(3); view(-30,70);
7 u' Z. O* T, I! Z+ o% v5 f1 dtitle('Crossed','FontSize',20);

+ ^, C  v. N+ [3 T, d; a8 Q( T. C结果如下：
* D5 v: c9 e; u! r+ H3 ?  S

7 ]- Y7 S: G5 h8 g( Y" N


, n$ L; ~% w2 q, X  j. B5 ~2 t

NingW

MatlabTools(1)mesh2tri