一直在写c++程序,不过对于一些作图程序来说,还是MATLAB比较实在。
de Casteljau算法是作贝塞尔曲线的一种高效的算法,其思想就是对[0,1]区间中所有的t,通过n个控制顶点不断递推得到一个顶点:下面是我的代码实现:
function deCasteljau(P,Q)
%P is 1*n matrix for X
%Q is 1*n matrix for Y
m=size(P);
n=m(2);
x=zeros(1,101);
y=zeros(1,101);
p=zeros(n);
q=zeros(n);
plot(P,Q,'r');
hold on
for i=1:n
p(i,1)=P(i);
q(i,1)=Q(i);
end
k=1;
for t=0:0.01:1
for j=2:n
for i=j:n
p(i,j)=t*p(i,j-1)+(1-t)*p(i-1,j-1);
q(i,j)=t*q(i,j-1)+(1-t)*q(i-1,j-1);
end
end
x(k)=p(n,n);
y(k)=q(n,n);
k=k+1;
end
plot(x,y,'b-')
下面是测试数据
x=[1,2,3];
y=[3,1,6];
deCasteljau(x,y)
结果为: