B样条曲线.ppt

《B样条曲线.ppt》由会员分享，可在线阅读，更多相关《B样条曲线.ppt（57页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期二、二、B B样条曲线样条曲线. .从从 Bezier Bezier 曲线到样条曲线曲线到样条曲线(1) Bezier (1) Bezier 曲线在应用中的不足：曲线在应用中的不足： 缺乏灵活性一旦确定了特征多缺乏灵活性一旦确定了特征多边形的顶点数边形的顶点数(m(m个个) )，也就决定了曲，也就决定了曲线的阶次线的阶次(m-1(m-1次次) )，无法更改；，无法更改； 控制性差当顶点数较多时，曲控制性差当顶点数较多时，曲线的阶次将较高，此时，特征多边形线的阶次将较高，此时，特征多边形对曲线形状的控

2、制将明显减弱；对曲线形状的控制将明显减弱；赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 不易修改由曲线的混合函数可不易修改由曲线的混合函数可看出，其值在开区间看出，其值在开区间 ( 0 , 1 ) ( 0 , 1 ) 内均不为内均不为零。因此，所定义之曲线在零。因此，所定义之曲线在 ( 0 t 1)( 0 t 1)的区间内的任何一点均要受到全部顶的区间内的任何一点均要受到全部顶点的影响，这使得对曲线进行局部修点的影响，这使得对曲线进行局部修改成为不可能。改成为不可能。（而在外形设计中，（而在外形设计中， 局部修改是随时要进行的）局部修改是随时

3、要进行的）赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期,!( )(1)!()!in ii nnBttti ni01tB0,3B1,3B2,3B3,330,3( )(1) Btt21,3( )3 (1) Bttt22,3( )3 (1) Bttt33,3( )Btti,n(t) i,n(t) 伯恩斯坦基函数伯恩斯坦基函数赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期为了克服为了克服 Bezier Bezier 曲线存在的问题，曲线存在的问题，Gordon Gordon 等人拓展了等人拓展了 Be

4、zier Bezier曲线，就曲线，就外形设计的需求出发，希望新的曲线外形设计的需求出发，希望新的曲线要：要：易于进行局部修改；易于进行局部修改； 更逼近特征多边形；更逼近特征多边形； 是低阶次曲线。是低阶次曲线。于是，用于是，用 n n次样条基函数替换了伯次样条基函数替换了伯恩斯坦基函数，构造了称之为样条恩斯坦基函数，构造了称之为样条曲线的新型曲线。曲线的新型曲线。赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期2.2.样条曲线的数学表达式样条曲线的数学表达式样条曲线的数学表达式为：样条曲线的数学表达式为：nknkkinitFPtP0,)()(

5、控制顶点控制顶点在上式中在上式中t1 , i= 0, 1, 2, , mt1 , i= 0, 1, 2, , m所以可以看出：所以可以看出：样条曲线是分段定样条曲线是分段定义的。如果给定义的。如果给定 m+n+1 m+n+1 个顶点个顶点 Pi ( i= Pi ( i=0, 1, 2, m+n)0, 1, 2, m+n)，则可定义，则可定义 m+1 m+1 段段 n n 次的参数曲线。次的参数曲线。 B样条基函数样条基函数赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期)!( !rnrnCrnknjnjnjnkjkntCntF01,)() 1(!

6、1)(式中： t 1 k = 0, 1, 2, , n 在以上表达式中：在以上表达式中：F Fk,nk,n( t )( t )为为n n次次B B样条基函数，也称样条基函数，也称样条分段混合函数。其表达式为：样条分段混合函数。其表达式为：赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期连接全部曲线段所组成的整条曲线称连接全部曲线段所组成的整条曲线称为为 n n 次样条曲线。依次用线段连接次样条曲线。依次用线段连接点点 P Pi+ki+k (k=0,1,n) (k=0,1,n)所组成的多边折所组成的多边折线称为样条曲线在第线称为样条曲线在第i i段

7、的特征多段的特征多边形。边形。 赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期. .二次二次样条曲线样条曲线在二次样条曲线中，在二次样条曲线中，n=2,k=0,1,2n=2,k=0,1,2故其基函数形式为：故其基函数形式为：22,222,1222220232,021)()122(21)()1(21!2!3)1(!2!3)2(!3!321)2()1(!21)(ttFtttFttttjtCtFjjjknjnjnjnkjkntCntF01,)() 1(!1)(赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二

8、学期 有了基函数，因此可写出二次样条有了基函数，因此可写出二次样条曲线的分段表达式为：曲线的分段表达式为：22, 212, 12, 0)()()()(iiiiPtFPtFPtFtP( i= 0,1,2,m )m+1段赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期写成一般的矩阵形式为：写成一般的矩阵形式为：2021022,011022121211)()(kkkBBBttBtFtP 式中，式中，k k为分段曲线的特征多边形为分段曲线的特征多边形的顶点：的顶点：B B0 0,B,B1 1,B,B2 2。对于第。对于第i i段曲线的段曲线的B Bk k

9、 即为：即为：P Pi i,P,Pi+1i+1,P,Pi+2i+2 连续的三个顶连续的三个顶点。（见下图）点。（见下图）赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期P3B:P0P0,P1,P2P2P1P1,P2,P3B:P4n=2,n=2,二次二次B B样条曲线样条曲线m+n+1m+n+1个顶点，三个顶点，三点一段，共点一段，共m+1m+1段。段。i=0P0,2(t)i=1P1,2(t)nknkkinitFPtP0,)()(赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期二次样条曲线的性质二次样

10、条曲线的性质先对 P(t)求导得：2100111211)(BBBttP然后分别将 t=0,t=0.5,t=1 代入 P(t)和 P(t)，可得：P(0)=1/2(B0+B1), P(1)=1/2(B1+B2); P(0)=B1-B0, P(1)=B2-B1; P(1/2)=1/21/2P(0)+P(1)+B1 P(1/2)=1/2(B2-B0)=P(1)- P(0)赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期与以上这些式子所表达的性质相符的曲线是何种形状：（见下图）B0P(0)P(1)MB2P(1/2)B1P(1/2)是什么曲线？与Bezie

11、r曲线有何差别？赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期结论：分段二次结论：分段二次B B样条曲线是一条抛样条曲线是一条抛物线；有物线；有n n个顶点定义的二次个顶点定义的二次B B样条曲样条曲线，其实质上是线，其实质上是n-2n-2段抛物线（相邻三段抛物线（相邻三点定义）的连接，并在接点处达到一点定义）的连接，并在接点处达到一阶连续。（见下图）阶连续。（见下图）P3P0P2P1P4F282.c 二次二次 B-样条曲线样条曲线 赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期. .三次样条曲

12、线三次样条曲线分段三次样条曲线由相邻四个顶点分段三次样条曲线由相邻四个顶点定义，其表达式为：定义，其表达式为：P(t)=FP(t)=F0,30,3(t)(t)B B0 0+F+F1,31,3(t)(t)B B1 1+F+F2,32,3(t)(t)B B2 2+F+F3,33,3(t)(t)B B3 3( (0 0t1)t1)可见，由可见，由 n n 个顶点定义的完整的三次个顶点定义的完整的三次样条曲线是由样条曲线是由 n-3 n-3 段分段曲线连接段分段曲线连接而成的。很容易证明，三次样条曲而成的。很容易证明，三次样条曲线在连接处达到二阶连续。线在连接处达到二阶连续。 *F281.c 三次三次

13、 B-样条曲线样条曲线 赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期样条曲线是一种非常灵活的曲线，曲线样条曲线是一种非常灵活的曲线，曲线的局部形状受相应顶点的控制很直观。这的局部形状受相应顶点的控制很直观。这些顶点控制技术如果运用得好，可以使整些顶点控制技术如果运用得好，可以使整个样条曲线在某些部位满足一些特殊的技个样条曲线在某些部位满足一些特殊的技术要求。如：术要求。如：可以在可以在曲线中构造一段直线；曲线中构造一段直线； 使曲线与特征多边形相切；使曲线与特征多边形相切； 使曲线通过指定点；使曲线通过指定点； 指定曲线的端点；指定曲线的端点

14、； 指定曲线端点的约束条件。指定曲线端点的约束条件。赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 对于特征多边形的逼近性对于特征多边形的逼近性 二次二次B B样条曲线优于三次样条曲线优于三次B B样条曲线样条曲线 三次三次BezierBezier曲线优于二次曲线优于二次BezierBezier曲线曲线 相邻曲线段之间的连续性相邻曲线段之间的连续性 二次二次B B样条曲线只达到一阶导数连续样条曲线只达到一阶导数连续 三次三次B B样条曲线则达到二阶导数连续样条曲线则达到二阶导数连续 角点的修改对曲线形状的影响角点的修改对曲线形状的影响 Bezi

15、erBezier曲线：修改一个角点将影响整条曲线的形状。曲线：修改一个角点将影响整条曲线的形状。 B B样条曲线：修改一个角点只影响该角点所在位置前后样条曲线：修改一个角点只影响该角点所在位置前后 三段曲线的形状。三段曲线的形状。B样条曲线的适用范围样条曲线的适用范围赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期角点重叠和角点共线(*) 二重角点二重角点 若要使B样条曲线与特征多边形相切，可运用二重角 点的方法。Q0Q1Q3Q2Q4P1P2Q（0-1）P0Q0Q1Q3Q2Q4P1P2Q（0-1）赤峰学院计算机系赤峰学院计算机系 计算机图形学计算

16、机图形学 08-0908-09第二学期第二学期 三重角点 若要使若要使B B样条曲线产生一个尖点，可运用三重角样条曲线产生一个尖点，可运用三重角 点的方法。点的方法。Q1Q4Q3Q2Q5P1P2Q0P0Q6P3P4赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 三角点共线 若要使B样条曲线产生反向弧切接的效果，可运用三角点共线的方法。Q1Q4Q3Q2P0P2Q0P1。赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 四角点共线 若要使B样条曲线段之间切接入一段直线，可运用四角点共线的方法。Q1Q4

17、Q3Q2P0P2Q0P1。Q5P3赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期6.5 Bzier曲面和B样条曲面 6.5.1 Bzier曲面赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 曲面的形状、位置由边界上的四个角点决定。中间角点只曲面的形状、位置由边界上的四个角点决定。中间角点只反映曲面的凹凸程度。反映曲面的凹凸程度。BzierBzier曲面及控制网格演示动画曲面及控制网格演示动画赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 B Be

18、zier曲面是Bezier曲线的扩展， B Bezier曲面的边界线就是由四条Bezier曲线构成的。三次Bezier曲线段由四个控制点确定，三次Bezier曲面片则由 控制点确定。16个控制点组成一个矩阵：B=Q00 Q 10 Q20 Q30 Q01 Q 11 Q21 Q31 Q02 Q12 Q22 Q32Q03 Q13 Q23 Q33Q00Q10Q20Q30Q01Q31Q32Q02Q03Q33Q13Q23Q11Q21Q12Q22wu赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期给定空间给定空间16个位置点个位置点bij，可以确定，可以确定一

19、张三次一张三次Bezier曲面片。曲面片。3 , 2 , 1 , 0)()(3,30ivBbvbjjiji303,*)()()(iiiuBvbubq 由曲线拓展为Bzier曲面rijuvuvV*uv 根据根据“线动成面线动成面”的思想，按设定间的思想，按设定间隔取隔取 ，在四条，在四条v v线上取点，沿线上取点，沿u u向生成三次向生成三次BezierBezier曲线：曲线： 1 ,0*v将将u,v向曲线方程合并得：向曲线方程合并得： 30303 ,3 ,303 ,),()()()()(),(ijjiiiijibvBuBuBvbvub首先生成四条首先生成四条v向的三次向的三次Bezier曲线：

20、曲线：赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期Bzier曲面的定义 在空间中给定在空间中给定( (n+1)(m+1)个点，称以个点，称以下张量积形式的下张量积形式的参数多项式曲面为参数多项式曲面为nm次的次的Bzierzier曲面曲面: : 贝塞尔曲面表达式如下：贝塞尔曲面表达式如下： n m P(u,v)= bi,jBi,n(u)Bj,m(v) 0u,v1 i=0 j=0 贝塞尔曲面中应用最广泛的是双贝塞尔曲面中应用最广泛的是双3次贝塞尔曲面，次贝塞尔曲面，它由给出的它由给出的4*4个网格点唯一决定个网格点唯一决定.赤峰学院计算机系赤峰

21、学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期P04P03P02P40P30P20P10P00P01P11P21P31P41P14P(u,0)P(0,v)图9.15 Bzier曲面的控制网格一般称一般称 为为 的控制顶点，把由的控制顶点，把由 和和 组成的网格组成的网格称为称为 的控制网格，记为的控制网格，记为 , ,如图如图9.159.15所示。所示。 ijP01,iiimPPP(0,1, )in01,jjnjPPP(0,1,)jmijP( , )P u v( , )P u v赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学

22、期第二学期BzierBzier曲面的矩阵表示是：曲面的矩阵表示是：0,000101,101110,1,01( )( )( , )( ) ( ) ( )( )mmmmnnn nm mnnnmJvJvu vJu JuJuJvPPPPPPPPPP赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期控制网格的四个角点控制网格的四个角点 是曲面是曲面 的四个端点。决定了曲线的形状，位置。的四个端点。决定了曲线的形状，位置。P03P02P30P20P10P01P11P21P31P41P14P(u,0)P(0,v)Bzier曲面的端点和边界线(1) (1) 端点位

23、置端点位置0000,mnnmPPPP0000(0,0),(0,1)(1,0),(1,1)mnnmPPPPPPPPBzierBzier曲面具有以下性质：曲面具有以下性质：P04P40P0044P( , )P u v赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期(2) (2) 边界线的位置边界线的位置(0, ) vP( ,0)uP(1, )vP( ,1)uP00010mP PP00100nP PP01nnnmP PP01mmnmP PP 的4条边界线 ， ， ， 是Bzier曲线，它们分别以 ， ， ， 为控制多边形。 P03P02P30P20P1

24、0P01P11P21P31P41P14P(u,0)P(0,v)Bzier曲面的端点和边界线P04P40P0044P( , )P u v赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 双一次双一次Bezier曲面：取曲面：取m=n=1。这是。这是一张双曲抛物面（马鞍面）一张双曲抛物面（马鞍面） ； 双二次双二次Bezier曲面：取曲面：取m=n=2。该曲。该曲面的四条边界是抛物线。面的四条边界是抛物线。 ； 双三次双三次Bezier曲面：取曲面：取m=n=3； 1 , 0,)()(),(00,vupuBuBvuPnimjijmjni1 mn注：矩

25、阵表示见课本注：矩阵表示见课本赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期端点 的u向切矢和v向切矢分别为 和 ，所以三角形 所在的平面在P P00点和曲面相切。同理，三角形 , , 所在的平面分别在点 , , 处与曲面相切。1000()nPP0100()mPP001001P P P010,1mmmP P P1,1nmnmn mP PP01,01nnnP PP0mPnmP0nP00P(3) (3) 端点的切平面端点的切平面10P20P01P11P21P31P02P22P12P32P13P23PBezier曲面的端点切平面00)0 , 0(PP

26、30)0 , 1 (PP03) 1 , 0(PP33) 1 , 1 (PP)0 ,(uP) 1 ,(uP), 0(vP), 1 ( vP赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期ijP(0,1, ,0,1, )in jm00P( , )u vP00010010 P PP P0mPnmP0nP由端点的切平面知， 是 在点 处的法向；其余各端点 ， ， 处法向的情况也类似。(4) (4) 端点的法向端点的法向(5) (5) 凸包性凸包性曲面 位于其控制顶点 的凸包内。( , )P u v赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学

27、 08-0908-09第二学期第二学期(6) (6) 仿射不变性仿射不变性曲面 的形状仅与点 的位置有关，而与坐标系的选择无关。(7) (7) 拟局部性拟局部性修改一个控制顶点时，曲面上距离它较近的点受影响较大。要改变曲面某部分的形状，只要交互调节相应的控制顶点即可。 ( , )u vP0001,nmPPP赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期曲面片1曲面片2 Beziet曲面片的拼接边界线P0,0Q3,0P3,3(Q0,3)P0,3Q3,3P3,1(Q0,1)P3,0(Q0,0)P3,2(Q0,2)P(u,v)Q(u,v)赤峰学院计算

28、机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期d243244434241343332312423222114131211321610002121216121103261212161612121610211210210210613261)1 () , (vvvdddddddddddddddduuuvupuvV1kd42d43d44d11d12d13d14d21d23d31d32d33d34C4C3C1C2V2kV3kV4kd41P(u,vK)d22q 双三次均匀B样条曲面P(u,v)的矩阵表示iiu U赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形

29、学 08-0908-09第二学期第二学期设节点向量设节点向量 , ( , ), ( , )分别是分别是对参数对参数 平面的平面的 轴和轴和 轴的分割，如图轴的分割，如图1 1所示。称下列张量积形所示。称下列张量积形式的参数曲面为式的参数曲面为 ( k( kn n,h,hm m) )阶的阶的B B样条曲面样条曲面 , ,u uk k-1-1u uu un n1 1, ,v vh h-1-1v vv vm m1 1其中其中 是空间中给定的是空间中给定的( (n+n+1)1)( (m m+1)+1)个网格点，通常称为个网格点，通常称为 的的控制顶点。控制顶点。 ， 分别是关于节点向量分别是关于节点向

30、量U,VU,V的的k k阶和阶和h h阶的阶的B B样条样条基函数。基函数。 jjvViu1iujv1jvuvuvkh ,0 0(,)() ()n mij ikjhi juvBuBv PPijP,()ikBu,( )j hBv( , )u vP6.5.2 B样条曲面的定义和性质vm 图1 uv平面的分割u0u1uiunv0v1vju2v20 1 2ii iimPPP P赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期由两组多边形由两组多边形 ( )和和 ( )组成的网格（如图组成的网格（如图2）称为）称为 的的控制网格，简记为控制网格，简记为 。

31、0,1,2,in012jjjnjPP PP0,1,2,jmi jP0 1 2ii iimPPP P00P10P20P30P01P11P21P31P02P22P12P32P03P23P33P图3.1.33 双三次B样条曲面片P P00P P01P P10P P04P P20P P40图2 B样条曲面及其控制网格赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期d24uvd42d43d44d11d12d13d14d21d23d31d32d33d34C4C3C1C2d41d22q 均匀均匀B样条曲面样条曲面给定给定16个顶点个顶点dij(i=1,2,3,

32、4 j=1,2,3,4)构成的特征网格，可以定义一张曲面片。构成的特征网格，可以定义一张曲面片。用用di1、di2、di3、di4(i=1,2,3,4 )构建四条构建四条V向曲线向曲线C1、C2、C3和和C4(图中虚线图中虚线);参数参数v在在0,1 之间取之间取值值vk ，对应于，对应于vk曲线曲线C1、C2、C3和和C4上可得到上可得到v1k、v2k、v3k和和v4k四个点，四个点，该四点构成该四点构成u向的一个向的一个特征多边形，定义一条特征多边形，定义一条新的曲线新的曲线P(u,vk)；uvC4C3C1C2V1kV3kV4kV2k当参数当参数vk在在0,1 之间取不同值时，之间取不同值

33、时，P(u,vk)沿箭头方向扫描，即得到沿箭头方向扫描，即得到由给定特征网格由给定特征网格dij(i=1,2,3,4 j=1,2,3,4)定义的双三次均匀定义的双三次均匀B样条样条曲面片曲面片P(u,v)。赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期与与BzierBzier曲面一样，曲面一样， 是对曲面是对曲面 的大致形状的勾画，的大致形状的勾画， 是对是对 的逼的逼近。近。B B样条曲面也具有局部调整性、凸包样条曲面也具有局部调整性、凸包性、几何不变性等，它的控制网格也是性、几何不变性等，它的控制网格也是人机交互的手段，也可以通过某些算法人

34、机交互的手段，也可以通过某些算法对其进行计算，这些都与对其进行计算，这些都与B B样条曲线的情样条曲线的情况类似。况类似。 00P10P20P30P01P11P21P31P02P22P12P32P03P23P33P图3.1.33 双三次B样条曲面片0 1 2ii iimPPP P0 1 2ii iimPPP P00P10P20P30P01P11P21P31P02P22P12P32P03P23P33P图3.1.33 双三次B样条曲面片赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期双三次B样条曲面片b样条曲面赤峰学院计算机系赤峰学院计算机系 计算机

35、图形学计算机图形学 08-0908-09第二学期第二学期6.6实验曲线的绘制方法实验曲线的绘制方法1 1、最小二乘法、最小二乘法 在科学研究中，通过实验或测量，可以获得大量的在科学研究中，通过实验或测量，可以获得大量的实验数据。一般在获得数据之后，对这些数据进行某实验数据。一般在获得数据之后，对这些数据进行某种处理，然后绘成图形。种处理，然后绘成图形。 但由于实验本身会受到各种具体因素的影响，使得但由于实验本身会受到各种具体因素的影响，使得通过实验测得的数据或多或少地带有误差。也就是说通过实验测得的数据或多或少地带有误差。也就是说，这些实验数据本身并不准确。因此如果仅仅是简单，这些实验数据本身

36、并不准确。因此如果仅仅是简单地将这些数据点连成曲线，那么这种看起来似乎很精地将这些数据点连成曲线，那么这种看起来似乎很精确的方法恰恰是不符合实际情况的，也是不可取的。确的方法恰恰是不符合实际情况的，也是不可取的。 正确的方法应该是用一条平滑的曲线以适当的方式正确的方法应该是用一条平滑的曲线以适当的方式来尽可能地靠近这些数据点，以弥补由于误差造成的来尽可能地靠近这些数据点，以弥补由于误差造成的数据点的跳动。数据点的跳动。赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 那么对于一系列的数据点（那么对于一系列的数据点（x xi i,y,yi i）（）

37、（i=1,2ni=1,2n），），所要绘制的曲线所要绘制的曲线y=f(x)y=f(x)，用什么样的标准来评价这条，用什么样的标准来评价这条曲线是否处于较为合理的状态呢？曲线是否处于较为合理的状态呢？ 通常把数据点的坐标值与曲线上对应的坐标值之通常把数据点的坐标值与曲线上对应的坐标值之差差作为评判的标准：作为评判的标准：iiiyxf)(i：称为残差f(xi)：为理论值yi：为相应的实测值 常用的评判方法是：使残差的平方和即常用的评判方法是：使残差的平方和即 达到最达到最小。这也就是所谓的最小二乘法。小。这也就是所谓的最小二乘法。2i赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08

38、-0908-09第二学期第二学期 用最小二乘法来绘制实验曲线，其实质也就是要找用最小二乘法来绘制实验曲线，其实质也就是要找出一个经验方程出一个经验方程y=f(x)y=f(x)用来描述这些点。所以，工作用来描述这些点。所以，工作的第一步首先要根据已知数据点的分布情况进行一下的第一步首先要根据已知数据点的分布情况进行一下预测，经验方程可能属于什么类型，比如说是线形函预测，经验方程可能属于什么类型，比如说是线形函数，还是其它阶次的多项式函数。数，还是其它阶次的多项式函数。2 2、用最小二乘法拟合直线、用最小二乘法拟合直线 设有测得的一组数据点（设有测得的一组数据点（x xi i,y,yi i）（）（

39、i=1,2ni=1,2n）。）。根据这些数据点的分布情况，预测到它们之间呈线根据这些数据点的分布情况，预测到它们之间呈线形关系：形关系：21axay那末，可按最小二乘法的原理建立起下面的式子：那末，可按最小二乘法的原理建立起下面的式子：niniiiiyaxa112212)(赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期（x xi i,y,yi i）为测得的已知数据点的值，故这个方程）为测得的已知数据点的值，故这个方程是关于是关于a a1 1和和a a2 2的函数，即：的函数，即：niiiyaxaaaf122121)(),( 当当a a1 1和和

40、a a2 2为何值时，该函数为何值时，该函数f(af(a1 1,a,a2 2) )能取得极小值能取得极小值。这是一个二元函数的求极值问题，其条件为：。这是一个二元函数的求极值问题，其条件为：0),(0),(121221aaafaaaf0)(20)(22121iiiiixyaxayaxa赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期展开整理后得：展开整理后得：iiiiiiyxxaxaynaxa22121写成矩阵形式为：写成矩阵形式为：iiiiiiyxyaaxxnx212 求解该方程组可解出未知数求解该方程组可解出未知数a a1 1和和a a2 2

41、的值，从而使线的值，从而使线形函数表达式：形函数表达式：21axay成为已知，并可根据该表达式绘出图形。成为已知，并可根据该表达式绘出图形。赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 1. 双二次Bezier曲面的4条边界都是抛物线，其特征网格有（ ）个顶点。 A）8B）9C）10D）16赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期贝塞尔曲面表达式如下： m n P(u,v)=Pi,jBi,m(v)Bj,n(u) 0u,v1 i=0 j=0贝塞尔曲面中应用最广泛的是双3次贝塞尔曲面，它由给

42、出的4*4个网格点唯一决定，P(u,v)=(u3 u2 u 1) 两个贝塞尔曲面块的连接交于一条公共边，为了保证在这条边上的光滑性，必须满足下面条件：(1)定义公共边的四点相等(2)公共边上四点与其在二块曲面上相邻的两个点共线贝塞尔曲面贝塞尔曲面贝塞尔曲面贝塞尔曲面1010-1 3 -3 1 3 -6 3 0-3 3 0 0 1 0 0 0P00P01P02P03P10P11P12P13P20P21P22P23P30P31P32P33v3v2v1-1 3 -3 1 3 -6 3 0-3 3 0 0 1 0 0 0赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第

43、二学期第二学期 曲面的形状、位置由边界上的四个角点决定。中间四个角点只反映曲面的凹凸程度。v(u,w)= U N B NT WT (0 u1 , 0 w1 )Bezier曲面得表达式：式中 U=u3 u2 u1 1 WT=w3 w2 w1 1 T 3 -6 3 0 -3 3 0 0 1 0 0 0-1 3 -3 1N=NT（与Bezier曲线相同）赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期写成，三个方向分量得形式：X(u,w)= U N Bx NT WT Y(u,w)= U N By NT WT (0 u1 , 0 w1 )Z(u,w)=

44、U N Bz NT WT赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期样条曲面也是样条曲线的推广，与三次Bezier曲面一样，三次样条曲面片也是由4 4控制点确定的。同样，16个控制点写成如下矩阵形式：B=Q00 Q 10 Q20 Q30 Q01 Q 11 Q21 Q31 Q02 Q12 Q22 Q32Q03 Q13 Q23 Q33与三次样条曲线一样，三次样条曲面也很好地解决了曲面片之间的连接问题。Q00Q10Q20Q30Q01Q31Q32Q02Q03Q33Q13Q23Q11Q21Q12Q22wu赤峰学院计算机系赤峰学院计算机系 计算机图形学计

45、算机图形学 08-0908-09第二学期第二学期样条曲面的表达式为：v(u,w)= U N B NT WT (0 u1 , 0 w1 )U、W、B矩阵与Bezier曲面是一样的。 3 -6 3 0 -3 0 3 0 1 4 1 0-1 3 -3 11/6N= （与样条曲线相同）写成X，Y，Z三个方向的分量形式:x(u,w)=U N Bx NT WTy(u,w)= U N By NT WT (0 u1 , 0 w1 )z(u,w)= U N Bz NT WT赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期Void B_face( float *Bx

46、, *By, *Bz, *projection_matrix44; float Fx, Fy)/绘制B样条曲面的程序，16个控制点的x,y,z坐标分别/在4x4的矩阵Bx,By,Bz中,projection_matrix44为投影变换矩阵,/Fx,Fy为投影后平移量 float N44, Nt44, U14, Wt41, temp_product1414; float product44x44, product14x14; ; float x, y, z, x1, y1,; construct_matrix44(N, -1, 3, -3, 1, 3, -6, 0, 4, -3, 3, 3, 1

47、, 1, 0, 0, 0); construct_matrix44(Nt,); /先求出product44x = N * Bx * Nt matrix44_matrix44(temp_product44, N, Bx); matrix44_matrix44(product44x, temp_prodct44, Nt); 同样算出product44y，product44z；赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期for ( float u = 0; u 1; u +=0.1) construct_matrix14( U, u*u*u, u*

48、u, u, 1); /product14x = U * product44x = U * N * Bx * Nt; matrix14_matrix44( product14x, U, product44x); matrix14_matrix44( product14y, U, product44x); matrix14_matrix44( product14z, U, product44x);赤峰学院计算机系赤峰学院计算机系 计算机图形学计算机图形学 08-0908-09第二学期第二学期 for ( float w = 0; w 1; w +=0.1) construct_matrix41(

49、Wt, w*w*w, w*w, w, 1); matrix14_matrix41( x, prouduct14x, Wt); matrix14_matrix41( y, prouduct14y, Wt); matrix14_matrix41( z, prouduct14z, Wt); /再做投影变换,Fx,Fy为在投影面上的平移量 projection_transform( x, y, z, x1, y1, projection_matrix44, Fx, Fy); if (w = 0) moveto( x1, y1); lineto( x1, y1); / for ( float w = 0; w 1; w +=0.1) / for ( float u = 0; u 1; u +=0.1)