非均匀三次B样条曲线插值的Jacobi-PIA算法.pdfVIP

第 27卷 第 3期 计算机辅助设计与图形学学报 V_01．27NO．3 2015年 3月 JournalofComputer-AidedDesign ComputerGraphics Mar．2O15 非均匀三次B样条曲线插值的Jacobi．PIA算法 刘晓艳，邓重阳 (杭州电子科技大学理学院 杭州 310018) (dey@hdu．edu．cn) 摘 要：为了求解非均匀三次 B样条曲线插值问题，基于解线性方程组的Jacobi迭代方法提出一种渐进迭代插值算 法—— Jacobi．PIA算法．该算法 以待插值点为初始控制多边形得到第 0层的三次B样条曲线，递归地求得插值给定 点集的三次B样条曲线；在每个迭代过程中，定义待插值点与第 k层的三次B样条 曲线上对应点的差向量乘以该点 对应的B样条系数的倒数为偏移向量，第 k层的控制顶点加上对应的偏移向量得到第 k+1层的三次B样条 曲线的控 制顶点．由于Jacobi—PIA算法在更新控制顶点时减少了一个减法运算，因而运算量更少．理论分析表明该算法是收敛 的．数值算例结果表明，Jacobi．P1A算法的收敛速度优于经典的渐进迭代插值算法，与最优权因子对应的带权渐进迭 代插值算法基本相 同． 关键词：迭代算法；曲线插值；非均匀三次B样条 中图法分类号：0245 Jacobi--PIAAlgorithm forNon--uniform CubicB--SplineCurveInterpolation LiuXiaoyanandDengChongyang (SchoolofScience,HangzhouDianziUniversity,Hangzhou 310018) Abstract：BasedontheJacobiiterativemethodforsolvingthesystem oflinearequations，weproposea progressiveiterativeapproximationmethodforinterpolatingasetofpointsbynon-·uniform cubicB--spline curves，(abbr．Jacobi—PIA)．InJacobi—PIA，thecontrolpointsoftheinitialcubicB—splinecurvearesetasthe pointstobeinterpolated，thencontrolpointsoftheinterpolationcubicB-splinecurvearederivediniteration manner．Ineachiteration，wedefinethedisplacementvectorasthedifferenceofthepointtobeinteprolated anditscorrespondingpointonthecubicB—splinecurveinthepreviouslevel，andthenthecontrolpointsof thecurrentlevelarederivedby thoseofpreviousleveladdingthecorresponding displacementvector． Jacobi—PIA algorithm haslesscomputationthanexistedPIA algorithm becauseofcuttingdownasubtraction in updating the controlpoints．Theoretical analysis shows thatJacobi—PIA algorithm is convergent． Numericalexamplesshow thattherateofconvergenceofJacobi—PIA isfasterthanthatofPIA algorithm， andisalmostequaltothatoftheweightedPIA algorithm withthebestweight． Keywords：iterativealgorithm；curveinterpolation；non--uniform cubicB--splinecurve 用 B样条曲线曲面插值给定点集是 CA