引en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

优化张力参数与边界条件的平面三次Cardinal样条

  • 李军成1

    机 构：湖南人文科技学院 数学与金融学院，湖南 娄底 417000

    邮 箱：lijuncheng82@126.com.

    作者简介：李军成(1982—)，ORCID: http://orcid.org/ 0000-0002-1904-4068，男，博士，副教授，主要从事计算机辅助几何设计及其应用研究, E-mail：lijuncheng82@126.com.

    ×
  • 刘成志1

    机 构：湖南人文科技学院 数学与金融学院，湖南 娄底 417000

    ×
  • 易叶青2

    机 构：湖南人文科技学院 信息学院，湖南 娄底 417000

    ×

1. 湖南人文科技学院 数学与金融学院，湖南 娄底 417000； 2. 湖南人文科技学院 信息学院，湖南 娄底 417000

中图分类号: TP391

DOI:10.3785/j.issn.1008-9497.2019.02.004

Planar cubic Cardinal spline with tension parameter and boundary condition optimization

  • Juncheng LI1

    Affiliation：College of Mathematics and Finance, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China

    Email：lijuncheng82@126.com.

    ×
  • Chengzhi LIU1

    Affiliation：College of Mathematics and Finance, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China

    ×
  • Yeqing YI2

    Affiliation：College of Information, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China

    ×

1. College of Mathematics and Finance, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China； 2. College of Information, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China

CLC: TP391

DOI:10.3785/j.issn.1008-9497.2019.02.004

参考文献 1
HEARND, BAKERM P.Computer Graphics: C Version [M]. 2nd ed. Upper Saddle River: Prentice Hall, 1996.
查找原文
参考文献 2
YOONH S, PARKT H.Smooth path planning method for autonomous mobile robots using cardinal spline [J]. Transactions of the Korean Institute of Electrical Engineers, 2010, 59 (4): 803-808.
查找原文
参考文献 3
FANANIA, YUNIARTIA, SUCIATIN.Geometric feature extraction of batik image using Cardinal spline curve representation [J]. Telkomnika, 2014, 12(2): 397-404.DOI:10.12928/telkomnika.v12i2.1915
查找原文
参考文献 4
吴晓勤, 严秀坤.三次C-Cardinal样条曲线及曲面[J]. 计算机工程与科学, 2006, 28(2): 48-50.DOI:10.3969/j.issn.1007-130X.2006.02.015
WUX Q, YANX K.The cubic C-Cardinal spline curve and surface [J]. Computer Engineering & Science, 2006, 28(2): 48-50.DOI:10.3969/j.issn.1007-130X.2006.02.015
查找原文
参考文献 5
宋爱平, 陶建明, 易旦萍, 等.可调形三次三角Cardinal插值样条曲线[J]. 计算数学, 2015, 37(1): 34-41.
SONGA P, TAOJ M, YID P, et al.Cubic trigonometric Cardinal interpolation spline with adjustable shape [J]. Mathematica Numerica Sinica, 2015, 37(1): 34-41.
查找原文
参考文献 6
吴晓勤, 韩旭里, 罗善明.三次H-Cardinal样条曲线及曲面[J]. 工程图学学报, 2008, 29(5): 83-88.DOI:10.3969/j.issn.1003-0158.2008.05.016
WUX Q, HANX L, LUOS M.Cubic H-Cardinal spline curves and surfaces [J]. Journal of Engineering Graphics, 2008, 29(5): 83-88.DOI:10.3969/j.issn.1003-0158.2008.05.016
查找原文
参考文献 7
李军成, 刘成志, 易叶青.带形状因子的C2连续五次 Cardinal 样条与Catmull-Rom 样条[J]. 计算机辅助设计与图形学学报, 2016, 28(11): 1821-1831.
LIJ C, LIUC Z, YIY Q.C2 continuous quintic Cardinal spline and Catmull-Rom spline with shape factors [J]. Journal of Computer-Aided Design & Computer Graphics, 2016, 28(11): 1821-1831.
查找原文
参考文献 8
LIJ C.Curvature variation minimizing Cardinal spline curves [J]. Journal of Modeling and Optimization, 2018, 10(1): 31-36.DOI:10.32732/jmo.2018.10.1.31
查找原文
参考文献 9
GONZALEZR C , WOODSR E, EDDINSS L.Digital Image Processing Using MATLAB [M]. Upper Saddle River: Prentice Hall, 2004.DOI:http://dx.doi.org/10.9774/GLEAF.978-1-909493-38-4_2
查找原文
参考文献 10
DEROSET D, BARSKYB A.Geometric continuity, shape parameters, and geometric constructions for Catmull-Rom splines [J]. ACM Transactions on Graphics, 1988(7): 1-41.DOI:10.1145/42188.42265
查找原文
参考文献 11
YUKSELC, SCHAEFERS, KEYSERJ.Parameterization and applications of Catmull–Rom curves [J]. Computer-Aided Design, 2011, 43(7): 747-755.DOI:10.1016/j.cad.2010.08.008
查找原文
参考文献 12
SCARPINITIM, COMMINIELLOD, PARISIR, et al.Nonlinear spline adaptive filtering [J]. Signal Processing, 2013, 93(4): 772-783.DOI:10.1016/j.sigpro.2012.09.021
查找原文
参考文献 13
HOFERM, POTTMANNH.Energy-minimizing splines in manifolds [J]. ACM Transaction on Graphics, 2004, 23(3): 284-293.DOI:10.1145/1015706.1015716
查找原文
参考文献 14
FANW, LEE C H, CHENJ H. A realtime curvature-smooth interpolation scheme and motion planning for CNC machining of short line segments [J]. International Journal of Machine Tools and Manufacture, 2015, 96: 27-46.DOI:10.1016/j.ijmachtools.2015.04.009
查找原文
参考文献 15
LUL Z, JIANGC K, HUQ Q.Planar cubic G1 and quintic G2 Hermite interpolations via curvature variation minimization [J]. Computers & Graphics, 2017, 70: 92-98.DOI:10.1016/j.cag.2017.07.007
查找原文
参考文献 16
JAKLIČG, ŽAGARE.Curvature variation minimizing cubic Hermite interpolants [J]. Applied Mathematics and Computation, 2011, 218(7): 3918-3924. DOI:10.1016/j.amc.2011.09.039
查找原文
参考文献 17
JAKLIČG, ŽAGARE.Planar cubic G1 interpolatory splines with small strain energy [J]. Journal of Computational and Applied Mathematics, 2011, 235(8): 2758-2765.DOI: 10.1016/j.cam.2010.11.025
查找原文
参考文献 18
AHN Y J, HOFFMANNC, ROSENP.Geometric constraints on quadratic Bézier curves using minimal length and energy [J]. Journal of Computational and Application Mathematics, 2014, 255: 887-897.DOI:10.1016/j.cam.2013.07.005
查找原文
参考文献 19
LIJ C.Planar T-Bézier curve with approximate minimum curvature variation [J]. Journal of Advanced Mechanical Design, Systems, and Manufacturing, 2018, 12(1): 29.DOI:10.1299/jamdsm.2018jamdsm0029
查找原文
参考文献 20
FARING.Geometric Hermite interpolation with circular precision [J]. Computer-Aided Geometric Design, 2008, 40(4): 476-479.DOI:10.1016/j.cad.2008.01.003
查找原文
参考文献 21
LUL Z.A note on curvature variation minimizing cubic Hermite interpolants [J]. Applied Mathematics and Computation, 2015, 259: 596-599.DOI:10.1016/j.amc.2014.11.113
查找原文
目录 contents

    摘要

    当数据点给定时,三次Cardinal样条的张力参数和边界条件均为自由变量，因此可对这些自由变量进行优化，以得到满足某种特定要求的最佳三次Cardinal样条。讨论了如何通过优化张力参数与边界条件使得构造的平面三次Cardinal样条尽可能光顺。首先,分析了三次Cardinal参数样条曲线形状的影响因素；然后，利用曲率变化能极小对三次Cardinal参数样条曲线的张力参数与边界条件进行优化，获得张力参数与边界条件的唯一解； 最后，给出了对应三次Cardinal样条函数的张力参数与边界条件的优化方法。实例表明，经曲率变化极小优化张力参数与边界条件后的三次Cardinal样条比三次Catmull-Rom样条更为光顺，插值效果更好。

    Abstract

    The tension parameter and boundary condition of the cubic Cardinal spline can be optimized to obtain the best cubic Cardinal spline satisfying some specific requirements, because they are free variables when the data points are fixed. How to construct the smooth planar cubic Cardinal spline by optimizing the tension parameter and boundary condition is discussed in this paper. Firstly, the influence factors on the shape of the parametric cubic Cardinal spline curve are analyzed. Then, the tension parameter and boundary condition of the parametric cubic Cardinal parameter spline curve are optimized by the curvature variation minimization, and the unique solution of the tension parameter and the boundary condition is obtained. Finally, the corresponding method of the cubic Cardinal spline function is given. Some examples show that the cubic Cardinal spline obtained by optimizing the tension parameter and boundary condition is smoother, and the interpolation effect is better than that of the cubic Catmull-Rom spline.

    在计算机辅助设计及其相关研究领域中， 插值曲线的构造一直都是重要的研究课题。 作为一种分段插值样条， 三次Cardinal样条［1］是构造插值曲线的一种重要方法。 三次Cardinal样条不仅满足C1连续， 而且无须求解方程系统即可直接插值于给定的数据点。 另外， 当数据点保持不变时， 三次Cardinal样条的插值效果还可通过自带的张力参数进行调控。 这些优点使得三次Cardinal样条被广泛应用于多工程领域［2,3］。 虽然， 三角Cardinal样条［4,5］、双曲Cardinal样条［6］以及C2连续的Cardinal样条［7］在某些方面较传统三次Cardinal样条更优，但它们的方程结构更为复杂， 因此不可否认，传统三次Cardinal样条仍具有较高的应用价值。

    注意到， 张力参数对三次Cardinal样条的形状有明显的影响， 为了使三次Cardinal样条能插值于给定的首、末数据点， 往往需要在数据点的两端添加2个边界条件， 故边界条件对首、末两段样条曲线的形状也有较大影响。 因此， 在利用三次Cardinal样条构造插值曲线时， 可通过修改张力参数和边界条件来实现对插值效果的调控， 这显然为插值曲线的构造提供了便利。 当然， 也可以在某种准则下确定张力参数和边界条件的最佳取值，使得构造的三次Cardinal样条满足某种特定的要求。 为使三次Cardinal样条曲线尽可能光顺， 文献［8］提出了一种利用曲率变化极小来确定张力参数和边界条件最佳取值的方法。 但文献［8］仅简单给出了确定张力参数和边界条件最佳取值的方法， 对其具体过程以及相应的三次Cardinal样条函数并未展开详细讨论。 为此， 本文进一步讨论如何通过优化张力参数与边界条件以使得构造的平面三次Cardinal样条曲线及其函数形式尽可能光顺的方法。 通过对平面三次Cardinal样条曲线及其函数形式的近似曲率变化进行极小化， 易获得张力参数与边界条件的唯一解， 从而使构造的三次Cardinal样条曲线及其函数形式尽可能光顺， 且具有更好的插值效果。

  • 1　问题的提出

    给定平面上n(n≥3)个数据点pi(i=0,1,⋯,n)， 三次Cardinal参数样条曲线［1］可表示为

    Ri(t)=b0(t)pi-1+b1(t)pi+b2(t)pi+1+b3(t)pi+2,
    （1）

    式中，i=1,2,⋯,n-2，0≤t≤1，bj(t)(j=0,1,2,3)为

    b0(t)=1-T-t3+2t2-t/2,b1(t)=3+Tt3-5+Tt2+2/2,b2(t)=b1(1-t),b3(t)=b0(1-t),
    （2）

    其中， 实数T称为张力参数。

    特别地， 当三次Cardinal样条用于插值用途时， 张力参数的最佳值取为T=0［9］， 此时对应的三次Cardinal样条又称为三次Catmull-Rom样条［10］。 三次Catmull-Rom样条在许多工程领域已得到了较广泛的应用［11,12］。

    由式（1）与式（2）不难发现， 当T=1时， 式（1）为

    Ri(t)=2t3-3t2+1pi+-2t3+3t2pi+1，

    此时三次Cardinal参数样条曲线将退化为一条直线段。 因此， 下面仅考虑T≠1的情形。

    式（1）与式（2）经计算可得，

    Ri(0)=pi，Ri(1)=pi+1；
    （3）
    R'i(0)=1-T2pi+1-pi-1,
    R'i(1)=1-T2pi+2-pi。
    （4）

    式（3）表明，每段三次Cardinal参数样条曲线插值于第2与第3个数据点，意味着整条样条曲线插值于除首、末数据点以外的其他数据点。由式（4）可得R'i(1)=R'i+1(0)，表明三次Cardinal参数样条曲线满足C1连续。另外，由式（4）不难发现，当所有数据点保持不变时，三次Cardinal参数样条曲线的形状将由张力参数T决定。例如，给定数据点pi(i=0,1,⋯,5)的坐标(xi,yi)为

    xi=0,-2,-2,0,0,-1，
    yi=0,0,8,8,13,12，

    张力参数T取不同值时的三次Cardinal参数样条曲线如图1所示。

    图1 张力参数取不同值时的三次Cardinal参数样条曲线

    图1 张力参数取不同值时的三次Cardinal参数样条曲线

    Fig.1 Parametric cubic Cardinal spline curve with different tensor parameters

    另一方面， 在许多实际问题中往往要求三次Cardinal参数样条曲线也能插值于首、末数据点， 此时，须在数据点pi(i=1,2,⋯,n)的两端补充2个边界条件：p-1与pn+1， 常用的边界条件取为p-1=p0，pn+1=pn。 而事实上， 当张力参数T固定， 2个边界条件取不同点时， 首、末两段三次Cardinal参数样条曲线的形状也将不同。 例如， 对于图1中给定的数据点， 取张力参数T=0， 2个边界条件取不同点时的三次Cardinal参数样条曲线如图2所示， 图中虚线对应的边界条件为p-1=p0，p6=p5； 实线对应的边界条件为p-1=p2，p6=p3。

    图2 边界条件取不同点时的三次Cardinal参数样条曲线

    图2 边界条件取不同点时的三次Cardinal参数样条曲线

    Fig.2 Parametric cubic Cardinal spline curve with different boundary condition

    综上， 当数据点给定时， 三次Cardinal参数样条曲线的形状由张力参数T以及2个边界条件p-1与pn+1决定。理论上， 可根据不同的需求将张力参数和边界条件取为任意值。 当然， 也可对张力参数与边界条件进行优化， 得到满足特定要求的三次Cardinal参数样条曲线。

  • 2 张力参数与边界条件的优化

    光顺性是曲线的重要几何特征之一。 在CAD及其相关研究领域， 构造光顺的曲线是一项重要的研究课题［13,14］。虽然曲线的光顺性无法定量描述， 但人们通常通过极小化曲线的能量函数来实现平面光顺曲线的构造［15］。应变能量和曲率变化能量成为研究人员描述平面曲线光顺性的2种常见方法［16,17,18,19］。本文利用曲率变化能极小来优化平面三次Cardinal参数样条的张力参数和边界条件， 使得构造的样条曲线尽可能光顺。

    平面参数曲线r(t)的曲率变化能定义为［20］

    E=∫01k'(t)2dt，
    （5）

    式中，κ(t)=r'(t)×r″(t)r'(t)3。

    由于式（5）高度非线性， 往往需要利用近似形式将其线性化。 假设曲线r(t)被弧长近似参数化， 则其曲率变化能可近似表示为［21］

    Eˆ=∫01r‴(t)2dt。
    （6）

    为讨论方便， 令T=1-2α， 则式（2）可改写为

    Ri(t)=Li(t)α+Mi(t)，
    （7）

    式中，

    Li(t):=-t3+2t2-tpi-1-t3-t2pi+t3-2t2+tpi+1+t3-t2pi+2，
    Mi(t):=2t3-3t2+1pi-2t3-3t2pi+1。

    于是， 据式（6）与（7）， 补充边界条件p-1与pn+1后的三次Cardinal参数样条曲线的应变能可近似表示为

    Eˆ=∑i=0n-1∫01R‴i(t)2dt=∑i=0n-1∫01L‴i(t)α+M‴i(t)2dt=α2∑i=0n-1ai+2α∑i=0n-1bi+∑i=0n-1ci，
    （8）

    式中，

    ai:=∫01L‴i(t)2dt=36pi-1+pi-pi+1-pi+22，
    bi:=∫01L‴i(t)M‴i(t)dt=72pi+1-pi×pi-1+pi-pi+1-pi+2,
    ci:=∫01M‴i(t)2dt=144pi+1-pi2。

    显然， 当数据点pi(i=0,1,⋯,n)固定时， 式（8）中的参数α(α=1-T/2)与边界条件p-1，pn+1为自由变量。因此，为了确定最优的张力参数与边界条件，使得三次Cardinal参数样条曲线的应变能近似极小化， 可得无约束优化问题

    minEˆ(α,p-1,pn+1)。
    （9）

    为了求解问题（9）， 首先给出以下引理。

    引理1给定平面上的向量a=(ax,ay)，b=(bx,by)， 记∂∂a:=(∂∂ax,∂∂ay)， 则有

    ∂a2∂a=2a，∂(a⋅b)∂a=b。

    证明由于a2=ax2+ay2，a⋅b=axbx+ayby， 故

    ∂a2∂ax=2ax，∂a2∂ay=2ay，
    ∂(a⋅b)∂ax=bx，∂(a⋅b)∂ay=by。

    于是，

    ∂a2∂a=∂a2∂ax,∂a2∂ay=(2ax,2ay)=2a，
    ∂(a⋅b)∂a=∂(a⋅b)∂ax,∂(a⋅b)∂ay=(bx,by)=b。

    证毕。

    若记ui:=pi-1+pi-pi+1-pi+2，vi:=pi+1-pi， 则由引理1,有

    ∂a0∂p-1=72u0，∂aj∂p-1=0，j=1,2,⋯,n-1，
    ∂an-1∂pn+1=-72un-1，∂aj∂pn+1=0，j=0,1,⋯,n-2，
    ∂b0∂p-1=72v0，∂bj∂p-1=0，j=1,2,⋯,n-1，
    ∂bn-1∂pn+1=-72vn-1，∂bj∂pn+1=0，j=0,1,⋯,n-2，
    ∂cj∂p-1=∂cj∂pn+1=0，j=0,1,⋯,n-1。

    于是，

    ∂Eˆ∂α=2α∑i=0n-1ai+2∑i=0n-1bi=72αu02+α∑i=1n-2ui2+αun-12+144u0⋅v0+∑i=1n-2ui⋅vi+un-1⋅vn-1，
    （10）
    ∂Eˆ∂p-1=∂a0∂p-1α2+2∂b0∂p-1α=72u0α2+144v0α，
    （11）
    ∂Eˆ∂pn+1=∂an-1∂pn+1α2+2∂bn-1∂pn+1α=-72un-1α2-144vn-1α。
    （12）

    于是有

    定理1给定平面上一列数据点pi(i=0,1,⋯,n)， 要使所构造的三次Cardinal参数样条曲线具有极小近似的曲率变化能， 则张力参数T(T≠1)与2个边界条件p-1与pn+1应分别取为

    T=1+4∑i=1n-2pi-1+pi-pi+1-pi+2pi+1-pi∑i=1n-2pi-1+pi-pi+1-pi+22，
    p-1=p2-3+T1-Tp1-p0，
    pn+1=pn-2+3+T1-Tpn-pn-1。

    证明极小化Eˆ(α,p-1,pn+1)， 可得方程组

    ∂Eˆ∂α=0,∂Eˆ∂p-1=0,∂Eˆ∂pn+1=0。
    （13）

    由式（11）与式（12）， 当α≠0时， 令∂Eˆ∂p-1=0与∂Eˆ∂pn+1=0，可得

    u0=-2αv0，un-1=-2αvn-1。
    （14）

    将式（14）代入式（10）， 并令∂Eˆ∂α=0， 经计算整理得

    α=-2∑i=1n-2ui⋅vi∑i=1n-2ui2。
    （15）

    将ui=pi-1+pi-pi+1-pi+2，vi=pi+1-pi，α=1-T/2分别代入式（14）与式（15），并经简单推导即得定理成立。证毕。

    下面通过数值实例来说明本文方法的有效性。由于三次Catmull-Rom样条被视为三次Cardinal样条中插值效果最好的［9］， 因此， 本文将与三次Catmull-Rom参数样条曲线进行比较， 且三次Catmull-Rom参数样条曲线的边界条件取常用情形， 即p-1=p0，pn+1=pn。

    例1对于图1中给定的数据点， 要使所构造的三次Cardinal参数样条曲线具有极小的近似曲率变化能， 经计算所得张力参数、2个边界条件及近似曲率变化能与三次Catmull-Rom参数样条曲线的对比如表1所示。

    表1 本文方法与Catmull-Rom的对比(例1)

    Table 1 The comparison between the proposed method and Catmull-Rom (example 1)

    样条曲线Catmull-Rom本文方法
    张力参数T=0T=-0.3844
    边界条件

    p-1=(0,0)

    p6=(-1,12)

    p-1=(1.7790,8.0000)

    p6=(-1.8895,6.1105)

    近似曲率

    变化能

    		6.2640×1034.3718×103

    构造的具有极小近似曲率变化能的三次Cardinal参数样条曲线（实线）与三次Catmull-Rom参数样条曲线（点虚线）如图3所示。

    图3 具有极小近似曲率变化能的三次Cardinal参数样条曲线(例1)

    图3 具有极小近似曲率变化能的三次Cardinal参数样条曲线(例1)

    Fig. 3 Parametric cubic Cardinal spline curve with minimum approximate curvature variation energy (example 1)

    例2从单位半圆上取定数据点pi(i=0,1,⋯,5)的坐标(xi,yi)为

    xi=1,1/2,0,-1/2,-1，
    yi=0,3/2,1,-3/2,0。

    要使所构造的三次Cardinal参数样条曲线具有极小的近似曲率变化能，经计算可得张力参数、2个边界条件及近似曲率变化能与三次Catmull-Rom参数样条曲线的对比如表2所示。

    表2 本文方法与Catmull-Rom的对比(例2)

    Table 2 The comparison between the proposed method and Catmull-Rom (example 2)

    样条曲线Catmull-Rom本文方法
    张力参数T=0T=0.0928
    边界条件

    p-1=(1,0)

    p6=(-1,0)

    p-1=(1.7046,-1.9525)

    p6=(-1.7046,-1.9525)

    近似曲率变化能54.346 33.101 6

    构造的具有极小近似曲率变化能的三次Cardinal参数样条曲线（实线）与三次Catmull-Rom参数样条曲线（长虚线）如图4所示。

    图4 具有极小近似曲率变化能的三次Cardinal参数样条曲线(例2)

    图4 具有极小近似曲率变化能的三次Cardinal参数样条曲线(例2)

    Fig. 4 Parametric cubic Cardinal spline curve with minimum approximate curvature variation energy (example 2)

    由例1与例2可知， 张力参数与边界条件优化后的三次Cardinal参数样条曲线较三次Catmull-Rom参数样条曲线更为光顺， 插值效果更好。

  • 3　三次Cardinal样条函数

    与三次Cardinal参数样条曲线对应， 本节讨论三次Cardinal样条函数。

    给定平面上n(n≥3)个数据点(xi,yi)(i=0,1,⋯,n)， 设xi为等距， 即xi+1-xi=h（h为常数）， 三次Cardinal样条函数可表示为

    Si(t)=b0(t)yi-1+b1(t)yi+b2(t)yi+1+b3(t)yi+2，
    （16）

    式中，i=1,2,⋯,n-2，xi≤x≤xi+1，t=(x-xi)/h，bj(t)(j=0,1,2,3)为式（2）表示的函数。

    由式（2）不难验证，

    Si(xi)=yi，Si(xi+1)=yi+1，
    （17）
    S'i(xi)=1-T2hyi+1-yi-1，
    S'i(xi+1)=1-T2hyi+2-yi。
    （18）

    由式（17）与式（18）可知，三次Cardinal样条函数插值于除首、末数据点外的其他数据点， 且函数满足C1连续。同样，可以补充2个边界条件(x-1,y-1)与(xn+1,yn+1)，使得三次Cardinal样条函数插值于包括首、末数据点在内的所有数据点。由于要求xi为等距，故有x-1=2x0-x1，xn+1=2xn-xn-1。显然，当所有数据点保持不变时，三次Cardinal样条函数的形状由张力参数T以及2个边界条件y-1与yn+1决定。

    与三次Cardinal参数样条曲线的情形类似， 三次Cardinal样条函数满足

    定理2给定平面上一列数据点(xi,yi)(i=0,1,⋯,n)， 要使所构造的三次Cardinal样条函数具有极小的曲率变化能， 则张力参数T(T≠1)与2个边界条件y-1与yn+1应分别取

    T=1+4∑i=1n-2yi-1+yi-yi+1-yi+2yi+1-yi∑i=1n-2yi-1+yi-yi+1-yi+22，
    y-1=y2-3+T1-Ty1-y0，
    yn+1=yn-2+3+T1-Tyn-yn-1。

    下面通过数值实例来说明本文方法的有效性， 同上节，将本文方法与三次Catmull-Rom样条函数进行比较， 且三次Catmull-Rom样条函数的边界条件为常用情形， 即y-1=y0，yn+1=yn。

    例3设y=sinx(0≤x≤2π)， 取xi=π4i，yi=sinxi(i=0,1,⋯,8)。 要使所构造的三次Cardinal参数样条曲线具有极小的近似曲率变化能， 经计算得张力参数与2个边界条件分别为T=-0.1716，y-1=-0.7071，y9=0.7071。 构造的具有极小近似曲率变化能的三次Cardinal样条插值函数（实线）与三次Catmull-Rom样条插值函数（点虚线）及对应的绝对误差曲线如图5所示。

    图5 样条插值函数及其绝对误差曲线(例3)

    图5 样条插值函数及其绝对误差曲线(例3)

    Fig.5 Spline interpolation functions and the absolute error curves (example 3)

    例4设y=11+x2(-1≤x≤1)， 取xi=-1+27i，yi=11+xi2(i=0,1,⋯,7)。 要使所构造的三次Cardinal参数样条曲线具有极小的近似曲率变化能，经计算所得张力参数与2个边界条件分别为T=-0.1316，y-1=0.4338，y9=0.4338。构造的具有极小近似曲率变化能的三次Cardinal样条插值函数（实线）与三次Catmull-Rom样条插值函数（点虚线）及对应的绝对误差曲线如图6所示。

    图6 样条插值函数及其绝对误差曲线(例4)

    图6 样条插值函数及其绝对误差曲线(例4)

    Fig.6 Spline interpolation functions and the absolute error curves (example 4)

    由例3与例4可知， 优化张力参数与边界条件后的三次Cardinal样条函数较三次Catmull-Rom样条函数的插值效果更好。

  • 4　结 语

    当数据点给定时，三次Cardinal样条的形状完全由张力参数和边界条件决定。虽然通过修改张力参数与边界条件可实现对三次Cardinal样条形状的调控，但有时只要求得张力参数和边界条件的最佳值，就能使构造的三次Cardinal样条满足某种特定的要求。为此，本文给出了一种利用曲率变化能极小来优化三次Cardinal样条的张力参数与边界条件的方法，该方法容易获得张力参数与边界条件的唯一解。相对于三次Catmull-Rom样条，优化张力参数与边界条件后的三次Cardinal样条不仅更为光顺，而且插值效果更好。

  • 参考文献(References)

    • 1

      HEARN D, BAKER M P.Computer Graphics: C Version [M]. 2nd ed. Upper Saddle River: Prentice Hall, 1996.

    • 2

      YOON H S, PARK T H.Smooth path planning method for autonomous mobile robots using cardinal spline [J]. Transactions of the Korean Institute of Electrical Engineers, 2010, 59 (4): 803-808.

    • 3

      FANANI A, YUNIARTI A, SUCIATI N.Geometric feature extraction of batik image using Cardinal spline curve representation [J]. Telkomnika, 2014, 12(2): 397-404.DOI:10.12928/telkomnika.v12i2.1915

    • 4

      吴晓勤, 严秀坤.三次C-Cardinal样条曲线及曲面[J]. 计算机工程与科学, 2006, 28(2): 48-50.DOI:10.3969/j.issn.1007-130X.2006.02.015

      WU X Q, YAN X K.The cubic C-Cardinal spline curve and surface [J]. Computer Engineering & Science, 2006, 28(2): 48-50.DOI:10.3969/j.issn.1007-130X.2006.02.015

    • 5

      宋爱平, 陶建明, 易旦萍, 等.可调形三次三角Cardinal插值样条曲线[J]. 计算数学, 2015, 37(1): 34-41.

      SONG A P, TAO J M, YI D P, et al.Cubic trigonometric Cardinal interpolation spline with adjustable shape [J]. Mathematica Numerica Sinica, 2015, 37(1): 34-41.

    • 6

      吴晓勤, 韩旭里, 罗善明.三次H-Cardinal样条曲线及曲面[J]. 工程图学学报, 2008, 29(5): 83-88.DOI:10.3969/j.issn.1003-0158.2008.05.016

      WU X Q, HAN X L, LUO S M.Cubic H-Cardinal spline curves and surfaces [J]. Journal of Engineering Graphics, 2008, 29(5): 83-88.DOI:10.3969/j.issn.1003-0158.2008.05.016

    • 7

      李军成, 刘成志, 易叶青.带形状因子的C2连续五次 Cardinal 样条与Catmull-Rom 样条[J]. 计算机辅助设计与图形学学报, 2016, 28(11): 1821-1831.

      LI J C, LIU C Z, YI Y Q.C2 continuous quintic Cardinal spline and Catmull-Rom spline with shape factors [J]. Journal of Computer-Aided Design & Computer Graphics, 2016, 28(11): 1821-1831.

    • 8

      LI J C.Curvature variation minimizing Cardinal spline curves [J]. Journal of Modeling and Optimization, 2018, 10(1): 31-36.DOI:10.32732/jmo.2018.10.1.31

    • 9

      GONZALEZ R C , WOODS R E, EDDINS S L.Digital Image Processing Using MATLAB [M]. Upper Saddle River: Prentice Hall, 2004.DOI:http://dx.doi.org/10.9774/GLEAF.978-1-909493-38-4_2

    • 10

      DEROSE T D, BARSKY B A.Geometric continuity, shape parameters, and geometric constructions for Catmull-Rom splines [J]. ACM Transactions on Graphics, 1988(7): 1-41.DOI:10.1145/42188.42265

    • 11

      YUKSEL C, SCHAEFER S, KEYSER J.Parameterization and applications of Catmull–Rom curves [J]. Computer-Aided Design, 2011, 43(7): 747-755.DOI:10.1016/j.cad.2010.08.008

    • 12

      SCARPINITI M, COMMINIELLO D, PARISI R, et al.Nonlinear spline adaptive filtering [J]. Signal Processing, 2013, 93(4): 772-783.DOI:10.1016/j.sigpro.2012.09.021

    • 13

      HOFER M, POTTMANN H.Energy-minimizing splines in manifolds [J]. ACM Transaction on Graphics, 2004, 23(3): 284-293.DOI:10.1145/1015706.1015716

    • 14

      FAN W, LEE C H, CHEN J H. A realtime curvature-smooth interpolation scheme and motion planning for CNC machining of short line segments [J]. International Journal of Machine Tools and Manufacture, 2015, 96: 27-46.DOI:10.1016/j.ijmachtools.2015.04.009

    • 15

      LU L Z, JIANG C K, HU Q Q.Planar cubic G1 and quintic G2 Hermite interpolations via curvature variation minimization [J]. Computers & Graphics, 2017, 70: 92-98.DOI:10.1016/j.cag.2017.07.007

    • 16

      JAKLIČ G, ŽAGAR E.Curvature variation minimizing cubic Hermite interpolants [J]. Applied Mathematics and Computation, 2011, 218(7): 3918-3924. DOI:10.1016/j.amc.2011.09.039

    • 17

      JAKLIČ G, ŽAGAR E.Planar cubic G1 interpolatory splines with small strain energy [J]. Journal of Computational and Applied Mathematics, 2011, 235(8): 2758-2765.DOI: 10.1016/j.cam.2010.11.025

    • 18

      AHN Y J, HOFFMANN C, ROSEN P.Geometric constraints on quadratic Bézier curves using minimal length and energy [J]. Journal of Computational and Application Mathematics, 2014, 255: 887-897.DOI:10.1016/j.cam.2013.07.005

    • 19

      LI J C.Planar T-Bézier curve with approximate minimum curvature variation [J]. Journal of Advanced Mechanical Design, Systems, and Manufacturing, 2018, 12(1): 29.DOI:10.1299/jamdsm.2018jamdsm0029

    • 20

      FARIN G.Geometric Hermite interpolation with circular precision [J]. Computer-Aided Geometric Design, 2008, 40(4): 476-479.DOI:10.1016/j.cad.2008.01.003

    • 21

      LU L Z.A note on curvature variation minimizing cubic Hermite interpolants [J]. Applied Mathematics and Computation, 2015, 259: 596-599.DOI:10.1016/j.amc.2014.11.113

  • 基本信息

    DOI:10.3785/j.issn.1008-9497.2019.02.004

    中图分类号: TP391

    文献标识码: A

    引用信息

    李军成,刘成志,易叶青.优化张力参数与边界条件的平面三次Cardinal样条[J].浙江大学学报(理学版),2019,46(02):164-171.

    JunchengLI,Chengzhi LIU,Yeqing YI.Planar cubic Cardinal spline with tension parameter and boundary condition optimization[J].Journal of Zhejiang University(Science Edition),2019,46(02):164-171.

    基金信息

    国家自然科学基金资助项目(61472135); 湖南省自然科学基金资助项目(2017JJ3124)； 湖南省教育厅资助科研项目（18A415）.

    稿件历史

    纸刊出版 : 2019-03-25
    收稿日期 : 2018-08-14

    • 李军成

    机 构：湖南人文科技学院 数学与金融学院，湖南 娄底 417000

    Affiliation：College of Mathematics and Finance, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China

    邮 箱：lijuncheng82@126.com.

    作者简介：李军成(1982—)，ORCID: http://orcid.org/ 0000-0002-1904-4068，男，博士，副教授，主要从事计算机辅助几何设计及其应用研究, E-mail：lijuncheng82@126.com.

    刘成志

    机 构：湖南人文科技学院 数学与金融学院，湖南 娄底 417000

    Affiliation：College of Mathematics and Finance, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China

    易叶青

    机 构：湖南人文科技学院 信息学院，湖南 娄底 417000

    Affiliation：College of Information, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China

    1008⁃9497-2019-46-2-164/alternativeImage/6a605db4-601d-49d5-82b1-1d3c50f9a458-F001.jpg
    1008⁃9497-2019-46-2-164/alternativeImage/6a605db4-601d-49d5-82b1-1d3c50f9a458-F002.jpg
    样条曲线Catmull-Rom本文方法
    张力参数T=0T=-0.3844
    边界条件

    p-1=(0,0)

    p6=(-1,12)

    p-1=(1.7790,8.0000)

    p6=(-1.8895,6.1105)

    近似曲率

    变化能

    		6.2640×1034.3718×103
    1008⁃9497-2019-46-2-164/alternativeImage/6a605db4-601d-49d5-82b1-1d3c50f9a458-F003.jpg
    样条曲线Catmull-Rom本文方法
    张力参数T=0T=0.0928
    边界条件

    p-1=(1,0)

    p6=(-1,0)

    p-1=(1.7046,-1.9525)

    p6=(-1.7046,-1.9525)

    近似曲率变化能54.346 33.101 6
    1008⁃9497-2019-46-2-164/alternativeImage/6a605db4-601d-49d5-82b1-1d3c50f9a458-F004.jpg
    1008⁃9497-2019-46-2-164/alternativeImage/6a605db4-601d-49d5-82b1-1d3c50f9a458-F005.jpg
    1008⁃9497-2019-46-2-164/alternativeImage/6a605db4-601d-49d5-82b1-1d3c50f9a458-F006.jpg

    图1 张力参数取不同值时的三次Cardinal参数样条曲线

    Fig.1 Parametric cubic Cardinal spline curve with different tensor parameters

    图2 边界条件取不同点时的三次Cardinal参数样条曲线

    Fig.2 Parametric cubic Cardinal spline curve with different boundary condition

    表1 本文方法与Catmull-Rom的对比(例1)

    Table 1 The comparison between the proposed method and Catmull-Rom (example 1)

    图3 具有极小近似曲率变化能的三次Cardinal参数样条曲线(例1)

    Fig. 3 Parametric cubic Cardinal spline curve with minimum approximate curvature variation energy (example 1)

    表2 本文方法与Catmull-Rom的对比(例2)

    Table 2 The comparison between the proposed method and Catmull-Rom (example 2)

    图4 具有极小近似曲率变化能的三次Cardinal参数样条曲线(例2)

    Fig. 4 Parametric cubic Cardinal spline curve with minimum approximate curvature variation energy (example 2)

    图5 样条插值函数及其绝对误差曲线(例3)

    Fig.5 Spline interpolation functions and the absolute error curves (example 3)

    图6 样条插值函数及其绝对误差曲线(例4)

    Fig.6 Spline interpolation functions and the absolute error curves (example 4)

    image /

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

    无注解

  • 参考文献(References)

    • 1

      HEARN D, BAKER M P.Computer Graphics: C Version [M]. 2nd ed. Upper Saddle River: Prentice Hall, 1996.

    • 2

      YOON H S, PARK T H.Smooth path planning method for autonomous mobile robots using cardinal spline [J]. Transactions of the Korean Institute of Electrical Engineers, 2010, 59 (4): 803-808.

    • 3

      FANANI A, YUNIARTI A, SUCIATI N.Geometric feature extraction of batik image using Cardinal spline curve representation [J]. Telkomnika, 2014, 12(2): 397-404.DOI:10.12928/telkomnika.v12i2.1915

    • 4

      吴晓勤, 严秀坤.三次C-Cardinal样条曲线及曲面[J]. 计算机工程与科学, 2006, 28(2): 48-50.DOI:10.3969/j.issn.1007-130X.2006.02.015

      WU X Q, YAN X K.The cubic C-Cardinal spline curve and surface [J]. Computer Engineering & Science, 2006, 28(2): 48-50.DOI:10.3969/j.issn.1007-130X.2006.02.015

    • 5

      宋爱平, 陶建明, 易旦萍, 等.可调形三次三角Cardinal插值样条曲线[J]. 计算数学, 2015, 37(1): 34-41.

      SONG A P, TAO J M, YI D P, et al.Cubic trigonometric Cardinal interpolation spline with adjustable shape [J]. Mathematica Numerica Sinica, 2015, 37(1): 34-41.

    • 6

      吴晓勤, 韩旭里, 罗善明.三次H-Cardinal样条曲线及曲面[J]. 工程图学学报, 2008, 29(5): 83-88.DOI:10.3969/j.issn.1003-0158.2008.05.016

      WU X Q, HAN X L, LUO S M.Cubic H-Cardinal spline curves and surfaces [J]. Journal of Engineering Graphics, 2008, 29(5): 83-88.DOI:10.3969/j.issn.1003-0158.2008.05.016

    • 7

      李军成, 刘成志, 易叶青.带形状因子的C2连续五次 Cardinal 样条与Catmull-Rom 样条[J]. 计算机辅助设计与图形学学报, 2016, 28(11): 1821-1831.

      LI J C, LIU C Z, YI Y Q.C2 continuous quintic Cardinal spline and Catmull-Rom spline with shape factors [J]. Journal of Computer-Aided Design & Computer Graphics, 2016, 28(11): 1821-1831.

    • 8

      LI J C.Curvature variation minimizing Cardinal spline curves [J]. Journal of Modeling and Optimization, 2018, 10(1): 31-36.DOI:10.32732/jmo.2018.10.1.31

    • 9

      GONZALEZ R C , WOODS R E, EDDINS S L.Digital Image Processing Using MATLAB [M]. Upper Saddle River: Prentice Hall, 2004.DOI:http://dx.doi.org/10.9774/GLEAF.978-1-909493-38-4_2

    • 10

      DEROSE T D, BARSKY B A.Geometric continuity, shape parameters, and geometric constructions for Catmull-Rom splines [J]. ACM Transactions on Graphics, 1988(7): 1-41.DOI:10.1145/42188.42265

    • 11

      YUKSEL C, SCHAEFER S, KEYSER J.Parameterization and applications of Catmull–Rom curves [J]. Computer-Aided Design, 2011, 43(7): 747-755.DOI:10.1016/j.cad.2010.08.008

    • 12

      SCARPINITI M, COMMINIELLO D, PARISI R, et al.Nonlinear spline adaptive filtering [J]. Signal Processing, 2013, 93(4): 772-783.DOI:10.1016/j.sigpro.2012.09.021

    • 13

      HOFER M, POTTMANN H.Energy-minimizing splines in manifolds [J]. ACM Transaction on Graphics, 2004, 23(3): 284-293.DOI:10.1145/1015706.1015716

    • 14

      FAN W, LEE C H, CHEN J H. A realtime curvature-smooth interpolation scheme and motion planning for CNC machining of short line segments [J]. International Journal of Machine Tools and Manufacture, 2015, 96: 27-46.DOI:10.1016/j.ijmachtools.2015.04.009

    • 15

      LU L Z, JIANG C K, HU Q Q.Planar cubic G1 and quintic G2 Hermite interpolations via curvature variation minimization [J]. Computers & Graphics, 2017, 70: 92-98.DOI:10.1016/j.cag.2017.07.007

    • 16

      JAKLIČ G, ŽAGAR E.Curvature variation minimizing cubic Hermite interpolants [J]. Applied Mathematics and Computation, 2011, 218(7): 3918-3924. DOI:10.1016/j.amc.2011.09.039

    • 17

      JAKLIČ G, ŽAGAR E.Planar cubic G1 interpolatory splines with small strain energy [J]. Journal of Computational and Applied Mathematics, 2011, 235(8): 2758-2765.DOI: 10.1016/j.cam.2010.11.025

    • 18

      AHN Y J, HOFFMANN C, ROSEN P.Geometric constraints on quadratic Bézier curves using minimal length and energy [J]. Journal of Computational and Application Mathematics, 2014, 255: 887-897.DOI:10.1016/j.cam.2013.07.005

    • 19

      LI J C.Planar T-Bézier curve with approximate minimum curvature variation [J]. Journal of Advanced Mechanical Design, Systems, and Manufacturing, 2018, 12(1): 29.DOI:10.1299/jamdsm.2018jamdsm0029

    • 20

      FARIN G.Geometric Hermite interpolation with circular precision [J]. Computer-Aided Geometric Design, 2008, 40(4): 476-479.DOI:10.1016/j.cad.2008.01.003

    • 21

      LU L Z.A note on curvature variation minimizing cubic Hermite interpolants [J]. Applied Mathematics and Computation, 2015, 259: 596-599.DOI:10.1016/j.amc.2014.11.113