MatlabTools(1)mesh2tri

mutougeda · 发表于 2020-8-7 13:24

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

将结构的四边形网格转换为三角形网格。
Demo:
clear all; close all; clc;

% 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');

% 2. Plot the original rectangular mesh
subplot(2,2,1);
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);
title('Forward slash','FontSize',20);

%%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);

%%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);
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);
title('Crossed','FontSize',20);

结果如下：

NingW · 发表于 2020-8-7 14:06

MatlabTools(1)mesh2tri

帐号		自动登录	找回密码
密码			注册

MatlabTools(1)mesh2tri

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1