设为首页
加入收藏
期刊导航
网站地图
首页
期刊
数学与物理
地球与环境
信息通讯
经济与管理
生命科学
工程技术
医药卫生
人文社科
化学与材料
会议
合作
新闻
我们
招聘
千人智库
我要投搞
办刊
期刊菜单
●领域
●编委
●投稿须知
●最新文章
●检索
●投稿
文章导航
●Abstract
●Full-Text PDF
●Full-Text HTML
●Full-Text ePUB
●Linked References
●How to Cite this Article
Pure Mathematics
理论数学
, 2014
, 4, 14-20
http://dx.doi.org/10.12677/pm.2014.41003
Published Online
January 2014 (http://www.hanspub.org/journal/
pm.html
)
OPEN ACCESS
14
Application of a Four Quadratic Rational Interpolation
Spline Curve
Lin Fu
Anhui University of Science & Technology, Huainan
Email:
fulin07211208@163.com
Received: Dec. 16
th
, 2013; revised: Dec
.
27
th
, 2013; accepted: Dec
.
29
th
, 2013
Copyright © 201
4
Lin Fu. This is an open access article d istributed und er the Creative Common s Attribution Li cense, which permits un-
restricted use, distribu tion, and reproduction in any medium, provided the original work is properly cited. In accordance of the Creativ e
Commons Attribution Licen se all Copyrights © 201
4
are reserved for Hans and the owner of the intellectual pr operty Lin Fu
.
All Copy-
rig
ht © 2014 are guarded by law and by Hans as a guardian.
Abstract:
A four
q
uadratic
rational spline interpolation for non closed curves w
as
well described. This paper
presents an approach of four
quadratic rational interpolation spline
and
explores the interpolation function
monotonicity
and
continuity.
Error
estimates confirm the conformality
and
numerical practical examples illu
-
strate the effectiveness of the method.
Keywords:
C u rve ; Shape P reservi ng; Effective
一元四次有理插值样条曲线的应用
符
琳
安徽理工大学,淮南
Email: fulin07211208@163.com
收稿日期:
2013
年
12
月
16
日;修回日期:
2013
年
12
月
27
日;录用日期:
2013
年
12
月
29
日
摘
要:
一元四次有理插值样条对于非封闭曲线进行了很好地数学描述。本文提出了一元四次有理插
值样条的方法,探究了这种插值函数的单调性,连续性,误差估计证实其保形性,最后用实际的数值
实例来说明该方法的有效性。
关键词:
曲线
;保形;有效
1.
引言
随着航空、造船、机械设计和制造等现代工业的蓬勃发展,计算机辅助几何设计,简称
CAGD
(Computer
Aided Geometric Design)
,逐步成为的一门新兴的交叉学科与边缘学科。作为
CAGD
系统基本几何元素,自由曲
线、曲面的表示、设计、显示、分析以及规格、处理
(
包括数据结构、数据库、图形的信息形式和调整方式等
)
问题,是
CAGD
的主要研究对象和内容,而用插值与逼近方法解决曲线、曲面造型问题是
CAGD
最基础的研究
课题。所谓插值问题就是从给定的离散点的值,去构造一个连续定义的
(
简单
)
函数,使得它与被逼近的函数在给
定点的值完全一致。插值法是数值逼近的一种最简单的重要方法,利用插值法可以通过函数在有限个点处的取
值情况估算出函数在其它点处的值。插值法是整个数值逼近的基础。它被广泛应用于方程求根、函数逼近、数
值微分、数值积分、积分和微分方程数值解等。有理样条由于其极好的保形性,光滑性已经成为一种广泛应用
的曲线曲面设计方法。
符琳
|
一元四次有理插值样条曲线的应用
OPEN ACCESS
15
近些年来,不少学者在有理二次,三次有理样条以及他们的性质和应用上做了研究
[1
-8]
。如:
Gregory
,
Delbourgo
在文献
[1]
构造了
2/2
型有理插值样条,探究了它的保单调问题。王强在文献
[7]
构造了双参数
3/1
型的
有理插值样条,
探究了其光滑性。
Delbourgo
在文献
[2]
构造了
2/1
型有理插值样条,并探究了它的保凸问题。王
仁宏,吴顺唐在文献
[3]
从实际的课题着手,构造了几种具有线性结构的有理插值样条并讨论了它们的解析。王
强,
段奇、
Hussain
等在文献
[4
-7]
也是用含参数的分段三次有理函数构造了满足约束条件插值样条。
本文构造四次分段有理插值曲线即分子、分母分别为四次、一次多项式的插值函数,并讨论了该函数的一
阶、二阶连续性和保单调的条件。
该函数含有二族参数
α
,
β
可自由调节以改变曲线的形状,并可通过导数
d
赋
值,无需求解线性方程组使得插值函数二阶连续。从理论上证实插值样条函数的保形性。
它具有以下优点:
1)
局部性,每个区间的函数表达式只与相邻的若干节点取值有关;
2)
保形性,无需增加
中间节点,通过参数选取可实现保形要求;
3)
具有显式表达式,便于计算分析;
4)
含有参数,便于交互式修改,
可以通过参数的改变而不是点的变动来修改插值函数。
本文结构如下:第一部分的引言介绍了该插值方法存在的背景,目前的研究领域以及本方法的优点之处。
第二部分具体讲述了这种方法。第三部分从单调性,连续性和误差估计三方面探究了其的保形性。第四部分用
数值例子证实了这种方法的可行性。第五部分是对文章的小结。
2.
一元四次有理插值样条曲线的构造
给定数据
( )
,,1, 2,,,
ii
xf in
=
其中
i
f
为被插函数在分划点
i
x
上的函数值,此处
12
n
xx x
<<<
。记
1
ii i
hx x
+
= −
,
( )
1
ii i
ff
+
∆= −
,当
[ ]
1
,
ii
x xx
+
∈
时,
[ ]
0,1 ,
t
∈
1, 2,,1
in
= −
,定义:
( )
( )
( )
( )( )
if 0
if 0
ii i
ii
ii
Qt Pt
wxwx htWt
f
∆≠
= +≡=
∆=
(1)
()( )( )( )( )
43 2
2 34
12 31
11 11
i iiii
Q tftZttZttZttft
αβ
+
=−+ −+−+−+
( )()
1
ii i
Pttt
αβ
= −+
且
( )
( )
1
21
3 11
3
33
3
ii iiii
ii ii
iiiii i
Zf dh
Zf f
Zf dh
αβ α
αβ
αβ β
+
++
=++
= +
=+−
i
α
,
i
β
称为形状参数都取正数,
i
d
表示函数在点
i
x
处的导数值,对
(1)
定义的函数求导可得:
( )
( )
( )
( )
( )
( )
4
4
0
11
1
2
1
1
j
j
ij
j
i
ii
A tt
w x Wth
tt
αβ
−
=
−
−
==
−+
∑
(2)
其中
( )
( )
( )
2
0
2
1
21
2
31
2
41
23
34
23
i ii
ii ii
iii iii
ii ii
i ii
Ad
Ad
A dd
Ad
Ad
α
α
αβ
β
β
+
+
+
=
= ∆−
=∆− −
= ∆−
=
由
(1)
、
(2)
容易验证函数满足下列插值性质:
符琳
|
一元四次有理插值样条曲线的应用
OPEN ACCESS
16
( )
ii
wx f
=
,
( )
11
ii
wx f
++
=
;
( )
( )
1
ii
wx d
=
,
( )
( )
1
11
ii
wx d
++
=
因此,函数
( )
[ ]
1
1
,
n
wxCx x
∈
是在
i
x
和
1
i
x
+
处的
Hermite
插值。
3.
一元四次有理插值曲线的保形性
下面从插值函数单调性、插值函数二阶连续性和插值函数误差估计来说明其保形性。
3.1.
插值函数单调性
下面我们讨论
(1)
式定义的一元四次有理样条插值函数的单调性。假设
( )
fx
在区间
[ ]
,
ab
单调递增,因此设
12
n
ff f
≤≤≤
,或
0
i
∆≥
。
选取导数值
i
d
满足:
0,1,2, ,
i
di n
≥=
在区间
[ ]
1
,
ii
xx
+
,
( )
wx
单调递增的充要条件为:
( )
( )
1
0
wx
≥
注意到由
(1)
式给出的插值函数,
( )
wx
的导数
( )
( )
1
wx
的分母恒正,因此只需
( )
( )
1
wx
的分子大于零。当
0
i
d
≥
时,如果
1
0
i
A
≥
,
2
0
i
A
≥
,
3
0
i
A
≥
成立,则
( )
( )
1
0
wx
≥
,因此,我们得到插值函数
( )
wx
单调的充分条件为:
0, 1,2,3
ij
Aj
≥=
即
3
ii
d
≤∆
(3)
1
3
ii
d
+
≤∆
(4)
1
4
ii i
dd
+
+ ≤∆
(5)
因此,我们就得到下述定理:
定理
1
:已知数据
( )
12
n
ff f
≤≤≤
单调递增,导数值
0
i
d
≥
,当
i
d
满足
(3)
、
(4)
、
(5)
时,则存在二族含有参
数
i
α
,
i
β
的四次有理插值函数
( )
[ ]
1
,
w xCab
∈
并且是单调递增的。
3.2.
插值函数二阶连续性
对任意
[ ]
1
,
ii
x xx
+
∈
,由
(1)
式定义的一元四次有理样条插值函数
( )
wx
求二阶导数可得:
( )
( )
( )
( )
4
4
0
2
3
1
,
1
j
j
ij
j
ii i
B tt
wx
h tt
αβ
−
=
−
=
−+
∑
其中
3 32
0
642
i iiiiiii
B dd
αα αβ
= ∆−−
232 23
11
186 862
iiiiiiiii iiiii
Bd dd
αβααβ αβα
+
=∆− ∆−−+
22 22
21
18186 6
iiiiiiii iiiii
B dd
αβαβα βαβ
+
=∆− ∆+−
3 2223
31 1
6 86182
iiiii iii iiiiii
B ddd
β αβαβαββ
++
= ∆++−∆−
符琳
|
一元四次有理插值样条曲线的应用
OPEN ACCESS
17
3 23
41 1
426
ii iii iii
Bd d
β αββ
++
= +−∆
因此有
( )
( )
3 23
2
111 11
3
11
1 1 11
11
42 6
426
iiiiii i
i
ii
iiiii i
ii
dd
wx
h
dd
h
β αββ
β
βαβ
β
−−− −−
−
−−
−− −−
−−
+ −∆
=
+ −∆
=
及
( )
( )
332
2
3
642
642
.
iiiiiii
i
ii
i iiiii
ii
dd
wx
h
dd
h
αα αβ
α
ααβ
α
+
∆− −
=
∆− −
=
由函数
( )
wx
在点
i
x
二阶连续,可得
( )
11
1
11
1
21
33
iii i
i
ii
iiii
ii
hh
d
hhhh
βα
αβ
−−
−
−−
−
∆+ ∆
=
++ +
(6)
由以上推导可得下述定理
定理
2
:
当
i
d
满足
(6)
式时,含有参数
i
α
,
i
β
的四次有理插值函数族
( )
,,
ii
wx
αβ
在区间
[ ]
,
ab
上二阶连续。
3.3.
插值函数误差估计
对于分片函数
(1)
只需考虑在区间
[ ]
1
,
ii
xx
+
上的情形。
定理
3
:
假设
( )
[ ]
4
,
f xCab
∈
,
( )
wx
是由
(1)
式定义的在区间
[ ]
1
,
n
xx
上的分片有理插值样条函数,则对任意
[ ]
1
,
ii
xx
+
下式成立:
( )( )
( )
4
4
3844 16
iii iii
i
hkh ch
wxf xf
λ
µ
− ≤++
其中
{ }
max ,
i ii
λ αβ
=
,
{ }
min ,
i ii
µ αβ
=
( )
( )
{ }
11
11
max ,
i iiii
k fdfd
++
= −−
,
1
i ii
cdd
+
= +
证明:
对于
[ ]
1
,
ii
x xx
+
∈
,令
( )
ii
txx h
= −
,
( )( )
( )
ii
Ft fxt
=
,记
( )( )
( )
ii
W twx t
=
,其中
( )
i ii
x txht
= +
。
将由
(1)
式定义的函数
( )
i
Wt
改写成下列形式:
( )( )( )
i ii
WtHt Kt
=+
其中
( )()
( )
( )
( )
( )
32
23
111
13131
iiiiiii ii
Htft fdhttfdhttft
++ +
=−++−+ −−+
令
( )()
( )
( )
( )
( )
( )
( )
32
11
* 23
11 1
131 31
iii iiiiii
Htftffh ttffh ttft
++ +
=−++−+ −−+
( )
()
( )
2
2
1
1
1
ii iiii
i
ii
d hdh
Kttt
tt
αβ
αβ
+
−
= −
−+
符琳
|
一元四次有理插值样条曲线的应用
OPEN ACCESS
18
则有
( )( )( )( )( )( )( )( )( )( )( )
**
iii iiiiiii
FtQtPtFtHtKtFtHt HtHt Kt
−=−−≤−+−+
(7)
因为函数
( )
*
i
Ht
是关于
( )
i
Ft
在
[ ]
1
,
ii
xx
+
上的两点
Hermite
插值,所以有:
( )( )
[ ]
( )
( )
4
4
4
*
0,1
4
1d
Max
384 384
d
i
ii i
t
h
FtH tFtf
t
∈
−≤≤
(8)
由函数
( )
*
i
Ht
,
( )
Ht
,
( )
i
Kt
的定义并将其带入化简可得:
( )
( )
( )
( )
( )
( )
( )
11
*
11
11
4
ii
i iiiiii
kh
HtHhttfd tfdt
++
− =−−−−−≤
(9)
及
( )
( )
( )
2
2
1
1 16
i iiii
i ii
i
ii i
hd d
hc
Kttt
tt
αβ
λ
αβ µ
−
= −≤
−+
(10)
将
(8)
式,
(9)
式和
(10)
式带入
(7)
式,即可得定理的结论。
4.
数值实例
下面通过一个例子来验证本文所述的方法。
表
1
、表
2
的数据分别取自
Akima (1970)
和
Sarfraz (2000)
。
Table 1. Akima data
表
1.
Akima
数据
i
1 2 3 4 5 6 7 8 9 10
11
i
x
0 2 3 5 6 8 9 11
12
14
15
i
y
10
10
10
10
10
10
10.5
15
30
60
85
Table 2.
Sarfraz data
表
2. Sarfraz
数据
i
1 2 3 4 5
i
x
0 6
10
29.5
30
i
y
1
15
15
25
30
导数值
{ }
i
d
通过已知数据
( )
,,1,2,,
ii
xf in
=
,由下式确定:
( )
11
1
11
1
,2,3, ,
21
33
iii i
i
i
i
iii i
i
i
hh
d in
hhh h
β
α
β
α
−−
−
−−
−
∆+ ∆
= =
++ +
() ()
1
1
1121 12
0if 0
otherwise
d
hh h
∆=
=
∆+∆−∆ +
用
(1)
定义的保形样条可绘制出相应的插值曲线,形状参数
i
α
,
i
β
可根据需要交互式选取,以调整曲线的形
状。图
1
和图
2
为取
1
i
α
=
,
2
i
β
=
的插值曲线,曲线保持了被插数据的单调性且具有较好的视觉效果。
符琳
|
一元四次有理插值样条曲线的应用
OPEN ACCESS
19
Figure 1. Shape preserving interpolat ion curves of Akima data
图
1.
Akima
数据的保形插值曲线
Figure 2. Shape preserving interpolation curves of Sarfraz data
图
2. Sarfraz
数据的保形插值曲线
5.
小结
本文构造了四次分段有理插值曲线即分子、分母分别为四次、一次多项式的插值函数,并讨论了该函数的
一阶、二阶连续性和保单调的条件。该函数含有二族参数
α
i
,
β
i
可自由调节以改变曲线的形状,并可通过导数
d
i
赋值,无需求解线性方程组使得插值函数二阶连续。
致谢
本课题在选题和研究过程得到了王强老师的亲切关怀和悉心指导。在此向王老师致以诚挚的谢意和崇高的
敬意。同时感谢编者对文章提出宝贵意见,使文章更加的完整。
参考文献
(References)
[1]
Gregory
, J.A. and Delbourgo.
R.
(1982
)
Piecewise rational quadratic inte rpolation to monotonic data
.
IMA Journal
of Numerical Analysis
,
2
,
123
-130.
[2]
Delbourgo
,
R. (1989
)
Shape preserving interpolation to convex data by rational functions with quadratic numerator and linear denominator.
IMA Journal of Numerical Analysis
,
9
, 123-136.
[3]
王仁宏
,
吴顺唐
(1978
)
关于有理
spline
函数
.
吉林大学自然科学学报
,
1
, 58
-70.
符琳
|
一元四次有理插值样条曲线的应用
OPEN ACCESS
20
[4]
Duan, Q
, Djidjeli
, K
, Price
,
W.G.
and Twizell,
E.H.
(1998
)
A rational cubic spline based on function values
.
Computers & Graphics
,
22
,
479
-486.
[5]
Duan, Q.
, Wang
, L. and Twizell,
E.H.
(2005
)
A rational interpolation based on function values and constrained control of the interpolant curves.
Applied Mathematics and Computation
,
1
, 311
-322.
[6]
Hussain,
M.Z.
and Hussain, M. (2006
)
Visualization of data subject to positive constraints.
Information and Computi ng Science
,
1
, 149
-160.
[7]
王强
(2005
)
三次保形有理插值
.
合肥工业大学学报
(
自然科学版
),
11
, 1461
-1464.
[8]
方逵
(2010
)
一种新的二元有理插值及其性质
.
工程图学学报
,
4
, 117-122.