1个回答
2010-12-19
展开全部
close all
[X1,Y1]=meshgrid(-1:0.05:2,-1:0.05:2);
X=[ linspace(0,1,10).' linspace(0,1,10).'
linspace(1,2,10).' 1+zeros(10,1)
2+zeros(10,1) linspace(1,2,10).'
linspace(2,-1,30).' 2+zeros(30,1)
-1+zeros(30,1) linspace(2,-1,30).'
linspace(-1,2,30).' -1+zeros(30,1)
2+zeros(10,1) linspace(-1,0,10).'
linspace(2,0,20).' zeros(20,1)
X1(:) Y1(:)];
C=[1:150;[2:150 1]]';
dt=DelaunayTri(X,C);
IO = inOutStatus(dt);
trisurf(dt(IO, :), dt.X(:,1), dt.X(:,2), zeros(size(dt.X(:,1))) )
hold on
trisurf(dt(~IO, :), dt.X(:,1), dt.X(:,2), exp(-dt.X(:,1)+dt.X(:,2)) )
hold off
shading interp
[X1,Y1]=meshgrid(-1:0.05:2,-1:0.05:2);
X=[ linspace(0,1,10).' linspace(0,1,10).'
linspace(1,2,10).' 1+zeros(10,1)
2+zeros(10,1) linspace(1,2,10).'
linspace(2,-1,30).' 2+zeros(30,1)
-1+zeros(30,1) linspace(2,-1,30).'
linspace(-1,2,30).' -1+zeros(30,1)
2+zeros(10,1) linspace(-1,0,10).'
linspace(2,0,20).' zeros(20,1)
X1(:) Y1(:)];
C=[1:150;[2:150 1]]';
dt=DelaunayTri(X,C);
IO = inOutStatus(dt);
trisurf(dt(IO, :), dt.X(:,1), dt.X(:,2), zeros(size(dt.X(:,1))) )
hold on
trisurf(dt(~IO, :), dt.X(:,1), dt.X(:,2), exp(-dt.X(:,1)+dt.X(:,2)) )
hold off
shading interp
参考资料: ilovematlab wanbaocheng
本回答被提问者采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询