datapara.m

malab语言实现bspline曲线拟合

function t=datapara(Q,m,e)

% 插值数据Q:1-m+1的参数化( chord length)

t(1) = 0;
d = 0;
for i = 2:m+1
    w = Q(i,:) - Q(i-1,:);
    s(i) = (norm(w,2))^e;
    d = d + s(i);
end
for i = 2:m
    t(i) = t(i-1) + s(i)/d;
end
t(m+1) = 1;