MatlabTools(1)mesh2tri

mutougeda

将结构的四边形网格转换为三角形网格。
; G( k- O* h$ ^5 c- }Demo:
clear all; close all; clc;

; _$ }7 E. M. H6 ?- _% 1. Create mesh grid
6 X7 V2 ~5 y6 S* e[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');
1 f$ b+ m  u% p7 B
% 2. Plot the original rectangular mesh
2 k" b/ k4 }' M- J  H" R" Fsubplot(2,2,1);
5 [$ D3 [( U- v9 z/ |; ^, s2 z& ^  bsuRF(X,Y,Z); hold on;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
+ w" K, D; a) W7 ztitle('Meshgrid','FontSize',20);
  R/ H% m7 e5 w$ J# x) H( a  J; F
3 y' o7 J9 N  Z, M%3.1 forward slash division of quadrilateral

) q* S- r% I5 s8 b[F,V]=mesh2tri(X,Y,Z,'f');
C=V(:,3); C=mean(C(F),2);
subplot(2,2,2);
patch('Faces',F,'Vertices',V,'FaceColor','flat','CData',C); hold on;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
title('Forward slash','FontSize',20);

1 Q6 L; A% A7 n4 J%%3.2  back slash division of quadrilateral
7 g$ P+ I6 [- l' @! X2 D[F,V]=mesh2tri(X,Y,Z,'b');
C=V(:,3); C=mean(C(F),2);
! j6 G. X: Z# @' M9 h# C- ^9 Tsubplot(2,2,3);
%Example of using original meshgrid coordinates instead
trisurf(F,X,Y,Z);
- P5 P5 x4 _. Yaxis 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
[F,V]=mesh2tri(X,Y,Z,'x');
% Replace Z-coordinates of added points by interpolated values if desired
, j( t+ A# |6 }2 K+ @3 `/ [4 FIND=(numel(X)+1):size(V,1);
0 V8 Q) N7 M. Z, V( ?& qZI = interp2(X,Y,Z,V(IND,1),V(IND,2),'cubic');
3 s0 Z8 j: Q. Y$ UV(IND,3)=ZI;
C=V(:,3); C=mean(C(F),2);
$ x; {0 J: |$ F; b2 U" \( N+ {subplot(2,2,4);
patch('Faces',F,'Vertices',V,'FaceColor','flat','CData',C); hold on;
( M* v# A2 j2 ^9 U  q& P5 E8 Paxis tight; axis square; grid on; axis off; view(3); view(-30,70);
title('Crossed','FontSize',20);

7 E  B$ d9 E' q, m2 [结果如下：

4 }2 ^* u$ G% E, n6 w- r

. u% W# n- u0 v- v3 Z/ F) r& o+ d
( q& A0 D9 U% t! W* Y9 q9 T" }
1 J) t; G/ c+ J

NingW

MatlabTools(1)mesh2tri