shuju.m

采用德布尔算法进行B样条的生成

%m为数据点个数减1，n为控制点个数减1，k为曲线次数，dim为维数
dim=3;

fid=fopen('H.dat');
q=fscanf(fid,'%g',[3,inf]);%%%[2,inf]表示返回的矩阵A的形状

fclose(fid);



m=length(q(1,:))-1;
for i=1:m+1
    
  plot3(q(1,i),q(2,i),q(3,i),'rd');
  hold on

end

k=3;
n=m+k-1;
[u,l,t,r]=fqxjdsl(k,m,q,dim,n);

d=kqx(n,u,q,m);
%plot(d(1,:),d(2,:));
m=1;
for u=0:0.001:1;
     p=doboor(d,u,l,t,r,k,dim,n);
     x(m)=p(1);
     y(m)=p(2);
     z(m)=p(3);
     m=m+1;
    
end
plot3(x,y,z,'r');
n=length(x);


fid=fopen('H_chazhi.dat','wt');
for i=1:n
fprintf(fid,'%9.6f  %9.6f  %1.0f',x(i),y(i),z(i));
fprintf(fid,'\n');
end
fclose(fid);