一直在写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)
结果为:
de Casteljau算法的matlab实现
时间: 2024-10-25 09:07:05