圓域有理q-Bézier曲線

呂雁燕, 劉 植, 劉曉雁

(1.合肥工業(yè)大學(xué) 數(shù)學(xué)學(xué)院,安徽 合肥 230009; 2.拉文大學(xué) 數(shù)學(xué)系,加利福利亞 拉文 91750)

圓域有理q-Bézier曲線

呂雁燕1, 劉 植1, 劉曉雁2

(1.合肥工業(yè)大學(xué) 數(shù)學(xué)學(xué)院,安徽 合肥 230009; 2.拉文大學(xué) 數(shù)學(xué)系,加利福利亞 拉文 91750)

文章基于一類廣義Bernstein基函數(shù)定義了圓域有理q-Bézier曲線,通過改變參數(shù)q的取值,可以得到一類有理q-Bézier曲線簇,并研究了該類曲線的基本性質(zhì)及De Casteljau型算法,用二次有理q-Bézier曲線可精確表示圓錐曲線。該方法比現(xiàn)有方法更加靈活,且表示范圍更大。數(shù)值實(shí)例表明,圓域有理q-Bézier曲線的研究具有一定的理論意義與應(yīng)用價(jià)值。

圓域算術(shù);圓域有理q-Bézier曲線;中心曲線;De Casteljau型算法;圓錐曲線

曲線和曲面的表示是CAD及CAM重要的研究方向之一。經(jīng)典Bernstein多項(xiàng)式在函數(shù)逼近論、計(jì)算幾何以及概率論等領(lǐng)域均有著重要的應(yīng)用。近20年來關(guān)于Bernstein多項(xiàng)式的推廣研究已取得豐富成果[1]。文獻(xiàn)[2]首次提出一類基于q-整數(shù)的廣義Bernstein多項(xiàng)式,它是經(jīng)典Bernstein多項(xiàng)式的一種推廣形式。文獻(xiàn)[3]基于該多項(xiàng)式提出了一類新的多項(xiàng)式基函數(shù),稱之為q-Bernstein基,進(jìn)而得到新的參數(shù)多項(xiàng)式曲線——q-Bézier曲線。q-Bézier曲線包括了Bézier曲線,且具有與Bézier曲線類似的性質(zhì)。參數(shù)q的引入賦予了q-Bézier曲線更加豐富的形狀控制能力。然而和經(jīng)典Bézier曲線一樣,q-Bézier曲線也不能精確表示除拋物線外的圓錐曲線。文獻(xiàn)[4]基于De Casteljau算法研究了有理q-Bézier曲線,并在一定范圍內(nèi)給出了圓錐曲線的有理q-Bézier表示，但其結(jié)果并不完整,文獻(xiàn)[4]中圓錐曲線的構(gòu)造只對參數(shù)q進(jìn)行討論,而忽略了權(quán)因子ω。

幾何造型中大量運(yùn)算都是基于浮點(diǎn)運(yùn)算環(huán)境[5],因此幾何體的表示不夠精確,幾何計(jì)算也是近似的。為了保證數(shù)值運(yùn)算的穩(wěn)定性以及計(jì)算結(jié)果的精確性,人們常采用區(qū)間方法[6-7]。 用二維區(qū)間代替點(diǎn)進(jìn)行運(yùn)算,從而能保證理論上的精確結(jié)果含于計(jì)算結(jié)果中,在一定程度上避免重要信息的丟失。然而,文獻(xiàn)[8]指出區(qū)間算法有2個(gè)缺點(diǎn):① 區(qū)間在計(jì)算過程中會逐漸擴(kuò)大;② 矩形區(qū)間在二維空間中不具有旋轉(zhuǎn)對稱性,導(dǎo)致區(qū)間曲線不再具有仿射不變性。為了克服這些缺點(diǎn),文獻(xiàn)[9]用圓盤代替控制頂點(diǎn)引入了圓域Bézier曲線,該曲線可看作由Bézier曲線上的每一點(diǎn)與相應(yīng)半徑構(gòu)成的動(dòng)圓掃過的平面區(qū)域。圓域Bézier曲線也可看成一種帶誤差的Bézier曲線,當(dāng)所有的圓域半徑都相同時(shí)即為等距曲線。相對于傳統(tǒng)的NURBS曲線,上述各種Bézier曲線克服了基本算法中存在的數(shù)值不穩(wěn)定問題。 有理q-Bézier曲線能精確表示圓錐曲線,因此將圓域Bézier曲線向有理圓域q-Bézier曲線的推廣是一件有意義的工作。

本文提出圓域有理q-Bézier曲線的概念,并研究了其基本性質(zhì),得到q-Bézier曲線的De Casteljau算法。重點(diǎn)討論了圓錐曲線帶的圓域有理二次q-Bézier表示,通過改變q和ω的取值可以得到更加豐富的圓錐曲線帶,并給出參數(shù)q以及權(quán)因子取不同值時(shí)的圓域有理q-Bézier曲線表示圓錐曲線的分類圖,拓展了現(xiàn)有方法的使用范圍。

1 預(yù)備知識

1.1 圓域運(yùn)算

設(shè)R為全體實(shí)數(shù)集,R+為全體非負(fù)實(shí)數(shù)集,N為全體自然數(shù)集。平面圓域

可表示為〈P0〉=〈P0;r0〉。其中,(x0,y0)∈R2，為圓域的中心;r0∈R+，為圓域的半徑;P0表示以原點(diǎn)為起點(diǎn),點(diǎn)(x0,y0)為終點(diǎn)的向量。

對于任意2個(gè)圓域〈P1;r1〉、〈P2;r2〉和實(shí)數(shù)λ∈R,定義圓域運(yùn)算如下：

λ〈P1;r1〉=〈λP1;|λ|r1〉,

〈P1;r1〉+〈P2;r2〉=〈P1+P2;r1+r2〉。

故n+1個(gè)圓域的線性組合可表示為：

(1)

1.2q-Bernstein基

在參數(shù)曲線的基表示方法中,基函數(shù)具有重要的作用?；趒-Bernstein基的有理q-Bézier曲線是有理Bézier曲線的推廣。

對任意n∈N,q∈R+,n次q-Bernstein基函數(shù)[3]定義為：

(2)

其中,連乘符號中沒有乘積項(xiàng)時(shí)視為1;[i]定義為：

q-二項(xiàng)式系數(shù)定義為：

特別地,當(dāng)i=0時(shí),上式為1。

對任意q∈(0,1],q-Bernstein基在區(qū)間[0,1]上是一組全正基。四次q-Bernstein基函數(shù)曲線(q=0.5)與經(jīng)典四次Bernstein基函數(shù)曲線(q=1)的對比如圖1所示。

圖1 四次q-Bernstein基與四次Bernstein基

2 圓域有理q-Bézier曲線

2.1 概念

定義1給定n+1個(gè)平面圓域〈Pi〉=〈Pi;ri〉和ωi∈R+,i=0,…,n。對t∈[0,1],定義平面n次圓域有理q-Bézier曲線為：

(3)

其中,Bi,n(t)為(2)式定義的n次q-Bernstein基函數(shù);〈Pi〉 (i=0,1,…,n)為控制圓域;ωi(i=0,1,…,n)為權(quán)因子。特別地,當(dāng)ωi為常數(shù)時(shí),〈R〉(t)退化為圓域q-Bézier曲線，即

若記

并稱之為n次有理q-Bernstein基函數(shù),則由(1)式知,(3)式也可表示為：

〈R〉(t)=〈R(t);r(t)〉,t∈[0,1],

〈Pi〉w={(wix,wiy,wi)∈R3|

在三維空間定義非有理圓域q-Bézier曲線為：

或

其中,t∈[0,1];Bi,n(t)為(2)式定義的n次q-Bernstein基函數(shù)。對〈R〉w(t)使用中心投影變換,即得相應(yīng)的圓域有理q-Bézier曲線。

2.2 性質(zhì)

圓域有理q-Bézier曲線具有與有理Bézier曲線相似的性質(zhì)，具體如下:

(1) 端點(diǎn)插值。〈R〉(0)=〈P0〉,〈R〉(1)=〈Pn〉，這是由于

(2) 凸包性。當(dāng)所有的權(quán)因子ωi>0,參數(shù)q∈(0,1]時(shí),圓域有理q-Bézier曲線位于控制圓域的凸包內(nèi)。事實(shí)上,

即〈R〉(t)是圓域〈Pi;ri〉(i=0,…,n)的凸組合。

(3) 仿射不變性。圓域有理q-Bézier曲線的中心曲線是一條有理參數(shù)q-Bézier曲線,因此它具有仿射不變性。對于圓域半徑r(t)，其在平移、旋轉(zhuǎn)、反射條件下依舊成立。

(4) 形狀可調(diào)性。給定控制圓域和權(quán)因子,可以通過改變參數(shù)q的取值來調(diào)整圓域有理q-Bézier曲線的形狀。如給定如下控制圓域：

〈P0;r0〉=〈(9,10);2.3〉,〈P1;r1〉=〈(10,24);3〉,

〈P2;r2〉=〈(18,44);1.5〉,〈P3;r3〉=〈(30,46);2〉,

〈P4;r4〉=〈(35,30);1.5〉;

權(quán)因子

ω0=1,ω1=2,ω2=3,ω3=4,ω4=5。

形狀參數(shù)q=0.1和q=0.9時(shí)對應(yīng)的圓域有理四次q-Bézier曲線如圖2所示。

圖2 q取不同值時(shí)的圓域有理四次q-Bézier曲線

由圖2可以看出,參數(shù)q的增大使圓域有理q-Bézier曲線更接近控制圓域。

2.3DeCasteljau型算法

(2) 對s=1,2,…,n;i=0,1,…,n-s,計(jì)算

類似經(jīng)典De Casteljau算法,上述算法的中間圓域可表示為:

3 圓錐曲線的精確表示

圓錐曲線是CAD系統(tǒng)中重要的幾何造型工具。圓域有理q-Bézier曲線不但可以生成自由形態(tài)的曲線曲面,而且能精確表示圓錐曲線。由(3)式可知,二次圓域有理q-Bézier曲線可表示為:

〈R〉(t)=

(4)

其中,B0,2(t)=(1-t)(1-qt)、B1,2(t)=(1+q)×t(1-t)、B2,2(t)=t2為二次q-Bernstein基函數(shù)。

為了便于幾何描述,本文取

ω0=ω2=1,ω1=ω,

且所有的圓域半徑都相等。

首先,當(dāng)(1-q)+ω(1+q)=0時(shí),(4)式的分母B0,2(t)+ωB1,2(t)+B2,2(t)=(1-t)2+t2恒不為0,二次圓域有理q-Bézier曲線(4)表示橢圓。

其次,當(dāng)(1-q)+ω(1+q)≠0時(shí),因?yàn)?/p>

所以可以通過變換控制圓域和權(quán)因子,把圓域有理二次q-Bézier曲線表示為圓域有理Bézier曲線的形式。新的控制圓域和權(quán)因子分別為:

〈P0′〉=〈P0〉,〈P1′〉=

此時(shí)圓錐曲線可以按q和ω的取值分類如下:

(1) 若(ω1′)2<1,即-4<(q+1)(ω-1)<0,則〈R〉(t)表示橢圓。(1-q)+ω(1+q)=0恰好也包含在該結(jié)論中。

(2) 若(ω1′)2=1,即(q+1)(ω-1)=-4或者(q+1)(ω-1)=0,則〈R〉(t)表示拋物線。

(3) 若(ω1′)2>1,即(q+1)(ω-1)<-4或者(q+1)(ω-1)>0,則〈R〉(t)表示雙曲線。

注意到,當(dāng)q=-1時(shí),〈P1′〉=〈P0〉,得到一條從〈P0〉到〈P2〉的直線段帶。對任意實(shí)數(shù)ω，ω1′<0,即(1-q)+ω(1+q)<0,得到互補(bǔ)弧線,以相反的順序遍歷。當(dāng)ω1′<0或q<-1或q>1時(shí),曲線不滿足凸包性質(zhì)。

在幾何造型中,上述結(jié)果可采用分布圖的形式表示,如圖3所示。在由q和ω構(gòu)成的坐標(biāo)系中,“E” 表示橢圓區(qū)域,“H” 表示雙曲線區(qū)域,位于E和H之間的邊界曲線是拋物線區(qū)域。 特別地,直線q=1即為經(jīng)典有理二次Bézier曲線。

圖3 圓域有理二次q-Bézier曲線的形狀分布圖

因此,與有理二次Bézier曲線相比,有理二次q-Bézier曲線表示圓錐曲線的范圍更廣,除了可以利用ω控制圓錐曲線的形狀外,也可以利用q的取值實(shí)現(xiàn)。

對于給定的控制圓域〈P0〉=〈(-3,0);0.15〉、 〈P1〉=〈(2,4);0.15〉、〈P2〉=〈(3,0);0.15〉和權(quán)因子ω0=ω2=1、ω1=ω,圓域有理二次q-Bézier曲線生成的各種圓錐曲線帶如圖4所示。圖4a為固定q、調(diào)整ω得到的圓錐曲線帶;圖4b為固定ω、調(diào)整q得到的圓錐曲線帶。

圖4 q和ω的取值不同得到各種圓錐曲線帶

4 結(jié) 論

本文提出圓域有理q-Bézier曲線的概念,改變參數(shù)q的取值可以得到一簇圓域有理q-Bézier曲線。給出了圓域有理q-Bézier曲線的基本性質(zhì)及相應(yīng)的De Casteljau型算法。與有理Bézier曲線一樣,圓域有理q-Bézier曲線也可精確表示各種圓錐曲線。具體給出了圓錐曲線帶的圓域有理二次q-Bézier曲線表示方法以及形狀分布圖,經(jīng)典有理Bézier方法以及文獻(xiàn)[4]的結(jié)果均包含在形狀分布圖中,故本文的方法在幾何造型中增加了構(gòu)造的靈活性,同時(shí)拓展了構(gòu)造的范圍。數(shù)值實(shí)例表明,通過改變參數(shù)q和ω的值可以得到更加豐富的圓錐曲線帶。

[1] 達(dá)買力汗·胡爾曼哈吉,郭清偉.Bézier曲線的擴(kuò)展及其應(yīng)用[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,37(6):764-768.

[2] PHILLIPS G M.Bernstein polynomials based on theq-integers[J].Annals of Numerical Mathematics,1997,4:511-518.

[3] ORUC H,PHILLIPS G M.q-Bernstein polynomials and Bézier curves[J].Journal of Computational and Applied Mathematics,2003,151(1):1-12.

[4] DISIBUYUK C,ORUC H.A generalization of rational Bernstein-Bézier curves[J].BIT Numerical Mathematics,2007,47(2):313-323.

[5] PATRIKALAKIS N M.Robustness issues in geometric and solid modeling[J].Computer-Aided Design,2000,32(11):629.

[6] MUDUR S P,KOPARKAR P A.Interval methods for processing geometric objects[J].IEEE Computer Graphics and Applications,1984,4(2):7-17.

[7] SEDERBERG T W,FAROUKI R T.Approximation by interval Bézier curves[J].IEEE Computer Graphics and Applications,1992,15(2):87-95.

[8] CHEN F L,WU Y.Degree reduction of disk Bézier curves[J].Computer Aided Geometric Design,2004,21(3):263-280.

[9] LIN Q,ROKNE J G.Disk Bézier curves[J].Computer Aided Geometric Design,1998,15(7):712-737.

Diskrationalq-Béziercurves

(1.School of Mathematics, Hefei University of Technology, Hefei 230009, China; 2.Dept. of Mathematics, University of La Verne, La Verne 91750, USA)

Disk rationalq-Bézier curves are presented using a class of generalized Bernstein basis. A family of disk rationalq-Bézier curves can be obtained by changing the value of parameterq. The basic properties of this kind of rational curves are discussed. De Casteljau type algorithm of disk rationalq-Bézier curves is also considered. The quadratic disk rationalq-Bézier curve can be used to represent conic section accurately. The proposed method is more flexible than the existing methods, and it indicates a larger range. Some numerical examples demonstrate that the study of disk rationalq-Bézier curves is of theoretical importance and practical significance.

disk arithmetic; disk rationalq-Bézier curve; center curve; De Casteljau type algorithm; conic section

2016-07-07；

2016-10-25

國家自然科學(xué)基金資助項(xiàng)目(11471093);安徽省教育廳自然科學(xué)重大研究資助項(xiàng)目(KJ2014ZD30)；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)經(jīng)費(fèi)資助項(xiàng)目(JZ2015HGXJ0175);安徽省省級質(zhì)量工程專業(yè)綜合改革試點(diǎn)資助項(xiàng)目(2012zy007)和名師工作室資助項(xiàng)目(2015msgzs126)

呂雁燕(1993-),女,安徽東至人,合肥工業(yè)大學(xué)碩士生;

劉 植(1976-),男,安徽金寨人,博士,合肥工業(yè)大學(xué)副教授,碩士生導(dǎo)師，通訊作者,E-mail:liuzhi314@126.com;

劉曉雁(1962-),女,安徽合肥人,博士,拉文大學(xué)教授.

10.3969/j.issn.1003-5060.2017.09.026

TP391.41

A

1003-5060(2017)09-1289-05

(責(zé)任編輯 朱曉臨)