• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    矩形域上的參數(shù)曲面自由變形算法

    2016-06-16 07:13:57葛先玉檀結(jié)慶

    張 莉 葛先玉 檀結(jié)慶,2

    1(合肥工業(yè)大學(xué)數(shù)學(xué)學(xué)院 合肥 230009)2(合肥工業(yè)大學(xué)計(jì)算機(jī)與信息學(xué)院 合肥 230009)(lizhang@hfut.edu.cn)

    矩形域上的參數(shù)曲面自由變形算法

    張莉1葛先玉1檀結(jié)慶1,2

    1(合肥工業(yè)大學(xué)數(shù)學(xué)學(xué)院合肥230009)2(合肥工業(yè)大學(xué)計(jì)算機(jī)與信息學(xué)院合肥230009)(lizhang@hfut.edu.cn)

    摘要針對(duì)參數(shù)曲面的自由變形,構(gòu)造了矩形域上一類帶平臺(tái)的、分片多項(xiàng)式形式的伸縮函數(shù).新的伸縮函數(shù)具備了對(duì)稱性、單點(diǎn)峰值性、線性峰值性和區(qū)域峰值性等良好的性質(zhì);且具備了簡(jiǎn)單的多項(xiàng)式形式,易于構(gòu)造和操控;由此構(gòu)造的伸縮因子中的各個(gè)參數(shù)具備明顯的幾何意義,特別適用于交互設(shè)計(jì).將基于伸縮因子構(gòu)造的變形矩陣作用于待變形的曲面,可以獲得諸如凹凸、剪切、伸縮和改變變形中心等各類變形效果,借助于疊加變形可進(jìn)一步實(shí)現(xiàn)豐富多樣的復(fù)雜變形效果.該文付出了大量的數(shù)值實(shí)例,展示了各類變形效果以及伸縮函數(shù)中各個(gè)不同參數(shù)對(duì)形狀的調(diào)控.

    關(guān)鍵詞參數(shù)曲面;伸縮因子;自由變形;矩形區(qū)域;變形效果

    在計(jì)算機(jī)圖形學(xué)領(lǐng)域中,物體的變形一直是一個(gè)熱點(diǎn)問題.從眾多的相關(guān)研究來看,我們可以將物體變形大致分為基于物理的變形[1]和基于幾何的變形.基于物理的變形方法是將物體看作物理實(shí)體而非單純的幾何體,根據(jù)客觀規(guī)律對(duì)物體進(jìn)行變形,得到較為逼真的效果.而基于幾何的變形方法則是將物體抽象為空間中的幾何體,運(yùn)用CADCAM中的數(shù)字幾何處理技術(shù)對(duì)物體形狀進(jìn)行修改、編輯.其數(shù)學(xué)本質(zhì)就是對(duì)其逐點(diǎn)映射以生成滿足特定條件的新的幾何形狀,與基于物理的變形方法相比較,該類方法計(jì)算量小、簡(jiǎn)單易行、但逼真度欠佳.

    有鑒于此,本文針對(duì)參數(shù)曲面的自由變形,通過構(gòu)造雙變量分片多項(xiàng)式型的伸縮函數(shù),給出了矩形域上一種交互良好,參數(shù)幾何意義明確,且能在點(diǎn)、線、矩形區(qū)域上達(dá)到峰值的參數(shù)曲面變形算法.這種算法無需依賴新的輔助工具,表達(dá)式中的每個(gè)參數(shù)都具備各自明顯的幾何屬性,可以在交互設(shè)計(jì)中靈活調(diào)用.

    1伸縮因子的定義和性質(zhì)

    1.1伸縮函數(shù)

    首先給出二元伸縮函數(shù)的定義.

    定義1. 設(shè)n為正整數(shù),r1,r2,s1,s2均為非負(fù)實(shí)

    數(shù),其中0≤r1

    D1={(x,y)||x-x0|≤r1,|y-y0|≤s1};

    D2={(x,y)|r1<|x-x0|≤r2,|y-y0|≤s1};

    D3={(x,y)||x-x0|≤r1,s1<|y-y0|≤s2};

    D4={(x,y)|r1<|x-x0|≤r2,s1<|y-y0|≤s2};

    記D=D1∪D2∪D3∪D4,定義2上伸縮函數(shù)如下:

    (1)

    為表述簡(jiǎn)便起見,這里的r=|x-x0|,s=|y-y0|.

    由定義1可見,二元伸縮函數(shù)g(x,y)的支撐區(qū)域D共有9個(gè)部分,其中D1為峰值區(qū)域,區(qū)域D2,D3,D4環(huán)繞在峰值區(qū)域D1的四周,具體分布如圖1所示:

    Fig. 1 Support area and peak region.圖1 支撐區(qū)域和峰值區(qū)域圖

    二元伸縮函數(shù)g(x,y)具有5個(gè)性質(zhì):

    Fig. 2  Peak properties of extension function g(x,y).圖2 伸縮函數(shù)g(x,y)的峰值性

    1) 區(qū)域峰值.(x,y)∈D1時(shí),g(x,y)取最大值1;(x,y)?D時(shí),g(x,y)取最小值0,取r1>0,s1>0,g(x,y)在D1上達(dá)到峰值,如圖2(a)所示.

    2) 線型峰值.任取r1和s1中一個(gè)為0但不同時(shí)為0,這時(shí)峰值區(qū)域D1退化為線段,g(x,y)在線段上達(dá)到峰值,如圖2(b)所示.

    3) 單點(diǎn)峰值.取r1=0,s1=0,這時(shí)峰值區(qū)域退化為單點(diǎn),g(x,y)具有單峰值,如圖2(c)所示.

    容易證明性質(zhì)1~3是成立的,性質(zhì)4,5的證明見附錄A.

    1.2伸縮因子

    定義2. 設(shè)h為實(shí)數(shù),定義帶參數(shù)n,h的二元伸縮因子:

    (2)

    其中,n為光滑參數(shù),h為伸縮參數(shù),D為支撐區(qū)域,最值區(qū)域?yàn)镈1.

    在式(2)的定義基礎(chǔ)上,可以得到伸縮因子具有5個(gè)與伸縮函數(shù)類似的性質(zhì):

    1) 區(qū)域最值.取r1和s1均大于0,當(dāng)h≥0時(shí),E(x,y)有最大值1+h,h<0時(shí),E(x,y)有最小值1+h.

    2) 線型最值.任取r1和s1中一個(gè)為0但不同時(shí)為0,E(x,y)在線段上達(dá)到最值.

    3) 單點(diǎn)最值.r1=0,s1=0,E(x,y)在單點(diǎn)取最值.

    4) 支撐區(qū)域D的邊界?D處的偏導(dǎo)值:

    5) 內(nèi)部交界處的偏導(dǎo)值:

    2空間曲面的變形與控制

    2.1變形的數(shù)學(xué)模型

    給定一張空間中的Cu(u≥1)類曲面P(x,y)=f(x,y)(亦可采用參數(shù)曲面形式),由定義2,伸縮因子為

    Eij(x,y)=1+hij×g(x,y),i,j=1,2,3,

    這組伸縮因子具有相同的支撐區(qū)域及光滑指數(shù),其中光滑指數(shù)n≤u+1.我們稱矩陣

    則稱F=(e1,e2,e3)MN(e1,e2,e3)T為變形矩陣,特殊地,取e1=(1,0,0)T,e2=(0,1,0)T,e3=(0,0,1)T時(shí),F=M.

    定義3. 在曲面定義域上選取中心點(diǎn)O′,變形后曲面Pd(x,y)與變形前曲面P(x,y)滿足關(guān)系:

    Pd(x,y)=F×(P(x,y)-O′)+O′,

    (3)

    P(x,y)-O′=(l1,l2,l3)(c1,c2,c3)T,

    則式(3)可以表示為

    (4)

    由此可以看出定義3中的變形有明確的幾何意義:在仿射坐標(biāo)系[O′,l1,l2,l3]下,對(duì)P(x,y)-O′的坐標(biāo)(c1,c2,c3)T作仿射變換,變換矩陣就是變形矩陣F.

    定理1. 待變形曲面為Bézier或NURBS曲面時(shí),變形后曲面為分片Bézier或分片NURBS曲面,即形式上仍為Bézier或NURBS曲面.

    證明. 作變換:

    (5)

    則定義2中伸縮因子可表示為

    (6)

    若待變形的曲面為K×S次Bézier曲面:

    它的控制網(wǎng)格為bi,j(i=0,1,…,K;j=0,1,…,S),在區(qū)域D4={(x,y)|r1<|x-x0|≤r2,s1<|y-y0|≤s2},此時(shí),0≤u<1,0≤v<1,作形如式(5)的變換可得:

    (7)

    (1+hun(2-u)nvn(2-v)n)×

    (8)

    注意到式(8)中的簡(jiǎn)單變形,最后2步借助于冪基和Bernstein基的相互轉(zhuǎn)換[20],由于伸縮因子為分片多項(xiàng)式,故而變形后的曲面為分片Bézier曲面.復(fù)雜的變形可類似處理.關(guān)于NURBS曲面的討論類似.

    2.2疊加變形

    通常一次變形不能得到理想的曲面,我們需要對(duì)同一個(gè)曲面進(jìn)行多次變形,這里假定為k次,令P0,d(x,y)=P(x,y). 由定義3可知,相鄰2次變形Pi,d(x,y),Pi-1,d(x,y),i=1,2,…,k,滿足遞推關(guān)系:

    (9)

    Fi=(ei1,ei2,ei3)MiNi(ei1,ei2,ei3)T,

    其中ei1,ei2,ei3和Mi分別為第i變形的主方向及伸縮矩陣,這里:

    特別地,當(dāng)變形中心均取同一點(diǎn)O′時(shí),經(jīng)過k次變形后曲面Pd(x,y)與初始曲面P(x,y)滿足關(guān)系:

    (10)

    2.3變形與控制

    在變形過程中,通過改變控制變形中心、峰值區(qū)域、光滑度、主方向等各個(gè)參數(shù),可以對(duì)待變形的初始曲面形狀進(jìn)行靈活的控制,以得到理想的曲面.

    1) 要控制曲面的變形中心,只需改變O′.

    2) 要控制曲面的峰值區(qū)域,可通過改變r(jià)1,s1.另外,r2和s2控制曲面變形的支撐區(qū)域.

    3) 要控制變形曲面在支撐區(qū)域邊界的光滑度,可通過改變n的大小,n越大,變形后曲面與待變曲面越靠近.

    4) 要控制曲面的變形主方向,只需改變e1,e2,e3.

    5) 改變hi i的符號(hào),可以控制變形曲面沿主方向ei的正向或負(fù)向變形,改變|hi i|的大小,控制曲面沿ei變形的幅度,其中i=1,2,3.改變h12,h13的符號(hào),可以控制曲面沿e1的正向或負(fù)向剪切,改變|h12|, |h13|的大小,可以得到不同的剪切效果,值越大,剪切效果越明顯.類似地,可以通過h21,h23控制曲面沿e2的剪切效果,通過h31,h32控制曲面沿e3的剪切效果.在實(shí)際應(yīng)用中,通常同時(shí)改變多個(gè)參數(shù),以達(dá)到理想的變形效果.

    3數(shù)值實(shí)例

    例1. 圖3、圖4以拋物面z=(2-x2-y2)8+2為例,取光滑指數(shù)n=4,依次展示了改變控制參數(shù)引起的實(shí)際效果.為了方便起見,在本例中取主方向?yàn)樽鴺?biāo)軸方向.圖3(a)為原始拋物面.取變形中心1=(0,0,0),參數(shù)r1=1.5,s1=1.5,r2=3.5,s2=3.5來確定變形的峰值區(qū)域及支撐區(qū)域,取參數(shù)hi j(i,j=1,2,3)中h33=2.5,其余為0,獲得圖3(b)所示變形效果,原始拋物面在支撐區(qū)域內(nèi)沿主方向e3(本例中即為z軸)正向伸展變形.上述參數(shù)不變,僅改變變形中心為2=(1,-1,0),所得伸展變形效果見圖3(c).取變形中心O′=(0,0,0),通過對(duì)參數(shù)r1,s1的調(diào)控可以讓初始拋物面達(dá)到不同類型峰值區(qū)域的效果,如圖4(a)(b)(c)所示,依次取得區(qū)域峰值、線型峰值、單點(diǎn)峰值.給定相同的支撐區(qū)域和峰值區(qū)域,取h11=2,得到曲面沿x軸凸出的效果圖(如圖4(d)所示),分別取h12,h13=1,得到如圖4(e)和圖4(f)兩種沿x軸方向不同的剪切效果.圖4中各個(gè)子圖具體參數(shù)賦值情況如表1所示.

    Fig. 3 Effects of changing deformation center O′.圖3 改變變形中心O′的效果

    Iconr1s1r2s2hij(i,j=1,2,3)Fig.4(a)1.51.533h33=-1,otherhij=0Fig.4(b)1.5042h33=2.5,otherhij=0Fig.4(c)0044h33=2.5,otherhij=0Fig.4(d)1122h33=2,h11=2,otherhij=0Fig.4(e)1122h33=2,h12=1,otherhij=0Fig.4(f)1122h33=2,h13=1,otherhij=0

    Fig. 4 Deformation effects after changing different parameters at the deformation center O′=(0,0,0)(see Table 1).圖4 O′=(0,0,0)處不同參數(shù)(如表1所示)變化得到的曲面變形效果

    例2. 以平面z=1(如圖5(a)所示)為例,展示了復(fù)雜變形效果.圖5(b)是在原平面基礎(chǔ)上沿著z軸連續(xù)4次上凸變形得到的.圖5(c)通過2次變形得到十字架模型.在圖5(d)中,對(duì)平面進(jìn)行4次疊加變形,用h33控制每次伸縮的長(zhǎng)度,h13控制其搖擺的角度.圖5(e)則展示了對(duì)平面連續(xù)5次伸縮得到類似于peaks函數(shù)曲面效果.在圖5(f)中我們同時(shí)控制h33,h11,使得變形沿z軸上凸且沿x軸正負(fù)方向拉伸,得到類似扇子的形狀.上述各圖表明,通過調(diào)控各參數(shù),變形操作可以在點(diǎn)、線和區(qū)域達(dá)到峰值,將變形復(fù)合和疊加,可以靈活地控制和調(diào)整變形曲面的形狀,得到豐富多彩的效果.

    Fig. 5 Complex deformation effects of the plane.圖5 平面上復(fù)雜變形效果

    Fig. 6 Complex deformation effects of cylindrical surface.圖6 圓柱面的復(fù)雜變形效果

    Fig. 7 Complex deformation effects of paraboloid.圖7 拋物面的復(fù)雜變形效果

    4結(jié)論

    本文針對(duì)參數(shù)曲面的自由變形,構(gòu)造了分片的雙變量多項(xiàng)式類型的伸縮函數(shù),考察了它在變形區(qū)域的邊界及內(nèi)部各個(gè)矩形區(qū)域交界處的偏導(dǎo)值,給出了一種交互性良好、幾何意義明確的新型參數(shù)曲面變形算法.該算法無需依賴任何輔助工具,不僅可以在單點(diǎn)處達(dá)到峰值,也可以在線、矩形區(qū)域上達(dá)到峰值,且表達(dá)式中的每個(gè)參數(shù)都具備各自明顯的幾何屬性.對(duì)于任意參數(shù)形式的曲面,我們可以通過調(diào)控各個(gè)變形參數(shù)得到期望的變形效果.

    本文所給出的參數(shù)曲面自由變形方法,本質(zhì)上是通過建立二元的伸縮函數(shù),借助于伸縮因子和變形矩陣來實(shí)現(xiàn)參數(shù)曲面的自由變形,盡管具有交互簡(jiǎn)單、幾何意義直觀等特點(diǎn),但變形效果也受伸縮函數(shù)的約束,無法實(shí)現(xiàn)現(xiàn)有自由變形算法中一些效果,如基于體的變形方法中可實(shí)現(xiàn)的任意拓?fù)涞淖冃巍⒍喾直媛首冃蝃10],基于曲面的變形方法中可實(shí)現(xiàn)的基于細(xì)分曲面的多分辨率變形[21]、基于掃掠曲面的多層次的聯(lián)動(dòng)變形[22]等.

    參考文獻(xiàn)

    [1]Che Yinghui, Liang Xiaohui, Zhao Qinping. A physically-based detail-preserving algorithm for real-time deformation[J]. Journal of Computer Research and Development, 2008, 45(3): 503-509 (in Chinese)(車英慧, 梁曉輝, 趙沁平. 一種基于物理的實(shí)時(shí)細(xì)節(jié)保持變形算法[J]. 計(jì)算機(jī)研究與發(fā)展, 2008, 45(3): 503-509)

    [2]Xu Gang, Wang Guozhao, Chen Xiaodiao. Free form deformation and its application[J]. Journal of Computer Research and Development, 2010, 47(2): 344-352 (in Chinese)(徐崗, 汪國(guó)昭, 陳小雕. 自由變形技術(shù)及其應(yīng)用[J]. 計(jì)算機(jī)研究與發(fā)展, 2010, 47(2): 344-352)

    [3]Barr A H. Global and local of deformation of solid primitive[J]. Computer Graphics, 1984, 18(3): 21-30

    [4]Ugüdükbay U, ōzgüc B. Free-form solid modeling using deformation[J]. Computer & Graphics, 1990, 14(34): 491-500

    [5]Sederberg T W, Parry R. Free-form deformations of solid geometric models[J]. Computer Graphics, 1986, 20(4): 151-160

    [6]Griessmair J, Purgathofer W. Deformation of solids with trivariate B-splines[C]Proc of Eurographics’89. New York: Wiley, 1989: 137-148

    [7]MacCracken R, Joy K I. Free-form deformations with lattices of arbitrary topology[J]. Computer Graphics, 1996, 30(4): 181-188

    [8]Huang Haiyun, Qi Feihu, Yao Zhihong. A leg motion modeling method based on NURBS free-form deformation[J]. Journal of Computer Research and Development, 2000, 37(6): 697-702 (in Chinese)(黃海赟, 戚飛虎, 姚志洪. 基于NURBS自由變形的腿部運(yùn)動(dòng)建模方法[J]. 計(jì)算機(jī)研究與發(fā)展, 2000, 37(6): 697-702)

    [9]Song Wenhao, Yang Xunnian. Free-form deformation with weighted T-spline[J]. The Visual Computer, 2005, 21(3): 139-151

    [10]Xu Gang, Wang Guozhao, Chen Xiaodiao. Free-form deformation with rational DMS-spline volumes[J]. Journal of Computer Science and Technology, 2008, 23(5): 862-873

    [11]Lazarus F, Coquillart S, Jancene P. Axial deformations: An intuitive deformation technique[J]. Computer Aided Design, 1994, 26(8): 607-613

    [12]Feng Jieqing, Ma Lizhuang, Peng Qunsheng. A new free-form deformation through the control of parametric surfaces[J]. Computer & Graphics, 1996, 20(4): 531-539

    [13]Feng Jieqing, Ma Lizhuang, Peng Qunsheng. Free-form deformation through embedding objects into the parametric domain of surface[J]. Journal of Computer-Aided Design & Computer Graphics, 1998, 10(3): 208-215 (in Chinese)(馮結(jié)青, 馬利莊, 彭群生. 嵌入?yún)?shù)空間的曲面控制自由變形方法[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào), 1998, 10(3): 208-215)

    [14]Zhao Yong, Xiao Chunxia, Shi Feng, et al. Rigid mesh deformation with detail preserving[J]. Journal of Computer Research and Development, 2010, 47(1): 1-7(in Chinese)(趙勇, 肖春霞, 石峰, 等. 保細(xì)節(jié)的網(wǎng)格剛性變形算法[J]. 計(jì)算機(jī)研究與發(fā)展, 2010, 47(1): 1-7)

    [15]Xu Gang, Hui Kinchuen, Ge Wenbing, et al. Direct manipulation of free-form deformation using curve-pairs[J]. Computer-Aided Design, 2013, 45(3): 605-614

    [16]Wang Xiaoping, Ye Zhenglin, Meng Yaqin, et al. Free-form deformation based on extension factor for parametric curve[J]. Journal of Computer-Aided Design & Computer Graphics, 2002, 14(1): 66-69 (in Chinese)(王小平, 葉正麟, 孫雅琴, 等. 基于伸縮因子的參數(shù)曲線自由變形[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào), 2002, 14(1): 66-69)

    [17]Wang Xiaoping, Ye Zhenglin, Meng Yaqin, et al. Space deformation of parametric surface based on extension function[J]. International Journal of CADCAM, 2002, 1(1): 23-32

    [18]Liu Zhi, Wu Hongyi, Zhang Li, et al. Method of polynomial factor for free-form deformation of parametric curves and surfaces[J]. Journal of Computer-Aided Design & Computer Graphics, 2009, 21(3): 412-418 (in Chinese)(劉植, 鄔弘毅, 張莉, 等. 參數(shù)曲線曲面自由變形的多項(xiàng)式因子方法[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào), 2009, 21(3): 412-418)

    [19]Lu Youtai, Zhou Laishui, Wang Zhiguo, et al. Free-form deformation based on power extension factor for parametric curves[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2012, 43(6): 793-798 (in Chinese)(陸友太, 周來水, 王志國(guó), 等. 基于乘冪伸縮因子的參數(shù)曲線自由變形[J]. 南京航空航天大學(xué)學(xué)報(bào), 2012, 43(6): 793-798) [20]Wang Guojin, Wang Guozhao, Zheng Jianmin. Computer Aided Geometric Design[M]. Beijing: China Higher Education Press, 2001: 14-15 (in Chinese)(王國(guó)瑾, 汪國(guó)昭, 鄭建民. 計(jì)算機(jī)輔助幾何設(shè)計(jì)[M]. 北京: 高等教育出版社, 2001: 14-15)

    [21]Feng Jieqing, Shao Jin, Jin Xiaogang, et al. Multiresolution free-form deformation with subdivision surface of arbitrary topology[J]. The Visual Computer, 2006, 22(1): 28-42

    [22]Yoon S H, Kim M S. Sweep-based free form deformation[J]. Computer GraphicsForum, 2006, 25(3): 487-496

    Zhang Li, born in 1976. PhD and professor in Hefei University of Technology. Member of China Computer Federation. Her main research interests include computer-aided geometric design & computer graphics.

    Ge Xianyu, born in 1991.Master candidate in Hefei University of Technology. His main research interests include computer-aided geometric design & computer graphics (shuxuegxy@126.com).

    Tan Jieqing, born in 1962. Professor and PhD supervisor in Hefei University of Technology. His main research interests include non-linear numerical approximation theory, computer-aided geometric design and computer graphics. (jieqingtan@hfut.edu.cn).

    附錄A.

    正文定義1中二元伸縮函數(shù)g(x,y)性質(zhì)4和性質(zhì)5的證明.

    證明.

    1) 性質(zhì)4證明.由正文圖1可知,支撐區(qū)域D的邊界是由x=x0±r2,y=y0±s2, 4條線段圍成,按照二元伸縮函數(shù)在各個(gè)區(qū)域上的表達(dá)式,將D的邊界如圖A1所示分成3類:虛線段、點(diǎn)劃線段和點(diǎn)線段.

    Fig. A1 Classified support area and peak region.圖A1 分類后的支撐區(qū)域和峰值區(qū)域圖

    (1) 求解虛線段處的偏導(dǎo)值

    二元伸縮函數(shù)g(x,y)在x=x0±r2作為區(qū)域D2邊界處的偏導(dǎo)值,即圖A1中虛線段.

    ①x=x0+r2時(shí),記:

    L1={(x,y)|x=x0+r2,y0-s1≤y≤y0+s1},

    由二元伸縮函數(shù)的定義知,(x,y)∈L1時(shí),

    g(x,y)=g1(x)=

    g(x,y)實(shí)際上是關(guān)于x的一元函數(shù),不妨記為g1(x),它關(guān)于y的偏導(dǎo)數(shù)值為0.從而當(dāng)l≠0時(shí),

    0≤k+l=i≤n-1,

    l=0時(shí),相當(dāng)于對(duì)函數(shù)g1(x)求解在x=x0+r2處直到n-1階的導(dǎo)數(shù)值,即:

    (x-x0-2r1+r2)n-k,

    其中:

    c=Cki(-1)i-kn(n-1)…

    (n-i+k+1)n(n-1)…(n-k+1),

    由于0≤i≤n-1,且n-i+k≥1,代入邊界信息,有:

    ②x=x0-r2時(shí),記:

    L2={(x,y)|x=x0-r2,y0-s1≤y≤y0+s1},

    由二元伸縮函數(shù)的定義知,當(dāng)(x,y)∈L2時(shí),

    g(x,y)=g2(x)=

    g(x,y)實(shí)際上是關(guān)于x的一元函數(shù),不妨記為g2(x),它關(guān)于y的偏導(dǎo)數(shù)值為零,即當(dāng)l≠0時(shí),

    l=0時(shí),求g2(x)在x=x0-r2處直到n-1階的導(dǎo)數(shù)值,類似于g1(x)的證明,我們有:

    (2) 求解點(diǎn)劃線段處的偏導(dǎo)值.

    從式(1)可以看出,g(x,y)在區(qū)域D2是x的單變量函數(shù),而在D3上恰好是y的單變量函數(shù),輪換變量x和y,上述對(duì)于D2邊界處偏導(dǎo)值的證明同樣適用于D3的邊界,從而可證明g(x,y)在點(diǎn)劃線段上偏導(dǎo)值為0.

    (3) 求解點(diǎn)線段處的偏導(dǎo)值

    圖A1中點(diǎn)線段為區(qū)域D4的邊界.由二元伸縮函數(shù)的對(duì)稱性,我們僅需考慮右上角的一部分,記:

    L3={(x,y)|x=x0+r2,y0+s1≤y≤y0+s2};

    L4={(x,y)|x0+r1≤x≤x0+r2,y=y0+s2};

    當(dāng)(x,y)∈L3,

    g(x,y)=g1(x)×g1(y)=

    g(x,y)為變量可分離的二元伸縮函數(shù),因此有:

    類似本證明步驟①②中關(guān)于x或者關(guān)于y的偏導(dǎo)值求解過程,可得:

    從而:

    這樣整個(gè)邊界?D的偏導(dǎo)值即證.

    2) 同理可以證明性質(zhì)5.

    證畢.

    Free Form Deformation Method of Parametric Surfaces on Rectangular Region

    Zhang Li1, Ge Xianyu1, and Tan Jieqing1,2

    1(SchoolofMathematics,HefeiUniversityofTechnology,Hefei230009)2(SchoolofComputerandInformation,HefeiUniversityofTechnology,Hefei230009)

    AbstractAccording to free form deformation of parametric surfaces, a new method based on extension function is proposed. It is made by piecewise polynomials and defined on rectangular region. Based on these, the new extension factor we constructed not only possesses perfect properties such as symmetry, single peak, linear peak and region peak, but also holds some parameters which have obvious geometric meanings. In real applications, the extension factor is very suitable for interactive design due to these properties and extraordinary parameters. Furthermore, the new extension factor is easy to construct and convenient to control because of its simple polynomial forms. Applying the deformation matrix constructed by the new extension factor to the original surfaces’ equations, plentiful deformation surfaces can be achieved. Deformation effects such as shearing, concave & convex, expand & contract and changing of deformation center can be obtained. With the help of superimposing, more complex deformation effects can be realized. Lots of numerical experiments are given at the end of paper which illustrate different kinds of deformation effects and plentiful adjusting effects of different parameters.

    Key wordsparametric surface; extension factor; free form deformation; rectangular region; deformation effects

    收稿日期:2014-11-28;修回日期:2015-08-11

    基金項(xiàng)目:國(guó)家自然科學(xué)基金-廣東聯(lián)合基金重點(diǎn)項(xiàng)目(U1135003);國(guó)家自然科學(xué)基金項(xiàng)目(61472466,61100126);安徽省自然科學(xué)基金項(xiàng)目(1508085QF116);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(JZ2015HGXJ0175);中國(guó)博士后科學(xué)基金項(xiàng)目(2015M571926)

    中圖法分類號(hào)TP391.41

    This work was supported by the National Natural Science Foundation of China-Guangdong Joint Foundation Key Project (U1135003), the National Natural Science Foundation of China (61472466,61100126), the Natural Science Fundation of Anhui Province of China (1508085QF116), and the Fundamental Research Funds for the Central Universities (JZ2015HGXJ0175), and the China Postdoctoral Science Foundation (2015M571926).

    波多野结衣高清作品| av专区在线播放| 女人被狂操c到高潮| 国产亚洲精品久久久com| 成人特级黄色片久久久久久久| 一区二区三区高清视频在线| 变态另类成人亚洲欧美熟女| 国产高清视频在线播放一区| 97碰自拍视频| 婷婷色综合大香蕉| 91麻豆精品激情在线观看国产| 国产欧美日韩精品亚洲av| 国产精品综合久久久久久久免费| 色哟哟哟哟哟哟| 女同久久另类99精品国产91| 免费观看在线日韩| 老熟妇仑乱视频hdxx| 国产高清视频在线播放一区| 男人舔奶头视频| 五月伊人婷婷丁香| 国产高清激情床上av| 丰满人妻一区二区三区视频av| 99热网站在线观看| 国产精品免费一区二区三区在线| 露出奶头的视频| 偷拍熟女少妇极品色| 男女之事视频高清在线观看| 亚洲精品日韩av片在线观看| 不卡一级毛片| av天堂在线播放| 无遮挡黄片免费观看| www日本黄色视频网| 午夜a级毛片| 国产国拍精品亚洲av在线观看| 最近视频中文字幕2019在线8| 一级毛片我不卡| 男人的好看免费观看在线视频| 最近2019中文字幕mv第一页| 亚洲第一电影网av| 亚洲成av人片在线播放无| 一个人看视频在线观看www免费| 亚洲欧美成人精品一区二区| 亚洲国产精品合色在线| 中文字幕av在线有码专区| 国产亚洲精品av在线| 亚洲最大成人中文| 看十八女毛片水多多多| 免费观看在线日韩| 成人毛片a级毛片在线播放| 午夜福利在线在线| 日韩精品有码人妻一区| 亚洲18禁久久av| 国产亚洲91精品色在线| 国国产精品蜜臀av免费| 男女之事视频高清在线观看| 最近手机中文字幕大全| 韩国av在线不卡| 女的被弄到高潮叫床怎么办| 网址你懂的国产日韩在线| 日韩中字成人| 成人午夜高清在线视频| 国产久久久一区二区三区| 婷婷亚洲欧美| 成人国产麻豆网| 狠狠狠狠99中文字幕| 又爽又黄无遮挡网站| av在线亚洲专区| 日韩国内少妇激情av| 国产乱人视频| 免费看美女性在线毛片视频| 又爽又黄无遮挡网站| 午夜福利成人在线免费观看| 国产黄色小视频在线观看| 成人性生交大片免费视频hd| 日韩欧美在线乱码| 日本免费一区二区三区高清不卡| 我要看日韩黄色一级片| 国产视频内射| 国产精品乱码一区二三区的特点| 久久这里只有精品中国| 男插女下体视频免费在线播放| 久久精品国产鲁丝片午夜精品| 国产白丝娇喘喷水9色精品| 少妇丰满av| 日日摸夜夜添夜夜爱| 中出人妻视频一区二区| 欧美日韩国产亚洲二区| 国产不卡一卡二| 天美传媒精品一区二区| 亚洲国产精品国产精品| 12—13女人毛片做爰片一| 免费av毛片视频| 男女之事视频高清在线观看| 欧美一区二区国产精品久久精品| 少妇高潮的动态图| 人人妻,人人澡人人爽秒播| 少妇的逼好多水| 日本一二三区视频观看| 久久人人爽人人片av| 国产色爽女视频免费观看| 午夜日韩欧美国产| 看免费成人av毛片| 白带黄色成豆腐渣| 国产aⅴ精品一区二区三区波| 可以在线观看毛片的网站| 国产精品国产高清国产av| 神马国产精品三级电影在线观看| 丰满的人妻完整版| 99在线人妻在线中文字幕| 变态另类丝袜制服| 一区二区三区四区激情视频 | 国产伦精品一区二区三区视频9| 变态另类丝袜制服| 嫩草影院入口| 亚洲一区二区三区色噜噜| 嫩草影视91久久| 久久精品夜夜夜夜夜久久蜜豆| 国产一区亚洲一区在线观看| 久久99热这里只有精品18| 国产一区亚洲一区在线观看| 久久综合国产亚洲精品| 国产一区二区亚洲精品在线观看| 高清日韩中文字幕在线| 婷婷色综合大香蕉| 久久精品人妻少妇| 日本黄色视频三级网站网址| 又黄又爽又免费观看的视频| 人人妻人人看人人澡| 亚洲在线观看片| av黄色大香蕉| 别揉我奶头 嗯啊视频| 美女 人体艺术 gogo| 欧美另类亚洲清纯唯美| 精品久久久噜噜| 18+在线观看网站| 久久久久国产精品人妻aⅴ院| 亚洲人成网站在线播放欧美日韩| 国产一区二区在线观看日韩| 国产精品精品国产色婷婷| 色播亚洲综合网| 国产精华一区二区三区| 国产精品久久电影中文字幕| 日韩在线高清观看一区二区三区| 我的老师免费观看完整版| 神马国产精品三级电影在线观看| 久久久久久久午夜电影| 美女高潮的动态| 看免费成人av毛片| 免费电影在线观看免费观看| 中文字幕精品亚洲无线码一区| 97超碰精品成人国产| 搡老岳熟女国产| 99久久无色码亚洲精品果冻| 人妻丰满熟妇av一区二区三区| 日本免费一区二区三区高清不卡| av在线天堂中文字幕| 男人狂女人下面高潮的视频| 全区人妻精品视频| 乱人视频在线观看| 亚洲av不卡在线观看| 婷婷精品国产亚洲av| 亚洲av美国av| 久久久久性生活片| 国产高清不卡午夜福利| 我要看日韩黄色一级片| 美女 人体艺术 gogo| 日本撒尿小便嘘嘘汇集6| 高清日韩中文字幕在线| 久久久精品欧美日韩精品| 狠狠狠狠99中文字幕| 麻豆国产av国片精品| 免费观看精品视频网站| 免费av不卡在线播放| 国产精品久久久久久亚洲av鲁大| 菩萨蛮人人尽说江南好唐韦庄 | 国产一区二区亚洲精品在线观看| а√天堂www在线а√下载| 亚洲最大成人手机在线| 级片在线观看| 久久久欧美国产精品| 变态另类成人亚洲欧美熟女| 欧美日韩在线观看h| АⅤ资源中文在线天堂| 欧美成人一区二区免费高清观看| 免费在线观看成人毛片| 在线播放国产精品三级| 美女被艹到高潮喷水动态| 尾随美女入室| 麻豆av噜噜一区二区三区| 久久亚洲国产成人精品v| 亚洲中文字幕一区二区三区有码在线看| 日韩三级伦理在线观看| 俺也久久电影网| 亚洲第一电影网av| 在线天堂最新版资源| 日韩三级伦理在线观看| 日韩高清综合在线| 欧美国产日韩亚洲一区| av天堂中文字幕网| 亚洲国产精品sss在线观看| 91在线精品国自产拍蜜月| 蜜桃久久精品国产亚洲av| 精品少妇黑人巨大在线播放 | 国产精品免费一区二区三区在线| 久久中文看片网| 国产亚洲91精品色在线| 俄罗斯特黄特色一大片| 毛片一级片免费看久久久久| 国内精品宾馆在线| 欧美中文日本在线观看视频| 亚洲图色成人| 精品久久久噜噜| 亚洲,欧美,日韩| av福利片在线观看| 久久久久久九九精品二区国产| 97碰自拍视频| 久久婷婷人人爽人人干人人爱| 人人妻人人澡欧美一区二区| 亚洲欧美日韩无卡精品| 最新中文字幕久久久久| 91麻豆精品激情在线观看国产| 国产精品爽爽va在线观看网站| 免费人成视频x8x8入口观看| 国产日本99.免费观看| 欧美潮喷喷水| 永久网站在线| 久久久午夜欧美精品| 色5月婷婷丁香| 日本-黄色视频高清免费观看| 精品一区二区免费观看| 嫩草影院新地址| 少妇熟女欧美另类| 国产探花在线观看一区二区| 国产成人精品久久久久久| 国产精品一及| 性插视频无遮挡在线免费观看| 亚洲欧美精品综合久久99| 不卡视频在线观看欧美| 人人妻人人澡欧美一区二区| 黄色视频,在线免费观看| 在线免费观看不下载黄p国产| 可以在线观看毛片的网站| 成人毛片a级毛片在线播放| 99在线人妻在线中文字幕| 国产成人aa在线观看| 欧美潮喷喷水| 身体一侧抽搐| 99热网站在线观看| 俺也久久电影网| 99热这里只有是精品50| 12—13女人毛片做爰片一| 舔av片在线| 成人二区视频| 国产精品亚洲一级av第二区| 亚洲国产精品成人综合色| 免费人成在线观看视频色| www日本黄色视频网| 能在线免费观看的黄片| 成人特级黄色片久久久久久久| 日本免费a在线| 人妻丰满熟妇av一区二区三区| 色噜噜av男人的天堂激情| 三级毛片av免费| 九九在线视频观看精品| 免费大片18禁| 午夜福利成人在线免费观看| 国语自产精品视频在线第100页| 少妇人妻一区二区三区视频| av视频在线观看入口| 日本三级黄在线观看| 国产一级毛片七仙女欲春2| 国产成人福利小说| 一级av片app| 国产视频一区二区在线看| 床上黄色一级片| 99热全是精品| 久久久久久久久久成人| 18禁黄网站禁片免费观看直播| 99九九线精品视频在线观看视频| 久久草成人影院| 亚洲第一区二区三区不卡| 日韩三级伦理在线观看| 国产中年淑女户外野战色| 免费av不卡在线播放| 国产一区二区激情短视频| 久久精品国产99精品国产亚洲性色| 国产精品一区二区三区四区免费观看 | 国产免费男女视频| 最新中文字幕久久久久| 久久久精品欧美日韩精品| 女的被弄到高潮叫床怎么办| 精品不卡国产一区二区三区| 97超视频在线观看视频| 能在线免费观看的黄片| 亚洲性久久影院| 国产精品av视频在线免费观看| 成人午夜高清在线视频| 校园春色视频在线观看| 欧美最新免费一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 看黄色毛片网站| 蜜臀久久99精品久久宅男| av视频在线观看入口| 乱码一卡2卡4卡精品| ponron亚洲| 乱人视频在线观看| 日本成人三级电影网站| 色吧在线观看| 一进一出抽搐动态| 成人永久免费在线观看视频| 热99在线观看视频| 亚洲成人久久爱视频| 伦精品一区二区三区| 男人舔女人下体高潮全视频| 最后的刺客免费高清国语| 亚洲无线观看免费| 黄色视频,在线免费观看| 一级毛片aaaaaa免费看小| 99热只有精品国产| 男女那种视频在线观看| 久久中文看片网| 亚洲人成网站在线播| 亚洲成人精品中文字幕电影| 可以在线观看的亚洲视频| 亚洲第一区二区三区不卡| 一区二区三区免费毛片| 三级男女做爰猛烈吃奶摸视频| 嫩草影院新地址| 成人精品一区二区免费| 欧美性猛交黑人性爽| 欧美一级a爱片免费观看看| 久久久成人免费电影| 午夜福利成人在线免费观看| 看非洲黑人一级黄片| 男人狂女人下面高潮的视频| 少妇人妻一区二区三区视频| 又爽又黄无遮挡网站| 日日啪夜夜撸| 日本撒尿小便嘘嘘汇集6| 久久热精品热| 别揉我奶头~嗯~啊~动态视频| 在线天堂最新版资源| 美女内射精品一级片tv| 亚洲成人av在线免费| 免费搜索国产男女视频| 欧美不卡视频在线免费观看| 免费人成视频x8x8入口观看| 亚洲av五月六月丁香网| 久久韩国三级中文字幕| 欧美人与善性xxx| av在线天堂中文字幕| 小蜜桃在线观看免费完整版高清| 亚洲精品国产成人久久av| 国产精品野战在线观看| 久久久精品94久久精品| 欧美高清成人免费视频www| 久久精品影院6| 淫妇啪啪啪对白视频| 我的女老师完整版在线观看| 老女人水多毛片| 国产精品一及| 精品久久久久久久久久免费视频| 免费电影在线观看免费观看| 精品一区二区三区av网在线观看| 国产欧美日韩精品亚洲av| 国产成人a区在线观看| 久久久久久久久大av| 国产久久久一区二区三区| 最新中文字幕久久久久| 色尼玛亚洲综合影院| 亚洲欧美成人综合另类久久久 | 午夜爱爱视频在线播放| 日韩一本色道免费dvd| 国产精品一区www在线观看| 两个人的视频大全免费| av女优亚洲男人天堂| 1000部很黄的大片| 最近2019中文字幕mv第一页| av中文乱码字幕在线| 热99re8久久精品国产| 嫩草影视91久久| 午夜激情欧美在线| 欧美日本亚洲视频在线播放| 久久草成人影院| 精华霜和精华液先用哪个| 男人狂女人下面高潮的视频| 男人舔女人下体高潮全视频| 亚洲国产精品国产精品| 日韩欧美三级三区| 成人一区二区视频在线观看| 综合色av麻豆| 特级一级黄色大片| 一区二区三区四区激情视频 | 国产不卡一卡二| 亚洲人成网站在线播| 国内揄拍国产精品人妻在线| 久久久久久久久久久丰满| 国产蜜桃级精品一区二区三区| 久久久久精品国产欧美久久久| 99热网站在线观看| 欧美一级a爱片免费观看看| 国产欧美日韩精品亚洲av| 一级毛片aaaaaa免费看小| 国产一区二区亚洲精品在线观看| 91在线观看av| 嫩草影院新地址| 日韩精品中文字幕看吧| 午夜福利在线观看吧| 国产精品亚洲一级av第二区| 1024手机看黄色片| 男人舔奶头视频| a级毛色黄片| 欧美日韩在线观看h| 99久国产av精品| 黄片wwwwww| 亚洲精品乱码久久久v下载方式| 男女做爰动态图高潮gif福利片| 内地一区二区视频在线| 亚洲av.av天堂| 观看免费一级毛片| 女生性感内裤真人,穿戴方法视频| 看黄色毛片网站| 深爱激情五月婷婷| 少妇高潮的动态图| 美女cb高潮喷水在线观看| 精品一区二区三区视频在线观看免费| 麻豆成人午夜福利视频| 又粗又爽又猛毛片免费看| 免费在线观看成人毛片| 精品久久久久久久久亚洲| 国产精品久久电影中文字幕| 丰满人妻一区二区三区视频av| 男女那种视频在线观看| 97在线视频观看| 欧美日本视频| 久久6这里有精品| 日本成人三级电影网站| 亚洲国产色片| 亚洲精品一区av在线观看| 欧美中文日本在线观看视频| 久久婷婷人人爽人人干人人爱| 老师上课跳d突然被开到最大视频| 国产男人的电影天堂91| 波多野结衣巨乳人妻| 久久综合国产亚洲精品| 国语自产精品视频在线第100页| 一区二区三区免费毛片| 日韩欧美 国产精品| 免费在线观看成人毛片| 久久精品久久久久久噜噜老黄 | 国产黄色视频一区二区在线观看 | 久久精品夜色国产| 高清午夜精品一区二区三区 | 亚洲av第一区精品v没综合| 插逼视频在线观看| 男女做爰动态图高潮gif福利片| 又爽又黄无遮挡网站| 亚洲精品日韩在线中文字幕 | 久久久久性生活片| 天天躁夜夜躁狠狠久久av| a级毛色黄片| 少妇被粗大猛烈的视频| 观看免费一级毛片| a级毛片a级免费在线| 久久久久久久久久成人| 97热精品久久久久久| 我的老师免费观看完整版| 免费在线观看成人毛片| 久久精品久久久久久噜噜老黄 | 色av中文字幕| 我的女老师完整版在线观看| 麻豆一二三区av精品| 午夜亚洲福利在线播放| 国产老妇女一区| 欧美色视频一区免费| 亚洲精品一区av在线观看| 精品99又大又爽又粗少妇毛片| 国产麻豆成人av免费视频| 国产三级在线视频| 欧洲精品卡2卡3卡4卡5卡区| АⅤ资源中文在线天堂| 国产精品99久久久久久久久| 国产精品久久视频播放| 国内少妇人妻偷人精品xxx网站| 免费看av在线观看网站| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看 | 永久网站在线| 国产日本99.免费观看| 一区福利在线观看| 青春草视频在线免费观看| 听说在线观看完整版免费高清| 老熟妇乱子伦视频在线观看| 久久久国产成人免费| 无遮挡黄片免费观看| 亚洲图色成人| 久久久久久久午夜电影| 欧美不卡视频在线免费观看| 国产欧美日韩一区二区精品| 亚洲乱码一区二区免费版| 亚洲美女搞黄在线观看 | 熟妇人妻久久中文字幕3abv| 97热精品久久久久久| 精品人妻偷拍中文字幕| 欧美另类亚洲清纯唯美| 波多野结衣高清作品| 搡老熟女国产l中国老女人| 欧美成人a在线观看| 日韩欧美 国产精品| 免费看日本二区| 免费不卡的大黄色大毛片视频在线观看 | 看非洲黑人一级黄片| 久久久久久久亚洲中文字幕| 亚洲性夜色夜夜综合| 波野结衣二区三区在线| 国产麻豆成人av免费视频| 久久精品国产自在天天线| 亚洲av成人精品一区久久| 亚洲最大成人av| 亚洲av五月六月丁香网| 精品人妻熟女av久视频| 色噜噜av男人的天堂激情| 成人综合一区亚洲| 亚洲av第一区精品v没综合| 日本黄色视频三级网站网址| 色视频www国产| 高清日韩中文字幕在线| 搡老熟女国产l中国老女人| 最后的刺客免费高清国语| 亚洲精品日韩在线中文字幕 | 三级男女做爰猛烈吃奶摸视频| 日韩 亚洲 欧美在线| 97人妻精品一区二区三区麻豆| 日本在线视频免费播放| 久久久久久久午夜电影| 免费av不卡在线播放| 最近的中文字幕免费完整| 少妇丰满av| av专区在线播放| 卡戴珊不雅视频在线播放| 久久久久久九九精品二区国产| 午夜精品国产一区二区电影 | 午夜亚洲福利在线播放| 国产一级毛片七仙女欲春2| 麻豆av噜噜一区二区三区| 亚洲一级一片aⅴ在线观看| 在线播放无遮挡| 插阴视频在线观看视频| 蜜臀久久99精品久久宅男| 精品人妻一区二区三区麻豆 | 精品国内亚洲2022精品成人| 两个人视频免费观看高清| 国产女主播在线喷水免费视频网站 | 日本-黄色视频高清免费观看| 成人鲁丝片一二三区免费| 国产一区亚洲一区在线观看| 不卡视频在线观看欧美| 欧美色视频一区免费| 久久精品国产自在天天线| 日日摸夜夜添夜夜添小说| 天天躁日日操中文字幕| 久久久午夜欧美精品| 老师上课跳d突然被开到最大视频| 美女xxoo啪啪120秒动态图| 99视频精品全部免费 在线| 久久久久性生活片| 99热只有精品国产| 亚洲最大成人中文| 精品一区二区三区视频在线| 色噜噜av男人的天堂激情| 真人做人爱边吃奶动态| 国模一区二区三区四区视频| 自拍偷自拍亚洲精品老妇| 国产在线精品亚洲第一网站| 欧美日韩乱码在线| 日韩欧美 国产精品| 高清毛片免费看| 嫩草影院精品99| 亚洲电影在线观看av| 欧美成人精品欧美一级黄| 女生性感内裤真人,穿戴方法视频| 人人妻人人澡人人爽人人夜夜 | 亚洲电影在线观看av| 成年女人毛片免费观看观看9| 国产国拍精品亚洲av在线观看| 亚洲国产精品合色在线| 热99re8久久精品国产| av免费在线看不卡| 欧美+日韩+精品| 亚洲国产色片| 久久久成人免费电影| 国产伦一二天堂av在线观看| 国产亚洲精品久久久com| 黄片wwwwww| 久久久久久久久中文| 日本一二三区视频观看| www日本黄色视频网| 赤兔流量卡办理| 春色校园在线视频观看| 亚洲成人精品中文字幕电影| 成人一区二区视频在线观看| 欧美xxxx黑人xx丫x性爽| 国产精品福利在线免费观看| 日韩大尺度精品在线看网址| 免费搜索国产男女视频| 免费看a级黄色片| 国产探花在线观看一区二区| 亚洲精品影视一区二区三区av| 人人妻人人澡欧美一区二区| 啦啦啦观看免费观看视频高清| 亚洲五月天丁香| av在线观看视频网站免费| 成人综合一区亚洲| 一区二区三区免费毛片| 成人三级黄色视频|