自由曲线和曲面的绘制.docVIP

计算机图形学课程实验 报 告 实验题目 自由曲线和曲面的绘制 班 级 计算091 姓 名 张 阳 学 号 3090811019 指导教师 胡 钢 日 期 2012-6-17 西安理工大学理学院应用数学系 二0一二年春季学期 (实验说明 实验目的：掌握自由曲线和曲面(包括Bezier曲线Direction): 1.每个学生单独完成； 2.开发语言为TurboC或C++，也可使用其它语言； 3.请在自己的实验报告上写明姓名、学号、班级； 4.每次交的实验报告内容包括：题目、试验目的和意义、程序制作步骤、主程序、运行结果图以及参考文件； 5. 自己保留一份可执行程序,考试前统一检查和上交。 (实验内容 (实验题一 1.1实验题目 上机编写一个能绘制Bezier曲线的程序，并调试成功Bezier曲线B样条曲线B样条曲线 Step1:选择Bezier或B样条曲线绘图； Step2:若选择用Bezier曲线绘图，输入控制顶点个数，依次输入控制顶点坐标； 若选择用B样条曲线绘图，输入B样条曲线段次数，输入控制顶点个数，依次输入控制顶点坐标； Step3:则将t [0,1]区间剖分成m等分，对于每一个ti，i=0,1,…,m，根据式（1）或（2）都可以计算出一个P(ti)，计算所有的P(ti)。 Step4:依次连接每个P(ti)，就得到Bezier或B样条曲线 1.4主程序 #includeiostream.h #includemath.h #includemalloc.h #includestdio.h #includegraphics.h //using namespace std; #define pi 3.1415926 /************** 二维 点类 **********************/ class Point { public: double x; double y; Point(int x=0,int y=0) //构造函数 { this-x = x; this-y = y; } void operator=(Point a) //重载=运算符 { x=a.x; y=a.y; } Point operator*(double a) //重载*运算符 { return Point(a*x,a*y); } Point operator+(Point a) //重载+运算符 { return Point(x+a.x,y+a.y); } void operator+=(Point a) //重载+=运算符 { x+=a.x; y+=a.y; }} /* ******************** 基本的函数 ************************* */ //求阶乘 long int Factorial(int n) { int i,sum=1; if(n==0) return 1; for(i=2;i=n;i++) sum *= i; return sum; } //求Bernstein基函数 double Bernstein(int i,int n,double t) { return (double)Factorial(n)/Factorial(i)/Factorial(n-i)*pow(t,i)*pow(1-t,n-i); } //求B_SplineBase基函数 double B_SplineBase(int i,int n,double t) { double sum = 0; int j; for(j=0;j=n-i;j++) { sum+=pow(-1,j)*Factorial(n+1)/Factorial(j)/Factorial(n+1-j)*pow(t+n-i-j,n); } return sum/Factorial(n); } /* ******************** 曲线 ************************* */ //求Bezier曲线 void Bezier(Point p[],int pn,int m) { //将t【0，1】分成多少等分 m double h = 1.0/m; double t = 0; int i,j; int *pB = (int *)malloc(sizeof(int)*(2*m+2));