计算几何Bézier曲面绘制MATLAB实现

问题描述

给定控制顶点

绘制Bézier曲面

图像绘制

%----------画线----------
xs=[0 1 2;
    0 1 2;
    0 1 2];
ys=[0 0 0;
    1 1 1;
    2 2 2];
zs=[1 1.5 1;
    1.5 1.7 1.5;
    1 1.5 1];
hold on
x=[xs(1,1),xs(1,2)];
y=[ys(1,1),ys(1,2)];
z=[zs(1,1),zs(1,2)];
for i=1:3
   for j=1:2
        x=[xs(i,j),xs(i,j+1)];
        y=[ys(i,j),ys(i,j+1)];
        z=[zs(i,j),zs(i,j+1)];
        plot3(x,y,z,'r','LineWidth',2.5);
        x=[xs(j,i),xs(j+1,i)];
        y=[ys(j,i),ys(j+1,i)];
        z=[zs(j,i),zs(j+1,i)];
        plot3(x,y,z,'r','LineWidth',2.5);
   end
end
%----------画点----------
xp=[0 1 2 0 1 2 0 1 2];
yp=[0 0 0 1 1 1 2 2 2];
zp=[1 1.5 1 1.5 1.7 1.5 1 1.5 1];
scatter3(xp,yp,zp,'b','filled')

%----------画面----------
U=0:0.05:1;V=0:0.05:1;    %这里可以通过改变步长来调整控制网格的疏密
[u,v]=meshgrid(U,V);
d=length(u);
Bu0=(1-u).^2;  Bu1=2.*u.*(1-u);  Bu2=u.^2;
Bv0=(1-v).^2;  Bv1=2.*v.*(1-v);  Bv2=v.^2;
X=zeros(d);Y=zeros(d);Z=zeros(d);
Bv=[Bv0;Bv1;Bv2];Bu=[Bu0;Bu1;Bu2];
for i=1:3
    for j=1:3
        nj=d*(j-1)+1;
        ni=d*(i-1)+1;
        X=X+xs(i,j).*Bu(nj:nj+(d-1),:).*Bv(ni:ni+(d-1),:);
        Y=Y+ys(i,j).*Bu(nj:nj+(d-1),:).*Bv(ni:ni+(d-1),:);
        Z=Z+zs(i,j).*Bu(nj:nj+(d-1),:).*Bv(ni:ni+(d-1),:);
    end
end
surf(X,Y,Z)
axis off

运行结果如下

（这里的动图是用ScreenToGif制作，十分简洁快速，强烈安利）

总结

这里的函数写的很傻……而且当控制顶点的数量改变，不是9个的时候，改动就十分的大，如果以后改良的话再写……