function u=LJD(n,k,x,y,z)
%使用累加弦長法確定節(jié)點(diǎn)矢量
%該算法取兩端點(diǎn)節(jié)點(diǎn)重度r=k，即首末控制點(diǎn)就是首末數(shù)據(jù)點(diǎn)
u=zeros(1,n+2*k-2);
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
    u(i)=0;
end
for i=k+1:n+k-1
    u(i)=u(i-1)+L(i-k);
end
for i=n+k:n+2*k-2
    u(i)=u(i-1);
end