实验五基本图形处理技术.docVIP

实验七 曲线与曲面 一、实验目的 根据曲线和曲面的基础知识和常用曲线的数学基础，对其算法进行程序设计，验证算法的正确性，并通过程序结果加深对常用曲线数学模型的理解。 二、实验任务 1．抛物线程序设计； 2．Hermite 曲线程序设计； 3．Bezier曲线的算法实现； 4．B样条曲线的程序设计 三、实验内容和实验步骤 。 任务一：Hermite 曲线程序设计 Hermite 曲线方程为： P(t)=FB=TMB=[ t3 t2 t 1 ] 程序设计时只考虑二维图形的显示，其代数形式为： x(t)=TMBx ， Bx =[ P0x P1x R0x R1x]T y(t)= TMBy ， By =[ P0y P1y R0y R1y]T 所以，只要给出Hermite曲线的起点坐标（P0x，P0y），终点坐标（P1x，P1y），以及起点处的切矢量（R0x，R0y）和终点处的切矢量（R1x，R1y），参数变量t在[0，1]的范围内分别取0.01，0.02，…，1，步长为0.01，取100个点，分别求出P(t)=[ x(t)，y(t)]，在计算机屏幕上显示出每个坐标点，即可绘出Hermite曲线。 任务二：Bezier曲线的算法实现 （1）算法描述 根据 Q(t)=∑Pi?Bi,n(t) = ∑Pi? C in ti(1-t)n-i 并且 C in =C(n,k)= C(n,k-1) Bezier曲线的C语言算法描述如下： #include math.h #include graphics.h void computeCoefficients(int n, int *c) { int k, i; for (k=0;k=n;k++) /*求C kn =n!/(k!(n-k)!= )*/ { int c[k]=1; for (i=n; i=k+1; i--) /*求c[k]=n*(n-1)…(k+1) */ c[k]*=i; for (i=n-k; i=2; i--) /*求c[k]/(n-k)!*/ c[k]/=i; } } void computepoint(float t,wcPt3 *pt,int ncontrols,wcPt3 *controls,int *c) { int i, n=ncontrols-1; float blend; pt-x=0.0; pt-y=0.0; pt-z=0.0; for (i=0; incontrols; i++) { blend = c[k]*powf(t,i)*powf(1-t,n-i); /*求 C kn ti(1-t)n-i */ pt-x+=controls[i].x*blend; /*求x(t)*/ pt-y+=controls[i].y*blend; /*求y(t)*/ pt-z+=controls[i].z*blend; /*求z(t)*/ } } void Bezier(wcPt3 *controls, int ncontrols, int m, wcPt3 *curve) { int *c=(int *) malloc(ncontrols *sizeof(int)); int i; computecoefficients(ncontrols-1,c); for (i=0; i=m; i++) computepoint(i/(float)m,curve[i],ncontrols,controls,c); free(c); } 在主程序中提供特征多边形的各个顶点坐标放入controls[]数组中，ncontrols为顶点的个数，m为曲线上取的样点数，比如m=100表示取100个样点。计算出曲线上的各个样点坐标放入curve[]数组中，这样可以通过相邻点连线绘出生成的Bezier曲线。 程序实现步骤：（工程名:BezierCurve） （ⅰ）C m n的函数实现，定义成成员函数，命名为Multiply_n。 C mn = = int Multiply_n(int m,