fqxjdsl.m

来自「采用德布尔算法进行B样条的生成」· M 代码 · 共 26 行

function [u,l,t,r]=fqxjdsl(k,m,q,dim,n)
%定义定义域内节点值，开曲线闭曲线通用
uu(1)=0;
for i=2:m+1
    for j=1:dim
        a(j)=(q(j,i)-q(j,i-1))^2;
    end
    l=sqrt(sum(a(1:j)));
uu(i)=uu(i-1)+l;
end

for i=1:m+1
    uu(i)=uu(i)/uu(m+1);
    u(k+i)=uu(i);
end



    u(1:k+1)=0;
    u(n+2:n+k+2)=1;
    t(1:m+1)=u(k+1:n+2);
    r(2:m)=1;
    r(1)=k+1;
    r(m+1)=k+1;
    l=m+1;