aline.m

matlab 闭曲线B样条的拟合

l=length(u1);
clear qq;clear u4;
global q1;
q1=2;
close all;
qq=zeros(l*q1,l);
for i=1:l
    if i>5
        for j=(i-5):i
            for t=0:q1-1
                i1=mod(i,l)+1;
               u4((i-1)*q1+t+1)=u1(i)+t*(u1(i1)-u1(i))/q1;
            end
        end
    else
        for j=(i-5)+l:l
            for t=0:q1-1
               i1=mod(i,l)+1;
               u4((i-1)*q1+t+1)=u1(i)+t*(u1(i1)-u1(i))/q1;
            end
        end
        for j=1:i
            for t=0:q1-1
              
               i1=mod(i,l)+1;
               u4((i-1)*q1+t+1)=u1(i)+t*(u1(i1)-u1(i))/q1;
            end
        end
    end
end
for i=1:l
    if i>5
        for j=(i-5):i
            for t=0:q1-1
               qq((i-1)*q1+t+1,j)=Ni(i,j,t,u1,u4);
               i1=mod(i,l)+1;
               
            end
        end
    else
        for j=(i-5)+l:l
            for t=0:q1-1
               qq((i-1)*q1+t+1,j)=Ni(i,j,t,u1,u4);
               i1=mod(i,l)+1;
               
            end
        end
        for j=1:i
            for t=0:q1-1
              qq((i-1)*q1+t+1,j)=Ni(i,j,t,u1,u4); 
               i1=mod(i,l)+1;
               
            end
        end
    end
end
% qq=qq';
u4=u4';
z=qq*P;
plot(u1,u2,'*',u4',z)