第六章 曲与曲面

?∑====-∞→∞→t n i i i n n dt dt t dP P P n L c 01

1)

(lim )(lim T dt dc dt dp dt dp dt dc dt dp

dt dp T dc dp c T dt

dp

dt dp dt dp if t dc

dp

T c P dc dp c P t P c

P t C r dt

dp

t r if P t P t t P P c ?=?==

=±==?≠→=??=→?=??→???→??=?-?+=?→?对比上两式：对于参数对于一般参数＝单位切矢量，则：为曲线参数，即如选择设弧长为点切线方向的方向为点有切线弦长

，:1

0:1lim )

()(C 00)()(0曲线过于平坦

如果切矢量远小于弦长曲线过顶点或回转

倍如果切矢量是弦长的：切矢量：单位切矢量明确概念：??n dt

dp

dc dp

)()()()(0)()(0

c P P t P P t c c t t c c dt t dP dt dc dt dt t dP c t ==?=?=?>=?=?可以用弧长参数表示曲线存在反函数的单调函数是关于参数dc z

d dc y d dc x d k c p dc

p d k c p dc dp T dc dT

T T T T c T c k T T T

T T T T T c T T T c T T c T T T T T c c c 1

)()()()()()lim ()lim (lim 1

lim ,2/1222

222222''2

2'

21210002

12102

12121212121??????++=?==?===???=??=∴=?

??=???=???=?=??

→?→?→??

→??

??????又又：ΘΘ第六章 曲线与曲面

一、 曲线、曲面参数表示的基础知识

1、 参数曲线的定义：切矢量、法矢量、曲率、挠率 §切矢量：坐标变量关于参数的变化率；

弧长：对正则曲线P （t ）参数从0到T 的弧长；

§曲率：曲线的弯曲变化率；

为单位主法线矢量点的法线）与主法线（通过曲率中心的法线

平行垂直的平面）法平面（通过该点与在同一平面点为中心向外辐射），以曲线某点有一束法线（为单位法矢量

为法矢量，法矢量的矢量垂直单位切矢量对于空间的参数曲线：

为曲率矢量，模为＝＝＝平行的单位矢量记为与垂直与线的切线方向单位切矢量，方向为曲N R N T R N T N 1

KN N N T :????????????KN K KN dc dT dc dT

dc

dT dc dT T ρ

??

?

???????=?=?=???=化直平面决定的平面法平面决定的平面密切平面决定的平面通过定点标系，下列关系成立：

组成互相垂直的直角坐为单位副法线矢量其中副法线的法线和垂直于设BT NB TN R

T B N B N T N T B B N T B N T N T B ,,,,0

,,R C R B 0

意点处的挠率等于的充要条件：曲线上任定理：曲线的平面曲线平均挠率为弧长处密切平面的夹角为参数邻域内取曲线上点，点的弧参数设曲线挠率反映曲线的钮挠性质不是常数非平面曲线副法矢量不变平面密切平面就是曲线所在平面曲线=????

?????+???=???c RQ c Q R C C Q dc

dB

dc

dB θ

θ设给定函数f(x)在两个点的值：y1=f(x1),y2=f(x2); 要求：线性函数b ax x y +==)(?近似代替)(x f y =；

如选择a, b ，使2211)(,)(y x y x ==??则)(x ?为f(x)的线性插值函数

两点式

点斜式2

1

1211

21

21)

()(y x x y x x x x x x y y y x -+-=

---+=?§法矢量

§挠率

2、 插值、逼近、拟合与光顺 -函数逼近的重要方法；

函数逼近问题与插值问题； 插值函数；

常用方法：线性插值，抛物线插值 线性插值：

抛物线插值（二次插值）：

§设已知f(x)在三个互异点x1,x2,x3的函数值为y1,y2,y3;

§要求：构造函数c bx ax x ++=2

)(?使该函数在节点Xi 处与f(x)在该点处的值相等； §求解：构造线性方程组，求参数a, b, c ，即构造了插值函数

§逼近：-插值的问题？

§型值点太多时→构造插值函数困难； §型值点多→误差大；

-解决：选择低阶函数,在某种意义上逼近型值点→最佳逼近 §常用方法：最小二乘法

-最小二乘法：逼近的效果由各点偏差的平方和最小或加权的方差最小；

§拟合：曲线、曲面的设计过程中，用插值或逼近方法是生成的曲线、曲面达到某些设计要求。

§光顺：拐点不要太多-曲线的拐点太多→视觉效果差； -平面曲线的相对光顺条件： §1）具有二阶几何连续；

§2）不存在多余拐点和奇异点； §3）曲率变化较小；

例：平面上的三次参数曲线段： 103

322103

32210≤≤+++=+++=t t b t b t b b y t a t a t a a x §相应拐点方程：

),,(),,(),,(,0223

213212b b b a a a r q p r qt pt ?==+-，式中

§若p ≠0，则可以构造表达式：p r p q I /2)/(2-= §当I>0时，相应曲线有两个实拐点； §当I ＝0时，曲线上出现一个尖点； §当I<0时，曲线上会出现一个二重点；

考虑三次参数曲线的代数形式：]1,0[)()()(012

233012233012233∈?????+++=+++=+++=t a t a t a t a t z a

t a t a t a t y a t a t a t a t x z

z z z y y y y x

x x x

矢量形式：式一0

12

233)(a t a t a t a t P +++=

令：[

]

123

t t t

T = []T

x x

a a a a C 0123=

x C t t t x ?=]0123[)(2'

给定边界条件：x x x

x R x R x P x P x 1'0'10)1()0()1()0(====

代入上两式，可得：

x x x

x x x x

x C R C R C P C P ?=?=?=?=]0123[]0100[]1111[]1000[1010 矩阵表示为：x C R R P P ??

?

???

?

?

??

???=?????

???????0123010011111000

101

0 对上述方程两端乘以4*4的矩阵，可得：

x x R R P P C ?

?

??

??????????????????

???----=101000010100123311

22 令：?????

???????=????????????----=1010;0001010012331122R R P P G M h x h ? 则称Mh 为Hermite 矩阵（常数阵），Gh 为几何矢量； ? X(t)=T*Mh*Ghx →P(t)=T*Mh*Gh

? 所以，只要给定Gh ，就可在0≤t ≤1范围内求出P(t)； ?

对于不同初始条件，Gh 不同，T ，Mh 均是相同的；

连续性条件：假定参数曲线段pi 以参数形式进行描述：]t ,[t t )(i1i0∈=t p p i i 参数连续性 几何连续性

? 参数连续性（传统的、严格的连续性）

称曲线P=P(t)在 0t t =处n 阶参数连续，如果它在0t 处n 阶左右导数存在，并且满足

n k dt t P d dt t P d t t k

k t t k

k Λ,1,0,)()(0

==

+

-

== 记号: n C

0阶参数连续性，记作C0连续性，是指曲线的几何位置连接，即)()(0)1()1(1++=i i i i t p t p

3、 参数曲线的代数形式和几何形式

4、 连续性定义

? 参数连续性

1阶参数连续性：记作C1连续性，指代表两个相邻曲线段的方程在相交点处有相同的一阶导数：

)

()()()(0)1()1(10)1()1(1++++'='=i i i i i i i i t p t p t p t p

2阶参数连续性：记作C2连续性，指两个相邻曲线段的方程在相交点处具有相同的一阶和二阶导数； ? 几何连续性

0阶几何连续性：记作G0连续性，与0阶参数连续性的定义相同，满足： )()(0)1()1(1++=i i i i t p t p

1阶几何连续性：记作G1连续性，指一阶导数在相邻段的交点处成比例

2阶几何连续性：记作G2连续性，指相邻曲线段在交点处其一阶和二阶导数均成比例

5、重新参数化、有理参数多项式曲线 重新参数化

目的：改变生成曲线的参数间隔，不改变曲线的形状与位置； 最简单的形式：曲线的走向变

重新参数化的一般形式：1）端点位置矢量不变；2）切矢量：重新参数化后的曲线与原来的曲线的几何系数之间存在一定的比例关系；

参数曲线的截断

参数曲线的分割：参数曲线被分割成具有任意长度的n 条新的参数曲线，求第i 条曲线的几何系数；参数曲线被等分成n 段曲线，即参数变量的间隔相等；

参数曲线的复合：把几条参数曲线段连接在一起，形成一条复合的参数曲线；构造新曲线的几何系数B 有理参数多项式曲线

? 目的：为更方便地控制曲线的形状；

? 原理：基于齐次坐标的概念，产生了用有理参数多项式构造曲线、曲面； ? 有理参数的优点：

- (1) 具有几何和透视投影变换不变性；

? 无理多项式表示的曲线：

? 生成曲线的离散点；

? 对这些离散点做透视投影变换，得到要求的曲线；

? 有理多项式表示的曲线：

? 对定义的曲线的控制点做透视投影变换； ? 用变换后的控制点生成要求的曲线；

- (2) 可精确的表示圆锥曲线、二次曲面，进而可统一几何造型算法；

? 研究内容包括：

- 参数多项式曲线的代数形式与几何形式； - 参数多项式曲线的矩阵表示； - 参数多项式曲线的生成；

10),()(0,≤≤=∑

=t t B P t C n

i n i i ∑

==n

i i,k i (u)N P C(u)0

11i 1

,)(

0t 1)(+++--???<≤=k i i i i u

t t u t u u N 其它若二、 常用的参数曲线

1、 Bezier 曲线 ? 定义：

- 一种以逼近为基础的参数曲线；

- 由一组折线集，或Bezier 特征多边形定义； - 曲线的起点、终点与多边形起点、终点重合；

- 多边形的第一个边与最后一个边表示了曲线在起点和终点的切矢量方向； - 形状趋于特征多边形的形状；

- 给定空间n+1个点的位置矢量：Pi ，则Bezier 曲线各点坐标的插值公式：

? Bezier 曲线的性质： 1）端点性质：

A)端点位置矢量：Bezier 曲线的起点、终点与其相应的特征多边形的起点、终点重合；

B)切矢量：Bezier 曲线的起点、终点的切线方向与其相应的特征多边形的第一条边及最后一条边的走向

一致；

C)曲率：Bezier 曲线在端点处的r 阶导数，只与（r+1）个相邻点有关，与更远的点无关； D)r 阶导函数的差分表示：N 次Bezier 曲线的r 阶导函数可用差分公式表示为： 2）对称性：Bezier 曲线及其特征多边形在起点处的几何性质与终点处相同；

若保持原Bezier 曲线的全部定点位置不变，仅把次序颠倒，形成新的顶点； 则新Bezier 曲线形状不变，只是走向相反；

3）凸包性：1）说明当t 在0与1区间变化时，对某个t 值，C （t ）是特征多边形各项点Pi 的加权平均，权

因子依次是Bi,n(t);

2）在几何图形上，Bezier 曲线是Pi 各点的凸线性组合，并且各点均落在特征多边形的凸包之中；

4）几何不变性：几何特性不随一定的坐标变换而变化的性质

Bezier 曲线的位置与形状仅与特征多边形的定点位置有关，不依赖坐标系的选择；

5）变差缩减性：如Bezier 曲线的特征多边形是一个平面图形，则直线与曲线的交点个数 ≤ 该直线和特征

多边形的交点个数→变差缩减性；

说明Bezier 曲线比特征多边形的波动小→Bezier 曲线比特征多边形所在的折线更光顺；

2、 B 样条曲线

目的：解决Bezier 曲线的不足（1972年，Gordon,Riesenfeld 扩展Bezier 曲线）；

1)控制多边形的顶点个数决定了Bezier 曲线的阶次,n 较大时特征多边形对曲线的控制减弱; 2)调和函数在整个区间内均不为零→不能作局部修改； 方法：用B 样条函数代替Bernstein 函数，从而：

1）改进了Bezier 特征多边形与Bernstein 多项式次数相关的问题； 2）克服了Bezier 曲线整体逼近的缺点； 均匀B 样条函数的定义：

已知有n+1个控制点的特征多边形，其顶点为： 则K 次（K ＋1阶）的B 样条曲线的表达式：

),,1,0(n i P i Λ=

],[)1()1()(1

,1,1,++++-∈+=-=i k i

k k i k i k i t t u u N u N u N ()()1

k 1,11

111,,1i 1,t )()(

0t

1)(+-++++++-++≤≤--+--=

???<≤=n k i i k i k i k i i k i i k i i i t u u N t t u

t u N t t t u u N t u u N 其它若()()()

()()()()()()

()()()

u N u u N u u N u N u

u N u u N u N u

u N u u N u N u u N u u N 2,42,33,32,32,23,22,22,13,12,12,03,03

73)3(363)2(353)1(3

43-+-=-+-=-+-=-+=)(2

1)(),(21)(3

210P P end p P P start p +=+=2

301)(,)(P P end p P P start p -='-='参数说明：k+1是曲线的阶数，k 为B 样条曲线的次数，曲线在连接点处具有(k-1)阶连续； 是节点值，

构成了k 次B 样条曲线的节点矢量，节点是非减序列； 且：

节点矢量：分为三种类型：均匀的，均匀非周期的和非均匀的；

- 节点沿参数轴均匀等距分布，即 =常数时，→均匀B 样条函数； - 节点沿参数轴分布不等距，即 ≠常数时，→非均匀B 样条函数。

均匀周期性B 样条曲线

- T=(-2,-1.5,-1,-0.5,0,0.5,1,1.5,2) - T=(0,1,2,3,4,5,6,7)

- 均匀B 样条的基函数呈周期性：

均匀周期性二次（三阶）B 样条曲线

取k=2，n=3，则k+n ＋1=6，不妨设节点矢量为：T=(0,1,2,3,4,5,6)：

曲线的起点和终点值：

均匀二次B 样条曲线起点和终点处的导数：

结论：

对于由任意数目的控制点构造的二次周期性B 样条曲线，曲线的起点位于头两个控制点之间，终点位于最后两个控制点之间；

对于高次多项式，起点和终点是m-1个控制点的加权平均值点。若某一控制点出现多次，样条曲线会更加接近该点。

由不同节点矢量构成的均匀B 样条函数所描绘的形状相同，可以看成是一个B 样条函数的简单平移； 在构造每个线段时，采用均匀B 样条函数比用非均匀B 样条函数工作量小，且外形设计的效果差别不大。

i t ],,,[1210++=k L t t t T Λi i t t -+1i i t t -+1k n L -=

;

)()(/)()(0,0,0,∑

∑∑=====n

i k i i n i k i i n i k i i i u R P u N W u N P W u C 1r 0 11

120<≤-===r r w w w ? B 样条曲线的性质 1．局部可控性

2. B 样条的凸组合性质

B 样条的凸组合性和B 样条基函数的数值均大于或等于0保证了B 样条曲线的凸包性，即B 样条曲线必处在控制多边形所形成的凸包之内。 3. 连续性 4. 导数；

B 样条曲线的导数可用其低阶的B 样条基函数和定点矢量的差商序列的线性组合表示。 可以证明，k 次B 样条曲线段之间达到k-1次的连续性； 5. 几何不变性：形状和位置与坐标系的选择无关；

6. 变差缩减性：平面内任一直线与曲线的交点个数不多于该直线与特征多边形的交点个数；

7. 造型的灵活性：

用B 样条曲线可构造直线段、尖点、切线等特殊情况：通过选择控制点的位置和节点的重复数。

3、非均匀有理B 样条（NURBS ）曲线

NURBS 方法:非均匀有理B 样条（Nonuniform Rational B-Spline ）方法

定义：由分段有理B 样条多项式基函数定义：

例：假定用定义在三个控制顶点和开放均匀的节点矢量上的二次（三阶）B 样条函数来拟合，于是，T=(0,0,0,1,1,1)，取权函数为：

有理B 样条的表达式为：

取不同的r 值得到各种二次曲线

NURBS 优点：

对标准的解析形式和自由曲线、曲面提供了统一的数学表达，无论始解析形状还是自由格式的形状均有统一的表示参数，便于工程数据库的存取和应用；

可通过控制点和权因子来灵活地改变形状；

对插入节点、修改、分割、集合插值等的处理工具比较有力； 具有透视投影变换和仿射变换的不变性；

非有理B 样条、有理及非有理Bezier 曲线、曲面是NURBS 的特例；

NURBS 存在问题：

比一般的曲线、曲面定义方法更费存储空间和处理时间； 权因子选择不当会造成形状畸变； 对搭接、重叠形状的处理相当麻烦；

象点的映射这类算法在NURBS 情况下会变得不稳定。

)(3,2)(3,1)(3,0)

(3,22)(3,11)(3,0011)(t t t t t t N N r

r

N N P N P r

r

N P t C +-++-+=

T H H H H H P P

P P B M B TM t P ]

[,000

10100

123311

22,)('1

'

1

=??

???

??

??

???----==T i i i i Z Z Z Z Z P P P P B M B TM t P ][,0001003303631331,)(3

21+++=????

?

?

??????----==Z B B B B B B B M B TM t P =?

?

????

??????----==,0141030303631331,)(;,)3;,)2;,)11

1111

1B B Z Z B B H H Z Z H H B B Z

H B B Z Z B

B Z Z H H B B H H Z Z B

H Z Z H H B B H H Z Z B B Z Z H H B

Z H H B M M B B M M B B M B M B M B B B M B M M B B M M B B M B M B M B B B M B M M B B M M B B M B M B M B B B M ------==?====?====?==ΘΘΘ，；求，已知：，；求，已知：，；求，已知：4、常用参数曲线的等价表示

Hermite ，Bezier ，B 样条三种非有理参数表示式的相互转换

对于同一条曲线，已知三种表示的一种形式，推导其它两种形式。 设已知三种曲线的表示形式都是几何形式： 则三次Hermite 曲线为：

三次Bezier 曲线为：

三次B 样条曲线为：

三、 常用的参数曲面：

1、 定义——参数曲面的表示：显式、隐式；参数式：从图形学角度，更便于用计算机表示与构造；

2、 重新参数化——参数方向变反

最简单形式：把参数变量u 或/和w 的方向变反；此方式不改变曲面片的形状；表示对一个曲面片改变参数方向的三种情况。

3、Bezier 曲面、B 样条曲面、非均匀有理B 样条（NURBS ）曲面

Bezier 曲面基于Bezier 曲线，可以方便的给出Bezier 曲面的定义及其性质，Bezier 曲线中的一些算法可以扩展到Bezier 曲面的情况。

B 样条曲面基于均匀B 样条的定义和性质，可得B 样条曲面的定义

非均匀有理B 样条（NURBS ）曲面由双参数分段有理多项式曲面定义的NURBS 曲面