空间曲线的高阶几何Hermite插值-应用数学专业毕业论文.docxVIP

摘要本文在给定空间曲线两个端点的位置、切方向、曲率法 摘要 本文在给定空间曲线两个端点的位置、切方向、曲率法 向量和挠率的情况下，用参数化五次B6zier曲线来对这条空 间曲线进行几何Hermite插值。我们证明了插值问题局部可 解，解有两个自由度，而且本文还给出了一种确定这两个自 由度的方法。并证明了在适当的假设下插值有6阶逼近度。 关键词：几何连续；插值；逼近度；曲率法向量；挠率；TGHI 中图法分类号TP391 4 4 bstract This paper uses parametric quinary B6zier curve to interpolate a given spacial curve，whose position，tangent direction】the curva— ture vector and torsion are prescribed at each knot．We prove that the interpolant exists locally with two degree of freedom，and we give a method to fix the two freedom．Moreover，we prove that Un— der appropriate assumptions the interpolant is 6th order accurate Keywords：Geometric Continuous；Interpolant；Accuracy；Curva— ture vector；Torsion；TGHI 3 致谢 致谢 在这里我要感谢我的导师曹沅老师，本文是在他的悉心 指导下完成的。并且，在研究生阶段的学习中，曹老师的谆 谆教诲让我获益匪浅，老师的无微关怀也让我将永远记住这 段求学的快乐时光。同时，也要感谢吴宗敏老师给这篇文章 的初稿提出宝贵意见，吴老师的帮助和教导也深刻的影响着 我。 我还要感谢在讨论班上，史建红、葛焰明、阎艳、李成、 王建华、陈荣华、熊正超、张立伟、张卫祥、杨艳、卜天奇等 同学提出的宝贵意见。我还想说的是，跟曹老师、吴老师以 及以上同学和朋友们在一起学习的时间我将永远记得，因为 这是我最难忘的时光。 预备知识§l 预备知识 §l Hermite插值 曲线插值是CAGD中一类基本而广泛的闽题，而多项式插值是所有插值中 最基本的、应用最广泛的方法(最早的多项式插值方法可能应该归功于I．New— ton)。曲线设计的一个最普通的问题就是点插值问题：给定一列相对于参数t。 的数据点Pi∈R4，i=0， ，n，需要找到一条曲线通过Pt。解决这个问题的 最早的方法就是找到一条多项式曲线通过这些点，这个多项式必须满足下面这些 插值条件： p(如)=P。；i=0，．．．，n． 对这个问题已经有多种解决方案，比如： Aitken’s Algorithm，Lagrange Polyno— mials，the Vandermonde Approach等等。在实际应用中，多项式播值不一定只能 插值点上的函数值，也可以插值别的一些信息，比如某些点上的导数信息，这样 就得到一个比一般的多项式插值更有用的插值方法，我们称之为Hermite插值。 在这里我们讨论两个节点的一般情形的Hermite插值：给定两点Po，P，，以及 两端的各阶导数玮，p{，i=1．m，J=1．m。(不失一股陛，我们设to=0，tl=h) 我们需要找到一个多项式P插值这些数据； p(o)=Po， p(h)=Pl， p(‘’(o)=pj，i=1．．m p(，)(^)=p{，J=1．．n 其中p(4】指P的i阶导数。我们记N=仇+n+1，上式共有N+1个插值条 件，它可以确定N+1个系数，所以一般情况下满足这些插值条件的多项式P的 最小次数为N，而且该多项式唯一。 2 我们把多项式P写成B6zier形式，则必须确定N+1个B6zier控制点bo， 我们把多项式P写成B6zier形式，则必须确定N+1个B6zier控制点bo， ，b_v。 其中两个能够很快确定： bo=Po，bN=P1 对余下的控制点，我们利用多项式Bgzier形式的导数特点： Ⅳ ，N-k p忙’(t)=(E b；B{。_Ⅳ(t))∽’=去∑△。b。B。，N一≈(t) i=0 7。i=0 其中，县蛹(t)=(?)(；)‘(1一；)”2，A％：=Ak-1b件l—Ak-1b：，△ob，=b‘。由 此，我们很容易解出b·， ，bⅣ一t，从而确定了这条多项式插值曲线。 设原睦线为，(t)，则上面的Hermite插值有如下的误差估计： ，。)一p(t)=币矿{可，(Ⅳ十u嬉)t“十1。～危)“+I，≤∈【0，叫 于是， l