de Casteljau算法的matlab实现

          一直在写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)

结果为：