单调数据的保形分段线性有理插值.pdfVIP

第二届全国几何设计与计算学术会议，2005合肥 单调数据的保形分段线性有理插值 王 强1，2檀结庆1 (1合肥工业大学计算机与信息学院，安徽合肥230009) (2安徽理工大学数理系，安徽淮南232001) 蔫要：本文构造了含参数的分段线性有理插值函数(分子、分母均为一次多项式)，通过适当选择形状参数，由 此函数产生的曲线一阶连续并且保单调。文中用张量积方法将此结果推广到二元矩形网格上的曲面插值，同时给 出了插值函数的误差估计及数值例子。 关键词。有理插值保形连续双线性曲面 (3)s(石)在【玉，‘】上是单调的， 1引言 则称J(工)是定义在【五，‘】上的分段线性保形有理插 有理曲线曲面作为一类重要的逼近工具，在计 值函数。下面我们给出上述函数的表达式，假设给 算机辅助设计与制造中有着广泛的应用。随着 定数据{(玉，五)，i=1，2，．．．，，ll，其中Z为被插函数在分 NURBS被确定为国际标准后，更奠定了有理函数划点毛上的函数值，此处而而…‘．记 在CAD中的主导地位．近年来，有理样条以及它们 ^=玉+．一‘，△。=(正¨一Z)，嚏，令t=(工一再)，吩。当 在形状控制中的应用已经引起人们的广泛关注【2弓】． 工∈ix,，玉+l】时，tE【0，1】，定义： 实际应用中通常不仅要求插值曲线过型值点，通常 还要求曲线保持由给定的数据所反映的原曲线的形 ％Ll—fJ十q+l‘ 状，如单调性、凹凸性等。 其中珥为正实数，容易验证由(2．1)定义的函数满足 在文献【1】讨论了三次有理样条的约束插值，文 下列插值性质： 献【2，3，4，5】中研究了三次有理样条及其保形性。本文 J(玉)=Z J(玉+1)=Z+l。 构造了分段线性有理插值函数(分子、分母均为一 次多项式)，通过适当选择形状参数，插值函数一阶 我们来讨论该函数的连续性和保形性。 连续且保单调。线性有理插值函数既克服了一次多 2．2连续性与保形-性 项式的不光滑性，同时又不会出现高次多项式产生 假设数据是严格单调递增的，即满足： 震荡。该函数形式简单，仅需给定函数值便可惟一 Z厶…￡或△f0 确定，便于推广到多元情形。文中给出了一元和二 我们知道由(2．1)式定义的j(工)是Co连续的，为了 元插值格式并研究它们的保形性，同时给出了插值 讨论J(工)的连续性，对(2．1)式求导，并化简： 函数的误差估计。 2一元插值函数 帕卜丽而a,af+lAi (2·2) 令 2．1插值函数的构造 砖D(再+o)=嘏(玉-o) 定义2．1如果函数s(工)满足条件： 得到 f=2，3，…n一1(2．3) ％+l=(△，一。／A，)q—I i=l，2，…，n： (I)j(五)=Z， 当参数嘭满足(2．