韓力文， 楚 瑛， 李 丁， 劉 鳳

（1. 河北師范大學(xué)數(shù)學(xué)與信息科學(xué)學(xué)院，河北 石家莊 050024；2. 河北省計算數(shù)學(xué)與應(yīng)用重點實驗室，河北 石家莊 050024）

基于Lupas q-模擬Bernstein算子的廣義Bézier曲線

韓力文1,2， 楚 瑛1， 李 丁1， 劉 鳳1

（1. 河北師范大學(xué)數(shù)學(xué)與信息科學(xué)學(xué)院，河北 石家莊 050024；2. 河北省計算數(shù)學(xué)與應(yīng)用重點實驗室，河北 石家莊 050024）

提出了一種全新的廣義Bézier曲線。首先，從Lupas q-模擬Bernstein算子出發(fā)，得到了一組有理函數(shù)，該函數(shù)帶有一個形狀參數(shù)，是經(jīng)典Bernstein基函數(shù)的自然推廣。然后，構(gòu)造了相應(yīng)的廣義Bézier曲線，本文稱之為Lupas q-Bézier曲線，并研究了其基本性質(zhì)。Lupas q-Bézier曲線具有與經(jīng)典Bézier曲線相類似的升階公式和de Casteljau算法。

計算機(jī)輔助幾何設(shè)計；Lupas q-模擬Bernstein算子；Lupas q-Bézier曲線；升階公式；de Casteljau算法

經(jīng)典的Bernstein算子具有很好的逼近性、收斂性和保形性，是逼近論中最重要的算子之一，被廣泛應(yīng)用于分析、幾何和計算機(jī)等領(lǐng)域。近年來，隨著q-微積分的發(fā)展，一類基于q-整數(shù)的廣義Bernstein算子得到迅速發(fā)展。1987年，Lupas首次提出包含q-整數(shù)的廣義Bernstein算子，即Lupas q-模擬Bernstein算子[1]。1996年P(guān)hillips提出的q-Bernstein算子[2]是目前研究比較廣泛的廣義Bernstein算子。

1972年，Bézier[3]采用經(jīng)典的Bernstein基函數(shù)構(gòu)造了Bézier曲線，為自由型曲線曲面的發(fā)展奠定了堅實的基礎(chǔ)。隨著廣義Bernstein算子的產(chǎn)生，經(jīng)典 Bézier曲線也得到進(jìn)一步推廣。2003年，Oruc和Phillips利用q-Bernstein算子的基函數(shù)構(gòu)造了q-Bézier曲線[4]。Disibuyuk等人分別于2007年和2008年定義了有理q-Bézier曲線[5]及張量積型的q-Bézier曲面[6]。近幾年，國內(nèi)學(xué)者也研究并構(gòu)造了含形狀參數(shù)的廣義 Bézier曲線曲面[7-8]。Simeonov等人更深入地為兩類廣義Bézier曲線建立了相應(yīng)的開花形式和細(xì)分過程[9-10]。

相比而言，Lupas q-模擬 Bernstein算子在CAGD中的研究較少。對該算子的研究工作主要集中于逼近論方面。1987年，Lupas研究了Lupas q-模擬Bernstein 算子的逼近性和保形性[1]。2006年，Ostrovska討論了該算子的一致收斂性[11]。Phillips于2010年指出Lupas q-模擬Bernstein 算子還沒有任何應(yīng)用性研究[12]。

首先，介紹Lupas q-模擬Bernstein 算子的表示形式，并從中提取出Lupas q-模擬Bernstein有理函數(shù)，進(jìn)而構(gòu)造了Lupas q-Bézier曲線。Lupas q-Bézier曲線具有仿射不變性、凸包性、變差縮減性、保凸性等。特別是，Lupas q-Bézier曲線的變差縮減性、升階性和Lupas q-de Casteljau算法為Lupas q-模擬Bernstein算子在CAGD中的應(yīng)用奠定了基礎(chǔ)。

1 Lupas q-模擬Bernstein算子

為了介紹Lupas q-模擬Bernstein算子，首先引入以下記號和定義：

定義1.1[13]對于給定的實數(shù) 0q＞ ，及任意i N∈ ，定義q-整數(shù)[i]如下：

事實上，對于給定實數(shù) 0q＞ 且 1q≠ ，q -整數(shù)是關(guān)于q的有理函數(shù)，而當(dāng) 1q= 時，q -整數(shù)為通常意義下的非負(fù)整數(shù)。

定義1.2[13]對于給定的實數(shù) 0q＞ ，及任意i ∈ N，定義q -階乘[i] !如下：

定義1.3[13]對于給定的實數(shù) q＞ 0，及任意整數(shù) n ≥ i≥ 0，定義q -二項式系數(shù)如下：。

特別地，q -二項式系數(shù)滿足帕斯卡型遞推關(guān)系式

定義1.4[1]（Lupas q -模擬Bernstein算子）令 f( x) ∈ C[0,1]，線性算子定義為：

2 Lupas q-模擬Bernstein有理函數(shù)及其遞推公式

從Lupas q-模擬Bernstein算子中，我們提取出n次Lupas q-模擬Bernstein有理函數(shù)

本文稱之為q-逆對稱性。而且，當(dāng) q= 1時，Lupas q-模擬 Bernstein有理函數(shù)退化為經(jīng)典 Bernstein基函數(shù)。 bn(t; q ), t ∈ [0,1], i = 0,1,… ,n 是線性i無關(guān)的。

如圖 1所，當(dāng) q=1.5時五次 Lupas q-模擬Bernstein有理函數(shù)的圖像。

圖1 五次Lupas q-模擬Bernstein有理函數(shù)(q=1.5)

Lupas q-模擬Bernstein有理函數(shù)具有與經(jīng)典Bernstein基函數(shù)類似的遞推公式。

定理 2.1 n次Lupas q-模擬Bernstein有理函數(shù)可由兩個n+1次Lupas q-模擬Bernstein有理函數(shù)遞推得到，即

證明：

定理 2.2 n次Lupas q-模擬Bernstein有理函數(shù)可由兩個n-1次Lupas q-模擬Bernstein有理函數(shù)遞推得到，即

和

利用式（1），得

利用式（2），得

3 Lupas q-Bézier曲線及其性質(zhì)

定義3.1（Lupas q-Bézier曲線） 給定 n+1個向量 Pi∈ R2(i =0,1,… ,n)及實數(shù) q＞ 0，稱n次參數(shù)曲線段為一條n次Lupas q-Bézier曲線。 Pi稱為控制頂點。依次用直線段連接相鄰兩個 Pi,i =0,1,2,…,n ，所得的n邊折線多邊形稱為Lupas q-Bézier曲線的控制多邊形。

Lupas q-Bézier曲線具有一個形狀參數(shù)，在控制多邊形不變的情況下，通過調(diào)整參數(shù)q可以調(diào)控曲線的形狀。如圖2所示，分別取q=1, 0.3, 5, 9時三次Lupas q-Bézier曲線的圖像，當(dāng)q=1時，Lupas q-Bézier曲線即為經(jīng)典的Bézier曲線。

圖2 取q=1, 0.3, 5, 9的三次Lupas q-Bézier曲線

定理 3.1 Lupas q-Bézier曲線具有如下基本性質(zhì)：

1） 曲線是幾何不變和仿射不變的。

2） 曲線位于控制多邊形的凸包內(nèi)。

3） 曲線插值于控制多邊形首尾兩端點，即P (0;q) =P0, P (1;q ) =Pn。

這表明Lupas q-Bézier曲線以控制多邊形的首尾兩邊為其起點和終點的切方向。

4)（q-逆對稱性）如果將一條以q0為參數(shù)的Lupas q-Bézier曲線的控制頂點逆序排列，得到的新的Lupas q-Bézier曲線和以1/q0為參數(shù)，原順序控制頂點為頂點的曲線是同一條曲線。

5）（退化性）當(dāng) q= 1時，Lupas q-Bézier曲線退化為經(jīng)典Bézier曲線。

證明：下面只對q-逆對稱性進(jìn)行證明，其他4條性質(zhì)可根據(jù)相應(yīng)的Lupas q-模擬Bernstein有理函數(shù)的性質(zhì)推導(dǎo)得到。

定理 3.2(變差縮減性) Lupas q-Bézier曲線具有變差縮減性，即Lupas q-Bézier曲線與所在平面內(nèi)的任一直線的交點個數(shù)不會超過它的控制頂點與該直線的交點個數(shù)。

記Lupas q-Bézier曲線為C，任取其所在平面內(nèi)的一直線L，記C與L的交點數(shù)為 I( C, L)。以L為橫軸建立直角坐標(biāo)系，由Lupas q-Bézier曲線的幾何不變性，記控制頂點的新坐標(biāo)為(xi, yi)(i = 0,1,… ,n)，其中 0≤ x0≤x1≤…≤xn=1記控制多邊形為P，則P與L的交點數(shù)為 I( P, L)。下面只需證明 I( C, L ) ≤ I( P, L)。

推論 Lupas q-Bézier曲線具有保凸性，即當(dāng)控制多邊形是凸的，則所定義的 Lupas q-Bézier曲線也是凸的。保凸性可看作是變差縮減性的特殊情況。

根據(jù)Lupas q-模擬Bernstein有理函數(shù)的遞推性質(zhì)，可以推導(dǎo)出Lupas q-Bézier曲線的升階公式與de Casteljau算法，它們均是經(jīng)典Bézier曲線的升階公式和de Casteljau算法的推廣形式。

定理3.3(升階公式) 一條n次Lupas q-Bé zier曲線可以形式上看作一條 n+ 1次的 Lupas q-Bézier曲線，即

其中

其中

說明：當(dāng)q=1時，Lupas q-Bézier曲線的升階公式退化為經(jīng)典Bézier曲線的升階公式。若記原n階曲線的控制頂點組成的向量為

升階后的 n+ 1次曲線的控制頂點組成的向量為

則可以把升階過程表示為：P(1)=TP.

稱為升階算子，它是一個(n + 2)× (n + 1)階矩陣。對r ∈ N不斷升階為n + r次 Lupas q-Bézier曲線，可得控制頂點為

當(dāng)r→∞時，控制多邊形 P(r)收斂到Lupas q-Bézier曲線。

定理 3.4 (de Casteljau算法) 一條 n次Lupas q-Bézier曲線可表示為分別由前后n個控制頂點決定的兩條 1n- 次Lupas q-Bézier曲線的線性組合，進(jìn)而得到Lupas q-Bézier曲線上某一點遞歸求值的de Casteljau算法：

或

說明：當(dāng)q=1時，Lupas q-Bézier曲線的Lupas q-de Casteljau算法退化為經(jīng)典Bézier曲線的de Casteljau算法。若記

則de Casteljau算法可表示為：

其中 Mr(t; q)是一個(n - r + 1)× (n - r+ 2)階矩陣，且

或

4 結(jié) 論

本文利用Lupas q-模擬Bernstein有理函數(shù)族構(gòu)造了一種全新的廣義Bézier曲線，即Lupas q-B ézier曲線。研究了Lupas q-Bézier曲線的仿射不變性、凸包性、插值端點性、q-逆對稱性、變差縮減性及保凸性。推導(dǎo)出Lupas q-Bézier曲線的升階公式和de Casteljau算法。

在該文研究的基礎(chǔ)上，將進(jìn)一步探索 Lupas q-Bézier曲線與有理Bézier曲線、q-Bézier曲線的關(guān)系，并對Lupas q-Bézier曲線的開花形式與細(xì)分過程的構(gòu)造進(jìn)行更為深入的研究。

[1] Lupas A. A q-analogue of the Bernstein operator [R]. University of Cluj-Napoca, Seminar on Numerical and Statistical Calculus, Preprint, 1987, (9): 85-92.

[2] Phillips G M. On generalized Bernstein polynomials [J]. Numerical Analysis: A. R. Mitchell 75thBirthday Volume, 1996: 263-269.

[3] Bézier P E. Numerical control-mathematics and applications [M]. London: John Wiley & Sons, 1972.

[4] Oruc H, Phillips G M. q-Bernstein polynomials and Bézier curves [J]. Journal of Computational and Applied Mathematics, 2003, 151: 1-12.

[5] Disibuyuk C, Oruc H. A generalization of rational Bernstein-Bézier curves [J]. BIT Numerical Mathematics, 2007, 47: 313-323.

[6] Disibuyuk C, Oruc H. Tensor product q-Bernstein polynomials [J]. BIT Numerical Mathematics, 2008, 48: 689-700.

[7] Han Xi’an, Ma Yichen, Huang Xili. A novel generation of Bézier curve and surface [J]. Journal of Computational and Applied Mathematics, 2008, 271: 180-193.

[8] Chen Jie, Wang Guojin, A new type of the generalized Bézier curves [J]. Applied Mathematics-A Journal of Chinese Universities, 2011, 26(1): 47-56.

[9] Simeonov P, Zafiris V, Goldman R. h-Blossoming: a new approach to algorithms and identities for h-Bernstein bases and h-Bézier curves [J]. Journal of Computer Aided Geometric Design, 2011, 28: 549-565.

[10] Simeonov P, Zafiris V, Goldman R. q-Blossoming: a new approach to algorithms and identities for q-Bernstein bases and q-Bézier curves [J]. Journal of Approximation Theory, 2012, 164: 77-104.

[11] Ostrovska S. On the Lupas q-analogue of the Bernstein operator [J]. Journal of Mathematics, 2006, 36(5): 1615-1629.

[12] Phillips G M. A survey of results on the q-Bernstein polynomials [J]. MA Journal of Numerical Analysis, 2010, 30: 277-288.

[13] Andrews G E, Askey R, Roy R. Special functions [M].London: Cambridge University Press, 1999.

Generalized Bézier Curves Based on Lupas q-analogue of Bernstein Operator

Han Liwen1,2, Chu Ying1, Li Ding1, Liu Feng1
( 1. College of Mathematics and Information Science, Hebei Normal University, Shijiazhuang Hebei 050024, China; 2. Hebei Province Key Laboratory of Computational Mathematics and Application, Shijiazhuang Hebei 050024, China )

This paper presents a novel generalization of Bézier curves. Firstly, a class of rational functions with one shape parameter is presented. It comes from the Lupas q-analogue of Bernstein operator and is a natural extension to classical Bernstein basis. Then, the corresponding generalized Bézier curves, the so-called Lupas q-Bézier curves, are also constructed and their properties are studied. The new generalized Bézier curves share the degree evaluation and de Casteljau algorithm of the classical Bézier curves.

computer aided geometric design; Lupas q-analogue of Bernstein operator; Lupas q-Bézier curves; degree elevation; de Casteljau algorithm

O 241.5

A

2095-302X (2013)04-0063-06

2012-09-02；定稿日期：2012-11-06

國家自然科學(xué)基金資助項目（61170107）；河北省教育廳自然科學(xué)研究項目（Q2012041）

韓力文（1974-），女，河北石家莊人，副教授，博士，主要研究方向為計算機(jī)輔助幾何設(shè)計，數(shù)字幾何處理。E-mail：hanliwen@sina.com