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

    在C語言中雙精度浮點(diǎn)數(shù)線性化相等比較的研究

    2017-10-13 04:23:24章國華
    船電技術(shù) 2017年1期
    關(guān)鍵詞:浮點(diǎn)數(shù)線性化位數(shù)

    章國華

    ?

    在C語言中雙精度浮點(diǎn)數(shù)線性化相等比較的研究

    章國華

    (武漢船舶職業(yè)技術(shù)學(xué)院,武漢430050 )

    以LCC軟件為設(shè)計基礎(chǔ),提出了雙精度浮點(diǎn)數(shù)比較大小算法的設(shè)計方法,通過程序驗證了算法的精確度,分別分析了絕對誤差和相對誤差方法設(shè)計算法,最后從雙精度浮點(diǎn)數(shù)在計算機(jī)中的表示方法角度提出了浮點(diǎn)數(shù)線性化后相等比較改進(jìn)的方法。

    比較 算法 雙精度 有效數(shù)字

    0 引言

    C程序設(shè)計的教科書中一般都不涉及浮點(diǎn)數(shù)的運(yùn)算,浮點(diǎn)數(shù)的運(yùn)算涉及到計算機(jī)的復(fù)雜的硬件結(jié)構(gòu)和理解起來有一定難度的表示方法[1]。然而浮點(diǎn)數(shù)的運(yùn)算在計算機(jī)軟件系統(tǒng)中可以說是無處不在的。大多數(shù)編程語言都將浮點(diǎn)數(shù)據(jù)類型作為基本類型,從計算機(jī)的CPU硬件到軟件編譯器、操作系統(tǒng)都涉及到浮點(diǎn)數(shù)的運(yùn)算。當(dāng)然,在用各種編程語言實現(xiàn)各種應(yīng)用程序設(shè)計的時候,即數(shù)據(jù)的處理,大都會碰到浮躁數(shù)的運(yùn)算[2]。例如,在用C語言寫應(yīng)用軟件的時候,各類整數(shù)的相等比較是直接使用“==”來判斷,正是由于浮點(diǎn)數(shù)在計算機(jī)中表示與整數(shù)在計算機(jī)中的表示具有完全不一樣的特殊性,浮點(diǎn)數(shù)的比較就完全不能簡單的用“==”來實現(xiàn)。通常會利用差值的絕對值的精度來判斷。

    假設(shè):f1和f2是兩個浮點(diǎn)數(shù),在C語言中規(guī)定有雙精度精度DBL_EPSILON。

    #defineDBL_EPSILON

    2.2204460492503131e-016

    /*smallest such that

    1.0+DBL_EPSILON!=1.0*/

    即可以用 fabs(f1-f2)<= DBL_EPSILON 來判斷f1和f2是否相等。如果要求更高的精度,也不能簡單地把DBL_EPSILON定得更小就行了。

    因為,DBL_EPSILON是C語言編譯器給出的一個不變的數(shù)據(jù),也就是誤差分析當(dāng)中的絕對誤差,使用一個完全固定的數(shù)值,對于雙精度類型的數(shù)據(jù)可以表達(dá)的整個數(shù)域來說是有缺陷的。例如對于f1和f2大小分別是1e-100 、2e-100附近這樣的數(shù)據(jù)的時候,它是不合適的,因為它們之差1e-100已經(jīng)遠(yuǎn)小于DBL_EPSILON,兩個不等的浮點(diǎn)數(shù)就可以判斷是相等的,出現(xiàn)了相等比較的錯誤。即使DBL_EPSILON已經(jīng)是最小的,而且它是由數(shù)據(jù)在計算機(jī)中表示時的位數(shù)決定的。例如,對于f1和f2大小是10和10.0000000000000000001這樣的數(shù)據(jù)的時候,它也不合適,因為10和10.0000000000000000001兩個數(shù)的絕對誤差為1e-19,小于DBL_EPSILON,判斷的結(jié)果當(dāng)然是相等,也就出錯了,因為兩個數(shù)據(jù)的有效數(shù)字位數(shù)超過了16位。適合浮點(diǎn)數(shù)相等比較的情況只是f1或者f2在數(shù)值1附近的時候,函數(shù)的注釋正好說明了這個問題,大范圍的浮點(diǎn)數(shù)的比較需要解決。

    1 雙精度浮點(diǎn)數(shù)比較相等方法的改進(jìn)

    以上用絕對誤差的方式比較浮點(diǎn)數(shù)的相等不合適,就用相對誤差方式嘗試解決浮點(diǎn)數(shù)相等比較的問題。相對誤差的算法如下:

    bool IsEqual(double a, double b, double relError )

    {

    return (( fabs ( (a-b)/a ) < relError )?TRUE :FALSE;)

    }

    這個浮點(diǎn)數(shù)比較相等的函數(shù)也是有問題的。用固定的函數(shù)的第一個形參做相對比較,在應(yīng)用過程中,調(diào)用IsEqual(a, b, relError ) 和 IsEqual(b, a, relError ) 的時候,由于參數(shù)a和b的數(shù)量級不一樣,得到的結(jié)果是不同的。顯然,當(dāng)實際調(diào)用的第一個參數(shù)是0的話,就有出現(xiàn)了除0溢出的錯誤,這樣的相對誤差比較顯然也是不合適的。

    顯而易見,可以把除數(shù)選取為a和b當(dāng)中絕對數(shù)值較大的就可以避免以上情況[3]:

    bool IsEqual(double a, double b, relError )

    {

    if (fabs(a)

    return ( fabs((a-b)/a) > relError ) ?TRUE : FALSE;

    else

    return (fabs( (a-b)/b) > relError ) ?TRUE : FALSE;

    };

    此算法就克服了絕對誤差算法對于太小的浮點(diǎn)數(shù)不能比較的缺陷。

    2 優(yōu)化浮點(diǎn)數(shù)比較大小的實現(xiàn)方法

    相對誤差算法能處理很小的浮點(diǎn)數(shù)相等的比較,但浮點(diǎn)數(shù)的運(yùn)算速度慢,能改用整數(shù)運(yùn)算就快了。為了從根本上解決問題,需要仔細(xì)研究浮點(diǎn)數(shù)在計算機(jī)中的表示。根據(jù)IEEE754浮點(diǎn)數(shù)的內(nèi)存結(jié)構(gòu),符號1位在最高位,指數(shù)11位在次高位,尾數(shù)52位在低位(不包括隱含的1位)。浮點(diǎn)數(shù)的大小對應(yīng)的是其指數(shù)和尾數(shù)的大小,人工比較時,先比較符號,再比較指數(shù),最后比較尾數(shù),就可知兩個浮點(diǎn)數(shù)是否相等。用C語言編程來比較浮點(diǎn)數(shù)的相等也應(yīng)該是同樣的道理。即根據(jù)浮點(diǎn)數(shù)的內(nèi)存結(jié)構(gòu),按照整數(shù)來理解進(jìn)行相等的比較,情況也是相同的。而且用這種方法對浮點(diǎn)數(shù)進(jìn)行比較的話,因為作為整數(shù)運(yùn)算效率比浮點(diǎn)數(shù)的高得多。比如

    double f1 = 1.58;

    double f2 = 1.57

    根據(jù)以上的定義,事實上f1>f2 是成立,能否用一種方法,證明(long long int)f1 > (long long int)f2 也是成立的。根據(jù)IEEE754的浮點(diǎn)結(jié)構(gòu)特點(diǎn),不是所有的浮點(diǎn)數(shù)都可以精確的表達(dá),可以精確表達(dá)的浮點(diǎn)數(shù)實際上是有限的,就是IEEE754枚舉的232個,事實上,絕大多數(shù)的浮點(diǎn)數(shù)數(shù)值是不能在計算機(jī)里精確表達(dá)的。

    目前雖然有80位的浮點(diǎn)數(shù)表示,這里只分析64位的情況(IEEE754)。尾數(shù)是53位的(暗含了第一個位數(shù)是1),對于可以精確表達(dá)的浮點(diǎn)數(shù)來說,對于IEEE754 單精度和雙精度浮點(diǎn)數(shù),能夠精確表示的整數(shù)的范圍為[4]:

    如果把這53位當(dāng)作整數(shù)來理解,對于隱含的1位,除浮點(diǎn)數(shù)零外,其它數(shù)都是1,不考慮也不影響比較的結(jié)果。對于任一浮點(diǎn)數(shù),把它當(dāng)作64位的整數(shù),根據(jù)二進(jìn)制計數(shù)原理,而且是線性分布的。這樣,將兩個浮點(diǎn)數(shù)看成是整數(shù),不需要轉(zhuǎn)換成整數(shù),將其對應(yīng)的整數(shù)做差值運(yùn)算,得到的整數(shù)表明的是兩個浮點(diǎn)數(shù)之間有多少個實際可以精確表達(dá)的浮點(diǎn)數(shù)的個數(shù)(對應(yīng)的指數(shù)相同的情況下是不言而喻的;指數(shù)不同的時候,也是同樣有效的,因為現(xiàn)在研究的是相等比較,不是大小比較,指數(shù)不同則涉及的是精確表達(dá)的浮點(diǎn)數(shù)的個數(shù)更多了,有利于相等的比較,而不會有相反的作用)。因此,對于兩個正的浮點(diǎn)數(shù),他們的大小比較就可以用 (long long int)f1 - (long long int)f2 來進(jìn)行比較了,即將兩個浮點(diǎn)數(shù)當(dāng)成整型數(shù)來相減,其差值的結(jié)果實際上就相當(dāng)于相對誤差了,這個相對誤差,不等同于普通意義上的相對誤差,它所表達(dá)的是,兩個浮點(diǎn)數(shù)之間可能還有多少個可以精確表達(dá)的浮點(diǎn)數(shù)。這樣通過指定這個閾值來控制兩個浮點(diǎn)數(shù)的比較就更有效了[5]。(long int)f1 - (long int)f2在這里只是算法的示意,C編譯器會報錯。具體實現(xiàn)見最后面的程序。

    對于兩個正的浮點(diǎn)數(shù)比較相等就有:

    bool IsEqual(double f1, double f2, int precisionFloatNum)

    {

    if ( abs ( (long long int)f1 - (long longint)f2 ) < precisionFloatNum ) return TRUE;

    }

    這里要用整數(shù)的絕對值比較函數(shù)而不能用浮點(diǎn)數(shù)的絕對值比較函數(shù),因為要將此時的浮點(diǎn)數(shù)看成是整數(shù),否則就無法實現(xiàn)以上分析的目的。但是負(fù)整數(shù)和正整數(shù)之間現(xiàn)在還不能進(jìn)行直接的比較,因為根據(jù)IEEE754的內(nèi)存結(jié)構(gòu),正整數(shù)和負(fù)整數(shù)是不同的,對應(yīng)的在計算機(jī)中表示的整數(shù)在整個區(qū)間是不連續(xù)的,正的最小的整數(shù)就是0,對應(yīng)的計算機(jī)表示的整數(shù)是0x0000000000000000,負(fù)的最大的整數(shù)就是-0,與之對應(yīng)的計算機(jī)表示的整數(shù)則是0x 8000000000000000,在IEEE754的表達(dá)當(dāng)中是有兩個0的,一個是 +0 一個是-0。而且按照 f1 == f2 的判斷 +0和-0是相等的,根據(jù)浮點(diǎn)數(shù)對應(yīng)的在計算機(jī)中表示的整數(shù)在整個區(qū)間的不連續(xù)性,具有以下特點(diǎn),+0 和正的計算機(jī)表示的浮點(diǎn)數(shù)可以當(dāng)成整數(shù)的方式直接進(jìn)行比較,-0 和負(fù)的計算機(jī)表示的浮點(diǎn)數(shù)可以當(dāng)成為整數(shù)的方式直接進(jìn)行比較,如果采取措施將把他們統(tǒng)一起來,計算機(jī)表示的浮點(diǎn)數(shù)可以當(dāng)成整數(shù)直接進(jìn)行整數(shù)比較方式算法就完備了。仔細(xì)分析負(fù)整數(shù)的結(jié)構(gòu),把負(fù)整數(shù)對應(yīng)的整數(shù)減去 -0 ,這樣,浮點(diǎn)數(shù)對應(yīng)的在計算機(jī)中表示的整數(shù)在整個區(qū)間就連續(xù)了。因此,所有的計算機(jī)表示用補(bǔ)碼的負(fù)整數(shù)經(jīng)過這次減法后,對應(yīng)的整數(shù)也都是原碼的負(fù)整數(shù)了,這樣整個整數(shù)比較就變得連續(xù)了就相當(dāng)于在整個浮點(diǎn)數(shù)范圍內(nèi)都是有效的了,只有被比較的浮點(diǎn)數(shù)連續(xù)或線性分布,將浮點(diǎn)數(shù)看成整數(shù)的比較算法才有意義。算法的關(guān)鍵是負(fù)數(shù)只需將最高位的1去掉,方法見程序的實現(xiàn)。注意,這里只是比較相等,不是比較大小,將被比較的兩個64位浮點(diǎn)數(shù)之間允許有多少個可以精確表達(dá)的浮點(diǎn)數(shù)作為比較相等的關(guān)鍵,從而實現(xiàn)快速和高精度的算法。

    最后的浮點(diǎn)數(shù)比較算法就是:

    /* 函數(shù): bool IsEqual(double f1, double f2, int precisionFloatNum) */

    /* 功能:兩個64位浮點(diǎn)數(shù)是否近似相等*/

    /* 輸入:兩個64位浮點(diǎn)數(shù)f1, f2 */

    /* precisionFloatNum 被比較的兩個64位浮點(diǎn)數(shù)之間允許有多少個可以精確表達(dá)的浮點(diǎn)數(shù) */

    /* 輸出: TRUE,兩個浮點(diǎn)數(shù)64位近似相等;FALSE 兩個64位浮點(diǎn)數(shù)不等 */

    bool IsEqual(double f1, double f2, long int precisionFloatNum)

    {

    void *p1=&f1;

    void *p2=&f2;

    long long int bits1=*((long long int*)p1);

    long long int bits2=*((long long int*)p2);

    if((bits1>0&&bits2>0)||(bits1<0&&bits2<0))

    { bits1=(bits1>0)?bits1:(bits1-0x8000000000000000);

    bits2=(bits2>0)?bits2:(bits2-0x8000000000000000);

    return ((abs(bits1-bits2))< precisionFloatNum) ? TRUE : FALSE;

    }

    else

    return FALSE;

    }

    通過編程實現(xiàn)浮點(diǎn)數(shù)相等的比較,相對誤差比較的算法和將浮點(diǎn)數(shù)當(dāng)在整數(shù)即線性化的比較的算法都能實現(xiàn)16位有效數(shù)字的精度,可實現(xiàn)最小到DBL_MIN數(shù)的比較。DBL_MIN的定義是:

    #defineDBL_MIN

    2.2250738585072014e-308

    /*min positive value*/

    在浮點(diǎn)數(shù)的絕對值大于DBL_MIN的基礎(chǔ)上,所有浮點(diǎn)數(shù)相等的比較現(xiàn)在只取決于有效數(shù)的位數(shù),由此完整實現(xiàn)了浮點(diǎn)數(shù)的相等比較。

    3 結(jié)論

    本文從浮點(diǎn)數(shù)在計算機(jī)內(nèi)部表示的角度出發(fā),介紹了雙精度浮點(diǎn)數(shù)比較大小的優(yōu)化算法,通過程序調(diào)試驗證了算法的有效性,并提出了改進(jìn)方法。對于兩個數(shù)據(jù)的有效數(shù)字位數(shù)超過了16位(253)的浮點(diǎn)數(shù)的相等比較,只有提高數(shù)據(jù)表示的位數(shù),才能實現(xiàn),例如,長雙精度浮點(diǎn)數(shù)有80位的數(shù)據(jù)寬度。從運(yùn)算速度上看,將浮點(diǎn)數(shù)當(dāng)成整數(shù)的浮點(diǎn)數(shù)線性化的算法在運(yùn)算速度和精度上優(yōu)于絕對和相對精度的算法,因為整數(shù)的減法和取絕對值運(yùn)算速度比浮點(diǎn)數(shù)的運(yùn)算要快。這里用LCC編譯器實現(xiàn)的雙精度的浮點(diǎn)數(shù)比較相等,其它編譯環(huán)境的算法則需要一定的修改。至于單精度的浮點(diǎn)數(shù)和長雙精度的浮點(diǎn)數(shù),修改程序的相關(guān)部分即可實現(xiàn)。

    [1] 杜叔強(qiáng). 淺析C語言中的浮點(diǎn)數(shù)[J]. 蘭州工業(yè)高等專科學(xué)校學(xué)報, 2010,17(5):27.

    [2] David Goldberg. What every computer scientist should know about doubling-point Arithmetic [M]. Computing Surveys, 1991.

    [3] 張宗杰,張明亮. C 語言中浮點(diǎn)數(shù)的存儲格式及其有效數(shù)字位數(shù)[J].計算機(jī)與數(shù)字工程, 2006,34(1):85.

    [4] http://blog.csdn.net/seizef/article/details/5571783.

    [5] http://cmdblock.blog.51cto.com/415170/600378.

    [6] 陳煒峰. 電磁脈沖模擬器及其應(yīng)用研究[J]. 南京: 東南大學(xué)博士學(xué)位論文, 2007.

    Approach to Compare to Equation of Double-precision Floating Point Numbers of Linearization in C Language

    Zhang Guohua

    (Wuhan Institute of Shipbuilding Technology, Wuhan430050, China)

    Based on LCC software for design, this paper puts forward the design method of the double precision doubling- point comparisons algorithm, the accuracy of the algorithm is verified in applications, and the algorithm is designed with the absolute error and relative error method respectively. Finally aimed at representation for double-precision doubling-point number, the paper presents the arithmetic for floating point linearization.

    compare; algorithms; double precision; significance digit

    TP 312.1.4

    A

    1003-4862(2017)01-0040-03

    2016-11-15

    章國華(1964-),男,副教授。研究方向:機(jī)電一體化技術(shù)教學(xué)與研究。Email:993468391@qq.com

    猜你喜歡
    浮點(diǎn)數(shù)線性化位數(shù)
    四種Python均勻浮點(diǎn)數(shù)生成方法
    電腦報(2021年11期)2021-07-01 08:26:31
    五次完全冪的少位數(shù)三進(jìn)制展開
    “線性化”在多元不等式證明與最值求解中的應(yīng)用
    基于反饋線性化的RLV氣動控制一體化設(shè)計
    非精確浮點(diǎn)數(shù)乘法器設(shè)計
    EHA反饋線性化最優(yōu)滑模面雙模糊滑??刂?/a>
    空間機(jī)械臂鎖緊機(jī)構(gòu)等效線性化分析及驗證
    遙感衛(wèi)星CCD相機(jī)量化位數(shù)的選擇
    “判斷整數(shù)的位數(shù)”的算法分析
    河南科技(2014年11期)2014-02-27 14:09:41
    基于分位數(shù)回歸的剪切波速變化規(guī)律
    亚洲va日本ⅴa欧美va伊人久久| 国产精品影院久久| 日韩欧美在线二视频| 午夜免费激情av| 精品久久久久久,| 三级国产精品欧美在线观看 | 麻豆av在线久日| 黄频高清免费视频| 宅男免费午夜| 男女下面进入的视频免费午夜| 午夜福利在线在线| 三级国产精品欧美在线观看 | 天天一区二区日本电影三级| 一区二区三区国产精品乱码| 亚洲欧美日韩高清在线视频| 久久久水蜜桃国产精品网| 国产熟女xx| 一级片免费观看大全| 欧美在线黄色| 亚洲欧美日韩高清在线视频| 国产成人精品无人区| av欧美777| 欧美+亚洲+日韩+国产| 欧美日韩黄片免| 在线观看www视频免费| 久久国产精品影院| 欧美 亚洲 国产 日韩一| 女人爽到高潮嗷嗷叫在线视频| 99精品在免费线老司机午夜| 成人一区二区视频在线观看| 成人一区二区视频在线观看| 欧美最黄视频在线播放免费| 日韩大尺度精品在线看网址| av片东京热男人的天堂| 午夜福利欧美成人| 午夜日韩欧美国产| 免费看日本二区| 国产精品av久久久久免费| 亚洲国产精品999在线| 成年版毛片免费区| 久久热在线av| 在线观看免费午夜福利视频| 亚洲av第一区精品v没综合| 少妇被粗大的猛进出69影院| 欧美日本亚洲视频在线播放| 不卡一级毛片| 亚洲欧美日韩东京热| 亚洲自偷自拍图片 自拍| 午夜日韩欧美国产| 午夜日韩欧美国产| 国产伦一二天堂av在线观看| 午夜日韩欧美国产| 精品久久久久久久人妻蜜臀av| 一级毛片精品| 国产成人啪精品午夜网站| 国内精品久久久久久久电影| 欧美色视频一区免费| 免费看十八禁软件| 国产精品久久久久久人妻精品电影| 十八禁网站免费在线| 午夜免费成人在线视频| 午夜福利高清视频| 免费在线观看影片大全网站| 妹子高潮喷水视频| 婷婷精品国产亚洲av| 欧美一级a爱片免费观看看 | 黑人操中国人逼视频| 精品国产乱子伦一区二区三区| 亚洲天堂国产精品一区在线| 亚洲aⅴ乱码一区二区在线播放 | 中出人妻视频一区二区| 欧美成人性av电影在线观看| 色噜噜av男人的天堂激情| 欧美精品啪啪一区二区三区| 亚洲国产欧美网| 可以在线观看毛片的网站| 精品国内亚洲2022精品成人| 欧美久久黑人一区二区| 精品一区二区三区av网在线观看| 国产单亲对白刺激| 亚洲精品粉嫩美女一区| 亚洲成人中文字幕在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 国产三级中文精品| 亚洲精品中文字幕在线视频| 久久久久久久久中文| 两个人看的免费小视频| 日韩欧美国产在线观看| 中文字幕最新亚洲高清| 一本精品99久久精品77| 两个人看的免费小视频| 免费观看精品视频网站| 九色国产91popny在线| 香蕉国产在线看| 日日干狠狠操夜夜爽| 久久久久久国产a免费观看| 久久香蕉国产精品| 日本熟妇午夜| 啦啦啦观看免费观看视频高清| 很黄的视频免费| 亚洲国产欧洲综合997久久,| 在线观看66精品国产| 欧美精品啪啪一区二区三区| 亚洲av电影不卡..在线观看| 亚洲欧美日韩高清专用| 欧美成人一区二区免费高清观看 | 免费高清视频大片| 视频区欧美日本亚洲| 我要搜黄色片| 色播亚洲综合网| 亚洲精华国产精华精| 精品欧美国产一区二区三| 久久久国产成人免费| 国产精品 国内视频| av在线天堂中文字幕| 亚洲av成人不卡在线观看播放网| 全区人妻精品视频| 欧美日韩瑟瑟在线播放| 50天的宝宝边吃奶边哭怎么回事| 国产日本99.免费观看| 嫁个100分男人电影在线观看| 天天添夜夜摸| 美女高潮喷水抽搐中文字幕| 日韩欧美精品v在线| 麻豆久久精品国产亚洲av| 岛国在线观看网站| 精品高清国产在线一区| 免费在线观看日本一区| 成人亚洲精品av一区二区| 精品无人区乱码1区二区| a在线观看视频网站| 欧美成人性av电影在线观看| 国产欧美日韩精品亚洲av| 国产高清激情床上av| 男女下面进入的视频免费午夜| 久久精品aⅴ一区二区三区四区| 村上凉子中文字幕在线| 成人三级黄色视频| 国内揄拍国产精品人妻在线| АⅤ资源中文在线天堂| or卡值多少钱| 51午夜福利影视在线观看| av超薄肉色丝袜交足视频| 久久精品成人免费网站| 免费观看人在逋| 中文字幕人妻丝袜一区二区| 色综合亚洲欧美另类图片| 国产精品久久久久久精品电影| 国产精品一区二区三区四区久久| 免费高清视频大片| 欧美一区二区国产精品久久精品 | 在线免费观看的www视频| 九色成人免费人妻av| 露出奶头的视频| 大型av网站在线播放| 女警被强在线播放| 亚洲免费av在线视频| 成人亚洲精品av一区二区| xxx96com| 欧美高清成人免费视频www| 无限看片的www在线观看| 黄色成人免费大全| 久久九九热精品免费| 亚洲精品美女久久久久99蜜臀| 99热只有精品国产| 久久精品人妻少妇| 欧美性长视频在线观看| 国产熟女xx| √禁漫天堂资源中文www| 无人区码免费观看不卡| 久9热在线精品视频| 99riav亚洲国产免费| 丝袜人妻中文字幕| 女警被强在线播放| 校园春色视频在线观看| 午夜精品一区二区三区免费看| 一本综合久久免费| 黄色成人免费大全| 身体一侧抽搐| 国产一区二区三区视频了| 黄色 视频免费看| 国产高清视频在线播放一区| 黄色 视频免费看| 99久久久亚洲精品蜜臀av| 91九色精品人成在线观看| 好看av亚洲va欧美ⅴa在| 欧美性长视频在线观看| 一级毛片女人18水好多| 精品久久久久久成人av| 国产成人欧美在线观看| 亚洲成人久久爱视频| 久久久国产成人免费| 久久精品91无色码中文字幕| 日韩av在线大香蕉| 免费人成视频x8x8入口观看| 亚洲七黄色美女视频| 日本免费a在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲中文av在线| 午夜免费观看网址| 长腿黑丝高跟| 国产私拍福利视频在线观看| 麻豆国产97在线/欧美 | 禁无遮挡网站| 精品久久蜜臀av无| 99久久综合精品五月天人人| 精品久久久久久久人妻蜜臀av| 国产精品乱码一区二三区的特点| 亚洲成人久久性| 色噜噜av男人的天堂激情| 男女做爰动态图高潮gif福利片| 欧美av亚洲av综合av国产av| 亚洲国产高清在线一区二区三| 国产精品1区2区在线观看.| 动漫黄色视频在线观看| 久久久久亚洲av毛片大全| avwww免费| 一级a爱片免费观看的视频| 亚洲18禁久久av| а√天堂www在线а√下载| 无人区码免费观看不卡| 欧美中文综合在线视频| 欧美高清成人免费视频www| 三级男女做爰猛烈吃奶摸视频| 亚洲av成人不卡在线观看播放网| 日韩精品免费视频一区二区三区| 亚洲中文av在线| 国产精品 国内视频| 在线观看日韩欧美| 国产高清videossex| 黑人操中国人逼视频| 老汉色∧v一级毛片| 757午夜福利合集在线观看| 国产午夜精品论理片| 国产伦在线观看视频一区| 亚洲专区国产一区二区| av福利片在线观看| 亚洲专区字幕在线| 在线十欧美十亚洲十日本专区| 日韩欧美 国产精品| 变态另类丝袜制服| 成人av在线播放网站| 啦啦啦观看免费观看视频高清| 精品国产超薄肉色丝袜足j| 免费在线观看完整版高清| 在线观看午夜福利视频| 特大巨黑吊av在线直播| 亚洲专区字幕在线| 999久久久国产精品视频| 巨乳人妻的诱惑在线观看| 国产免费av片在线观看野外av| 宅男免费午夜| 欧美性猛交黑人性爽| 日本a在线网址| 高清毛片免费观看视频网站| 18禁黄网站禁片午夜丰满| 国产69精品久久久久777片 | 国产黄色小视频在线观看| 亚洲人与动物交配视频| 亚洲欧洲精品一区二区精品久久久| 特级一级黄色大片| 亚洲成人久久性| 久久精品影院6| 一级作爱视频免费观看| 人妻丰满熟妇av一区二区三区| 精品福利观看| 久久国产精品影院| 黄片小视频在线播放| 欧美zozozo另类| 精品日产1卡2卡| 国产亚洲精品av在线| 男人舔女人的私密视频| 亚洲18禁久久av| 男插女下体视频免费在线播放| 精品欧美一区二区三区在线| 好看av亚洲va欧美ⅴa在| 神马国产精品三级电影在线观看 | 国产片内射在线| 欧美日韩国产亚洲二区| 久久久久久免费高清国产稀缺| 欧美黑人精品巨大| 亚洲狠狠婷婷综合久久图片| 亚洲人成伊人成综合网2020| 波多野结衣巨乳人妻| 成人av一区二区三区在线看| 日韩欧美一区二区三区在线观看| av视频在线观看入口| 国产单亲对白刺激| 一本综合久久免费| 色哟哟哟哟哟哟| 免费在线观看完整版高清| 日本免费a在线| 一本久久中文字幕| 夜夜夜夜夜久久久久| 国产精品野战在线观看| 老熟妇乱子伦视频在线观看| 此物有八面人人有两片| 亚洲18禁久久av| 母亲3免费完整高清在线观看| 一本精品99久久精品77| 欧美日韩精品网址| 变态另类成人亚洲欧美熟女| 久久久久久人人人人人| 男女床上黄色一级片免费看| 18禁裸乳无遮挡免费网站照片| 黄片大片在线免费观看| 成人亚洲精品av一区二区| 又黄又爽又免费观看的视频| 久久亚洲精品不卡| 久久天躁狠狠躁夜夜2o2o| 最新在线观看一区二区三区| 一区二区三区国产精品乱码| 黄色 视频免费看| 亚洲熟女毛片儿| 精品一区二区三区视频在线观看免费| 国产男靠女视频免费网站| 国产精品一及| 欧美日韩中文字幕国产精品一区二区三区| 99国产精品一区二区三区| 日本黄色视频三级网站网址| 村上凉子中文字幕在线| 又黄又爽又免费观看的视频| 国产乱人伦免费视频| 久久久久久九九精品二区国产 | 欧美人与性动交α欧美精品济南到| 国产av一区在线观看免费| 欧美乱色亚洲激情| 国产亚洲av高清不卡| 国产精品久久久av美女十八| а√天堂www在线а√下载| 亚洲第一欧美日韩一区二区三区| 老熟妇乱子伦视频在线观看| 精品国产美女av久久久久小说| av超薄肉色丝袜交足视频| 久久久久精品国产欧美久久久| 国语自产精品视频在线第100页| 亚洲片人在线观看| 久久天躁狠狠躁夜夜2o2o| 国产片内射在线| 国产成人精品久久二区二区免费| 麻豆成人av在线观看| 在线视频色国产色| 日日干狠狠操夜夜爽| 99热这里只有精品一区 | 国产视频内射| 中文字幕av在线有码专区| 精品国产超薄肉色丝袜足j| 亚洲片人在线观看| 久久精品综合一区二区三区| 国产成人一区二区三区免费视频网站| 丁香六月欧美| www.自偷自拍.com| 久久天堂一区二区三区四区| 一个人免费在线观看电影 | 国模一区二区三区四区视频 | 久久久久国内视频| 国产精品 欧美亚洲| 在线a可以看的网站| 免费在线观看影片大全网站| 在线永久观看黄色视频| 人成视频在线观看免费观看| 人妻久久中文字幕网| 老司机在亚洲福利影院| 久久精品国产清高在天天线| 国产精品久久视频播放| 国产精华一区二区三区| 中出人妻视频一区二区| 欧美丝袜亚洲另类 | 男女那种视频在线观看| 午夜激情av网站| 精品人妻1区二区| 琪琪午夜伦伦电影理论片6080| 久久九九热精品免费| 久久久久精品国产欧美久久久| 少妇被粗大的猛进出69影院| 99热这里只有精品一区 | 亚洲专区中文字幕在线| 在线视频色国产色| 天天添夜夜摸| 国产麻豆成人av免费视频| 午夜福利成人在线免费观看| netflix在线观看网站| 成年版毛片免费区| 男人舔女人的私密视频| 一二三四社区在线视频社区8| www.www免费av| 国产精品99久久99久久久不卡| 大型黄色视频在线免费观看| 中文字幕熟女人妻在线| aaaaa片日本免费| 91九色精品人成在线观看| 婷婷精品国产亚洲av在线| 天堂√8在线中文| bbb黄色大片| 亚洲欧美日韩东京热| 国产精品一及| 可以在线观看的亚洲视频| 特级一级黄色大片| 香蕉av资源在线| 久久久国产成人免费| 超碰成人久久| 午夜免费观看网址| 精品欧美国产一区二区三| 亚洲av成人一区二区三| 久久久久久人人人人人| 88av欧美| 午夜福利在线观看吧| 国产在线精品亚洲第一网站| 男人舔女人的私密视频| 成人午夜高清在线视频| 国产aⅴ精品一区二区三区波| 18禁观看日本| 夜夜爽天天搞| 精品国产乱子伦一区二区三区| 日韩中文字幕欧美一区二区| 久久精品aⅴ一区二区三区四区| 久久香蕉激情| 亚洲国产日韩欧美精品在线观看 | 69av精品久久久久久| 99久久综合精品五月天人人| 国产黄片美女视频| 欧美黄色片欧美黄色片| 男人的好看免费观看在线视频 | 最好的美女福利视频网| 18禁美女被吸乳视频| xxxwww97欧美| 国内精品久久久久精免费| 在线a可以看的网站| 国产午夜精品久久久久久| 国产v大片淫在线免费观看| 亚洲中文字幕日韩| 亚洲欧美一区二区三区黑人| 黑人欧美特级aaaaaa片| 亚洲一区二区三区色噜噜| 久久久久久久久久黄片| 亚洲五月婷婷丁香| 脱女人内裤的视频| 男人舔奶头视频| 精品国产乱码久久久久久男人| 久久精品国产99精品国产亚洲性色| 日本熟妇午夜| 黄色女人牲交| 黄片大片在线免费观看| 看免费av毛片| 欧美又色又爽又黄视频| 香蕉丝袜av| 国产精品久久电影中文字幕| 欧美激情久久久久久爽电影| 少妇被粗大的猛进出69影院| 久久九九热精品免费| 国产激情偷乱视频一区二区| 女人高潮潮喷娇喘18禁视频| 在线观看美女被高潮喷水网站 | 国产精品一区二区三区四区久久| 在线观看日韩欧美| 搡老熟女国产l中国老女人| www日本黄色视频网| 中亚洲国语对白在线视频| 日本在线视频免费播放| 不卡一级毛片| 日本免费a在线| 一区二区三区国产精品乱码| 亚洲中文字幕日韩| 香蕉丝袜av| 欧美国产日韩亚洲一区| 国产乱人伦免费视频| 国产黄a三级三级三级人| 亚洲精品国产一区二区精华液| www.自偷自拍.com| 亚洲国产精品成人综合色| 亚洲国产欧美网| 久久这里只有精品中国| 欧美日韩福利视频一区二区| 国产成+人综合+亚洲专区| 午夜a级毛片| 精品熟女少妇八av免费久了| 午夜激情av网站| 欧美性猛交黑人性爽| 国产高清videossex| 18禁观看日本| 免费人成视频x8x8入口观看| xxx96com| 亚洲美女视频黄频| 欧美丝袜亚洲另类 | 日韩欧美免费精品| 日韩欧美国产在线观看| 国产爱豆传媒在线观看 | 国产精品电影一区二区三区| 久久久国产欧美日韩av| 在线观看午夜福利视频| 亚洲中文日韩欧美视频| 久久精品影院6| 国内精品久久久久精免费| 久久婷婷人人爽人人干人人爱| 免费在线观看成人毛片| 日韩精品中文字幕看吧| 男女床上黄色一级片免费看| 桃色一区二区三区在线观看| 又紧又爽又黄一区二区| 欧美性猛交黑人性爽| 蜜桃久久精品国产亚洲av| 一区二区三区高清视频在线| 亚洲七黄色美女视频| 亚洲在线自拍视频| 女人高潮潮喷娇喘18禁视频| 亚洲片人在线观看| 欧美人与性动交α欧美精品济南到| 亚洲专区字幕在线| 88av欧美| 91麻豆av在线| 少妇粗大呻吟视频| cao死你这个sao货| 国产精品野战在线观看| 久久精品aⅴ一区二区三区四区| 日韩有码中文字幕| 亚洲午夜精品一区,二区,三区| 九色成人免费人妻av| 美女扒开内裤让男人捅视频| 悠悠久久av| 久久亚洲真实| 欧美高清成人免费视频www| 岛国视频午夜一区免费看| 欧美中文日本在线观看视频| 亚洲一区二区三区色噜噜| 国产高清有码在线观看视频 | 激情在线观看视频在线高清| 亚洲国产欧美人成| 黄色丝袜av网址大全| 2021天堂中文幕一二区在线观| 男插女下体视频免费在线播放| 久久精品亚洲精品国产色婷小说| 美女高潮喷水抽搐中文字幕| 久久草成人影院| 久久婷婷成人综合色麻豆| 国产成人影院久久av| 欧美不卡视频在线免费观看 | 观看免费一级毛片| 国产不卡一卡二| 久久久久性生活片| 欧美+亚洲+日韩+国产| 真人一进一出gif抽搐免费| 亚洲欧美精品综合一区二区三区| 欧美日韩亚洲综合一区二区三区_| 日日夜夜操网爽| 宅男免费午夜| 国产精品 欧美亚洲| 在线观看一区二区三区| 国产成人欧美在线观看| 人人妻,人人澡人人爽秒播| 一级作爱视频免费观看| 精品国产亚洲在线| cao死你这个sao货| 老司机深夜福利视频在线观看| 90打野战视频偷拍视频| 国产又黄又爽又无遮挡在线| 少妇人妻一区二区三区视频| bbb黄色大片| 一进一出抽搐gif免费好疼| 高清毛片免费观看视频网站| 国产精品亚洲一级av第二区| 亚洲黑人精品在线| 久久久久免费精品人妻一区二区| 天堂动漫精品| 91大片在线观看| 久久精品国产亚洲av香蕉五月| 国产成人啪精品午夜网站| 一a级毛片在线观看| 夜夜夜夜夜久久久久| 日本免费一区二区三区高清不卡| 在线观看免费日韩欧美大片| 国产精品久久久久久久电影 | 成年版毛片免费区| 免费看美女性在线毛片视频| 精品国内亚洲2022精品成人| 长腿黑丝高跟| 久久人妻av系列| 国产精品久久久久久精品电影| 成熟少妇高潮喷水视频| 亚洲乱码一区二区免费版| 热99re8久久精品国产| 久久久久久九九精品二区国产 | 国产亚洲av高清不卡| 午夜福利在线观看吧| 在线观看66精品国产| 日韩大码丰满熟妇| 亚洲成人免费电影在线观看| 黄片大片在线免费观看| 国产视频一区二区在线看| 亚洲黑人精品在线| 亚洲最大成人中文| 日韩欧美免费精品| 国产熟女xx| 久久久久久久久久黄片| 亚洲国产中文字幕在线视频| 国产在线观看jvid| 又爽又黄无遮挡网站| 国产三级在线视频| 非洲黑人性xxxx精品又粗又长| 老熟妇仑乱视频hdxx| 真人一进一出gif抽搐免费| 欧美三级亚洲精品| 桃红色精品国产亚洲av| 国产视频一区二区在线看| 午夜福利在线观看吧| av福利片在线观看| 亚洲国产欧美一区二区综合| 国产精品国产高清国产av| 岛国在线免费视频观看| 成人高潮视频无遮挡免费网站| 性色av乱码一区二区三区2| 一本大道久久a久久精品| 在线观看美女被高潮喷水网站 | 国模一区二区三区四区视频 | 久久精品亚洲精品国产色婷小说| 天天添夜夜摸|