function CK=CurveDerivsAlg2(n,p,U,P,u,d,CK)
% 计算B样条曲线上的点及其直到d阶导矢量;
% input: n,p,U,P,d,u,d
% output: CK (可认为是A(u)或者w(u),如果P用Pw取代的话！)
% clc;
% clear;
% p=3;
% u=2.5;
% P=[1 2 3 4 5 6 7];
% U=[0 0 0 1 2 3 4 4 5 5 5];
% m=size(U,2)-1;
% n=m-p-1;
% d=1;
%调试程序段1；
% U=[0 0 0 0 0.25 0.5 0.75 1 1 1 1];
% w=[1 0.6 0.85 1 0.85 0.6 1];
% Px=[0 -200 -200 0 200 200 0];
% Py=[0 -200 200 0 -200 200 0];
% p=3;
% n=6;
%
%
CK(:)=0.0;
N1(1:p+1,1:p+1)=0.0;
PK(:,:)=0.0;
du=min(d,p);
for k=0:1:d %%%%%
CK(k+1)=0.0;
end
span=FindSpan1(n,p,u,U);
N1=AllBasisFuns(span,u,p,U,N1);
PK=CurveDerivCpts(n,p,U,P,du,span-p,span,PK);%PK(k,i)即k阶导曲线的
%第i个控制点;
for k=0:1:du
CK(k+1)=0.0;
for j=0:1:p-k
CK(k+1)=CK(k+1)+N1(j+1,p-k+1)*PK(k+1,j+1);
end
end
nurbs.zip （20个子文件） 
AllBernstein.m 329B CurveDerivCpts.m 484B CurveDerivsAlg2.asv 873B AllBernstein2.asv 435B CurveDerivsAlg2.m 889B Bernstein.m 339B CurvePoint.asv 332B CurveDerivsAlg1.asv 491B AllBernstein2.m 442B AllBernstein.asv 304B CurveDerivsAlg1.m 545B CurvePoint.m 372B AllBasisFuns.asv 649B AllBasisFunsnew.asv 579B AllBasisFunsnew1.asv 368B BasisFuns.m 564B Curvature.m 586B AllBasisFuns.m 458B Bernstein.asv 343B CurveDerivCpts.asv 465B资源评论
