NURBS曲线的曲率计算

NURBS曲线的曲率计算

这两天做课题需要用到NURBS曲线，仔细探究了NURBS曲线的计算和求导等，主要看的书籍是《The NURBS book》，代码主要参考nurbs工具箱-M语言，下面将计算曲率的相关代码贴出并加以注释。主要算法参考《The NURBS book》P127,ALGORITHM A4.2

脚本函数，用半径为1的圆测试

ut = linspace(0.0,1.0,50);

crv = nrbcirc(); 
crv.order
crv.coefs
crv.knots

[p, w] = nrbeval(crv,ut);
p = p./repmat(w,[3,1])
tem = [];
for i=1:50
    tem = [tem;norm(p(:,i))];
end
plot(p(1,:),p(2,:))
 axis equal
 grid