MatlabTools(1)mesh2tri

mutougeda

将结构的四边形网格转换为三角形网格。
8 n& M8 @3 L( V$ V: A2 aDemo:
1 K3 c  {1 `1 A. Y1 }clear all; close all; clc;
( g) @. q& u! ^6 e$ w
% 1. Create mesh grid
' s. G  n$ I6 S! c/ w' _% X[X,Y] = meshgrid(linspace(-10,10,25));
5 ^6 U4 G) V2 N3 {: }$ E2 ?Z = sinc(sqrt((X/pi).^2+(Y/pi).^2));
figure('units','normalized','Position',[0 0 1 1],'Color','w'); colordef('white');
- [/ I6 K" C- M! x1 p
% 2. Plot the original rectangular mesh
# {8 R3 |; M, X! n0 g) Wsubplot(2,2,1);
/ `* a/ ~1 N$ z: q# W& MsuRF(X,Y,Z); hold on;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
title('Meshgrid','FontSize',20);

" }6 E5 Q# ]) a+ t' t3 {  W%3.1 forward slash division of quadrilateral

9 ^4 K" Y: I, I[F,V]=mesh2tri(X,Y,Z,'f');
C=V(:,3); C=mean(C(F),2);
subplot(2,2,2);
6 T0 n! O  }7 X7 ]3 t8 lpatch('Faces',F,'Vertices',V,'FaceColor','flat','CData',C); hold on;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
+ {( J: a9 p: }- p9 p. {title('Forward slash','FontSize',20);
! S1 ]8 P8 w$ P6 ]9 l6 I
%%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);
$ W3 c8 C. M3 d8 Q. Y; taxis tight; axis square; grid on; axis off; axis off; view(3); view(-30,70);
- x1 Q5 P. k. Z; q. P4 T( S& Vtitle('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
) Q4 `4 i. N7 v% N7 F% M3 kIND=(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;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
3 i  U9 m0 G  V$ P0 [title('Crossed','FontSize',20);

结果如下：

% P6 D3 y! z: n" O- c2 Z

NingW

MatlabTools(1)mesh2tri