双三次Bezier曲面

1、实验六 双三次Bezier曲面一、实验目的根据Bizer曲面的基础知识和数学基础，对其算法进行程序设计，验证算法的正确性，并通过程序结果加深对常用曲面数学模型的理解。二、实验任务（2学时）Bezier曲面算法及其程序设计。3、 实验内容和实验步骤1、算法描述Bezier曲面是由Bezier曲线拓广而来，以两组正交的Bezier曲线控制点构造空间网格来生成曲面。mn次张量积形式的 Bezier曲面的定义如下(参照教材P200式7-20)：（u，v）0，10，1 双三次Bezier曲面定义如下(参照教材P201式7-21)： （u，v）0，10，1展开上式，有代入得到：令则有：生成曲面时可以通过先

2、固定u, 变化v得到一簇Bezier曲线；然后固定v，变化u得到另一簇Bezier曲线，两簇曲线交织生成Bezier曲面。2、要求：根据给定的16个控制顶点：P00(200,20,0),P01(150,0,100),P02(50,-130,100),P03(0,-250,50)；P10(150,100,100),P11(100,30,100),P12(50,-40,100),P13(0,-110,100)；P20(140,280,90),P21(80,110,120),P22(30,30,130),P23(-50,-100,150)；P30(150,350,30),P31(50,200,150

3、),P32(0,50,200),P33(-70,0,100)；使用斜等测投影绘制双三次Bizer网格曲面。3、程序实现步骤：（工程名:BezierCurve2）步骤1：创建“BezierCurve2”工程文件；步骤2：创建类class：“ P2”及“P3”；注意：P2 为二维点，含有两个成员变量“x，y”，P3含有三个成员变量“x，y，z”。单击右键-“new class”-选择类型“Generic Class”名称为“ P2”及“P3”，添加成员变量（添加至“class EACH_ENTRY public:”之内）：class P2 public:P2();virtual P2();doub

4、le x;double y;class P3 public:P3();virtual P3();double x;double y;double z;步骤3：包含头文件并定义顶点对象。1）在“class CBezierCurve2View : public Cview”之前添加代码“#include P3.h”及“#include P2.h”；#include P3.h/包含三维坐标点类#include P2.h/包含二维坐标点类2）在“class CBezierCurve2View : public CView”内添加代码：P3 P344;/三维顶点P2 P244;/二维顶点步骤4：添加成员

5、函数。1）C m n的函数实现，定义成员函数，命名为Multiply_n。方法及过程参照“实验六 曲线及曲面生成算法（一）”；2）伯恩斯坦多项式Bm,n(t)的函数实现，添加CBezierCurve2View:bernstein。方法及过程参照“实验六 曲线及曲面生成算法（一）”；步骤5：在OnDraw（）中添加代码：注意：添加头文件#include math.h/数学头文件#define Round(d) int(floor(d+0.5)/四舍五入宏定义上述两行代码添加至：类CBezierCurve2View的所有成员函数代码之前。1）窗口坐标变换，设置客户区中心为原点CBezierCurv

6、e2Doc* pDoc = GetDocument();ASSERT_VALID(pDoc);/ TODO: add draw code for native data hereCRect rect;/定义 客户区矩形GetClientRect(&rect);/获得客户区的大小pDC-SetMapMode(MM_ANISOTROPIC);/pDC自定义坐标系pDC-SetWindowExt(rect.Width(),rect.Height();/设置窗口范围pDC-SetViewportExt(rect.Width(),-rect.Height();/设置视区范围,x轴水平向右，y轴垂直向上p

7、DC-SetViewportOrg(rect.Width()/2,rect.Height()/2);/客户区中心为原点rect.OffsetRect(-rect.Width()/2,-rect.Height()/2);2）初始化16个控制顶点P300.x=20; P300.y=0; P300.z=200;/P00P301.x=0; P301.y=100;P301.z=150;/P01P302.x=-130;P302.y=100;P302.z=50; /P02P303.x=-250;P303.y=50; P303.z=0; /P03P310.x=100; P310.y=100;P310.z=15

8、0;/P10P311.x=30; P311.y=100;P311.z=100;/p11P312.x=-40; P312.y=100;P312.z=50; /p12P313.x=-110;P313.y=100;P313.z=0; /p13P320.x=280; P320.y=90; P320.z=140;/P20P321.x=110; P321.y=120;P321.z=80; /P21P322.x=30; P322.y=130;P322.z=30; /P22P323.x=-100;P323.y=150;P323.z=-50;/P23P330.x=350; P330.y=30; P330.z=1

9、50;/P30P331.x=200; P331.y=150;P331.z=50; /P31P332.x=50; P332.y=200;P332.z=0; /P32P333.x=0; P333.y=100;P333.z=-70;/P33 3) 三维控制顶点投影（采用斜等测法）至二维点： for(int i=0;i4;i+)for(int j=0;jSelectObject(&NewPen);for(int i1=0;i1MoveTo(Round(P2i10.x),Round(P2i10.y);for(int j1=1;j1LineTo(Round(P2i1j1.x),Round(P2i1j1.y

10、);for(int j2=0;j2MoveTo(Round(P20j2.x),Round(P20j2.y);for(int i2=1;i2LineTo(Round(P2i2j2.x),Round(P2i2j2.y);pDC-SelectObject(pOldPen);NewPen.DeleteObject();/ 以上代码绘制控制多边形5）控制顶点标注序号CString str; pDC-SetTextColor(RGB(0,0,255); for(int i3=0;i34;i3+)for(int j3=0;j3TextOut(Round(P2i3j3.x+10),Round(P2i3j3.y

11、+2),str);/ 以上代码实现“控制顶点标注序号”6）“曲面”生成CPen PenRed(PS_SOLID,1,RGB(255,0,0);/定义红色笔 pDC-SelectObject(&PenRed);/选择红色笔绘制曲面double dt1=0.01,dt2=0.01;double x,y,u,v; double BU03,BU13,BU23,BU33;double BV03,BV13,BV23,BV33;/U，V两个方向，三次曲面，共8个基函数for(u=0;u=1;u=u+dt1)for(v=0;vMoveTo(Round(x),Round(y);elsepDC-LineTo(Round(x),Round(y);/以上双重循环程序u=0、u=dt、u=2dtu=1的共1/dt+1条“纵向”曲线段for(v=0;v=1;v=v+dt2)for(u=0;uMoveTo(Round(x),Round(y);elsepDC-LineTo(Round(x),Round(y);/此双重循环程序v=0、v=dt、v=2dtv=1的共1/dt+1条“横向”曲线段/上述两批“纵向”及“横向”曲线段，犹如“织布”般“纵横交错”,dt