nq.m

matlab 闭曲线B样条的拟合

clear u1;
clear u2;
clear q1;
clear l;
clear i1;
clear i;
clear qq;
clear t;
clear l;
clear u4;
clear z;
clear q1;
clear NN;
load ah.txt
u1=ah(:,1);
u2=ah(:,2);
global l;
l=length(u1)-1;
for i=0:l
    if i>5
        for j=(i-5):i
            NN(i+1,j+1)=ni3(i+1,u1,j,l);
        end
    else
        for j=(i-5)+l:l
            NN(i+1,j+1)=ni3(i+1,u1,j,l);
        end
        for j=0:i
            NN(i+1,j+1)=ni3(i+1,u1,j,l);
        end
    end
end
P=NN\u2;
%PZ=[PZ;P'];