递推算法 matlab,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)

结果为：

原文：http://blog.csdn.net/alop_daoyan/article/details/25053431