B样条曲线与曲面

1、四、B样条曲线与曲面Bezier曲线具有很多优越性，但有二点不足：1) 特征多边形顶点数决定了它的阶次数，当n较大时，不仅计算量增大，稳定性降低，且控制顶点对曲线的形状控制减弱；2) 不具有局部性，即修改一控制点对曲线产生全局性影响。1972年Gordon等用B样条基代替Bernstein基函数，从而改进上述缺点。B样条曲线的数学表达式为：Rn(U)PikWn(U)在上式中，0VuV1;i=0,1,2,-,m所以可以看出：B样条曲线是分段定义的。如果给定m+n+1个顶点Pi(i=0,1,2,-,m+n)，则可定义m+1段n次的参数曲线。在以上表达式中：M,n(u)为n次BW条基函数，也称B样条

2、分段混合函数。其表达式为：1nNk,n(U)n!(1)jCnj1(Unkj)n式中：0Vu<1k=0,1,2,-,n1.均匀B样条曲线1一次均匀B样条曲线的矩阵表示空间n+1个顶点P(i=0,1,，n)定义n段一次(k=0,1,n=1)均匀B样条曲线，即每相邻两个点可构造一曲线段Pi(u),其定义表达为：形。P(u)u(1u)1Pi0PiPi1+uPiN0,1(u)Pi1+第i段曲线端点位置矢量:Ni,1,.,n;0u(u)PiP(0)Pi1,P(1)Pi,且一次均匀B样条曲线就是控制多边空间n+1个顶点的位置矢量pi(i=0,1,n)定义n-1段二次(k=0,1,2,n=2)均匀B样条

3、曲线，每相邻三个点可构造-一曲线段Pi(u)(i=1,-,n1),其定义表达为：P(u)1u2u1121Pi1220Pii1,.,n1;0u12110Pi1111=2!(12u+u2)Pi1+2!(1+2u-2u2)P+2!u2Pi+12二次均匀B样条曲线的=N0,2(U)Pi1+N1,2(u)Pi+N2,2(U)R+1囹3-11二次B样条曲线端点位置矢重：P(0)0.5(Pi1Pi)P爵归")，即曲线的起点和终点分别位于控制多边形Pi-1Pi和PiPi+1的中点。化成Pi1PiPi1直线边上的一段直线。若Pi1、Pi、Pi1三个顶点位于同一条直线上，Pi(u)蜕端点一阶导娄攵矢量：

4、P(0)PiPi1Pi(1)Pi1PiP(0)Pi1Pip(1)Pi2R1,即曲线的起点切矢和终点切矢分别和二边重合，且相邻两曲线段在节点处具有一阶导数连续。二阶导数矢量：Pi(0)Pi12巳Pi1PiPi，即曲线段内任何点处二阶导数相等，且相邻两曲线段在节点处二阶导数不连续。三次均匀B样条曲线/端点位置矢量：P(0)16(P'14Pi巳1)P*6(P，4P'1Pi2),即起点位于三角形Pi-1PiPi+1中线PiM1的1/3处，终点位于三角形线的端点并不通过控制点。PiPi+iPi+2中线Pi+1M2的1/3处。可见B样条曲端点一阶导数矢量：Pi(0)(Pi1Pi1)/2,P

5、(1)(Pi2Pi)/2Pi1(0),即曲线起空间n+1个顶点的位置矢量P(i=0,1,ooo,n)构造n2段三次(k=0,1,2,3，四阶n=3)均匀B样条曲线段，每相邻四个点可定义-曲线段Pi(u)(i=1,。,n-2),其定义表达为：1Pi(u)1u3u2u13631=3!(1u)3Pi1+3!3 31Pi1630Pii1,.,n2;0u1030Pi14 10Pi2(46u2+3u3)Pi+3(1+3u+3u23u3)Pi+1+3!u3Pi+2=N0,3(u)Pi1+N13(u)Pi+N2,3(u)Pi+1+N3,3(u)Pi+2点的切矢平行于Pi-1PiPi+1的底边Pi-1Pi+1，

6、其模长为底边Pi-1Pi+1长的1/2,同样曲线终点的切矢平行于PiPi+1Pi+2的底边PiPi+2,其模长也为底边PiPi+2长的1/2。且相邻两曲线段具有一阶导数连续二阶导数矢量：Pi(°)Pi12PiPi1,Pi(1)Pi2Pi1Pi2Pi1(0),即曲线段在端点处的二阶导数矢量等于相邻两直线边所形成的平行四边形的对角线，且两曲线段在节点处具有二阶导数连续(因RP(°)。若R1、P、R1三个顶点位于同一条直线上，三次均匀B样条曲线将产生拐点；若Pi1、P、Pi1、Pi2四点共线，则Pi(U)变成一段直线；若Pi1、P、Pi1三点重合，则Pi(U)过Pi点。三顶点共线

7、(c)二直节点和三更节点图.L旋三决琳茶曲线的一些特例思考：用作图法绘制下图均匀三次B样条曲线。Bezier曲线段与段之间B样条曲线段与段之间具有天然的连续性，具有整体的光滑特性，而必须光滑拼接。因此在商用系统中B样条方法应用更为广泛。2. B样条曲线的性质1 局部性空间n+1个控制顶点Pi（i=0,1,，n）构造（n-k+1）段k次（k+1阶）B样条曲线段，且每一曲线段Pi（u）（i=，nk+1）由Pi1、Pi、Pik1等k+1个控制顶点确定，与其它控制点无关。2 整体性和连续性一般情况下（即无重节点、重顶点），n+1个控制顶点所构造的（nk+1）段k次（k+1阶）B样条曲线段组成一完整的B

8、样条曲线，曲线段与段之间具有Ck1阶函数连续性（或Gk1阶几何连续性），当有K重顶点时，将可能产生尖点（前面已介绍），虽然仍满足函数连续，但不满足几何连续。4 几何不变性改变坐标系不改变曲线形状。5 变差缩减性与Bezier曲线性质相同。（5）造型的灵活性由于其良好的局部特性，可以方便构造低次的复杂曲线，且编辑顶点对曲线形状的改变是局部的；由于其整体性和连续性，曲线具有整体的光滑性。正因如此，B样条曲线比Bezier应用更为广泛，为商用系统普遍采用。缺点：首末两端点不通过控制顶点，与其优点比较微不足道。3. 均匀双二次B样条曲面已知曲面的控制点Pij（i,j0,1,2）,参数u，w,且U，W0

9、,1,kl2,构造步骤是：a、沿w向构造均匀二次B样条曲线，即有：P0(w)Ww1-2P0P0P0WMP0P0P0经转置后:P0(w)P00PoiP02MIWT同上可得:Pi(w)P10P11P12MBWTP2(w)P20P21P22MBWT,0b、再沿u向构造均匀二次S(u,w)UMBPo(w)P00P01P02P1(w)UMbRoR1p12mBwtP2(w)P20P21P22B样条曲线，即可得到均匀二次B样条曲面:4. 简记为：S(u,w)uMbPmBwt均匀双三次B样条曲面已知曲面的控制点Pij(l,j0,1,2,3),参数u，w,且u，w01,kl3,构造双三次B样条曲面的步骤同上述。

10、a、沿w向构造均匀三次B样条曲线，有：P°(w)P00P01P02P03MTWTP1(w)P10P11P12P13MBWTP2(w)P20P21P22P23MBwTP3(w)P30P31P32P33mBwtb、再沿u向构造均匀三次B样条曲线，此时可认为顶点沿滑动，每组顶点对应相同的，当值由0到1连续变化，即形成均匀双三次B样条曲面。此时表达式为：Po(w)P1(w)tTS(u,w)UMBUMbPMBw1P2(w)P3(w)133声三次E律*it询片P00P01P02P031331P10P11P12P1313630MBP20P21P22P2363030P30P31P32P331410上

11、式也可表达为:S（U,w）=N0,3（u）N3,3（w）TNi,3（u）N0,3（u）N0,3（u）Pij4x4N0,3（w）N1,3（w）N2,3（w）Pi,jP,j1Pi,j2P,j3N°,3（w）Pi1,jPi1,j1P1,j2Pi1,j3N,3（w）Pi2,jP2,j1Pi2,j2Pi2,j3N2,3（w）P3,jP3,j1Pi3,j2Pi3,j3N3,3（w）B样条曲面其定义如下:即任意单张均匀双三次B样条曲面片Si,j(u,w)是由Pk,i(k=i,.,i+3,l=j,-,j+3）6个控制点定义而成。B样条曲面具有B样条曲线的多种性质，曲面片与片之间具有天然的连续性。仍以

12、均匀双三次曲面为例的说明B样条曲面的性质。(1)均匀双三次B样条曲面的顶点不经过任何特征网格顶点，且仅与各角点对应的9个特征网格顶点有关，如：Sij（0,0）=1/36（Pi,j+Pi,j+2+Pi+2,j+Pi+2,j+2）+1/9（Pi,j+1+Pi+1,j+Pi+2,j+1+Pi+1,j+2）+4/9Pi+i,j+i，同理可得Sj（0,1）、Sij（1,0）、Sij（1,1）。(2)均匀双三次B样条曲面的边界曲线仍为B样条曲线,该边界B样条曲线由对应的三条边对于由控制点Pij(l0,1,'",m,j01."n)组成的均匀双三次Si,j（u,w）N°,

13、3（u）N1,3（u）N2,3（u）N3,3（u）界特征网格顶点确定，由B样条曲面得定义可得:Si,j（u,0）N°,3（u）N1,3（U）N2,3（U）N3,3（U）Pi,jP，di,j1P,j2Pi,j31/6P1,jP1,j1Pi1,j2Pi1,j32/3P2,jPi2,j1Pi2,j2P2,j31/6P3,jPi3,j1Pi3,j2Pi3,j30N°,3（u）N1,3（U）N2,3（U）N3,3（U）1/6P,j2/3R,j11/6Pi,j21/6Pi1,j2/3R1,j11/6Pi1,j21/6Pi2,j2/3R2,j11/6Pi2,j21/6Pi3,j2/3P3

14、,j11/6Pi3,j2同理可得Sij(u,1）、Sj（0,w）、Sij（1,w）。推广之:沿B样条曲面任何等参数的截线均为一B样条曲线(读者证明)。(3)均匀双三次B羊条曲面边界的跨界一阶切矢只与定义该边界的顶点及相邻二排顶点(共三排顶点)有关，Sui,ji（u,1）No,3（u）Ni,3（u）N2,3（u）N3,3（u）Pi,jiP,jPi,jiP,j20Pii,jiPii,jPi,jiPii,j2i/6Pi2,jiP2,jPi2,jiPi2,j22/3P3,jiP3,jPi3,jiP3,j2i/6No,3（u）Ni,3（u）N2,3（u）N3,3（u）Sui,j（u,0）No,3（u）N

15、i,3（u）N2,3（u）N3,3（u）No,3（u）Ni,3（u）N2,3（u）N3,3（u）i/6Pi,j2/3P,jii/6Pi,j2i/6Pi,j2/3Pii,jii/6Pi,j2i/6Pi2,j2/3Pi2,jii/6P2,j2i/6Pi3,j2/3P3,jii/6P3,j2Pi,jP,jiPi,j2P,j3i/6Pii,jPi,jiPii,j2Pii,j32/3Pi2,jPi2,jiPi2,j2P2,j3i/6P3,jPi3,jiPi3,j2P3,j30i/6Pi,j2/3Pi/6Pi,j2i/6Pi,j2/3Pii,jii/6Pi,j2i/6Pi2,j2/3Pi2,jii/6P2

16、,j2i/6Pi3,j2/3P3,ji"6P3,j2Suii,j（u,1）依次可得Sui,j（u,i）,Swi,j（0,w）,Swi，j(i,w)。可见均匀三次B羊条曲面具有一阶函数连续性。同理可得Sui,j（u,0）,Sui,j（u,i）,Swij(0,w),Swi,j(i,w),其跨界二阶导矢只与定义该边界的及相邻两排顶点(共三排顶点)有关，且均匀三次B样条曲面具有二阶函数连续性。(4) 几何不变性。(5) 对称性。(6) 凸包性。B样条曲面的线框图绘制方法：先沿等参数方向离散成网格点，然后依次连线绘制。Bezier方法更灵活，因此应用更广由此可见，B样条方法能够很方便绘制复杂曲

17、面，显然比泛。5. B样条曲线与曲面的递推表达1)B样条曲线的定义定义：由前面的内容得知，三次均匀B样条曲线的基函数为：N0,3(u)=3!(1u)3Ni，3(u)=3！(46u2+3u3)N2，3(u)=3!(1+3u+3u23u3)N3,3(u)=3!u3上述基函数图形如下图所示:已知n+1个控制点P(i=0,1,ooo,n),也称为特征多边形的顶点,K次(k+1阶)B样条曲线的表达式是:C(u)nPiNi,k(u)i0k<=n其中Ni,k(u)是调和函数,也称为基函数，按照递归公式可定义为:Ni,°(u)若tiu其它tiNi,k(u)(uti)Ni,k1(u)tikti(

18、tik1u)Ni1,k1(u)tik1ti1式中L是节点值，且为非减序列,0/00Tt0,t1,.,tnk1构成了K次(k+1阶)B样条基函数的节点矢量，每一基函数由对应的K+2个节点决定;上式也表明，高次B样条函数可用低次的B样条函数来表示，由此式可得其递推计算方法。由基函数的示意图可知B样条基函数具有局部支撑特性，即Bi,k(u)xti,tixti,ti节点矢量所含节点数目由控制顶点Pi(i=0,1,ooo,n)和曲线次数k所确定(节点数=n+k+2),显然，基函数个数=控制点数。当节点沿参数轴是均匀等距分布的，则表示均匀B样条函数，其节点值ti-ti-1=常数；例如：当k=3,ti-ti

19、-1=1,u0,1,则上述基函数可表示为均匀三次B样条函数，并通过变换可得如下表达：Ni,3(u)=3!(1u)Ni+1,3(u)=3!(46u2+3u3)Ni+2,3(u)=3!(1+3u+3u23u3)Ni+3,3(u)=3!u3当节点沿参数轴的分布是不等距的，则表示非均匀B样条函数，即节点值tti-1丰常数。均匀B样条和非均匀B样条曲线一般不通过控制多边形首末两点。若需B样条曲线具有较好的端点性质(即通过端点)，实际应用中常引入准均匀B样条，即在节点矢量中两端节点具有k+1个重复度：t1=tk,tn+1=tn+2=tn+k+1这样构造的准均匀B样条曲线将通过控制多边形首末两点。下图为例：

20、构造n=5,k=2的准均匀B样条曲线的节点矢量为u=0,0,0,1,2,3,4,4,4。03.1.24准均也三次E样条曲线当门=7,k=3的准均匀三次B样条曲线的节点矢量可定义为u=0,0,0,0,1,2,3,4,5,5,5,5。若n=3,k=3的节点矢量u=0,0,0,0,1,1,1,1,此时三次B样条曲线转化为三次Bezier曲线。推而广之，若n=k,节点矢量为k1k1，此时K次B样条曲线为K次Bezier曲线。性质：局部性：K次B样条曲线仅在K+1个区间内非0。换句话说，每段k次B样条曲线只涉及k+1个基函数，并由k+1个顶点所定义。第i段K次B样条曲线仅由Pi,Pi+1,000,Pi+k共k+1个顶点所控制，与其它点无关；反之，修改一个控制顶点，其影响范围为与该顶点有关的k+1段。凸组合性质(凸包性)连续性：在无重节点的情况下，K次B样条在节点处具有k-1次连续性，如三次B样条具有二阶连续；若在某节点处具有m重节点，则K