ljd2.m

de Boor 求值算法求作二次和三次B-样条插值曲线

function u=LJD2(n,k,x,y,z)
%使用累加弦长法确定节点矢量
%%该算法取两端点节点重度r=k+1
u=zeros(1,n+2*k);
for i=1:n-1    
    L(i)=sqrt((x(i+1)-x(i))^2+(y(i+1)-y(i))^2+(z(i+1)-z(i))^2);
end
for i=1:k+1
    u(i)=0;
end
for i=k+2:n+k
    u(i)=u(i-1)+L(i-k-1);
end
for i=n+k+1:n+2*k
    u(i)=u(i-1);
end