EDA365电子论坛网

标题: MatlabTools(1)mesh2tri [打印本页]

---

作者: mutougeda    时间: 2020-8-7 13:24
标题: MatlabTools(1)mesh2tri
将结构的四边形网格转换为三角形网格。
Demo:
clear all; close all; clc;

6 E1 L5 X: y4 E5 ~3 e; q: _# T% 1. Create mesh grid
[X,Y] = meshgrid(linspace(-10,10,25));
8 d, i" l/ [9 R9 oZ = sinc(sqrt((X/pi).^2+(Y/pi).^2));
figure('units','normalized','Position',[0 0 1 1],'Color','w'); colordef('white');
; E. O) w! B6 ]3 D
% 2. Plot the original rectangular mesh
subplot(2,2,1);
& {: N* _( Q- ]4 m; s8 @9 N9 @surf(X,Y,Z); hold on;
axis tight; axis square; grid on; axis off; view(3); view(-30,70);
title('Meshgrid','FontSize',20);

%3.1 forward slash division of quadrilateral

[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);
+ h9 T- w; s0 B8 u6 ptitle('Forward slash','FontSize',20);

) \4 m$ ~/ C4 u* _' F%%3.2  back slash division of quadrilateral
! j" V9 P; ^" q' V1 f[F,V]=mesh2tri(X,Y,Z,'b');
7 q3 o: T# v+ u6 I. I; U3 E* uC=V(:,3); C=mean(C(F),2);
: C) x- l0 g# H2 t8 E* e# m* xsubplot(2,2,3);
) E8 Y: K- t5 }( R%Example of using original meshgrid coordinates instead
trisurf(F,X,Y,Z);
' w$ N0 @9 z0 Y3 B/ i9 L' V/ z9 @axis tight; axis square; grid on; axis off; axis off; view(3); view(-30,70);
, a, D# D; x- t* ~9 H) otitle('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
$ s- b6 K; Z% A: P) S) M+ T; AIND=(numel(X)+1):size(V,1);
) G7 q5 ~! d1 M2 ~8 GZI = interp2(X,Y,Z,V(IND,1),V(IND,2),'cubic');
+ D% j% w+ s. Z/ Q" ~5 X+ MV(IND,3)=ZI;
C=V(:,3); C=mean(C(F),2);
$ N$ Z, u5 H2 N0 f0 gsubplot(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);
title('Crossed','FontSize',20);
/ w0 J  S2 Z: V/ U3 J  G
结果如下：


/ i# M, j' E. }3 d

---

作者: NingW    时间: 2020-8-7 14:06
MatlabTools(1)mesh2tri

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2