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

    廣義模糊K調(diào)和均值聚類的近紅外光譜生菜儲藏時間鑒別

    2016-07-12 12:59:05武小紅潘明輝
    光譜學(xué)與光譜分析 2016年6期
    關(guān)鍵詞:調(diào)和生菜儲藏

    武小紅,潘明輝,武 斌,嵇 港,孫 俊

    1. 江蘇大學(xué)電氣信息工程學(xué)院, 江蘇 鎮(zhèn)江 212013 2. 江蘇大學(xué)機械工業(yè)設(shè)施農(nóng)業(yè)測控技術(shù)與裝備重點實驗室, 江蘇 鎮(zhèn)江 212013 3. 江蘇大學(xué)京江學(xué)院, 江蘇 鎮(zhèn)江 212013 4. 滁州職業(yè)技術(shù)學(xué)院信息工程系,安徽 滁州 239000

    廣義模糊K調(diào)和均值聚類的近紅外光譜生菜儲藏時間鑒別

    武小紅1,2,潘明輝3,武 斌4,嵇 港1,孫 俊1,2

    1. 江蘇大學(xué)電氣信息工程學(xué)院, 江蘇 鎮(zhèn)江 212013 2. 江蘇大學(xué)機械工業(yè)設(shè)施農(nóng)業(yè)測控技術(shù)與裝備重點實驗室, 江蘇 鎮(zhèn)江 212013 3. 江蘇大學(xué)京江學(xué)院, 江蘇 鎮(zhèn)江 212013 4. 滁州職業(yè)技術(shù)學(xué)院信息工程系,安徽 滁州 239000

    生菜的儲藏時間是影響生菜新鮮程度的重要因素。為了快速、無損和有效地鑒別生菜的儲藏時間,以歐式距離的p次方代替模糊K調(diào)和均值聚類(FKHM)中歐式距離的平方提出了一種廣義模糊K調(diào)和均值聚類(GFKHM)算法并將該算法應(yīng)用于鑒別生菜的儲藏時間。以60個新鮮生菜樣本為研究對象,采用Antaris Ⅱ近紅外光譜分析儀每隔12 h檢測生菜的近紅外漫反射光譜,共檢測三次,光譜掃描的波數(shù)范圍為10 000~4 000 cm-1。首先用主成分分析(PCA)對1 557維的生菜近紅外光譜進行降維處理以減少冗余信息,取前20個主成分,經(jīng)過PCA處理后得到20維的數(shù)據(jù)。然后用線性判別分析(LDA)提取光譜數(shù)據(jù)的鑒別信息以提高聚類的準確率,取鑒別向量數(shù)為2,則LDA將20維的數(shù)據(jù)轉(zhuǎn)換為2維數(shù)據(jù)。最后以模糊C-均值聚類(FCM)的類中心作為FKHM和GFKHM的初始聚類中心,分別運行FKHM和GFKHM計算模糊隸屬度以實現(xiàn)生菜儲藏時間的鑒別。結(jié)果表明,GFKHM的鑒別準確率能達到92.5%,F(xiàn)KHM的鑒別準確率為90.0%,GFKHM具有比FKHM更高的鑒別準確率。GFKHM的聚類中心比FKHM更逼近真實類中心。GFKHM的收斂速度明顯快于FKHM。采用近紅外光譜技術(shù)同時結(jié)合GFKHM,PCA和LDA為快速和無損地鑒別生菜儲藏時間提供了一種新的方法。

    近紅外光譜; 生菜; 儲藏時間; 線性判別分析; 模糊K調(diào)和均值聚類

    引 言

    生菜中含有豐富的營養(yǎng)物質(zhì)[1],它是人們經(jīng)常食用的蔬菜之一。隨著人們生活水平的不斷提高,對生菜的品質(zhì)要求越來越高。另一方面,隨著社會的進步,工作繁忙的人們將生菜儲藏在冰箱里留待日后食用,但是隨著儲藏時間的增加,生菜中的亞硝酸鹽和微生物會遞增破壞了生菜品質(zhì)而導(dǎo)致不能食用[2]。所以儲藏時間是影響生菜品質(zhì)的重要因素,研究一種簡單、快速、非破壞性的生菜貯藏時間的判定方法是非常必要的。

    近紅外光譜技術(shù)是根據(jù)樣品內(nèi)部的有機官能團(O—H,C—H,N—H,S—H)的自身振動吸收近紅外光譜區(qū)相應(yīng)波長的能量,從而在光譜中產(chǎn)生能量躍遷的表現(xiàn),是一種高效、快速的現(xiàn)代分析技術(shù)。其廣泛應(yīng)用于多個領(lǐng)域,特別是在農(nóng)產(chǎn)品和食品的無損檢測領(lǐng)域[3]。目前國內(nèi)外學(xué)者常用近紅外光譜技術(shù)建立農(nóng)產(chǎn)品和食品的定性和定量分析模型。定量分析是對樣品的組成成分及含量進行檢測,而樣品的品種鑒別[4],果實成熟鑒別[5],產(chǎn)地溯源[6], 儲藏時間鑒別[7]等均屬于定性分析。

    模糊K調(diào)和均值聚類[8]是將模糊的概念應(yīng)用到K調(diào)和均值聚類[9-10],保留了K調(diào)和均值聚類對初始值不敏感的優(yōu)點,提高了聚類準確度。模糊K調(diào)和均值聚類采用歐式距離作為度量方式,歐式距離適合簇形狀是超橢球[11]。在生菜檢測時獲取的近紅外光譜一般為高維數(shù)據(jù),經(jīng)過維數(shù)壓縮和特征提取后數(shù)據(jù)的簇拓撲結(jié)構(gòu)類似橢球形狀。以歐式距離的p次方代替模糊K調(diào)和均值聚類中歐式距離的平方,提出一種廣義模糊K調(diào)和均值聚類(GFKHM)算法以提高聚類的準確率。

    用近紅外光譜分析儀采集生菜不同儲藏時間的近紅外漫反射光譜,用主成分分析進行生菜近紅外光譜的壓縮,應(yīng)用線性判別分析提取的鑒別信息,并用廣義模糊K調(diào)和均值聚類(GFKHM)算法進行分類,同時與模糊K調(diào)和均值聚類(FKHM)算法進行對比,以便建立更準確的生菜儲藏時間鑒別模型。

    1 近紅外光譜采集

    采集新鮮生菜樣本60個,將生菜清理干凈后放入貼有標簽的保鮮袋中,再放入4 ℃冰箱里儲藏備用。每隔12 h取出生菜進行近紅外光譜檢測,共檢測三次,得到60×3=180個近紅外光譜數(shù)據(jù)。采集光譜時實驗室溫度保持在(15±2)℃,相對濕度在70%左右,Antaris Ⅱ近紅外光譜分析儀開機預(yù)熱1 h。采用反射積分球模式采集生菜近紅外光譜,近紅外光譜分析儀掃描每個樣品32次以獲取樣品的漫反射光譜均值。光譜掃描的波數(shù)范圍為10 000~4 000 cm-1,掃描間隔為3.856 cm-1,采集到每個樣品的近紅外光譜是1 557維的數(shù)據(jù)。為減少誤差,每個樣本采樣3次,取其平均值作為最終的實驗數(shù)據(jù)。生菜樣本的漫反射近紅外光譜圖如圖1所示。

    圖1 生菜樣本近紅外光譜

    2 廣義模糊K調(diào)和均值聚類

    給定一個無標記的含有n個樣本的數(shù)據(jù)集X={x1,x2, …,xn}?Rq,將數(shù)據(jù)集X劃分為k(2≤k

    (1)

    式(1)中n為樣本數(shù),k為樣本類別數(shù),dij=‖xi-cj‖為xi-cj的歐式距離,xi為第i個樣本,cj為第j類的類中心,m為權(quán)重指數(shù)。uji為第i(i=1,2,…,n)個樣本隸屬于第j(j=1,2,…,k)類的隸屬度值,uji的模糊隸屬度矩陣U=[uji]k×n,U?Mfk,Mfk是數(shù)據(jù)集X的模糊k劃分空間

    求解式(1)的優(yōu)化問題可得到如下方程

    (2)

    (3)

    則廣義模糊K調(diào)和均值聚類可通過以下迭代運算實現(xiàn):

    (1)初始化

    固定生菜近紅外光譜樣本數(shù)n(+∞>n>1),樣本類別數(shù)k(n>k>1),權(quán)重指數(shù)m(+∞>m>1)和p(+∞≥p≥2),最大迭代數(shù)rmax,誤差上限值ε,初始化類中心cj,0;

    (2) 計算第r(r=1,2,…,rmax)次迭代時的隸屬度值uji,r

    隸屬度值uji,r表示第r次迭代計算時第i(i=1,2,…,n)個樣本隸屬于第j(j=1,2,…,k)類的隸屬度值,dij,r=‖xi-cj,r-1‖,‖xi-cj,r-1‖為xi-cj,r-1的歐式距離,xi為第個樣本,cj,r-1為第r-1次迭代時第j類的類中心,dil,r=‖xi-cl,r-1‖,cl,r-1為第r-1次迭代時第l類的類中心;

    (3) 計算第r次迭代時的類中心cj,r(j=1,2,…,k)

    其中cj,r為第r次迭代時第j類的類中心;

    當(dāng)maxj‖cj,r-cj,r-1‖<ε或者r=rmax時,迭代結(jié)束; 否則,r=r+1,返回步驟(2)繼續(xù)迭代計算。

    3 結(jié)果與討論

    3.1 生菜近紅外光譜的壓縮和鑒別信息提取

    60個生菜樣本分為三個儲藏時間: “儲藏時間1”,“儲藏時間2”和“儲藏時間3”。 “儲藏時間1”表示第1次采集生菜近紅外光譜; “儲藏時間2”表示第1次采集后經(jīng)過12h第2次采集生菜近紅外光譜; “儲藏時間3” 表示第2次采集后經(jīng)過12h第3次采集生菜近紅外光譜。每個儲藏時間看作一類,則三個儲藏時間為三類,共有180個近紅外光譜樣本數(shù)據(jù)。由于這些數(shù)據(jù)的維數(shù)是1 557維,這給計算和分析帶來困難。另外,這些近紅外光譜數(shù)據(jù)含有一些無關(guān)的冗余信息。為解決這些問題,本文采用主成分分析將1 557維的生菜近紅外光譜樣本數(shù)據(jù)壓縮為20維數(shù)據(jù)以保留大部分的有用光譜信息和減少無用的冗余信息。為了提高分類的準確率,采用線性判別分析(LDA)提取數(shù)據(jù)的鑒別信息。取訓(xùn)練集樣本數(shù)為60個(即每類20個),測試集樣本數(shù)120(即每類40個),鑒別向量數(shù)為2,LDA將20維的測試集樣本數(shù)據(jù)投影到其鑒別向量(由訓(xùn)練樣本計算得到)上可得到LDA的得分圖(二維數(shù)據(jù))如圖2所示。由圖2可看出,“儲藏時間1”和“儲藏時間2”的數(shù)據(jù)存在重疊部分,而“儲藏時間3”與“儲藏時間1”和“儲藏時間2”的數(shù)據(jù)幾乎沒有重疊。重疊部分往往會導(dǎo)致誤分類,給分類器的設(shè)計帶來一定難度。

    圖2 LDA得分圖

    3.2 計算初始聚類中心

    設(shè)置模糊C均值聚類(FCM)[12]的初始參數(shù)值: 測試集樣本數(shù)n=120,樣本類別數(shù)k=3,權(quán)重指數(shù)m=2.0,最大迭代數(shù)rmax=100,誤差上限值ε=0.000 01。然后對圖2的測試數(shù)據(jù)運行模糊C-均值聚類直到迭代終止以得到三個類別的聚類中心

    c1,0=[0.014 5 0.005 8]

    c2,0=[-0.004 4 0.005 9]

    將上述聚類中心作為FKHM和GFKHM的初始聚類中心。

    3.3 FKHM和GFKHM的模糊聚類分析

    設(shè)置FKHM和GFKHM的初始參數(shù)值: 生菜近紅外光譜測試集樣本數(shù)n=120,樣本類別數(shù)k=3,權(quán)重指數(shù)m=2.0和p≥2,最大迭代數(shù)rmax=100,誤差上限值ε=0.000 01,確定初始聚類中心(3.2節(jié)得到的類中心)。

    3.3.1 聚類準確率

    分別對測試集樣本數(shù)據(jù)運行FKHM和GFKHM,GFKHM的權(quán)重指數(shù)p=10時運行到迭代計算終止時得到的隸屬度值,再根據(jù)GFKHM的隸屬度值得到的聚類準確率為92.5%,而FKHM迭代計算終止時得到的隸屬度值并依據(jù)該隸屬度值得到聚類準確率為90.0%。當(dāng)GFKHM的權(quán)重指數(shù)p取值不同時其準確率如表1所示,當(dāng)GFKHM的權(quán)重指數(shù)p=2時,GFKHM算法就是FKHM算法,此時兩者聚類準確率相同; 當(dāng)GFKHM的權(quán)重指數(shù)p=4時,GFKHM的聚類準確率低于FKHM的聚類準確率; 當(dāng)GFKHM的權(quán)重指數(shù)p=8,10,12,…,24時,GFKHM的聚類準確率高于FKHM的聚類準確率。

    表1 GFKHM的聚類準確率

    3.3.2 聚類中心分析

    對測試樣本集運行GFKHM和FKHM,可得到GFKHM的聚類中心矩陣CGFKHM為

    FKHM的聚類中心矩陣CFKHM為

    取測試樣本的每類均值作為真實的類中心矩陣CTrue[13]

    用矩陣的Frobenius范數(shù)[14]衡量聚類中心矩陣C*(*表示GFKHM/FKHM)逼近真實類中心矩陣CTure的程度[15]

    E*=‖C*-CTrue‖F(xiàn)

    E*越小則表示C*越逼近真實類中心。計算結(jié)果:EGFKHM=0.002 4,EFKHM=0.003 9。EGFKHM

    3.3.3 聚類收斂分析

    GFKHM和FKHM進行迭代計算時,第r次的聚類迭代時的收斂狀況用參數(shù)δr=maxj‖cj,r-cj,r-1‖來計算,當(dāng)δr<ε=0.000 01或者迭代次數(shù)r=rmax=100時,GFKHM和FKHM終止迭碟計算。GFKHM和FKHM的收斂狀況如圖3所示。由圖3可知,經(jīng)過有限次的迭代計算,GFKHM

    圖3 GFKHM和FKHM的收斂狀況

    和FKHM都能在最后實現(xiàn)算法收斂,但是,GFKHM經(jīng)過4次迭代實現(xiàn)了收斂而FKHM要經(jīng)過35次迭代才能最后收斂??梢?,GFKHM的收斂速度明顯快于FKHM。

    4 結(jié) 論

    為了進一步擴展模糊K調(diào)和均值聚類(FKHM)的應(yīng)用范圍,以歐式距離的p次方代替模糊K調(diào)和均值聚類(FKHM)中歐式距離的平方,提出了廣義模糊K調(diào)和均值聚類(GFKHM),并結(jié)合近紅外光譜技術(shù)建立生菜儲藏時間鑒別模型。與FKHM相比較,GFKHM建立的模型具有聚類準確率高,聚類中心更接近真實類中心,收斂速度快的優(yōu)點。實驗結(jié)果表明: GFKHM結(jié)合近紅外光譜技術(shù),主成分分析和線性判別分析能快速,無損和高效地實現(xiàn)生菜儲藏時間鑒別,為實現(xiàn)生菜近紅外光譜的聚類分析提供了一種新方法。

    [1] SUN Jun, JIN Xia-ming, MAO Han-ping, et al(孫 俊,金夏明,毛罕平,等). Chinese Journal of Analytical Chemistry(化學(xué)分析), 2014, 42(5): 672.

    [2] Peng X, Yang J, Cui P, et al. LWT-Food Science and Technology, 2015, 60(1): 300.

    [3] Simon A H, Stewart F G, Emmanuelle C, et al. Food Chemistry, 2013, 136(3-4): 1557.

    [4] Warawut S, Guangli N, Liu R, et al. Computers and Electronics in Agriculture, 2013, 91: 87.

    [5] Panmanas S, Munehiro T, Takayuki K, et al. Journal of Food Engineering, 2012, 112(3): 218.

    [6] He W, Zhou J, Cheng H, et al. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2012, 86: 399.

    [7] Sun J, Jin X M, Mao H P, et al. Journal of Pure and Applied Microbiology, 2013, 7: 747.

    [8] ZHAO Heng, YANG Wan-hai, ZHANG Gao-yu(趙 恒,楊萬海,張高煜). Journal of Xidian University(西安電子科技大學(xué)學(xué)報), 2005, 32(4): 603.

    [9] Abdeyazdan M. The Journal of Supercomputing, 2014, 68(2): 574.

    [10] Hung C H, Chiou H M, Yang W N. Applied Mathematical Modelling, 2013, 37(24): 10123.

    [11] Bobrowski L, Bezedek J C. IEEE Transaction on SMC, 1991, 21(3): 545.

    [12] Gou J, Hou F, Chen W, et al. Neurocomputing, 2015, 151(3): 1293.

    [13] Pal N R, Pal K, Bezdek J C. IEEE Trans. Fuzzy Systems, 2005, 13(4): 517.

    [14] D?sc?lescu S, Iovanov M C, PreduS. Linear Algebra and Its Applications, 2013, 439(10): 3166.

    [15] WU Xiao-hong, ZHOU Jian-jiang(武小紅,周建江). Acta Electronica Sinica(電子學(xué)報), 2008, 36(10): 1996.

    (Received Jan. 28, 2015; accepted Apr. 22, 2015)

    Discrimination of Lettuce Storage Time Using Near Infrared Spectroscopy Based on Generalized Fuzzy K-Harmonic Means Clustering

    WU Xiao-hong1, 2, PAN Ming-hui3, WU Bin4, JI Gang1, SUN Jun1, 2

    1. School of Electrical and Information Engineering, Jiangsu University, Zhenjiang 212013, China 2. Key Laboratory of Facility Agriculture Measurement and Control Technology and Equipment of Machinery Industry, Jiangsu University, Zhenjiang 212013, China 3. Jingjiang College, Jiangsu University, Zhenjiang 212013, China 4. Department of Information Engineering, Chuzhou Vocational Technology College, Chuzhou 239000, China

    Lettuce storage time is an important factor affecting the lettuce freshness. To realize the non-destructive, rapid and effective discrimination of lettuce storage time, generalized fuzzy K-harmonic means (GFKHM) clustering was proposed by introducing the pth power of Euclidean distance into fuzzy K-harmonic means (FKHM) clustering to replace the square of Euclidean distance in FKHM, and furthermore GFKHM was applied in the discrimination of lettuce storage time. Sixty fresh lettuce samples were prepared as the research object, and the near infrared reflectance (NIR) spectra of lettuces were collected by Antaris Ⅱ near infrared spectrometer with a spectral range of 10 000~4 000 cm-1for three 12-hour detections. Firstly, the 1 557-dimensional NIR spectra were reduced by principal component analysis (PCA) to decrease redundant information. After the first 20 principal components were selected, PCA translated the 1 557-dimensional NIR spectra into the 20-dimensional data. Secondly, linear discriminant analysis (LDA) was used to extract the discriminant information from the 20-dimensional data to improve the clustering accuracy. With the first two discriminant vectors, LDA translated the 20-dimensional data into the two-dimensional data. Finally, the cluster centers from fuzzy C-means clustering (FCM) were set as the initial cluster centers for FKHM and GFKHM and fuzzy membership values of FKHM and GFKHM were calculated to identify lettuce storage time. The experimental results showed that the discrimination accuracy of GFKHM has achieved 92.5% which was higher than that of FKHM. The cluster centers of GFKHM were much closer to the true cluster centers in comparison with FKHM. Furthermore, the convergence of the GFKHM was significantly faster than FKHM. Near infrared spectroscopy coupled with GFKHM, PCA and LDA could cluster NIR spectra of lettuce quickly and correctly, and this provided a fast and nondestructive method for identifying lettuce storage time.

    Near infrared spectroscopy; Lettuce; Storage time; Linear discriminant analysis; Fuzzy K-harmonic means clustering

    2015-01-28,

    2015-04-22

    國家自然科學(xué)基金項目(31101082),江蘇高校優(yōu)勢學(xué)科建設(shè)工程項目PAPD資助

    武小紅,1971年生,江蘇大學(xué)電氣信息工程學(xué)院副教授 e-mail: wxh_www@163.com

    O657.3

    A

    10.3964/j.issn.1000-0593(2016)06-1721-05

    猜你喜歡
    調(diào)和生菜儲藏
    脆嫩爽口的生菜
    中老年保健(2022年6期)2022-08-19 01:44:10
    五味調(diào)和醋當(dāng)先
    生菜怎么吃更健康
    從“調(diào)結(jié)”到“調(diào)和”:打造“人和”調(diào)解品牌
    西藥藥品的儲藏探討
    調(diào)和映照的雙Lipschitz性質(zhì)
    生菜?你愛吃圓的還是散葉兒的?
    食品與健康(2017年3期)2017-03-15 18:07:48
    生菜有道
    美食堂(2017年2期)2017-02-14 00:16:54
    便于儲藏的美味
    “蘊藏”和“儲藏”
    精品少妇一区二区三区视频日本电影 | 国产一卡二卡三卡精品 | 亚洲国产最新在线播放| 最新的欧美精品一区二区| 啦啦啦在线免费观看视频4| 大香蕉久久网| av在线app专区| 一级毛片我不卡| 欧美日韩亚洲国产一区二区在线观看 | 国产深夜福利视频在线观看| 亚洲综合色网址| 午夜福利,免费看| 丝袜在线中文字幕| 免费黄色在线免费观看| 国产精品亚洲av一区麻豆 | 久久性视频一级片| 午夜影院在线不卡| 高清欧美精品videossex| 少妇的丰满在线观看| 亚洲激情五月婷婷啪啪| 十八禁网站网址无遮挡| 国产一级毛片在线| 久久国产亚洲av麻豆专区| svipshipincom国产片| 国产激情久久老熟女| 久久久国产一区二区| 日韩av在线免费看完整版不卡| 男人添女人高潮全过程视频| 一二三四在线观看免费中文在| 一级毛片电影观看| 国产熟女午夜一区二区三区| 国产男女超爽视频在线观看| 秋霞伦理黄片| 免费黄网站久久成人精品| 中文乱码字字幕精品一区二区三区| 成人三级做爰电影| 十八禁人妻一区二区| 久久久欧美国产精品| 免费看av在线观看网站| 最近手机中文字幕大全| 国产亚洲最大av| 午夜老司机福利片| 伊人亚洲综合成人网| 久久亚洲国产成人精品v| 日韩一区二区三区影片| 丝袜美足系列| 亚洲欧洲精品一区二区精品久久久 | 另类亚洲欧美激情| 在线观看www视频免费| 国产免费一区二区三区四区乱码| 制服丝袜香蕉在线| 国产毛片在线视频| 日韩视频在线欧美| 狠狠精品人妻久久久久久综合| 2021少妇久久久久久久久久久| 午夜日韩欧美国产| 亚洲欧洲精品一区二区精品久久久 | 亚洲一码二码三码区别大吗| 女性生殖器流出的白浆| 精品人妻在线不人妻| 无遮挡黄片免费观看| 久久久久久久大尺度免费视频| 麻豆精品久久久久久蜜桃| 国产有黄有色有爽视频| 国产毛片在线视频| 国产色婷婷99| 亚洲av在线观看美女高潮| 免费看av在线观看网站| 丝袜脚勾引网站| www日本在线高清视频| 大片免费播放器 马上看| 日日爽夜夜爽网站| 亚洲,一卡二卡三卡| 丁香六月欧美| 国产1区2区3区精品| 中文字幕最新亚洲高清| 国产精品一区二区在线不卡| 搡老乐熟女国产| 国产人伦9x9x在线观看| 亚洲av电影在线观看一区二区三区| 女人精品久久久久毛片| 一级毛片 在线播放| 成人国产麻豆网| 日本一区二区免费在线视频| 又大又黄又爽视频免费| 在线观看免费午夜福利视频| 69精品国产乱码久久久| 国产 一区精品| 男人操女人黄网站| 男人操女人黄网站| 日本一区二区免费在线视频| 一区二区av电影网| 国产极品粉嫩免费观看在线| 国产精品无大码| 在线观看免费日韩欧美大片| 高清黄色对白视频在线免费看| 色网站视频免费| 岛国毛片在线播放| 亚洲成人国产一区在线观看 | 国产在视频线精品| 久久青草综合色| 两个人看的免费小视频| 免费av中文字幕在线| 777米奇影视久久| 精品久久蜜臀av无| 日韩一卡2卡3卡4卡2021年| 久久久久精品久久久久真实原创| 国产黄频视频在线观看| 国产成人精品久久二区二区91 | 汤姆久久久久久久影院中文字幕| 国产av码专区亚洲av| 中文字幕人妻丝袜制服| videosex国产| 亚洲成色77777| 亚洲av电影在线进入| 黄色一级大片看看| 亚洲精品美女久久av网站| 国产一区二区三区综合在线观看| 国产xxxxx性猛交| tube8黄色片| 国产亚洲av片在线观看秒播厂| 久久久国产欧美日韩av| 丝袜喷水一区| 巨乳人妻的诱惑在线观看| 1024香蕉在线观看| 老汉色av国产亚洲站长工具| 99国产综合亚洲精品| 国产精品久久久久久人妻精品电影 | 国产精品久久久久久久久免| av线在线观看网站| 亚洲第一av免费看| a级片在线免费高清观看视频| 精品人妻熟女毛片av久久网站| 午夜福利免费观看在线| 国产福利在线免费观看视频| 亚洲人成电影观看| 亚洲美女视频黄频| 国产成人欧美在线观看 | 999精品在线视频| 超碰97精品在线观看| 欧美乱码精品一区二区三区| 免费av中文字幕在线| 国产欧美日韩综合在线一区二区| 久久久久久久国产电影| a 毛片基地| 精品国产国语对白av| 欧美另类一区| 亚洲激情五月婷婷啪啪| 伊人久久国产一区二区| 男女高潮啪啪啪动态图| 久久久久精品久久久久真实原创| 亚洲人成77777在线视频| 黄片小视频在线播放| 卡戴珊不雅视频在线播放| 亚洲国产精品一区二区三区在线| 久久精品亚洲熟妇少妇任你| 超色免费av| 午夜福利乱码中文字幕| 18禁观看日本| 一二三四中文在线观看免费高清| 欧美最新免费一区二区三区| 满18在线观看网站| 国产精品二区激情视频| 国产探花极品一区二区| av一本久久久久| 精品一区二区免费观看| 亚洲av国产av综合av卡| 美女脱内裤让男人舔精品视频| 人体艺术视频欧美日本| 宅男免费午夜| tube8黄色片| 男人添女人高潮全过程视频| 精品一品国产午夜福利视频| 亚洲第一av免费看| 日韩电影二区| 亚洲视频免费观看视频| 少妇精品久久久久久久| 九九爱精品视频在线观看| 乱人伦中国视频| 亚洲图色成人| 国产欧美日韩一区二区三区在线| 久久久久久久大尺度免费视频| 国产精品蜜桃在线观看| 天天躁夜夜躁狠狠躁躁| av线在线观看网站| 精品亚洲乱码少妇综合久久| 国产日韩欧美亚洲二区| 国产又爽黄色视频| 久久久久久久国产电影| 精品亚洲乱码少妇综合久久| av线在线观看网站| 赤兔流量卡办理| 亚洲国产精品999| 亚洲男人天堂网一区| 欧美激情 高清一区二区三区| 18禁国产床啪视频网站| kizo精华| 欧美黑人欧美精品刺激| 久久久久精品人妻al黑| 亚洲国产日韩一区二区| 国产成人精品在线电影| 国产一区二区三区av在线| 最黄视频免费看| 久久久精品国产亚洲av高清涩受| 亚洲国产av影院在线观看| 99热网站在线观看| 亚洲av男天堂| 一本—道久久a久久精品蜜桃钙片| 国产欧美日韩综合在线一区二区| 国产精品 欧美亚洲| 国产成人av激情在线播放| 成年女人毛片免费观看观看9 | 人人妻人人爽人人添夜夜欢视频| 久久精品aⅴ一区二区三区四区| 美女中出高潮动态图| 亚洲精品一二三| 男女床上黄色一级片免费看| 亚洲国产日韩一区二区| 国产成人精品在线电影| 欧美在线一区亚洲| 亚洲一区二区三区欧美精品| 在线观看免费日韩欧美大片| 日本午夜av视频| 国产精品亚洲av一区麻豆 | 国产在视频线精品| 99久国产av精品国产电影| 王馨瑶露胸无遮挡在线观看| 国精品久久久久久国模美| 精品卡一卡二卡四卡免费| 卡戴珊不雅视频在线播放| 亚洲 欧美一区二区三区| 99久久精品国产亚洲精品| 欧美中文综合在线视频| 亚洲国产中文字幕在线视频| 久久久精品94久久精品| 国产成人av激情在线播放| xxxhd国产人妻xxx| 国产精品一区二区在线不卡| 校园人妻丝袜中文字幕| 国产成人a∨麻豆精品| 国产片内射在线| 国产精品一国产av| 亚洲色图 男人天堂 中文字幕| 日韩av不卡免费在线播放| 国产精品一区二区在线观看99| 国产伦人伦偷精品视频| 精品久久久久久电影网| 久久久国产精品麻豆| 亚洲人成77777在线视频| 肉色欧美久久久久久久蜜桃| 99国产综合亚洲精品| 宅男免费午夜| 性高湖久久久久久久久免费观看| 国产成人精品在线电影| videosex国产| 国产成人精品久久久久久| 老鸭窝网址在线观看| 久久天躁狠狠躁夜夜2o2o | 999久久久国产精品视频| 国产日韩欧美视频二区| 欧美成人午夜精品| 亚洲三区欧美一区| 999久久久国产精品视频| 天天躁夜夜躁狠狠躁躁| 99久久人妻综合| 久久国产精品男人的天堂亚洲| 99国产综合亚洲精品| 久久免费观看电影| 久久免费观看电影| www.精华液| 欧美乱码精品一区二区三区| 黄网站色视频无遮挡免费观看| 熟妇人妻不卡中文字幕| 亚洲欧美色中文字幕在线| 赤兔流量卡办理| 久久久久久久久久久免费av| 免费人妻精品一区二区三区视频| 亚洲一级一片aⅴ在线观看| 亚洲精品国产一区二区精华液| 99久久99久久久精品蜜桃| 国产精品国产三级国产专区5o| 丝袜喷水一区| 欧美在线黄色| 欧美变态另类bdsm刘玥| 国产男女超爽视频在线观看| 交换朋友夫妻互换小说| 男人操女人黄网站| 久久久久精品国产欧美久久久 | 黄色怎么调成土黄色| 亚洲综合精品二区| 一边摸一边做爽爽视频免费| 在线观看www视频免费| 国产欧美日韩一区二区三区在线| 国产精品香港三级国产av潘金莲 | 最近的中文字幕免费完整| 亚洲精品中文字幕在线视频| 国产乱人偷精品视频| 亚洲国产精品999| 欧美精品一区二区免费开放| 99久久精品国产亚洲精品| 国产人伦9x9x在线观看| 一区二区三区乱码不卡18| 亚洲精品成人av观看孕妇| 亚洲精品第二区| 亚洲图色成人| 国产av一区二区精品久久| 黄色视频不卡| 伊人亚洲综合成人网| 69精品国产乱码久久久| 国产亚洲午夜精品一区二区久久| 亚洲精品av麻豆狂野| 美国免费a级毛片| 一二三四在线观看免费中文在| av线在线观看网站| 久久天堂一区二区三区四区| www.av在线官网国产| 亚洲国产精品一区二区三区在线| 日本av免费视频播放| 精品久久蜜臀av无| 久久精品久久精品一区二区三区| 99久国产av精品国产电影| 欧美成人精品欧美一级黄| 国产精品香港三级国产av潘金莲 | 国产亚洲av片在线观看秒播厂| 久久久精品免费免费高清| a级片在线免费高清观看视频| 国产精品香港三级国产av潘金莲 | 精品亚洲成国产av| 国产在线一区二区三区精| 精品一品国产午夜福利视频| a级片在线免费高清观看视频| 老司机影院毛片| 成人黄色视频免费在线看| 精品久久蜜臀av无| 中文字幕高清在线视频| 欧美精品亚洲一区二区| 亚洲图色成人| 欧美老熟妇乱子伦牲交| 欧美日韩综合久久久久久| 91精品国产国语对白视频| 国产精品免费视频内射| 80岁老熟妇乱子伦牲交| 午夜福利在线免费观看网站| 两性夫妻黄色片| 成人国产av品久久久| 国产一区二区 视频在线| 国产精品久久久久成人av| 亚洲国产成人一精品久久久| 亚洲成av片中文字幕在线观看| 母亲3免费完整高清在线观看| 中国三级夫妇交换| 亚洲国产最新在线播放| 99久久人妻综合| 国产男女超爽视频在线观看| 亚洲精品,欧美精品| 久久99热这里只频精品6学生| 国产精品久久久av美女十八| 69精品国产乱码久久久| 九色亚洲精品在线播放| 欧美老熟妇乱子伦牲交| 中国国产av一级| 伊人久久国产一区二区| 国产精品国产三级专区第一集| 99久久综合免费| h视频一区二区三区| 亚洲美女视频黄频| 丝袜脚勾引网站| 国产欧美亚洲国产| 欧美国产精品一级二级三级| 免费少妇av软件| 这个男人来自地球电影免费观看 | 亚洲人成网站在线观看播放| 大片电影免费在线观看免费| 一本色道久久久久久精品综合| 国产成人啪精品午夜网站| 欧美日韩精品网址| 国产乱人偷精品视频| 99精国产麻豆久久婷婷| 国产精品秋霞免费鲁丝片| 国产欧美日韩综合在线一区二区| 久久久久久久精品精品| 9热在线视频观看99| 新久久久久国产一级毛片| 欧美黄色片欧美黄色片| 只有这里有精品99| 久久ye,这里只有精品| 午夜激情av网站| 日韩精品免费视频一区二区三区| 最近最新中文字幕大全免费视频 | 国产又爽黄色视频| 亚洲欧美一区二区三区黑人| av又黄又爽大尺度在线免费看| 午夜福利影视在线免费观看| 欧美av亚洲av综合av国产av | 视频在线观看一区二区三区| av片东京热男人的天堂| 欧美日韩国产mv在线观看视频| 亚洲欧美激情在线| 国产一卡二卡三卡精品 | 最近中文字幕2019免费版| 少妇人妻久久综合中文| 国产一区亚洲一区在线观看| 人人澡人人妻人| 在线 av 中文字幕| h视频一区二区三区| 最近最新中文字幕大全免费视频 | 伊人亚洲综合成人网| 亚洲av福利一区| 亚洲四区av| 亚洲伊人色综图| 一级,二级,三级黄色视频| 人妻人人澡人人爽人人| 精品人妻熟女毛片av久久网站| 我的亚洲天堂| 久久久精品免费免费高清| 极品人妻少妇av视频| 激情视频va一区二区三区| 中文天堂在线官网| 国产精品嫩草影院av在线观看| 人妻一区二区av| 看免费av毛片| 另类亚洲欧美激情| 久久久国产一区二区| av线在线观看网站| 一区二区三区精品91| 一边摸一边抽搐一进一出视频| 90打野战视频偷拍视频| 黄色怎么调成土黄色| 国产又色又爽无遮挡免| 2018国产大陆天天弄谢| 国产深夜福利视频在线观看| 看非洲黑人一级黄片| 人妻人人澡人人爽人人| 精品福利永久在线观看| 免费在线观看视频国产中文字幕亚洲 | 欧美变态另类bdsm刘玥| 日日撸夜夜添| 韩国精品一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 在线精品无人区一区二区三| 亚洲成人av在线免费| 国产免费现黄频在线看| 亚洲激情五月婷婷啪啪| av天堂久久9| 久久久久精品国产欧美久久久 | 欧美亚洲 丝袜 人妻 在线| 无限看片的www在线观看| 性高湖久久久久久久久免费观看| 在线观看三级黄色| 日本91视频免费播放| 精品国产乱码久久久久久小说| 一级黄片播放器| 51午夜福利影视在线观看| 嫩草影院入口| 久久精品aⅴ一区二区三区四区| 久久ye,这里只有精品| 久久精品aⅴ一区二区三区四区| 大码成人一级视频| 晚上一个人看的免费电影| 香蕉国产在线看| 久久久久久免费高清国产稀缺| 欧美日韩亚洲综合一区二区三区_| 亚洲,一卡二卡三卡| 99re6热这里在线精品视频| 亚洲伊人久久精品综合| 老熟女久久久| 亚洲欧洲国产日韩| 亚洲成人免费av在线播放| 女人精品久久久久毛片| 国产在线一区二区三区精| 777米奇影视久久| 午夜福利视频在线观看免费| 视频在线观看一区二区三区| 欧美成人精品欧美一级黄| 国产精品无大码| 一边摸一边做爽爽视频免费| 看非洲黑人一级黄片| 飞空精品影院首页| 亚洲av日韩在线播放| 欧美成人午夜精品| 这个男人来自地球电影免费观看 | 婷婷色麻豆天堂久久| 丝袜脚勾引网站| av.在线天堂| 亚洲欧美一区二区三区久久| 亚洲视频免费观看视频| 性少妇av在线| 国产亚洲一区二区精品| 欧美激情极品国产一区二区三区| 只有这里有精品99| 高清欧美精品videossex| 久久精品熟女亚洲av麻豆精品| 99久久精品国产亚洲精品| 免费在线观看视频国产中文字幕亚洲 | 亚洲成国产人片在线观看| 亚洲精品国产av成人精品| 成年动漫av网址| 两性夫妻黄色片| 我要看黄色一级片免费的| 亚洲人成网站在线观看播放| 欧美日韩一级在线毛片| 亚洲四区av| 99re6热这里在线精品视频| 最近最新中文字幕免费大全7| 秋霞伦理黄片| 国产精品三级大全| 亚洲在久久综合| 9色porny在线观看| 99re6热这里在线精品视频| 人妻人人澡人人爽人人| 中文字幕另类日韩欧美亚洲嫩草| netflix在线观看网站| 亚洲国产精品一区二区三区在线| 久久久久网色| av免费观看日本| 中文字幕制服av| 视频区图区小说| videos熟女内射| 各种免费的搞黄视频| a级片在线免费高清观看视频| 成人漫画全彩无遮挡| 亚洲精品久久成人aⅴ小说| 黄色怎么调成土黄色| 国产在线视频一区二区| 我的亚洲天堂| avwww免费| 搡老乐熟女国产| 黑人巨大精品欧美一区二区蜜桃| 99热网站在线观看| 久久精品国产a三级三级三级| 夫妻午夜视频| 亚洲人成网站在线观看播放| 国产成人啪精品午夜网站| 国产又色又爽无遮挡免| 黑人欧美特级aaaaaa片| 中文字幕另类日韩欧美亚洲嫩草| 色婷婷久久久亚洲欧美| 看免费成人av毛片| 卡戴珊不雅视频在线播放| 性高湖久久久久久久久免费观看| 下体分泌物呈黄色| 午夜日本视频在线| 韩国高清视频一区二区三区| 最近最新中文字幕大全免费视频 | 成年人午夜在线观看视频| 满18在线观看网站| 精品少妇内射三级| 亚洲av中文av极速乱| 久久久久久免费高清国产稀缺| 哪个播放器可以免费观看大片| 欧美成人午夜精品| 国产97色在线日韩免费| 美女午夜性视频免费| 免费久久久久久久精品成人欧美视频| 亚洲人成网站在线观看播放| 国产日韩一区二区三区精品不卡| 在线天堂最新版资源| 国产精品熟女久久久久浪| 国产成人精品久久久久久| 午夜福利影视在线免费观看| 天堂中文最新版在线下载| 亚洲精品av麻豆狂野| 精品国产乱码久久久久久男人| 老司机亚洲免费影院| av福利片在线| 80岁老熟妇乱子伦牲交| 午夜91福利影院| 大香蕉久久成人网| 波多野结衣av一区二区av| 两个人看的免费小视频| 中文字幕另类日韩欧美亚洲嫩草| 国产精品久久久久久精品古装| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产亚洲av涩爱| 一边摸一边抽搐一进一出视频| 久久久久精品久久久久真实原创| 赤兔流量卡办理| 一级毛片我不卡| 国产精品成人在线| 亚洲国产看品久久| 嫩草影院入口| 91国产中文字幕| 搡老岳熟女国产| 一二三四中文在线观看免费高清| 麻豆av在线久日| 在线观看人妻少妇| 黑人欧美特级aaaaaa片| 在线看a的网站| 国产欧美日韩一区二区三区在线| 欧美另类一区| 你懂的网址亚洲精品在线观看| 国产无遮挡羞羞视频在线观看| 一级片'在线观看视频| 成人亚洲精品一区在线观看| 亚洲三区欧美一区| 国产男人的电影天堂91| 久久久国产精品麻豆| 99热全是精品| 成人国产av品久久久| 最近中文字幕2019免费版| 中文字幕精品免费在线观看视频| 日韩人妻精品一区2区三区| 亚洲图色成人| 少妇的丰满在线观看| 欧美在线黄色| 国产精品女同一区二区软件| 大香蕉久久网| 观看美女的网站| 亚洲国产欧美一区二区综合| 久久久精品区二区三区| 观看av在线不卡| 一区二区日韩欧美中文字幕| 午夜免费鲁丝| 狠狠婷婷综合久久久久久88av| 国产精品国产三级国产专区5o|