Matlab问答day10

1.Matlab的delaunay三角剖分函数

clear
clc
%假设已知数据：
x = rand(5,1);
y = rand(5,1);
TRI = delaunay(x',y');
triplot(TRI,x,y);%绘图

即连续平面上的点集，使之完全由三角形构成，且满足一定条件（互不相交、包含）等；

clear
clc

%假设已知数据：
x = rand(5,1);
y = rand(5,1);
z = rand(5,1);
TRI = delaunay(x,y,z);
trisurf(TRI,x,y,z);%绘图

三维也可以完成，只是全部图形由四面体构成。还有trimesh函数等。

参考：Matlab论坛

2. 画三维隐函数的一个数值方法：

Matlab画3维图形的ezplot函数先前提过，但那是用参数方程的方式画图。

对一个隐函数f(x,y,z)来讲，其实它已经到四维空间了，其实x,y,z三点确定一个f。

当我回答论坛一个问题:sin(z)+z+x^2+sin(x)+y+sin(y)=1时，我考虑的确是这个隐函数在f(x,y,z)=1处的“切体”（对等于切面）

因此，我给出的方法是：

clear
clc
syms z
counter =1;
f = @(x,y,z)z+x^2+sin(x-y)-y+sin(y)-1;%给出隐函数
for x = 0:0.3:pi%指定x,y范围
    for y = 0:0.3:pi
        temp(counter) = solve(f(x,y,z)==0,z);%按已知求解方程
        counter = counter + 1;
    end
end
x = 0:0.3:pi;
y = 0:0.3:pi;
temp2 = double(temp);%数值化数据
[X,Y] = meshgrid(x,y);
temp3 = reshape(temp2,size(X));
surf(X,Y,temp3)

可是这个方法十分的不稳定，理论也不太成熟。

还好论坛上的大神们解决了这个问题：

点击打开链接：from Matlab论坛

用到了isosurface， patch， isonormals三个函数。