计算机图形学(Hermite-curve).doc

计算机图形学
课程名：计算机图形学 　　　主讲教师：陈学工
教材：《计算机图形学基础》（第1版）/唐泽圣 清华大学出版社，2001年。
参考书：《计算机图形学》（第1版）/倪明田 北京大学出版社，2000年。

第1讲二维图形
—抛物线样条曲线和Hermite曲线
基本点：自由曲线（曲面）；曲线拟合；曲线插值；抛物线样条曲线；Hermite曲线。
重点：Hermite曲线。
难点：向量和矩阵运算。
疑点：
形状比较复杂、不能用二次方程来表示的曲线（曲面）称为自由曲线（曲面），也称为复杂曲线(曲面)，通常以三次参数方程来表示。给定一个点列，用该点列来构造曲线的方法称为曲线拟合。已知曲线上的一个点列，求曲线上的其他点的方法称为曲线插值。
1．抛物线参数样条曲线
样条是一根富有弹性的细木条或类似物，其两端连接着起固定作用的压铁。通过调整样条两端的压铁可以改变样条的形态，它是手工绘制自由曲线的一种工具。沿着样条绘制的曲线称为样条曲线。样条曲线可以表示为参数多项式曲线或分段参数多项式曲线。
给定一个点列P1,P2,…,Pn，对相邻三个点用抛物线来拟合，相邻抛物线在公共区间内用权函数t进行调配，所得到的曲线称为抛物线参数样条曲线。该曲线向量形式可表示为：

其中，Si为Pi,Pi+1,Pi+2三个点决定的抛物线。
2．Hermite曲线
Hermite曲线是以曲线的两个端点P0、P1和端点处的切向矢量R0、R1为边界条件的三次参数曲线。空间自由曲线三次参数方程的一般形式可表示为：
或

其中Q(t)=(x(t),y(t),z(t)),a=(ax,ay,az),b=(bx,by,bz),c=(cx,cy,cz), d=(dx,dy,dz)。令T=[t3 t2 t 1], ,那么Q（t）可以表示为Q（t）=TMhGh。Mh是由初始条件确定的一个矩阵，该矩阵不唯一，通常取。Fh(t)=TMh确定的一组函数称为调和（或基）函数。
对于端点处坐标和切线方向都相同的Hermite曲线，它们的形状随着切向矢量的长度的变化而变化。此外，如果取与z相关的系数为0，则得到平面上的Hermite曲线。
Hermite曲线通过端点和端点处的切向矢量来控制曲线的形状，一条单独的Hermite曲线不适合于用来拟合一个空间点列，它主要用在样条曲线中用来表示其中的某一段曲线。
第2讲 二维图形—三次参数样条曲线
基本点：三次参数样条曲线；约束条件。
重点：三次参数样条曲线。
难点：三次参数样条曲线。
疑点：

给定参数节点{ti}(i=0,1,…,n)和相应的点列{Pi}(i=0,1,…,n),求二阶导数连续的分段三次参数多项式曲线P(t)( )，使得该曲线能够拟合已知的点列，即P（ti）=Pi(i=0,1,…,n)。P(t)称为三次样条曲线。
三次样条曲线P(t)在每一个区间[ti,ti+1]()上的分段曲线记为Pi(t)。如果用Hermite曲线表示Pi(t)，则可以用Hermite 三次参数曲线来描述传统的样条曲线。
1．Hermite曲线的二阶导数形式
由前一讲可知，Hermite曲线可以表示为：
Q(t)=Fh1(t)Q(0)+Fh2(t)Q(1)+Fh3(t)Q’(0)+Fh4(t)Q’(1)
如果Hermite曲线在端点处二阶导数已知，则端点处的一阶导数可以用端点处的二阶导数来表示，最后可以把Hermite曲线用二阶导数形式表示为：

2．三次参数样条曲线
设有点列{Pi}(i=0,1,…,n)，用Hermite三次参数曲线将相邻点连接起来，使得最终的曲线在已知点处具有连续的二阶导数，该曲线是一条三次样条曲线。
Pi+1是前后两条Hermite曲线的端点，根据这两曲线的方程分别求该点处的二阶导数得到

(t=1)

(t=0)
由于曲线具有二阶连续导数，由上面两式得到

对于已知的点列，可以得到n-2个类似于上式的方程，其中未知的一阶导数有n个。如果再给定三次样条曲线在起点和终点处的切向矢量或二阶导数，则可以求出各一阶导数和分段Hermite曲线，进而可以确定三次样条曲线。
三次样条曲线在两端点处的边界条件是根据问题的物理要求来确定的。常用的端点约束条件有如下三种：
（1）自由端
在两端点处的二阶导数为零，由此得到

（2）夹持端
在端点处的切向矢量已知，即，，其中E1，En为单位切向矢量。
（3）抛物端
假定分段曲线的最初和最未段为抛物线，即三次项的系数为0，则在这两段曲线上二阶导数为常数，由此可以得到

例：已知平面上三点P1(0,0)，P2(1,1)和P2(2,0)，用自由端约束条件计算三次Hermite样条曲线Q(t)使得它拟合点列{P1，P2，P3}，并将Q(t)表示为参数序列{0，1，2}上的分段参数多项式函数。
解：由Hermite样条曲线二阶导数的连续性及自由端约束条件得到

解方程组得到

由控制点的坐标得到
由Hermite曲线方程

得到

所求曲线方程为：

4