曲线与曲面

⎰∑====-∞→∞→t n i i i n n dt dt t dP P P n L c 011)(lim )(lim T dt dc dt dp dt dp dt dc dt dpdt dp T dc dp c T dtdpdt dp dt dp if t dcdpT c P dc dp c P t P cP t C r dtdpt r if P t P t t P P c ⋅=⇒===±==⇒≠→=∆∆=→∆=∆⇒→∆⇒⇒→∆⇒=∆-∆+=∆→∆对比上两式：对于参数对于一般参数＝单位切矢量，则：为曲线参数，即如选择设弧长为点切线方向的方向为点有切线弦长，:10:1lim )()(C 00)()(0曲线过于平坦如果切矢量远小于弦长曲线过顶点或回转倍如果切矢量是弦长的：切矢量：单位切矢量明确概念：⇒⇒n dtdpdc dp)()()()(0)()(0c P P t P P t c c t t c c dt t dP dt dc dt dt t dP c t ==⇒=⇒=⇒>=⇒=⎰可以用弧长参数表示曲线存在反函数的单调函数是关于参数k dc zd dc y d dc x d k c p dcp d k c p dc dp T dc dTT 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 1lim ,2/1222222222''22'212100021210212121212121=⎥⎦⎤⎢⎣⎡++=⇒==⇒===⋅∆∆=∆∆=∴=⋅∆∆=∆⋅∆=∆∆⇒=∆=∆⋂→∆→∆→∆⋂→∆⋂⋂ρϕϕϕϕϕ曲率半径：又又： 为单位主法线矢量点的法线）与主法线（通过曲率中心的法线平行垂直的平面）法平面（通过该点与在同一平面点为中心向外辐射），以曲线某点有一束法线（为单位法矢量为法矢量，法矢量的矢量垂直单位切矢量对于空间的参数曲线：为曲率矢量，模为＝＝＝平行的单位矢量记为与垂直与线的切线方向单位切矢量，方向为曲N R N T R N T N 1KN N NT :⇒⇒⇒⇒⇒⇒⇒⇒⋅⇒⇒⇒KN K KN dc dT dc dT dcdT dc dT T ρ⎪⎩⎪⎨⎧⇒⇒⇒⇒⨯=⨯=⨯=⇒⇒⋅=化直平面决定的平面法平面决定的平面密切平面决定的平面通过定点标系，下列关系成立：组成互相垂直的直角坐为单位副法线矢量其中副法线的法线和垂直于设BT NB TN RT B N B N T N T B B N T B N T N T B ,,,,第六章 曲线与曲面一、 曲线、曲面参数表示的基础知识1、 参数曲线的定义：切矢量、法矢量、曲率、挠率 §切矢量：坐标变量关于参数的变化率；弧长：对正则曲线P （t ）参数从0到T 的弧长；§曲率：曲线的弯曲变化率；§法矢量,,R C R B 0意点处的挠率等于的充要条件：曲线上任定理：曲线的平面曲线平均挠率为弧长处密切平面的夹角为参数邻域内取曲线上点，点的弧参数设曲线挠率反映曲线的钮挠性质不是常数非平面曲线副法矢量不变平面密切平面就是曲线所在平面曲线=∆∆⇒⎭⎬⎫∆∆∆+⇒⇒⇒=⇒⇒⇒c RQ c Q R C C Q dcdBdcdBθθ设给定函数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)的线性插值函数两点式点斜式21211212112121)()(y x x x x y x x 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）曲率变化较小；例：平面上的三次参数曲线段： 10332210332210≤≤+++=+++=t t b t b t b b y t a t a t a a x §相应拐点方程：),,(),,(),,(,022*******b 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[)()()(012233012233012233∈⎪⎩⎪⎨⎧+++=+++=+++=t a t a t a t a t z at a t a t a t y a t a t a t a t x z z z z y y y y xx x x矢量形式：式一012233)(a t a t a t a t P +++=令：[]123t t t T = []Tx xa a a a C 0123=x C t t t x ⋅=]0123[)(2'给定边界条件：x x xx R x R x P x P x 1'0'10)1()0()1()0(====代入上两式，可得：x x xx xx xx C R C R C P C P ⋅=⋅=⋅=⋅=]0123[]0100[]1111[]1000[1010 矩阵表示为：x C R R P P ⋅⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡01230100111110001010 对上述方程两端乘以4*4的矩阵，可得：x x R R P P C ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⋅⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=10100001010012331122 令：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=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 kk t t kk ,1,0,)()(0==+-== 记号: n C0阶参数连续性，记作C0连续性，是指曲线的几何位置连接，即)()(0)1()1(1++=i i i i t p t p3、 参数曲线的代数形式和几何形式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 p2阶参数连续性：记作C2连续性，指两个相邻曲线段的方程在相交点处具有相同的一阶和二阶导数； ▪ 几何连续性0阶几何连续性：记作G0连续性，与0阶参数连续性的定义相同，满足： )()(0)1()1(1++=i i i i t p t p1阶几何连续性：记作G1连续性，指一阶导数在相邻段的交点处成比例2阶几何连续性：记作G2连续性，指相邻曲线段在交点处其一阶和二阶导数均成比例5、重新参数化、有理参数多项式曲线 重新参数化目的：改变生成曲线的参数间隔，不改变曲线的形状与位置； 最简单的形式：曲线的走向变重新参数化的一般形式：1）端点位置矢量不变；2）切矢量：重新参数化后的曲线与原来的曲线的几何系数之间存在一定的比例关系；参数曲线的截断参数曲线的分割：参数曲线被分割成具有任意长度的n 条新的参数曲线，求第i 条曲线的几何系数；参数曲线被等分成n 段曲线，即参数变量的间隔相等；参数曲线的复合：把几条参数曲线段连接在一起，形成一条复合的参数曲线；构造新曲线的几何系数B 有理参数多项式曲线▪ 目的：为更方便地控制曲线的形状；▪ 原理：基于齐次坐标的概念，产生了用有理参数多项式构造曲线、曲面； ▪ 有理参数的优点：- (1) 具有几何和透视投影变换不变性；▪ 无理多项式表示的曲线：▪ 生成曲线的离散点；▪ 对这些离散点做透视投影变换，得到要求的曲线；▪ 有理多项式表示的曲线：▪ 对定义的曲线的控制点做透视投影变换； ▪ 用变换后的控制点生成要求的曲线；- (2) 可精确的表示圆锥曲线、二次曲面，进而可统一几何造型算法；▪ 研究内容包括：- 参数多项式曲线的代数形式与几何形式； - 参数多项式曲线的矩阵表示； - 参数多项式曲线的生成；10),()(0,≤≤=∑=t t B P t C ni n i i ∑==ni i,k i (u)N P C(u)0()()1k 1,11111,,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 ut u N t t t u u N 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 样条曲线的表达式：参数说明：k+1是曲线的阶数，k 为B样条曲线的次数，曲线在连接点处具有(k-1)阶连续；),,1,0(n i P i =],[)1()1()(1,1,1,++++-∈+=-=i k ik k i k i k i t t u u N u N u N ()()1k 1,11111,,1i 1,t )()(0t1)(+-++++++-++≤≤--+--=⎩⎨⎧<≤=n k i i k i k i k i i k i i k i i i t u u N t t ut u N t t t u u N t u u N 其它若()()()()()()()()()()()()u N u u N u u N u N uu 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,0373)3(363)2(353)1(343-+-=-+-=-+-=-+=)(21)(),(21)(3210P P end p P P start p +=+=2301)(,)(P P end p P P start p -='-=' 是节点值，构成了k 次B 样条曲线的节点矢量，节点是非减序列； 且：节点矢量：分为三种类型：均匀的，均匀非周期的和非均匀的；- 节点沿参数轴均匀等距分布，即 =常数时，→均匀B 样条函数； - 节点沿参数轴分布不等距，即 ≠常数时，→非均匀B 样条函数。