deboor.m

来自「关于样条插值均匀与非均匀的示例」· M 代码 · 共 12 行

%************递归,德布尔递推算法***************************************
function y=DeBoor(d,u,kn,i,L)
%控制点序列,u节点序列,kn插入的节点,u(i)<kn<u(i+1),L递推层次
% k B样条次数
k=3;
if L==0
     y=d(i,:);
     return;
end

temp=(kn-u(i))/(u(i+k+1-L)-u(i));
y=(1-temp)*DeBoor(d,u,kn,i-1,L-1)+temp*DeBoor(d,u,kn,i,L-1);