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

    改進(jìn)的超越函數(shù)分段線性逼近方法

    2016-07-19 19:17:09田征杜慧敏黃小康
    計(jì)算機(jī)應(yīng)用 2016年7期
    關(guān)鍵詞:斜率分段線性

    田征 杜慧敏 黃小康

    摘要:針對(duì)超越函數(shù)計(jì)算中所采用的分段線性逼近算法存在的無(wú)法提前確定精度及部分區(qū)間資源浪費(fèi)的問(wèn)題,提出一種改進(jìn)的分段線性逼近超越函數(shù)算法。該算法由預(yù)定義的逼近區(qū)間端點(diǎn)計(jì)算出用于逼近的線性函數(shù),根據(jù)被逼近函數(shù)的凹凸性對(duì)所計(jì)算線性函數(shù)進(jìn)行調(diào)整,在此基礎(chǔ)上計(jì)算出預(yù)定義逼近區(qū)間內(nèi)調(diào)整后函數(shù)與被逼近函數(shù)之間的最大誤差;按照所需精度的要求,自動(dòng)調(diào)整逼近區(qū)間,通過(guò)該過(guò)程的迭代,獲得了較少分段次數(shù)。算法結(jié)果在Matlab上進(jìn)行仿真,仿真結(jié)果表明,所提算法的分段數(shù)相比等分法減少了60%。所提算法在保證精度的前提下,降低了查找表(LUT)的資源消耗。

    關(guān)鍵詞:

    分段線性逼近;超越函數(shù);查找表;資源浪費(fèi);優(yōu)化分段方法

    中圖分類號(hào): TP391.75 文獻(xiàn)標(biāo)志碼:A

    0引言

    圖形處理器(Graphics Processing Unit, GPU)是各種嵌入式系統(tǒng)、個(gè)人機(jī)(Personal Computer, PC)、工作站和游戲機(jī)等不可缺少的重要部件。浮點(diǎn)超越函數(shù)單元是GPU數(shù)據(jù)通路中的重要部件,其性能直接影響圖形渲染效果[1-2]。

    目前在現(xiàn)場(chǎng)可編程門陣列(FieldProgrammable Gate Array, FPGA)中計(jì)算超越函數(shù)常用的方法有級(jí)數(shù)近似法、查表法(Look Up Table, LUT)、坐標(biāo)旋轉(zhuǎn)數(shù)字計(jì)算(Coordinate Rotation Digital Computer, CORDIC)算法和分段線性逼近法等。

    其中,級(jí)數(shù)近似法展開式較為復(fù)雜,硬件實(shí)現(xiàn)復(fù)雜,資源消耗較大。查表法雖然計(jì)算簡(jiǎn)單、易于實(shí)現(xiàn),但是所需存儲(chǔ)單元隨著計(jì)算精度的提高呈指數(shù)形式增加,資源消耗大[3]。CORDIC算法[4]作為一種便于FPGA實(shí)現(xiàn)的超越函數(shù)計(jì)算方法得到了廣泛關(guān)注,但其收斂速度與數(shù)據(jù)的表示精度成反比,當(dāng)精度要求較高時(shí),算法的迭代次數(shù)較多,計(jì)算延遲會(huì)增大。與之相比,分段線性逼近法[5-6]將低階多項(xiàng)式與較小的查找表相結(jié)合,資源消耗少,速度較快,被廣泛應(yīng)用于傳感網(wǎng)絡(luò)中的數(shù)據(jù)壓縮[7]、非線性模型到線性模型的轉(zhuǎn)換[8]以及圖形圖像處理[9]等領(lǐng)域,已成為在計(jì)算資源有限條件下超越函數(shù)計(jì)算的一種較為理想的選擇。

    為了提高函數(shù)的計(jì)算速度,在硬件設(shè)計(jì)中一般利用查找表實(shí)現(xiàn)分段線性逼近算法[6],段數(shù)越多,查找表越大,誤差則越小,因此,分段線性逼近算法需要在查找表、精度和分段數(shù)目之間尋求一個(gè)合理的平衡。其中有均方誤差法[10]、區(qū)間2k等分法[11]、面積法[12]等多種分段算法,但這些算法都是在計(jì)算完成后才可以知道分段精度,難以在精度已知的條件下完成分段數(shù)的計(jì)算。

    文獻(xiàn)[6]中提出了一種超越函數(shù)計(jì)算的最佳等距分段線性逼近(Optimal EquiDistant PieceWise Linear approximation,OED_PWL)方法,該方法可以通過(guò)調(diào)節(jié)分段數(shù)目來(lái)完成計(jì)算精度的靈活控制,從而可以在保證精度的條件下完成分段數(shù)的計(jì)算。相較區(qū)間2k等分法其算法性能、資源消耗等方面的表現(xiàn)更好;但是,根據(jù)文獻(xiàn)[6]提出的方法進(jìn)行tanh函數(shù)計(jì)算時(shí),部分區(qū)間的計(jì)算誤差較大,原因是tanh函數(shù)斜率在某些分段較平緩,而在另外一些分段斜率較陡峭,若采用均勻分段的方法,則存在分段平緩處的資源浪費(fèi)問(wèn)題。在后面的實(shí)驗(yàn)結(jié)果中,本文會(huì)與OED_PWL算法進(jìn)行對(duì)比。

    因此本文提出一種改進(jìn)的超越函數(shù)分段線性逼近方法,通過(guò)自動(dòng)識(shí)別逼近誤差進(jìn)行分段使其不斷細(xì)分,實(shí)現(xiàn)了一個(gè)精度已知條件下較優(yōu)的分段方法。

    1分段線性逼近算法基本原理

    分段線性逼近算法的原理是把非線性特性曲線分成若干個(gè)區(qū)段,在每個(gè)區(qū)段中用直線段近似地逼近特性曲線。算法原理如圖1所示。

    分段線性逼近算法是根據(jù)一定的分段方法,如面積法、等分法等,將曲線f(x)所在的一段區(qū)間[A,B]劃分為若干段,在每一段里面通過(guò)一個(gè)線性函數(shù)來(lái)進(jìn)行逼近:

    y(x)ax+b(1)

    其中:a表示線性函數(shù)的斜率;b表示線性函數(shù)的偏移量。其中,在每一段逼近區(qū)間[x1,x2]里面,區(qū)間的兩個(gè)端點(diǎn)需要滿足:

    {f(x1)=y(x1), f(x2)=y(x2)}(2)

    在通過(guò)某種方法得到[A,B]上超越函數(shù)曲線的分段情況以及各個(gè)段的逼近系數(shù)ai、bi后,即可完成超越函數(shù)的分段線性逼近。

    算法硬件實(shí)現(xiàn)分為3步:

    步驟1按照某種方法進(jìn)行分段區(qū)間的劃分,求取所有分段點(diǎn)。

    步驟2在每個(gè)分段區(qū)間,根據(jù)兩個(gè)端點(diǎn)的值求取直線段f(x)=ax+b,將每一段的ai和bi存放在查找表中以待使用。

    步驟3利用查找表,讀取每一段直線方程的系數(shù),通過(guò)乘法和加法計(jì)算每一段函數(shù)的逼近值,從而實(shí)現(xiàn)超越函數(shù)的分段線性逼近。

    2改進(jìn)的分段線性逼近算法

    2.1算法原理

    如何在給定的精度下,得到一個(gè)較好的分段,使得應(yīng)用分段計(jì)算值滿足給定的精度要求是分段線性逼近算法的一個(gè)難題。本文提出一種基于區(qū)間2k等分法,通過(guò)迭代來(lái)獲得分段區(qū)間的分段逼近方法,接下來(lái)對(duì)其原理進(jìn)行詳細(xì)介紹。

    區(qū)間2k等分法會(huì)將[A,B]上的曲線均勻地分為2k段,k的取值取決于查找表的大小及資源的分配情況。這是分段線性逼近中最常用的一種分段方法,具有實(shí)現(xiàn)簡(jiǎn)單、精度較高等優(yōu)點(diǎn);但是由于超越函數(shù)曲線的斜率在定義域內(nèi)是不均勻分布的,部分計(jì)算域內(nèi)曲線段斜率較大,因此采用直線段線性逼近誤差較大;而另外一部分計(jì)算域曲線段斜率較小,逼近誤差較小。為了滿足誤差較大曲線段的誤差要求,往往會(huì)導(dǎo)致分段數(shù)過(guò)多,導(dǎo)致斜率較小曲線段的資源浪費(fèi)。

    因此本文提出一種改進(jìn)型的超越函數(shù)分段線性逼近算法,可以在精度已知的前提下,完成對(duì)分段數(shù)的最優(yōu)選取,使得計(jì)算函數(shù)誤差在一定的范圍內(nèi)。算法的主要思想為:對(duì)斜率較大、誤差較大的函數(shù)區(qū)間進(jìn)行進(jìn)一步細(xì)分,避免了對(duì)曲線平滑部分分段的浪費(fèi),減少了查找表資源的消耗。

    以對(duì)數(shù)lb x在區(qū)間[1,5]的計(jì)算說(shuō)明算法的思想:首先整個(gè)逼近區(qū)間等分為[1,3]和[3,5]回復(fù):分段區(qū)間沒(méi)有問(wèn)題兩段,并對(duì)兩段分別進(jìn)行線性逼近,得到line1的兩條直線段,可以看出,在[1,3]區(qū)間的直線逼近誤差比較大,因此,對(duì)區(qū)間[1,3]繼續(xù)進(jìn)行細(xì)分,分為更小的兩個(gè)區(qū)間:[1,2]和[2,3]在等于2時(shí),兩個(gè)區(qū)間都包含在內(nèi),其中一個(gè)需用開區(qū)間來(lái)表示,請(qǐng)明確?;貜?fù):分段區(qū)間沒(méi)有問(wèn)題,并對(duì)其分別進(jìn)行線性逼近,得到兩條誤差較低的直線段line2,完成分段區(qū)間的細(xì)分過(guò)程。

    2.2算法實(shí)現(xiàn)流程

    本文實(shí)現(xiàn)的算法流程如圖3所示,通過(guò)迭代的方法來(lái)完成分段區(qū)間的計(jì)算。

    圖4說(shuō)明了算法的一次逼近過(guò)程,其中逼近區(qū)間段的中點(diǎn)處逼近誤差為d,假設(shè)原逼近直線y0=ax+b向曲線方向偏移d/2,得到新的逼近直線y1=ax+b+d/2,新的直線更加逼近曲線。

    本文在求取被逼近函數(shù)f(x)與新的逼近直線y1 (x)之間的最大誤差dmax=f(x)-y1(x),對(duì)dmax求導(dǎo)數(shù),令導(dǎo)數(shù)等于零,將得到x的值,代入dmax計(jì)算出最大誤差。

    在本文的改進(jìn)型分段線性逼近算法中,所有數(shù)據(jù)采用雙精度類型,在每個(gè)分段區(qū)間上,判斷被逼近函數(shù)的凹凸性以及該函數(shù)在區(qū)間內(nèi)是否存在拐點(diǎn),求取逼近直線以及被逼近函數(shù)之間的最大誤差,然后與給定誤差進(jìn)行比較:若小于給定誤差,則停止該段的繼續(xù)細(xì)分;若大于給定誤差,則對(duì)該段逼近區(qū)間繼續(xù)細(xì)分。

    若刪除,可從此處刪除偽代碼。

    下面是改進(jìn)型的分段線性逼近算法的C語(yǔ)言偽代碼。

    圖3的表述,與此處的“c語(yǔ)言偽代碼”?是否有重復(fù)表達(dá)現(xiàn)象,是否可以刪除某一種表述方式?請(qǐng)明確。

    回復(fù):其中分段區(qū)間沒(méi)有問(wèn)題,還有就是那個(gè)關(guān)于偽代碼和流程圖重復(fù)問(wèn)題,我和我們老師商量了,可以刪除也可以不刪除,如果要?jiǎng)h除,就把偽代碼部分刪除。(附件是刪除偽代碼的稿件)

    程序前

    void subsection(double x0, double x1)

    {

    double y0=x0坐標(biāo)處要逼近的函數(shù)值f(x0);

    double y1=x1坐標(biāo)處要逼近的函數(shù)值f(x1);

    y′=f ″(x)=sign/*對(duì)逼近函數(shù)求導(dǎo)判斷凹凸性,sign=1表示凸函數(shù),sign=0表示凹函數(shù)*/

    令f ′(c)=0求出c值,c則為函數(shù)逼近區(qū)間內(nèi)的拐點(diǎn);

    if(y′=0)/*函數(shù)有拐點(diǎn)*/

    {

    重新劃分逼近區(qū)間為[x0,c]及[c,x1],確保每段新的逼近區(qū)間都是單調(diào)的;

    }

    else;

    mid=(x0+x1)/2.0;

    a=(y1-y0)/(x1-x0);

    b=y1-a*(x1);

    max=d=f(mid)-(a*mid+b);/*mid坐標(biāo)處的逼近誤差*/

    if(sign==1)

    b1=b+d/2.0;

    else

    b1=b-d/2.0;

    y1(x)=ax+b1;//得到新的逼近直線y1

    d=f(x)-y//被逼近函數(shù)與逼近函數(shù)的差值

    d′=f ′(x)-y′//對(duì)函數(shù)d求導(dǎo)數(shù)

    d′=f(m)-y1′(m)=0//求函數(shù)的極值點(diǎn)m

    max=f(m)-y1(m);//逼近區(qū)間內(nèi)的最大誤差

    if (max<指定誤差){

    打印x0,x1坐標(biāo),并打印逼近直線參數(shù)m及n到指定文檔;

    }

    else {

    subsection(x0,mid);

    subsection(mid,x1);

    }

    }

    //說(shuō)明:x0與x1為要逼近的函數(shù)取值范圍;

    程序后

    3實(shí)驗(yàn)結(jié)果與對(duì)比

    為了更好地說(shuō)明問(wèn)題,采用1E-3和1E-4作為預(yù)設(shè)的精度,可以滿足移動(dòng)端圖形圖像處理等應(yīng)用的精度要求;取值范圍設(shè)定為[1,2],可以充分地說(shuō)明逼近的正確性。表1為各個(gè)超越函數(shù)在區(qū)間2k等分法和本文算法下得到的分段數(shù)量對(duì)比。

    可以看出,改進(jìn)后的分段線性逼近算法對(duì)斜率變化較大的超越函數(shù)作用會(huì)比較明顯。而且如果精度提升,2k分段數(shù)會(huì)劇烈增長(zhǎng),而本文提出的算法分段數(shù)增加較慢。

    表2為本文算法與其他兩種算法在相同實(shí)驗(yàn)條件下的分段結(jié)果比較,可以看出:在保證誤差精度的前提下,本文算法的分段數(shù)相比等分法減少了60%以上,相比OED_PWL法減少了20%。對(duì)于FPGA硬件實(shí)現(xiàn)而言,分段數(shù)越少,需要的查找表則越小,從而節(jié)省了硬件查找表資源。

    針對(duì)超越函數(shù)中的對(duì)數(shù)函數(shù)lb (1+x),文獻(xiàn)[9,13-14]等多篇文獻(xiàn)都對(duì)其進(jìn)行了分段區(qū)間的優(yōu)化。表3為本文算法對(duì)對(duì)數(shù)函數(shù)的分段逼近結(jié)果與其他文獻(xiàn)算法的對(duì)比,由于其他文獻(xiàn)的分段算法沒(méi)有做到精度已知下的分段劃分,因此實(shí)現(xiàn)精度各不相同。

    由表2與表3的實(shí)驗(yàn)結(jié)果可以得出:

    1)與其他已經(jīng)提出的算法相比,在相同精度下,本文的算法劃分的逼近區(qū)間更少,則查找到所屬分段的時(shí)間就會(huì)變少,所占據(jù)的存儲(chǔ)空間更小,更適合用硬件實(shí)現(xiàn)。

    2)根據(jù)表3,與文獻(xiàn)[9]相比,分段數(shù)目相差不多的情況下,本文的算法可以大大降低逼近誤差。

    綜上所述,改進(jìn)型的分段線性逼近算法同樣適用于tanh(x)函數(shù)、lb x函數(shù)等多種超越函數(shù)算法。

    4對(duì)數(shù)函數(shù)的FPGA實(shí)現(xiàn)

    完成算法的基本設(shè)計(jì)之后,本文將其應(yīng)用在某個(gè)超越函數(shù)運(yùn)算單元中,用于計(jì)算以2為底的對(duì)數(shù)轉(zhuǎn)換結(jié)果,其電路結(jié)構(gòu)如圖5所示。

    對(duì)數(shù)運(yùn)算的原理如下所示:

    lb x=k+lb (1+f)(3)

    其中,對(duì)于lb (1+f)進(jìn)行分段線性逼近,絕對(duì)值模塊計(jì)算a_in的絕對(duì)值。前導(dǎo)零檢測(cè)模塊通過(guò)檢測(cè)前導(dǎo)零的個(gè)數(shù)來(lái)得到首1的位置k′。然后通過(guò)取首1之后的剩余位數(shù)來(lái)得到尾數(shù)值f。 f的高6位被用作尋址查找表來(lái)逼近lb (1+f)的非線性部分。其中,n表示尾數(shù)的位數(shù),將首1的位置數(shù)減去尾數(shù)位數(shù),將會(huì)得到首1的實(shí)際位置k。最終的轉(zhuǎn)換結(jié)果會(huì)通過(guò)連接k值和逼近區(qū)域部分的值來(lái)得到。

    其中LUT存放查找表的系數(shù),直接關(guān)系到最終消耗的硬件資源的多少。文獻(xiàn)[9]采用了15段的分段區(qū)間劃分,而經(jīng)過(guò)算法的自動(dòng)優(yōu)化,在同等精度4.1E-003條件下,本文的分段區(qū)間為7段,大大降低了查找表的資源消耗。

    5結(jié)語(yǔ)

    本文實(shí)現(xiàn)的改進(jìn)后的超越函數(shù)分段線性逼近算法可以在精度可控條件下,實(shí)現(xiàn)一個(gè)在不同取值范圍內(nèi)分段逼近區(qū)間的較優(yōu)劃分,能夠有效地降低分段區(qū)間的數(shù)量,減少查找表資源消耗。計(jì)算了多種不同的超越函數(shù),并與其他文獻(xiàn)算法進(jìn)行對(duì)比,其結(jié)果充分表明:本文算法在保證計(jì)算精度的前提下,能夠有效地降低資源消耗,并且適用于不同的超越函數(shù)計(jì)算,而且其精度已知的特點(diǎn)也進(jìn)一步增加了算法的適用性和靈活性。同時(shí),本文提出了這種改進(jìn)的算法,已經(jīng)用于自主開發(fā)的GPU上,并經(jīng)過(guò)FPGA驗(yàn)證。接下來(lái)的研究方向是采用查找表與更高階的多項(xiàng)式相結(jié)合的方式對(duì)超越函數(shù)進(jìn)行逼近,研究分段區(qū)間中的函數(shù)最優(yōu)逼近方式。

    參考文獻(xiàn):

    [1]

    焦繼業(yè),穆榮,郝躍,等.面向移動(dòng)圖形頂點(diǎn)處理器的高性能低功耗定點(diǎn)特殊函數(shù)運(yùn)算單元[J].電子與信息學(xué)報(bào),2011,33(11):2764-2770.(JIAO J Y, MU R, HAO Y, et al. High performance and low power fixedpoint special function unit for mobile vertex processors [J]. Electronics & Information Technology, 2011, 33(11): 2764-277.)

    [2]

    NICKLOLLS J, DALLY W J. The GPU computing era [J]. IEEE Micro, 2010, 30(2): 56-69.

    [3]

    吳慶達(dá),何書專,潘紅兵,等.32位浮點(diǎn)數(shù)正余弦函數(shù)FPGA實(shí)現(xiàn)方法[J].微電子學(xué)與計(jì)算機(jī),2012,29(1):113-116.(WU Q D, HE S Z, PAN H B, et al. Implementations of 32bits fixed and floating point trigonometric functions with FPGA [J]. Microelectronics & Computer, 2012, 29(1): 113-116.)

    [4]

    鮑宜鵬.一種CORDIC算法優(yōu)化及32位浮點(diǎn)反正切函數(shù)FPGA實(shí)現(xiàn)[J].電子與封裝,2015,15(3):22-25.(BAO Y P.One improved CORDIC algorithm of calculating 32bit floating the arctangent functions with FPGA [J]. Electronics and Packaging, 2015, 15(3): 22-25.)

    [5]

    FLOREA C, VERTAN C. Piecewise linear approximation of logarithmic image processing models for dynamic range enhancement [J]. IEEE Transactions on Power Systems, 2011, 26(4): 2581-2583.

    [6]

    王少軍,張啟榮,彭宇,等.超越函數(shù)FPGA計(jì)算的最佳等距分段線性逼近方法[J].儀器儀表學(xué)報(bào),2014,35(6):1209-1216.(WANG S J, ZHANG Q R, PENG Y, et al. Optimal equidistant piecewise linear approximation algorithm for the computation of transcendental functions in FPGA [J]. Chinese Journal of Scientific Instrument, 2014, 35(6): 1209-1216.)

    [7]

    張建明,林亞平,傅明,等.傳感網(wǎng)絡(luò)中誤差有界的分段逼近數(shù)據(jù)壓縮算法[J].軟件學(xué)報(bào),2011,22(9):2149-2165.(ZHANG J M, LIN Y P, FU M, et al. Piecewise approximation based data compression algorithm with error bound in wireless sensor networks [J]. Journal of Software, 2011, 22(9): 2149-2165.)

    [8]

    江釗.分段線性逼近法在梯級(jí)水電站優(yōu)化調(diào)度中的應(yīng)用[D].武漢:華中科技大學(xué),2012:4-5.(JIANG Z. Application of piecewise linear approximation method to the optimal scheduling of casades hydroelectric plants [D]. Wuhan: Huazhong University of Science and Technology, 2012: 4-5.)

    [9]

    NAM B G, KIM H, YOO H J. Power and areaefficient unified computation of vector and elementary functions for handheld 3D graphics systems [J]. IEEE Transactions on Computers, 2008, 57(4): 490-504.

    [10]

    DUNHAM J G. Optimum uniform piecewise linear approximation of planar curves [J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1986, PAMI8(1): 67-75.

    [11]

    BAJGER M, OMONDI A. Lowerror, highspeed approximation of the sigmoid function for large FPGA implementation [J]. Journal of Signal Processing Systems for Signal Image and Video Technology, 2008, 52(2): 137-151.

    [12]

    周輝,李濤,邢啟江,等.數(shù)字曲線的線性逼近和分段識(shí)別[J].大連理工大學(xué)學(xué)報(bào),1997,37(5):576-580.(ZHOU H, LI T, XING Q J, et al. Linear approximation and piecewise identification of digital curves [J]. Journal of Dalian University of Technology, 1997, 37(5): 576-580.)

    [13]

    NAM B G, KIM H, YOO H J. A lowpower unified arithmetic unit for programmable handheld 3D graphics systems [J]. IEEE Journal of SolidState Circuits, 2007, 42(8): 1767-1778.

    猜你喜歡
    斜率分段線性
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    一類連續(xù)和不連續(xù)分段線性系統(tǒng)的周期解研究
    線性回歸方程的求解與應(yīng)用
    物理圖像斜率的變化探討
    物理之友(2020年12期)2020-07-16 05:39:16
    分段計(jì)算時(shí)間
    二階線性微分方程的解法
    求斜率型分式的取值范圍
    基于子孔徑斜率離散采樣的波前重構(gòu)
    3米2分段大力士“大”在哪兒?
    太空探索(2016年9期)2016-07-12 10:00:04
    MMC-MTDC輸電系統(tǒng)新型直流電壓斜率控制策略
    国产精品久久久久久亚洲av鲁大| 我要搜黄色片| 国产精品精品国产色婷婷| 国产精品久久久久久人妻精品电影| 久久九九热精品免费| 欧美中文日本在线观看视频| 国产激情偷乱视频一区二区| 身体一侧抽搐| 丰满人妻一区二区三区视频av | 女人被狂操c到高潮| 两个人看的免费小视频| 国内精品美女久久久久久| 国产真实伦视频高清在线观看 | 午夜福利欧美成人| 日本免费a在线| 国产亚洲av嫩草精品影院| 成熟少妇高潮喷水视频| 久久精品国产亚洲av涩爱 | 欧美性猛交黑人性爽| 一边摸一边抽搐一进一小说| 国产单亲对白刺激| 一个人免费在线观看电影| 精品久久久久久久末码| 国产高清激情床上av| 久久精品国产99精品国产亚洲性色| a级毛片a级免费在线| 亚洲国产精品sss在线观看| 国产午夜福利久久久久久| 亚洲中文字幕日韩| 亚洲成av人片在线播放无| 久久久久久久亚洲中文字幕 | 性色av乱码一区二区三区2| 亚洲真实伦在线观看| 久久久久国产精品人妻aⅴ院| 免费一级毛片在线播放高清视频| 国产午夜精品论理片| 欧美乱码精品一区二区三区| 亚洲avbb在线观看| 成人鲁丝片一二三区免费| 日本黄色片子视频| 亚洲欧美激情综合另类| 日日夜夜操网爽| 成年女人看的毛片在线观看| 香蕉av资源在线| 嫩草影院入口| 九色国产91popny在线| 黄色日韩在线| 免费电影在线观看免费观看| 国产亚洲av嫩草精品影院| 法律面前人人平等表现在哪些方面| 成人永久免费在线观看视频| 日韩人妻高清精品专区| 黄色成人免费大全| 国产精品久久电影中文字幕| 久久久久久人人人人人| 久久久久久九九精品二区国产| 日本免费一区二区三区高清不卡| 亚洲最大成人中文| 手机成人av网站| 手机成人av网站| 色老头精品视频在线观看| 精品不卡国产一区二区三区| 97碰自拍视频| 久久人人精品亚洲av| 一a级毛片在线观看| 国产成年人精品一区二区| 久久久久久久午夜电影| 久久6这里有精品| 精品欧美国产一区二区三| 少妇人妻精品综合一区二区 | 欧美日韩国产亚洲二区| 国产精品亚洲美女久久久| 亚洲无线观看免费| 国产精品久久久久久久久免 | 精品久久久久久久毛片微露脸| 亚洲人成电影免费在线| 国产精品99久久久久久久久| 99久久精品热视频| 母亲3免费完整高清在线观看| 国产成人啪精品午夜网站| 国产精品1区2区在线观看.| 久久久成人免费电影| 国模一区二区三区四区视频| 亚洲天堂国产精品一区在线| 亚洲精品在线美女| 亚洲av电影在线进入| 97超级碰碰碰精品色视频在线观看| 国产午夜精品论理片| 美女大奶头视频| 成人午夜高清在线视频| 久久久久精品国产欧美久久久| 亚洲中文字幕日韩| 亚洲狠狠婷婷综合久久图片| 亚洲av电影不卡..在线观看| 国产老妇女一区| 丰满人妻一区二区三区视频av | 精品不卡国产一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 成年女人毛片免费观看观看9| 手机成人av网站| 色综合亚洲欧美另类图片| 18禁裸乳无遮挡免费网站照片| 日本一本二区三区精品| 中文字幕av成人在线电影| 中文字幕熟女人妻在线| 久久精品国产99精品国产亚洲性色| 亚洲国产精品合色在线| 亚洲18禁久久av| 在线观看免费午夜福利视频| 国产精品亚洲美女久久久| eeuss影院久久| 欧美成人a在线观看| 九色国产91popny在线| 夜夜爽天天搞| 高清在线国产一区| 男人舔奶头视频| 性色av乱码一区二区三区2| 亚洲avbb在线观看| 全区人妻精品视频| 19禁男女啪啪无遮挡网站| 国产中年淑女户外野战色| 国产午夜精品论理片| 在线观看一区二区三区| 欧美日韩福利视频一区二区| 日韩精品中文字幕看吧| 老汉色av国产亚洲站长工具| 好男人电影高清在线观看| 久久精品综合一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩黄片免| a级毛片a级免费在线| 老司机午夜福利在线观看视频| 色综合婷婷激情| 精品人妻1区二区| 免费观看人在逋| x7x7x7水蜜桃| 丁香六月欧美| 免费看a级黄色片| 岛国视频午夜一区免费看| 桃红色精品国产亚洲av| 一本综合久久免费| 88av欧美| 成年女人毛片免费观看观看9| 国产精品1区2区在线观看.| 一二三四社区在线视频社区8| 狂野欧美白嫩少妇大欣赏| 熟妇人妻久久中文字幕3abv| 国产成年人精品一区二区| 国产中年淑女户外野战色| www日本在线高清视频| 他把我摸到了高潮在线观看| 无人区码免费观看不卡| av天堂中文字幕网| 久久精品综合一区二区三区| 国产av一区在线观看免费| 麻豆久久精品国产亚洲av| 最近视频中文字幕2019在线8| 亚洲欧美精品综合久久99| 在线观看免费视频日本深夜| 乱人视频在线观看| 国产 一区 欧美 日韩| 91久久精品电影网| 夜夜爽天天搞| 看黄色毛片网站| 一a级毛片在线观看| 午夜精品在线福利| 成人高潮视频无遮挡免费网站| 国产精品,欧美在线| 操出白浆在线播放| 欧美日韩综合久久久久久 | 日韩欧美 国产精品| 叶爱在线成人免费视频播放| xxx96com| 亚洲精品成人久久久久久| 性色avwww在线观看| 一级毛片高清免费大全| 国产精品久久久人人做人人爽| 最近在线观看免费完整版| 香蕉久久夜色| 美女黄网站色视频| 精品午夜福利视频在线观看一区| 国产精品综合久久久久久久免费| 母亲3免费完整高清在线观看| 免费人成在线观看视频色| 日韩成人在线观看一区二区三区| 午夜精品在线福利| 内射极品少妇av片p| 久久久久久久久大av| 一进一出抽搐动态| 国产精品久久久久久久电影 | 国产精品98久久久久久宅男小说| 噜噜噜噜噜久久久久久91| 亚洲男人的天堂狠狠| 色哟哟哟哟哟哟| 午夜亚洲福利在线播放| 啦啦啦韩国在线观看视频| av国产免费在线观看| 美女cb高潮喷水在线观看| 国产亚洲精品av在线| 久久精品国产亚洲av香蕉五月| 精品一区二区三区视频在线 | 亚洲成av人片免费观看| 757午夜福利合集在线观看| 男女床上黄色一级片免费看| 亚洲人成网站在线播| 成人av一区二区三区在线看| 中国美女看黄片| 午夜视频国产福利| 亚洲,欧美精品.| 欧美3d第一页| 亚洲专区中文字幕在线| 国产精华一区二区三区| 白带黄色成豆腐渣| 欧美一区二区亚洲| 免费在线观看日本一区| 免费av毛片视频| 国产视频内射| 国产精品综合久久久久久久免费| 日韩高清综合在线| www日本在线高清视频| 少妇的逼好多水| 婷婷丁香在线五月| av中文乱码字幕在线| 性色avwww在线观看| 91九色精品人成在线观看| 国产精品一及| 亚洲avbb在线观看| 精品久久久久久,| 久久九九热精品免费| 天天一区二区日本电影三级| 精品久久久久久久毛片微露脸| 国产av一区在线观看免费| 国产av在哪里看| 桃红色精品国产亚洲av| 欧美不卡视频在线免费观看| 国产精品精品国产色婷婷| 黑人欧美特级aaaaaa片| 午夜福利欧美成人| 性欧美人与动物交配| 国内精品美女久久久久久| 99国产精品一区二区蜜桃av| 国产伦在线观看视频一区| 国产高清视频在线观看网站| 人妻久久中文字幕网| 成人国产一区最新在线观看| 精品欧美国产一区二区三| 激情在线观看视频在线高清| 午夜福利免费观看在线| 脱女人内裤的视频| 亚洲,欧美精品.| 中出人妻视频一区二区| 国产av在哪里看| 欧美成狂野欧美在线观看| 在线十欧美十亚洲十日本专区| 变态另类丝袜制服| 国产成人a区在线观看| 午夜福利欧美成人| 少妇的丰满在线观看| 日韩免费av在线播放| 中亚洲国语对白在线视频| 亚洲国产日韩欧美精品在线观看 | 中文资源天堂在线| 99精品欧美一区二区三区四区| 99久久综合精品五月天人人| 亚洲精品一区av在线观看| 日本一二三区视频观看| 国产精品1区2区在线观看.| 欧美一区二区亚洲| 久久精品国产自在天天线| 午夜激情欧美在线| 99久久久亚洲精品蜜臀av| 最近最新免费中文字幕在线| 色老头精品视频在线观看| 国产精品爽爽va在线观看网站| 国产又黄又爽又无遮挡在线| 搡女人真爽免费视频火全软件 | 午夜福利高清视频| 欧美成狂野欧美在线观看| 国产黄片美女视频| 免费观看人在逋| 成人一区二区视频在线观看| 深夜精品福利| 两个人看的免费小视频| 宅男免费午夜| 欧美另类亚洲清纯唯美| 亚洲国产精品999在线| 日韩有码中文字幕| 99久久久亚洲精品蜜臀av| 免费看十八禁软件| 变态另类成人亚洲欧美熟女| 亚洲欧美日韩卡通动漫| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲无线在线观看| 麻豆一二三区av精品| 丰满人妻一区二区三区视频av | 久久国产乱子伦精品免费另类| 欧美日韩一级在线毛片| 免费在线观看影片大全网站| 一本一本综合久久| 伊人久久精品亚洲午夜| 最新中文字幕久久久久| 日韩 欧美 亚洲 中文字幕| 亚洲国产精品sss在线观看| 露出奶头的视频| 亚洲av不卡在线观看| 动漫黄色视频在线观看| 亚洲av美国av| 99久久精品国产亚洲精品| 国产亚洲精品av在线| 在线视频色国产色| 午夜福利在线观看吧| 成人欧美大片| 欧美+亚洲+日韩+国产| 乱人视频在线观看| 精品日产1卡2卡| 欧美不卡视频在线免费观看| 国产精品久久久人人做人人爽| 性色avwww在线观看| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 99热只有精品国产| 好男人在线观看高清免费视频| 免费在线观看日本一区| 日韩精品青青久久久久久| 18禁在线播放成人免费| 亚洲人成网站在线播| 一区二区三区高清视频在线| 老鸭窝网址在线观看| 成人精品一区二区免费| 啦啦啦免费观看视频1| 欧美黑人巨大hd| 成人三级黄色视频| 偷拍熟女少妇极品色| 99久国产av精品| 成人av一区二区三区在线看| 99国产综合亚洲精品| 欧美丝袜亚洲另类 | 亚洲国产精品sss在线观看| 99国产精品一区二区蜜桃av| 国产毛片a区久久久久| av在线天堂中文字幕| 亚洲国产欧美网| 欧美日韩中文字幕国产精品一区二区三区| 小蜜桃在线观看免费完整版高清| 一个人免费在线观看的高清视频| 看黄色毛片网站| 久久久久久国产a免费观看| 国产精品电影一区二区三区| 国产精品三级大全| 小蜜桃在线观看免费完整版高清| 亚洲狠狠婷婷综合久久图片| 久久精品亚洲精品国产色婷小说| 男女之事视频高清在线观看| 少妇高潮的动态图| 欧美成人性av电影在线观看| 少妇高潮的动态图| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| 99国产精品一区二区蜜桃av| 日韩精品中文字幕看吧| 色播亚洲综合网| 99久久九九国产精品国产免费| 精品免费久久久久久久清纯| 中文资源天堂在线| 午夜福利免费观看在线| 一本综合久久免费| 亚洲av熟女| 日本熟妇午夜| 欧美黑人巨大hd| 黄色成人免费大全| 久久性视频一级片| 成人av一区二区三区在线看| 热99re8久久精品国产| 久久精品人妻少妇| 成年免费大片在线观看| www日本黄色视频网| 中出人妻视频一区二区| 草草在线视频免费看| 国产成人啪精品午夜网站| 搡女人真爽免费视频火全软件 | 一级a爱片免费观看的视频| 性欧美人与动物交配| 91av网一区二区| 国产精品久久久久久久久免 | 国产蜜桃级精品一区二区三区| 日本在线视频免费播放| 99热只有精品国产| 国产欧美日韩精品一区二区| 亚洲av电影在线进入| 国产熟女xx| 日韩欧美免费精品| 国产精品国产高清国产av| 女人十人毛片免费观看3o分钟| 免费无遮挡裸体视频| 久久精品91无色码中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 18美女黄网站色大片免费观看| 一级黄色大片毛片| 给我免费播放毛片高清在线观看| 亚洲专区中文字幕在线| 亚洲一区高清亚洲精品| 一区福利在线观看| 国产 一区 欧美 日韩| 久久久久久久久久黄片| 亚洲av中文字字幕乱码综合| 最近视频中文字幕2019在线8| 好看av亚洲va欧美ⅴa在| 极品教师在线免费播放| 日韩国内少妇激情av| 欧美一区二区亚洲| 成人国产综合亚洲| 俺也久久电影网| 老汉色∧v一级毛片| 午夜精品久久久久久毛片777| 日韩免费av在线播放| 国产中年淑女户外野战色| 女人高潮潮喷娇喘18禁视频| 亚洲欧美精品综合久久99| 国产不卡一卡二| 欧美成人免费av一区二区三区| 日韩精品青青久久久久久| 亚洲成av人片免费观看| 亚洲国产精品合色在线| 99精品欧美一区二区三区四区| 亚洲avbb在线观看| av在线蜜桃| 桃红色精品国产亚洲av| 午夜两性在线视频| 亚洲av一区综合| 99久久精品一区二区三区| 久久精品国产清高在天天线| 欧美中文综合在线视频| 成人国产一区最新在线观看| 午夜精品在线福利| 亚洲精品亚洲一区二区| 久久精品国产99精品国产亚洲性色| 精华霜和精华液先用哪个| 欧美精品啪啪一区二区三区| 一边摸一边抽搐一进一小说| 国产爱豆传媒在线观看| 日本黄色片子视频| 精华霜和精华液先用哪个| 久久伊人香网站| 日韩欧美三级三区| 少妇的逼水好多| 国产伦人伦偷精品视频| 又紧又爽又黄一区二区| 麻豆国产av国片精品| 亚洲av成人不卡在线观看播放网| www国产在线视频色| 婷婷六月久久综合丁香| 精品一区二区三区av网在线观看| 淫秽高清视频在线观看| 成人高潮视频无遮挡免费网站| 岛国在线免费视频观看| 亚洲 欧美 日韩 在线 免费| 在线免费观看不下载黄p国产 | 国产乱人视频| 成人亚洲精品av一区二区| 免费电影在线观看免费观看| 亚洲中文日韩欧美视频| 亚洲国产精品成人综合色| 亚洲电影在线观看av| 美女 人体艺术 gogo| 3wmmmm亚洲av在线观看| 两个人看的免费小视频| 舔av片在线| 岛国在线观看网站| 嫩草影院入口| 婷婷精品国产亚洲av| 国产精品日韩av在线免费观看| 国产黄色小视频在线观看| 国产熟女xx| 俄罗斯特黄特色一大片| 午夜视频国产福利| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 哪里可以看免费的av片| 国产精品电影一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 国产精品电影一区二区三区| 日本 av在线| 蜜桃久久精品国产亚洲av| 国产精品香港三级国产av潘金莲| 日本三级黄在线观看| 热99在线观看视频| 国产精品 欧美亚洲| 宅男免费午夜| 欧美一级a爱片免费观看看| 亚洲熟妇中文字幕五十中出| 看免费av毛片| 国产一区二区激情短视频| 99精品在免费线老司机午夜| 哪里可以看免费的av片| 夜夜躁狠狠躁天天躁| 老汉色av国产亚洲站长工具| 久久久久久久亚洲中文字幕 | 亚洲激情在线av| 蜜桃久久精品国产亚洲av| 久久精品国产自在天天线| 日韩欧美在线二视频| 国产精品,欧美在线| 日韩欧美国产在线观看| 美女黄网站色视频| 日韩大尺度精品在线看网址| 亚洲内射少妇av| 精品久久久久久久人妻蜜臀av| 亚洲熟妇熟女久久| 啦啦啦观看免费观看视频高清| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 亚洲精品粉嫩美女一区| 在线国产一区二区在线| 舔av片在线| 国产伦在线观看视频一区| 久久久久久久精品吃奶| 亚洲精品一卡2卡三卡4卡5卡| 成人无遮挡网站| 亚洲中文字幕日韩| 好男人在线观看高清免费视频| 大型黄色视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 午夜福利18| 一卡2卡三卡四卡精品乱码亚洲| 丰满的人妻完整版| 国产aⅴ精品一区二区三区波| 精品乱码久久久久久99久播| 国产色婷婷99| 国产精品日韩av在线免费观看| 啦啦啦韩国在线观看视频| 麻豆国产97在线/欧美| 亚洲午夜理论影院| 国产精品女同一区二区软件 | 桃红色精品国产亚洲av| 成年女人看的毛片在线观看| av女优亚洲男人天堂| 老汉色∧v一级毛片| 国产亚洲精品久久久com| 亚洲中文字幕一区二区三区有码在线看| 美女cb高潮喷水在线观看| 午夜福利在线观看免费完整高清在 | xxxwww97欧美| 亚洲精品一区av在线观看| 性欧美人与动物交配| 精华霜和精华液先用哪个| 亚洲久久久久久中文字幕| 国产蜜桃级精品一区二区三区| 国内久久婷婷六月综合欲色啪| 久久中文看片网| 日韩欧美国产在线观看| 深夜精品福利| 国产aⅴ精品一区二区三区波| 日韩欧美精品v在线| 国产视频内射| 亚洲精品成人久久久久久| 黄色片一级片一级黄色片| 一级黄色大片毛片| 国产中年淑女户外野战色| 黄色女人牲交| 亚洲国产日韩欧美精品在线观看 | 12—13女人毛片做爰片一| 99久国产av精品| 国产一区二区在线av高清观看| 色老头精品视频在线观看| 成人av在线播放网站| 国产高清激情床上av| 在线播放无遮挡| 给我免费播放毛片高清在线观看| 村上凉子中文字幕在线| 淫妇啪啪啪对白视频| 欧美激情在线99| 国产三级中文精品| 精品久久久久久久久久久久久| 偷拍熟女少妇极品色| 精品一区二区三区人妻视频| 国产精品久久电影中文字幕| 亚洲狠狠婷婷综合久久图片| 亚洲精品在线美女| 成人亚洲精品av一区二区| 九九热线精品视视频播放| 熟女电影av网| 99久国产av精品| 久久精品综合一区二区三区| 国产精品亚洲一级av第二区| 少妇的丰满在线观看| 国产一区二区激情短视频| 在线观看免费视频日本深夜| 久久精品国产99精品国产亚洲性色| 亚洲熟妇熟女久久| 三级男女做爰猛烈吃奶摸视频| 亚洲成a人片在线一区二区| 天堂动漫精品| h日本视频在线播放| e午夜精品久久久久久久| 天堂动漫精品| 一个人看的www免费观看视频| av片东京热男人的天堂| 国产高清视频在线观看网站| 亚洲色图av天堂| 桃色一区二区三区在线观看| 亚洲精品美女久久久久99蜜臀| 少妇高潮的动态图| 97超级碰碰碰精品色视频在线观看| 国产精品一区二区三区四区久久| 国产精品国产高清国产av| 久久久久性生活片| 国产欧美日韩精品亚洲av| 变态另类成人亚洲欧美熟女| 免费看日本二区| 欧美日韩中文字幕国产精品一区二区三区| 中文亚洲av片在线观看爽| 成人亚洲精品av一区二区| 乱人视频在线观看| 99久久久亚洲精品蜜臀av| 中文字幕人妻熟人妻熟丝袜美 |