1. 引言
自由曲线曲面的研究一直以来受到学者们的广泛关注,Delgado和Peña [1] 提出了一类新型参数曲线,称为DP曲线,该曲线具有端点插值性、数值计算稳定性、线性计算复杂度、曲线保形性等性质。为了解决DP曲线的形状调节问题,陈杰等 [2] 提出带多个形状参数的DP-NTP曲线,形状参数可以较好地进行外形调节;陈福来 [3] 结合形状参数和DP基提出了另一种形式的广义三次DP基及曲线,并讨论了曲线参数连续和几何连续的拼接条件;彭兴璇等 [4] 通过在端点处的切矢引入形状参数,包括张力参数
和偏移参数
,构造了一类带形状参数的三次DP曲线;张迪等 [5] 将三次DP基函数的定义区间由
推广到
,重新参数化构造了α-DP曲线。
在几何外形设计和机器制造中,圆锥曲线占有重要的地位,由于圆锥曲线(抛物线段除外)不能用多项式曲线精确表示,许多学者将注意力集中在多项式曲线逼近包括圆弧在内的圆锥曲线上,Floater [6] 给出了任意奇数阶的Bézier曲线对圆锥曲线的逼近,有最佳逼近阶
;Ahn [7] 通过在端点处具有
或
连续性的四次Bézier曲线提供了圆锥曲线的两种逼近方法;胡倩倩 [8] 提出了一种任意阶的受约束的Bézier曲线逼近圆锥曲线的算法;Davood Bakhshesh [9] 提出了一种利用最小二乘法将端点插值的任意阶DP曲线进行圆锥曲线逼近的算法;张伟红等 [10] 人根据圆弧端点和端点切向量提出等弧长多项式逼近方法,给出了五次PH多项式逼近圆弧的精确表示。
在现有DP曲线及逼近算法基础上,本文构造了带形状参数的四次DP基函数,给出了带形状参数的四次DP曲线的定义和性质,明确了形状参数的几何意义,并在Davood Bakhshesh提出的算法基础上考虑带形状参数的四次DP曲线逼近圆锥曲线,改善了逼近曲线的灵活度,最后给出了数值例子。
2. 带形状参数的四次DP基函数
定义1对于
,
,称关于
的多项式
为带形状参数的四次DP基函数;其中,参数
为形状参数。
带形状参数的四次DP基函数具有如下性质:
1) 非负性,对
,
,有
。
2) 权性,对
,
,有
。
3) 端点性质,即
4) 退化性,当
且
时,带形状参数的四次DP基函数退化为四次DP基函数。
5) 最值性,
在
上均有一个最大值。
6) 拟对称性,当
且
时,
。
3. 带形状参数的四次DP曲线及性质
定义2给定5个控制顶点
,称曲线:
(1)
为带形状参数的四次DP曲线,其中
为带形状参数的四次DP基函数。
带形状参数的四次DP曲线具有如下性质:
1) 端点性质,即
2) 拟对称性,当
且
时,由控制多边形
和
生成的带形状参数的四次DP曲线相同,只是定向相反。
3) 凸包性,由带形状参数的四次DP基函数的非负性、权性知,曲线落在控制顶点生成的凸包内。
4) 几何不变性与仿射不变性,由带形状参数的四次DP基的规范性得到。
5) 退化性,当
且
时,带形状参数的四次DP曲线退化为四次DP曲线。
4. 形状参数对曲线形状的影响
(2)
其中,
为四次DP基函数,
,
,
,
,
。由式(2)知,带形状参数的四次DP曲线式(1)可看成以
为控制顶点的四次DP曲线。
当修改形状参数
的值时,相当于调整控制顶点
的位置。
由小变大时,
在
上向
靠近,
由小变大时,
在
上向
靠近。
当修改形状参数
的值时,相当于调整控制顶点
的位置。
、
、
是点
关于
的重心坐标,可以调节
在
里取值,从而灵活地调控生成曲线的形状。
5. 带形状参数的四次DP曲线逼近圆锥曲线
5.1. 圆锥曲线的矩阵表示
圆锥曲线可用标准的二次有理Bézier曲线表示为
其中,
是控制顶点。当
时,
表示椭圆;当
时,
表示抛物线;当
时,
表示双曲线;当
时,
表示原来圆锥曲线段的剩余部分。
令
,
,
,
,则有
。
5.2. 带形状参数的四次DP曲线的矩阵表示
将式(1)中带形状参数的四次DP曲线写成如下形式:
将上式改写为如下形式:
即
其中,
。
接下来令
,
,
,
,可以得到:
(3)
5.3. 带形状参数的四次DP曲线逼近圆锥曲线
在给定圆锥曲线
的情况下,利用带形状参数的四次DP曲线进行圆锥曲线逼近的关键问题是找到带形状参数的四次DP曲线的控制顶点
,Davood Bakhshesh [9] 在2014年提出了一种基于DP曲线的圆锥曲线逼近算法,该算法利用最小二乘法将端点插值的任意阶DP曲线进行圆锥曲线逼近,本文主要采用该算法求出带形状参数的四次DP曲线的控制顶点。
为了找到带形状参数的四次DP曲线的控制顶点
,令
,
,且在
范数下距离函数
最小化,
(4)
将式(4)表示为矩阵形式:
为了使
最小化,有
,即
化简为:
(5)
其中,
由于对任意的非零3维列向量
,有
,则H为正定矩阵,且H可逆。
将式(5)中
化简得到:
其中,
,
。
可按如下方法 [8] 计算:
,
设
,则
且
,
当
时,有
其中,
。
由此,可将式(5)化为:
(6)
故得,
(7)
因此,带形状参数的四次DP曲线的控制顶点
,可表示成如下形式:
(8)
5.4. 数值例子
例1:取圆锥曲线
的控制顶点为
,
,
,且
,圆锥曲线
为椭圆(如图1(a))。令
,利用上述算法求出带形状参数的四次DP曲线的控制点为
将带形状参数的四次DP曲线进行圆锥曲线逼近(如图1(b)),在
范数下计算误差距离函数,这里,我们引入最大误差
来评估近似结果,可得到
。
(a)(b)
Figure 1. (a) The conic section
; (b) The quartic DP curve with shape parameters when
,
,
,
图1. (a) 圆锥曲线
;(b)
,
,
,
时的带形状参数的四次DP曲线
例2:取圆锥曲线
的控制顶点为
,
,
,且
,圆锥曲线
为双曲线(如图2(a))。令
,
,
,
,
,求出带形状参数的四次DP曲线的控制点为
将带形状参数的四次DP曲线进行圆锥曲线逼近(如图2(b)),在
范数下计算误差距离函数,得到
。
(a)(b)
Figure 2. (a) The conic section
; (b) The quartic DP curve with shape parameters when
,
,
,
,
图2. (a) 圆锥曲线
;(b)
,
,
,
,
,时的带形状参数的四次DP曲线