问题描述
给定控制顶点
绘制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个的时候,改动就十分的大,如果以后改良的话再写……