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

实验七曲线与曲面 一、实验目的 根据曲线和曲面的基础知识和常用曲线的数学基础，对其算法进行程序设计，验证算法的正确性，并通过程序结果加深对常用曲线数学模型的理解。 二、实验任务 1．抛物线程序设计； 2．Hermite曲线程序设计； 3．Bezier曲线的算法实现； 4．B样条曲线的程序设计 三、实验内容和实验步骤 。 任务一：Hermite曲线程序设计 Hermite曲线方程为： P0 P1R P0 P1 R0 R1 -2 1 1 -3 3 -2 -1 0 0 1 0 1 0 0 0 P(t)=FB=TMB=[t3t2t1] 程序设计时只考虑二维图形的显示，其代数形式为： x(t)=TMBx，Bx=[P0xP1xR0xR1x]T y(t)=TMBy，By=[P0yP1yR0yR1y]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曲线的算法实现 n（1）算法描述n n n i=0i=0根据Q(t)=∑Pi?Bi,n(t)=∑Pi?Cinti(1-t)n-i i=0 i=0 n-k+1 n-k+1 n并且Cin=C(n,k)=C(n,k-1) n Bezier曲线的C语言算法描述如下： #includemath.h #includegraphics.h voidcomputeCoefficients(intn,int*c) { intk,i; (n-k)!n*(n-1)…(k+1)for(k=0;k=n;k++)/*求Ckn=n!/(k!(n-k)!=)*/ (n-k)! n*(n-1)…(k+1) { intc[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; } } voidcomputepoint(floatt,wcPt3*pt,intncontrols,wcPt3*controls,int*c) { inti,n=ncontrols-1; floatblend; 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);/*求Cknti(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)*/ } } voidBezier(wcPt3*controls,intncontrols,intm,wcPt3*curve) { int*c=(int*)malloc(ncontrols*sizeof(int)); inti; 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） （ⅰ）Cmn的函数实现，定义成成员函数，命名为Multiply_n。 （m+1 （m+1）(m+2)…(n-1).n n! (n-m)!m!(n-m)!Cmn== (n-m)! m!(n-m)! intMultiply_n(intm,n) { inti,j,a; if(m!=0) { a=1; for(i=m+1;i=n;i++)//求（m+1）(m+2)…(n-1).n a=a*i; for(j=1;j=n-m;j++)//求(n-m)!和Cmn a=a/j; returna; } else return1; } （ⅱ）伯恩斯坦多项式Bm,n(t)的函数实现 Bm,n(t)=Cmntm(1-t)n-m DoubleBernstein(intm,intn,doubl