draw.m

malab语言实现bspline曲线拟合

function draw(m,p,Chrom,Q)

% draw函数将误差最小的个体bestchrom输出图象

t = Chrom(1:m+1);
n = Chrom(m+2);
u = knotpara(t,p,n,m); 
for i = 1:m-1
    w = basicfun(n,p,t(i+1),u);
    N(i,:) = w(1:n+1)';
end
tt = Q(2:m,:);
QQ = N*pinv(N)*tt;

clf;
axis([0 1 -0.3 0.3]);
hold on;
plot(Q(:,1),Q(:,2),'b-',Q(:,1),Q(:,2),'bo');
hold on;
plot(QQ(:,1),QQ(:,2),'r-',QQ(:,1),QQ(:,2),'r+');