第四章 曲线和曲面.pptVIP

**第125页，共199页，星期日，2025年，2月5日设己知三次Bezier曲线P(t)的控制顶点是P0，P1，P2，P3，在P()处将曲线分为两段，求出前半段的控制顶点Q0，Ql，Q2，Q3和后半段的控制顶点R0，R1，R2，R3。有算法如下voidsplit_Bezier(PointP[]){ PointR[4],Q[4]; inti,j; for(i=0;i=3;i++) R[i]=P[i]; for(i=0;i=2;i++){ Q[i]=R[0]; for(j=0;j=2-i;j++){ //分别对相邻两控制点间的线段进行分裂 R[j].x=(R[j].x+R[j+1].x)/2; R[j].y=(R[j].y+R[j+1].y)/2; } } Q[3]=R[0];}**第126页，共199页，星期日，2025年，2月5日分裂算法的计算**第127页，共199页，星期日，2025年，2月5日根据Bezier曲线的凸包性质，知道曲线上任意一点到线段P0P3的距离，小于P1和P2到线段P0P3距离中的较大者，即有：任意事先给定的对画出曲线近似程度的要求ε0，可以取max(d(P1，P0P3)，d(P2，P0P3))ε为分裂停止的条件。**第128页，共199页，星期日，2025年，2月5日voidnew_split_Bezier(PointP[]){ PointR[4],Q[4]; inti,j; constdoubleepsilon=0.01; if(maxdistance(P)epsilon) /*maxdistance(P)为求max(d(P1,P0P3)，d(P2,P0P3))的函数*/ { MoveTo(P[0].x,P[0].y); LineTo(P[3].x,P[3].y); } else { for(i=0;i=3;i++)R[i]=P[i]; for(i=0;i=2;i++) { Q[i]=R[0]; for(j=0;j=2-i;j++) { R[j].x=(R[j].x+R[j+1].x)/2; R[j].y=(R[j].y+R[j+1].y)/2; } } Q[3]=R[0]; new_split_Bezier(Q);new_split_Bezier(R); }}**第129页，共199页，星期日，2025年，2月5日5.有理Bezier曲线**第130页，共199页，星期日，2025年，2月5日图中h0=h1=h3=1，当h2=0、1/2、1、2、4时曲线逐渐地靠近P2点**第131页，共199页，星期日，2025年，2月5日6.Bezier曲面若在空间给定(m+1)(n+1)个控制点，Vij，i=0，1，…，m，j=0，1，…，n，令上式曲面为m×n次的Bezier曲面**第132页，共199页，星期日，2025年，2月5日当m=n=1，公式成为：设v00，v01，v10，v11四点依次是(0，0，0)，(1，0，0)，(0，1，0)，(0，0，1)，则可得P1,1(u，w)的坐标形式的参数方程为：消去参数，就得马鞍面方程：**第133页，共199页，星期日，2025年，2月5日马鞍面**第134页，共199页，星期日，2025年，2月5日当m=n=3，曲面成为：**第135页，共199页，星期日，2025年，2月5日7.双三次bezier曲面和双三次Coons曲面之间的关系设同一个曲面片，用Coons曲面形式确定它，又用Bezier形式确定它，现在来考查、矩阵M与矩阵B有什么关系。求解上式得：**第136页，共199页，星期日，2025年，2月5日**第137页，共199页，星期日，2025年，2月5日**第138页，共199页，星期日，2025年，2月5日四个角点00，0l，10和11分别对应控制点V00，V03，V30和V33。切向量，例如00w=3(V01-V00)，表明00w，是沿从V00到V01边的方向，大小是两个位置向量差的3倍，其余切向量是类似的。V10，V00沿w方向的切向量之差的三倍V01,V00沿u方向的切向量之差的三倍**第139页，共199页，星期日，2025年，2月5日**第140页，共199页，星期日，2025年，2月5日扭曲向量只与中间四个控制点有关系**第141页，共199页，星期日