1. 引言
用Bernstein基函数构造Bezier曲线是计算机辅助几何设计中经典的方法 [1] - [8] ,其原因是Bernstein基函数具有良好的保形性质,而这些基函数与来自二项分布的概率密度函数密不可分,所以,概率密度函数在曲线曲面设计中得到重要运用。文献 [9] 中,Goldman研究了负二项分布所构建的Bezier曲线的性质,在文献 [10] 中,Goldman和Morin将泊松分布用于解析函数的几何重构。所以,CAGD领域的一些学者将二项分布进行推广,得到了新的基函数,从而得到了具有特别重要实用价值的算子。在众多广义分布中寻找合适的且计算量不大的概率密度函数具有重大意义。1993年,Madsen引入了一种具有用于CAGD几何造型潜力的密度函数 [11] 。2020年,杨军等人在此基础上又将密度函数的形式进行变形,进而得到了一种广义Bezier算子,该算子定义如下 [12] :
,
其中基函数定义为:
,
为参数,
,
,
。
注1. 当
时,规定
;当
或
时,规定
。
注2. 当
时,该算子退化为经典Bernstein算子
。
该算子也可以表示为:
。 (1)
在文献 [12] 中,杨军等人证明了该广义Bezier曲线具有经典Bezier曲线的几何不变性,通过分析导出了该类曲线的广义de Casteljau算法,从而得出可用于离散造型的细分算法。在此基础上,我们将研究该广义Bezier算子的保形性质。本文结构如下:首先在第二节中,我们介绍有关定理证明所需的一些引理,即引理2.1~2.3。第三节,证明该算子在单调性、凸性、星形性和半可加性等方面的保形性质。第四节,给出了有关光滑性定理的证明。第五节,通过具体函数给出数值实验,进一步验证新方法的有效性。
2. 所需引理
引理2.1.
1)
;
2)
。
证明:由(1)式,
1)
;
2)
。
引理2.2. [13] 对任何连续模
,存在上凸连续模
,使得对
,有
,其中常数2不能再减小,其中
定义如下:如果
是连续的,非减的,下半可加的,且
,则称
是连续模。
引理2.3. [1] 设
是
到自身的线性正算子列,I是有限或无穷区间,且
,
,
,若
是I上的上凸、单调递增、连续函数,则有
。
3. 有关保形的定理
首先,我们给出函数有关的几何性质的概念。
设函数
的定义域为
。
定义3.1. (单调性)
,
,若
,则称
在I上单调递增;若
,则称
在I上单调递减。
注3. 若
存在,则
单调递增
;
单调递减
。
函数的单调性在理论上可以解决不等式的证明,求最值,方程解的唯一性等问题。
定义3.2. (凸性)
,
,若
,则称
在I上为下凸函数;若
,则称
在I上为上凸函数。
注4. 若
存在,则
为I上为下凸函数
;
为I上为上凸函数
。
函数的凸性在解决极值,经济学,博弈论,运筹学,优化理论等领域的问题方面发挥了广泛而重要的作用。
定义3.3. (星形性)若
非负连续,
是
上不增(或不减)函数且
,则称
是关于原点的星形函数。
注5. 星形的定义也可为:设
非负连续,
,
,有
。
定义3.4. (半可加性)
,若
,则称
是下半可加的;若
,则称
是上半可加的。
函数这些保形性质,彼此之间相互独立,它们对曲线曲面保形重构拟合至关重要。
定理3.1. (单调性)若
在
上单调递增(或减),则对一切
,
在
上也单调递增(或减)。
证明:因为
,由(1)式,我们可以得到
(2)
若
在
上单调递增,由(2)式知
,所以
在
上单调递增;
若
在
上单调递减,由(2)式知
,所以
在
上单调递减。
定理3.2. (凸性)若
在
上下凸(或上凸),则对一切
,
在
上也下凸(或上凸)。
证明:由于
,
在(2)式基础上继续求导,
若
在
上下凸,则
,所以
在
上下凸;
若
在
上上凸,则
,所以
在
上上凸。
定理3.3. (星形性)若
是关于原点的星形函数,则对一切
,
也是关于原点的星形函数。
证明:
,
对上式求导,得
若
不增,则
,即
,所以
不增。
注6. 若
不减时,
也不减的情况同理可证。
定理3.4. (半可加性)若
是下半可加且
,则对一切
,
也是下半可加的。
证明:
,由(1)式,
令
,则
,
即
,下半可加性得证。
注7. 若
是上半可加且
时,
也是上半可加的情况同理可证。
4. 光滑性
对于
,记
为
的连续模 [13] 。
定理4.1. 对一切
,
,若
,则有
。
证明:由于
并注意到
,故
定理4.2. 若
,则
,
。
证明:
且
,
,类似于定理4.1证明过程中的讨论,
注8. 这里我们使用了引理2.1~2.3。
定理4.3. 若
是连续模函数,则对一切
,
也是
上的连续模函数。
证明:由
是连续模函数可知,
在原点函数值为零,且具有下半可加性、单调递增性,所以
也具有下半可加性、单调递增性、连续性。又
。综上所述,
也是
上的连续模函数。
5. 数值实验
由于本文讨论的广义Bezier算子中含有一个参数α,它可以对曲线的形状进行调控,从而使该类算子在CAGD领域具有应用潜力。对一些没有凸包性要求的情况,参数α有更大的取值范围,比如说,α可取负数,我们从模拟图中可以看到逼近的效果(见图1~4)。
下面分别取
(见图1、图2),
(见图3、图4),进行测试。
Figure 1. Estimate of the function
,
,
图1.
的逼近情况,
,
Figure 2. Estimate of the function
,
,#
图2.
的逼近情况,
,
Figure 3. Estimate of the function
,
,
图3.
的逼近情况,
,
Figure 4. Estimate of the function
,
,
图4.
的逼近情况,
,
基金项目
国家自然科学基金(11871191);河北省教育厅重点基金(ZD2019053);河北师范大学重点基金(L2020Z03)。
NOTES
*通讯作者。