非均匀有理B样条.docVIP

非均匀有理B样条 （学习记录和上机练习） 非均匀有理B样条，通常简称为NURBS(Non-Uniform Rational B-Splines)。NURBS是非有理B样条、有理以及非有理Bezier曲线曲面的推广。 一、Bezier曲线 Bezier曲线 一条n 次Bezier曲线可以表示为 0≤≤1 （1） 其中，基函数（也称为混合函数）是著名的n次Bernstein多项式，其定义为 （2） 式中的几何系数称为控制点。 //计算参数u的n次基函数值并存在B[]中 void Bernstein(int n,double u,float B[]) { double u1; int j,k; B[0]=1.0; u1=1.0-u; for(j=1;j=n;j++) { float saved=0.0; for(k=0;kj;k++) { float temp=B[k]; B[k]=saved+u1*temp; saved=u*temp; } B[j]=saved; } } 有理Bezier曲线 一条n次有理Bezier曲线的定义为 0≤≤1 （3） 其中是标量，称为权因子。 图1、3次Bezier曲线 图2、3次有理Bezier曲线 图3、n=3 Bernstein图形 图4、n=9 Bernstein图形 图5、Bezier双三次曲面片 图6、有理Bezier双三次曲面片 二、B样条基函数 1、B样条基函数的定义和性质 令是一个实数序列，即≤。其中，称为节点，称为节点矢量，用表示第i个p次（p+1阶）B样条基函数，其定义为： （4） 由此可知： （1）是一个阶梯函数，它在半开区间外都为零； （2）当p0时，是两个p-1次基函数的线性组合； （3）计算一组基函数时需要事先指定节点矢量和次数p； （4）上式中可能出现0/0，我们规定0/0=0； （5）是定义在整个实数轴上的分段多项式函数，它在区间上有意义； （6）半开区间称为第i个节点区间，它的长度可以为零，因为相邻节点可以是相同的； （7）计算p次基函数的过程生成一个三角阵列。 性质： 1 局部支撑性，如果 ，则＝0。 2 非负性，对于所有的i,p和u，有≥0。 3 规范性，对于任意的节点区间 ，当 时 。 4 可微性，在节点区间内部，是无限次可微的（在每个节点区间内部，它是一个多项式），在节点处是p-k次连续可微的，其中k是节点的重复度。 2、B样条基函数的导数 1 基函数的求导公式为： （5） 对基函数 求导得到一般的求导公式： （6） 为了完整性，另一个计算B样条基函数各阶导数的公式（参考[Butt76]）: （7） //------------------------------------------------------------------------- //计算所有非零B样条基函数并返回其值 //i为参数u所在的节点区间下标 void BasisFunction(int i,int p,float u,float