bsplinefit.m

采用德布尔算法进行B样条的生成

function [x,y]=BsplineFit(m,n,k,CPX,CPY)
%    [x,y]=Bspline(m,n,k,CPX,CPY)
%     m为型值点数目，比如m＝40
%     n＋2为控制点，比如n＝2
%     k－1为B样条曲线阶数，比如 k＝4
%     CPX,CCPY 控制点的横、纵坐标，为n+2阶行向量
%     x,y为型值点的坐标，均为m阶行向量

u=linspace(0,1,m);
for i=1:m    
    [tnu(:,i),dnu(:,i)]=nik(k,n,u(i));
end
x=CPX*tnu;
y=CPY*tnu;