jchzh.m

来自「de Boor 求值算法求作二次和三次B-样条插值曲线」· M 代码 · 共 22 行

function D=JCHZH(n,x,u)
%使用基于基的运算的插值法求控制顶点
%%该算法取两端点节点重度r=k+1
D=zeros(1,n+2);
DT=zeros(n+2,n+2);
for j=1:2
    DT(1,j)=(-1)^j;
end
    DT(n+2,n+2)=1;
for i=2:n+1
    for j=i-1:i
        DT(i,j)=Ji2(j-1,i,u);
    end
end

P=zeros(1,n+2);
P(1,1)=(u(4)-u(3))/2;
for i=2:n+2
    P(1,i)=x(i-1);
end

D=(inv(DT)*(P'))';