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

    二維相關(guān)可見-近紅外光譜結(jié)合支持向量機評價豬肉新鮮度

    2018-10-08 02:50:32王文秀彭彥昆孫宏偉魏文松鄭曉春
    食品科學 2018年18期
    關(guān)鍵詞:光譜信息新鮮度波段

    王文秀,彭彥昆*,孫宏偉,魏文松,鄭曉春

    (中國農(nóng)業(yè)大學工學院,國家農(nóng)產(chǎn)品加工技術(shù)裝備研發(fā)分中心,北京 100083)

    生鮮肉作為蛋白質(zhì)的主要來源,在人們的日常飲食結(jié)構(gòu)中占據(jù)非常重要的地位。2016年我國肉品的總產(chǎn)量達8 625萬 t,整個肉品行業(yè)發(fā)展?jié)摿薮?。然而,由于肉品本身富含營養(yǎng)物質(zhì)和水分,容易被其他微生物入侵而發(fā)生系列變化,導(dǎo)致肉的腐敗變質(zhì),給消費者的健康安全造成威脅。因此,對肉的新鮮度進行快速評價具有非常重要的現(xiàn)實意義。

    可見-近紅外光譜技術(shù)具有無損傷、操作簡單、無需樣品前處理的優(yōu)點,已被證實可以對生鮮肉的品質(zhì)安全參數(shù)進行快速分析[1-4]。目前,基于光譜技術(shù)對肉品新鮮度評價可以參照GB 2707—2016《鮮(凍)畜、禽產(chǎn)品》預(yù)測揮發(fā)性鹽基氮值(total volatile basic nitrogen,TVB-N)評價新鮮度。也可借助化學計量學方法,直接將光譜信息與新鮮度建立聯(lián)系作出評價。馬世榜等[5]利用近紅外光譜結(jié)合支持向量機(support vector machine,SVM)對牛肉的新鮮度進行定性判別,驗證集正確率為88.89%。蔡健榮等[6]利用近紅外光譜對豬肉中TVB-N進行定量分析,預(yù)測集相關(guān)系數(shù)為0.823 8。然而,由于近紅外光譜數(shù)據(jù)量大,光譜信息復(fù)雜且具有多重共線性,因此不少學者展開特征波長篩選的研究。如馬世榜[7]、劉飛[8]等利用無信息變量消除和連續(xù)投影算法提取有效變量,建立了更為簡化的TVB-N預(yù)測模型。這些方法大多通過多次采樣分析,構(gòu)建多個子模型,通過評價模型誤差確定最優(yōu)變量組合。然而,不同方法篩選的變量個數(shù)和分布常常有所不同,部分變量缺乏化學意義上的響應(yīng)機理解釋,因此變量優(yōu)選仍是目前近紅外光譜分析面臨的重要問題。

    二維相關(guān)光譜通過將光譜信號擴展到第二維上,可以使原來一維光譜上的弱峰和重疊峰更加清晰,從而提高光譜分辨率。同時,通過添加一定的外界擾動,對動態(tài)光譜進行二維相關(guān)計算,可以解析光譜細微特征的變化,尋找與微擾相關(guān)的敏感變量[9-11]。目前,已有學者將二維相關(guān)光譜應(yīng)用在牛奶摻雜鑒別[12]、食用油摻雜檢測[13]、中藥摻偽[14]、蛋白飼料原料判別[15]、紅棗營養(yǎng)成分變化解析[16]等的研究上,表明通過分析二維相關(guān)光譜特性,可以觀察到各吸收峰隨外部擾動的變化關(guān)系以及對外部擾動的敏感程度。結(jié)合本實驗的研究,將所有樣品(從新鮮到腐敗)整體作為待研究體系,以TVB-N為外部擾動,原始光譜可看作是外擾誘發(fā)的動態(tài)光譜。通過二維光譜計算,解析同步光譜,觀察外部擾動對不同波長下動態(tài)光譜的影響,可以從機理上確立與外部擾動相關(guān)的敏感波長。目前將二維相關(guān)近紅外光譜應(yīng)用在生鮮肉新鮮度評價的研究尚鮮見報道。

    為探究二維相關(guān)光譜優(yōu)選新鮮度特征變量的可行性,本研究以生鮮豬肉中TVB-N為外界干擾對光譜進行調(diào)制,通過研究在外部擾動存在下體系的變化,解析二維相關(guān)同步光譜和自相關(guān)譜,分析隨外部擾動變化的光譜特征信息,確定與TVB-N相關(guān)的特征變量,建立生鮮豬肉的快速判別模型。

    1 材料與方法

    1.1 材料

    冷鮮豬肉背最長肌購買于北京美廉美超市。將其修整為長×寬×高約8 cm×5 cm×2.5 cm的肉塊,放置在自封袋中運輸至無損檢測實驗室。將所有樣品按順序編號后無積壓放置在4 ℃冰箱中保存,每12 h取出2 個樣品進行光譜采集和理化值測定,共進行15 d,獲得有效樣品58 個。

    1.2 儀器及光譜采集

    自行搭建的實驗用可見-近紅外光譜采集系統(tǒng)主要包括鹵鎢燈光源、AvaSpec-2048光譜儀(荷蘭Avantes公司)、光纖、環(huán)形光導(dǎo)、計算機等硬件。光譜范圍為350~1 100 nm,分辨率為0.59 nm。采集光譜前,將儀器預(yù)熱約30 min,然后依次采集參比光譜、暗背景光譜以實現(xiàn)校正,最后采集樣品的光譜信息。整個實驗過程中,保證光纖探頭與待測物距離一致,每個樣品采集5 個不同位置處光譜信息并求取平均值,從而減少實驗誤差。光譜數(shù)據(jù)的獲取和保存通過光譜儀配套軟件Avasoft完成。

    1.3 理化指標測定

    光譜信息采集完成后,參照GB 5009.228—2016《食品中揮發(fā)性鹽基氮的測定》中的方法對豬肉樣品中的TVB-N進行測定,并參照文獻[17]將TVB-N小于15 mg/100 g的樣品定義為新鮮肉,將TVB-N介于15~25 mg/100 g之間的樣品定義為次新鮮肉,將TVB-N大于25 mg/100 g的樣品定義為腐敗肉,以此為標準對豬肉的新鮮度進行劃分。

    1.4 數(shù)據(jù)處理與分析

    首先基于全波段光譜信息建立豬肉新鮮度的判別模型。由于近紅外光譜采集過程中常常伴隨基線漂移、雜散光等噪音信號,因此本研究采用標準正態(tài)變量變換(standard normal variate,SNV)、歸一化的預(yù)處理方法來減少干擾信息,并進行分析比較。建模采用了SVM判別分析方法,該方法以結(jié)構(gòu)風險最小化為準則構(gòu)造決策超平面,使不同種類樣品的間隔最大[18]。

    以TVB-N為外部微擾,進行二維相關(guān)同步光譜解析,尋找與新鮮度評價相關(guān)的特征變量。第1步:根據(jù)TVB-N實測值,從最小值和最大值之間以均勻濃度梯度共選取10 個代表性樣品用于二維相關(guān)光譜分析。第2步:提取光譜特征,本研究采用包絡(luò)線去除方法來擴大較弱的特性信息,同時壓抑背景光譜[19]。第3步:根據(jù)包絡(luò)線去除后的光譜,選擇不同TVB-N對應(yīng)光譜具有明顯差異的敏感波段。第4步:對上述敏感波段分別進行二維相關(guān)同步光譜分析,獲取其二維同步光譜及自相關(guān)譜。第5步:通過解析同步譜和自相關(guān)譜,明確與TVB-N變化密切相關(guān)的特征波長。

    基于所選的特征波長,分別從原始、經(jīng)SNV預(yù)處理、歸一化預(yù)處理后的光譜中提取相應(yīng)的光譜數(shù)據(jù)組成新的光譜矩陣,建立優(yōu)化的判別模型。整個數(shù)據(jù)分析流程圖如圖1所示。上述預(yù)處理和建模過程在Matlab 2012a平臺下利用PLS-Toolbox完成,二維光譜解析在薩斯喀徹溫大學加拿大光源編寫的MIDAS軟件下完成。

    圖1 數(shù)據(jù)分析流程圖Fig. 1 Flow chart of data analysis

    2 結(jié)果與分析

    2.1 新鮮度統(tǒng)計結(jié)果

    根據(jù)測定的TVB-N對58 個豬肉樣品的新鮮度進行統(tǒng)計分析,共有14 個新鮮肉、34 個次新鮮肉以及10 個腐敗肉。將所有樣品按照3∶1的比例分為校正集和驗證集,校正集有43 個樣品,驗證集有15 個樣品。將3 類樣品分別賦值為1、2、3,各子集樣品的新鮮度情況如表1所示。

    表1 3 種新鮮度豬肉的類別賦值及其在校正集和驗證集中分布Table 1 Category assignment of three kinds of pork and their distribution in calibration and prediction set

    2.2 一維可見-近紅外光譜分析

    由于原始光譜信息兩端噪音較大,信噪比低,因此只截取400~1 000 nm波段范圍內(nèi)的光譜用于后續(xù)定性判別分析。從新鮮、次新鮮和腐敗豬肉中各取一個樣品進行光譜信息的對比分析,其原始一維近紅外光譜如圖2a所示,3 類樣品整體上趨勢一致,在545~585 nm之間具有明顯差異,該波段范圍與肉中色素的變化有關(guān)[20]。經(jīng)過SNV和歸一化預(yù)處理后的光譜曲線如圖2b和c所示,2 種預(yù)處理方法均能有效消除基線漂移、光散射等噪音信號,提高光譜分辨率和靈敏度。

    圖2 樣品光譜曲線Fig. 2 Spectral curves of samples

    2.3 二維相關(guān)同步光譜分析

    為尋找能對豬肉新鮮度快速評價的關(guān)鍵變量,以TVB-N為外部微擾條件,選擇10 條代表性光譜并進行包絡(luò)線去除,以增強對TVB-N變化敏感波段的光譜特征,同時也增強每條曲線上各個波段之間光譜信息的對比性。包絡(luò)線去除前后的光譜曲線對比如圖3所示,經(jīng)過包絡(luò)線去除后光譜集中在0~1之間,且在460、590 nm及960 nm波長處出現(xiàn)明顯的波谷,可見光波段范圍內(nèi)的波谷與肉中肌紅蛋白的濃度和狀態(tài)有關(guān),960 nm波長處與N—H鍵的二級倍頻有關(guān)。在進行二維相關(guān)光譜分析之前,為盡可能詳盡的挖掘有效變量,避免微弱的特征信息被隱藏,結(jié)合包絡(luò)線去除后的光譜曲線,將不同TVB-N光譜曲線之間具有差異的波段細分為7 個子區(qū)間,分別進行二維同步光譜解析。這7 個子區(qū)間分別為400.1~429.4、430.6~494.8、496.0~550.4、555.1~584.0、585.7~680.3、835.9~954.7、955.8~999.5 nm。

    圖3 代表性光譜包絡(luò)線去除前后對比圖Fig. 3 Comparison of representative spectra before and after continuum removal

    對7 個子區(qū)間分別進行二維光譜分析,其三維投影圖和自相關(guān)譜圖分別如圖4、5所示。三維立體圖形的X軸和Y軸為獨立的變量軸,一般為波長或波數(shù),本研究中的變量為波長;Z軸為因變量軸,本研究中以相關(guān)強度為Z軸構(gòu)成三維投影圖[21]。同步譜反映了動態(tài)光譜信號隨外部擾動發(fā)生變化的協(xié)同程度,它關(guān)于主對角線對稱分布,可以從三維同步圖中直觀地觀察到自相關(guān)峰的位置、強度。自相關(guān)峰位于主對角線上,它的強度反映了該變量處光譜信號隨外部擾動變化的程度[22-28]。結(jié)合自相關(guān)譜圖,在波峰位置處相關(guān)強度大,表明該變量處光譜信號對微擾較敏感,即為所要尋找的與新鮮度密切相關(guān)的特征信息。

    結(jié)合圖4a、5a,在400.1~429.4 nm波長范圍內(nèi),主對角線上411 nm和428 nm波長處出現(xiàn)微弱的自相關(guān)峰,表明這兩個波長下的光譜信號對外部微擾敏感,與TVB-N的變化有一定的關(guān)系。根據(jù)Benesch等[29]的研究,415 nm為氧合血紅蛋白的吸收峰,430 nm為脫氧血紅蛋白的吸收峰,本研究中的波長位置與之稍有偏移,這可能與所用光譜儀器的響應(yīng)有關(guān)。結(jié)合圖4b、5b,在430.6~494.8 nm波長范圍內(nèi),434、461 nm和490 nm波長處出現(xiàn)自相關(guān)峰,434 nm為脫氧肌紅蛋白的吸收峰,該處的相關(guān)強度高于461 nm和490 nm,從三維投影圖的顏色參考圖中也可以直觀地看出。結(jié)合圖4c、5c,在496.0~550.4 nm波段范圍內(nèi),508、521 nm和537 nm波長處出現(xiàn)自相關(guān)峰,且537 nm處的強度高于前二者,表明該位置處對應(yīng)的官能團隨TVB-N變化更加明顯。根據(jù)Liu Yang等[30]的研究,505 nm為高鐵肌紅蛋白的吸收峰,535 nm為氧合肌紅蛋白的吸收峰,本研究的特征變量雖稍有偏移,但仍可將其歸屬為二者的吸收峰。結(jié)合圖4d、5d,在555.1~584.0 nm波段范圍內(nèi),560、569 nm和580 nm處出現(xiàn)自相關(guān)峰,其中580 nm處的強度最高。參考Mancini等[20]研究可知,530~580 nm波段范圍為肌紅蛋白色素的特征波段。結(jié)合圖4e、5e,在585.7~680.3 nm波段范圍內(nèi),595 nm波長處出現(xiàn)較強的自相關(guān)峰,在620、630 nm和667 nm波長處出現(xiàn)3 個相對弱的自相關(guān)峰。其中,630 nm為硫化肌紅蛋白的吸收峰,這是由于肉在放置過程中細菌產(chǎn)生的H2S氣體與肌紅蛋白結(jié)合,形成了硫化肌紅蛋白。結(jié)合圖4f、5f,在835.9~954.7 nm波段范圍內(nèi),937 nm波長處出現(xiàn)一個較強的自相關(guān)峰,其與C—H鍵的三級倍頻有關(guān)。結(jié)合圖4g、5g,在955.8~999.5 nm波段范圍內(nèi),977 nm波長處出現(xiàn)較弱的自相關(guān)峰,與O—H鍵的一級倍頻有關(guān)。

    圖4 不同波段范圍內(nèi)豬肉樣品的二維相關(guān)同步光譜Fig. 4 2D-correlation spectra of pork samples in different spectral regions

    圖5 不同波段范圍內(nèi)的自相關(guān)譜Fig. 5 Autocorrelation spectra in different spectral regions

    綜合各個波段,可以發(fā)現(xiàn)通過解析二維相關(guān)同步光譜,得到的17 個特征變量多集中在可見波段。這與肉在貯藏放置過程中,由于微生物繁殖產(chǎn)生氫過氧化物,引起肉中色素發(fā)生變化有關(guān)。同時,由于厭氧菌分解蛋白質(zhì)和氨基酸等產(chǎn)生吲哚、甲胺和H2S等物質(zhì),脂肪氧化形成硫代巴比妥酸等物質(zhì),這些均會促進肌紅蛋白的轉(zhuǎn)化和血紅素的氧化[31]。

    2.4 定性判別分析模型建立

    表2 基于全波段光譜的SVM判別分析模型結(jié)果Table 2 Discrimination results of SVM models based on full-band spectra

    首先,基于全波段光譜數(shù)據(jù)建立豬肉新鮮度的SVM判別模型,并比較原始、經(jīng)SNV預(yù)處理和歸一化預(yù)處理光譜的建模效果,如表2所示??梢娀谶@3 類光譜信息,在校正集中分別有6、4、5 個樣品誤判,驗證集中均有1 個樣品誤判,且均為腐敗樣品被誤判,總體判別正確率為87.93%、91.38%和89.66%,同時也可看出預(yù)處理對模型的效果有所改善,預(yù)測總體正確率有所提高。

    然后,利用解析二維相關(guān)光譜得到的17 個特征變量,建立優(yōu)化的判別分析模型?;谒x特征波長,利用原始光譜信息建立的SVM模型,校正集中次新鮮樣品全部判斷正確,新鮮和腐敗樣品各有1和2 個誤判,驗證集樣品全部判斷正確,總體判別正確率為94.83%。利用SNV預(yù)處理后光譜建立的SVM模型,校正集中新鮮和次新鮮樣品全部判斷正確,腐敗樣品有1 個誤判,驗證集樣品全部判斷正確,總體判別正確率為98.28%。利用歸一化預(yù)處理后光譜信息建立的模型效果與SNV預(yù)處理后建模效果一致,總體判別正確率為98.28%,其校正集和驗證集樣品判別結(jié)果如圖6所示,可見除第40號腐敗樣品被誤判為次新鮮外,其余樣品均判斷正確。同時也進一步驗證了光譜的預(yù)處理有助于改善模型,提高模型預(yù)測效果。優(yōu)化后的判別模型,僅利用了總體變量1.61%(原始變量個數(shù)為1 059)的特征變量,取得了更佳的模型結(jié)果。這表明通過二維同步光譜分析選出的17 個變量,能夠反映與豬肉新鮮度密切相關(guān)的特征信息,表征不同TVB-N的光譜變化,簡化判別模型,縮短運算處理時間。

    圖6 基于歸一化預(yù)處理光譜特征變量的判別模型結(jié)果Fig. 6 Discrimination results using feature variables based on spectra after normalization

    3 結(jié) 論

    以TVB-N為外界干擾,通過解析二維相關(guān)同步光譜和自相關(guān)譜,共提取到17 個與新鮮度有關(guān)的特征變量?;谠?、SNV預(yù)處理及歸一化預(yù)處理后的光譜,利用特征波長對應(yīng)光譜建立的SVM模型效果優(yōu)于全波段光譜,總體判別正確率為94.83%、98.28%和98.28%。同時也說明,光譜的預(yù)處理有助于改進模型效果??梢?近紅外光譜結(jié)合二維相關(guān)分析,可以有效的解析肉在腐敗變質(zhì)過程中的光譜特征信息變化,確定與新鮮度相關(guān)的敏感變量。該方法簡單快速,為近紅外光譜特征變量篩選提供了一種新的方法。二維相關(guān)光譜包含信息豐富,選擇的外擾不同得到的譜圖不同,可進一步借助化學計量學方法,對各官能團的變化進行有效解析,提取特征光譜信息。

    猜你喜歡
    光譜信息新鮮度波段
    春日暖陽
    基于光譜和Gabor紋理信息融合的油桃品種識別
    基于植被光譜信息的龜裂堿土堿化程度預(yù)測研究
    綠色科技(2017年20期)2017-11-10 18:54:19
    傅立葉變換光譜儀的研究現(xiàn)狀與光譜信息分析原理
    M87的多波段輻射過程及其能譜擬合
    日常維護對L 波段雷達的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    一種基于區(qū)域生長的高分辨率遙感影像道路提取方法
    基于TTI的冷鮮羊肉新鮮度研究
    不同新鮮度金槍魚肉蒸煮品質(zhì)的研究
    幾種魚體新鮮度快速檢測方法的研究進展
    av免费观看日本| 中国美白少妇内射xxxbb| 日韩电影二区| 亚洲,欧美,日韩| 亚洲成人av在线免费| 国产精品国产三级专区第一集| 婷婷色综合大香蕉| 各种免费的搞黄视频| 成人18禁高潮啪啪吃奶动态图| 日本av手机在线免费观看| 春色校园在线视频观看| 国产一级毛片在线| 婷婷成人精品国产| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区国产| 少妇 在线观看| 18禁动态无遮挡网站| 少妇的逼水好多| 国产成人精品在线电影| 91国产中文字幕| 国产男女内射视频| 国产成人免费无遮挡视频| 一边摸一边做爽爽视频免费| 丝袜喷水一区| 男女无遮挡免费网站观看| 亚洲经典国产精华液单| 亚洲美女黄色视频免费看| 国产男女内射视频| 2021少妇久久久久久久久久久| 在线亚洲精品国产二区图片欧美| 国产不卡av网站在线观看| 亚洲欧美一区二区三区黑人 | 一级,二级,三级黄色视频| 七月丁香在线播放| 久久久久久久精品精品| 少妇人妻精品综合一区二区| 老女人水多毛片| h视频一区二区三区| 欧美激情 高清一区二区三区| 丁香六月天网| 人人妻人人爽人人添夜夜欢视频| 成人毛片a级毛片在线播放| 少妇人妻久久综合中文| 国产视频首页在线观看| 91精品国产国语对白视频| 日韩一区二区三区影片| 久久热在线av| av电影中文网址| 免费在线观看黄色视频的| 香蕉丝袜av| 久久人妻熟女aⅴ| 99热国产这里只有精品6| 制服诱惑二区| 91精品三级在线观看| 久久久久久久久久久久大奶| 亚洲一码二码三码区别大吗| 国产精品 国内视频| 街头女战士在线观看网站| 国产精品无大码| 国产精品秋霞免费鲁丝片| 欧美亚洲 丝袜 人妻 在线| 日本-黄色视频高清免费观看| 成人毛片a级毛片在线播放| 制服诱惑二区| 日韩制服骚丝袜av| 亚洲综合色网址| 国产高清国产精品国产三级| 亚洲 欧美一区二区三区| 亚洲成人手机| 女性被躁到高潮视频| 99热全是精品| 女的被弄到高潮叫床怎么办| 免费黄网站久久成人精品| 男女午夜视频在线观看 | 国产成人精品无人区| 男女边摸边吃奶| 91aial.com中文字幕在线观看| 黄片无遮挡物在线观看| 亚洲美女黄色视频免费看| 草草在线视频免费看| 五月伊人婷婷丁香| 免费高清在线观看视频在线观看| 日本色播在线视频| 一区二区av电影网| 日韩制服骚丝袜av| 国产欧美另类精品又又久久亚洲欧美| 97在线视频观看| 中国国产av一级| 观看美女的网站| 久久综合国产亚洲精品| 一二三四中文在线观看免费高清| 一级a做视频免费观看| 18禁在线无遮挡免费观看视频| 综合色丁香网| 久久99热这里只频精品6学生| av免费观看日本| 欧美人与性动交α欧美精品济南到 | 王馨瑶露胸无遮挡在线观看| 高清欧美精品videossex| 九色亚洲精品在线播放| 国产精品久久久久久久久免| 日韩成人伦理影院| 90打野战视频偷拍视频| 纵有疾风起免费观看全集完整版| 久久99蜜桃精品久久| 中文字幕最新亚洲高清| 丝袜人妻中文字幕| 久久久亚洲精品成人影院| 插逼视频在线观看| 免费大片黄手机在线观看| 日韩在线高清观看一区二区三区| 男女午夜视频在线观看 | 老熟女久久久| 在线观看免费视频网站a站| 捣出白浆h1v1| 国产av码专区亚洲av| 日韩av在线免费看完整版不卡| 国产一区二区在线观看日韩| 精品久久蜜臀av无| 精品少妇久久久久久888优播| 久久精品国产综合久久久 | 国产高清三级在线| 午夜免费观看性视频| 久热久热在线精品观看| 国产在线一区二区三区精| h视频一区二区三区| 看十八女毛片水多多多| 精品福利永久在线观看| 久久久久久久久久久久大奶| 交换朋友夫妻互换小说| 看非洲黑人一级黄片| 亚洲精品色激情综合| 亚洲欧洲日产国产| 国产成人免费无遮挡视频| 九色亚洲精品在线播放| 久久精品国产鲁丝片午夜精品| av视频免费观看在线观看| 性高湖久久久久久久久免费观看| av网站免费在线观看视频| 欧美日韩精品成人综合77777| 丝袜人妻中文字幕| 国产熟女午夜一区二区三区| av免费在线看不卡| 欧美日韩精品成人综合77777| 桃花免费在线播放| 有码 亚洲区| 一区二区三区四区激情视频| 日韩 亚洲 欧美在线| 天天躁夜夜躁狠狠躁躁| 日韩一区二区视频免费看| 热99久久久久精品小说推荐| 亚洲人成网站在线观看播放| 捣出白浆h1v1| 亚洲中文av在线| 人体艺术视频欧美日本| 婷婷色麻豆天堂久久| 国产视频首页在线观看| 春色校园在线视频观看| 黄色 视频免费看| 精品国产一区二区三区四区第35| 日韩,欧美,国产一区二区三区| 捣出白浆h1v1| 日韩精品有码人妻一区| www.色视频.com| 777米奇影视久久| av在线播放精品| 老司机亚洲免费影院| 女的被弄到高潮叫床怎么办| 欧美97在线视频| 熟女av电影| 欧美日本中文国产一区发布| 你懂的网址亚洲精品在线观看| 91aial.com中文字幕在线观看| 久久久国产一区二区| 亚洲色图综合在线观看| 一二三四在线观看免费中文在 | 成人综合一区亚洲| 精品少妇黑人巨大在线播放| 日本午夜av视频| 国产精品免费大片| av天堂久久9| 一级黄片播放器| 久久久国产欧美日韩av| 精品酒店卫生间| 午夜福利网站1000一区二区三区| 一级a做视频免费观看| 美女内射精品一级片tv| 少妇熟女欧美另类| 久久国产精品男人的天堂亚洲 | 伦理电影大哥的女人| 男人添女人高潮全过程视频| 亚洲精华国产精华液的使用体验| 亚洲精品美女久久av网站| 成年人免费黄色播放视频| 午夜福利在线观看免费完整高清在| 视频在线观看一区二区三区| 日本av免费视频播放| 亚洲伊人久久精品综合| 国产男女内射视频| 丰满少妇做爰视频| 久久人人爽人人片av| 国产色婷婷99| 国产av精品麻豆| av又黄又爽大尺度在线免费看| 日日摸夜夜添夜夜爱| 下体分泌物呈黄色| 午夜免费男女啪啪视频观看| 肉色欧美久久久久久久蜜桃| 最黄视频免费看| av一本久久久久| 亚洲色图综合在线观看| 久久久久久久国产电影| 国产乱来视频区| 亚洲成人一二三区av| 国产亚洲精品第一综合不卡 | 人人妻人人澡人人看| 丰满迷人的少妇在线观看| 国产精品一区www在线观看| 99久久人妻综合| 久久久久久久久久成人| 午夜福利网站1000一区二区三区| 精品一区在线观看国产| 亚洲欧美清纯卡通| 国产精品嫩草影院av在线观看| 亚洲,一卡二卡三卡| 亚洲精品国产av蜜桃| 亚洲成av片中文字幕在线观看 | 少妇熟女欧美另类| 成人综合一区亚洲| 久久久久精品人妻al黑| 久久久久国产网址| 日日撸夜夜添| 国产成人免费无遮挡视频| 超碰97精品在线观看| 久久亚洲国产成人精品v| 久久人妻熟女aⅴ| 两个人看的免费小视频| 午夜激情av网站| 一个人免费看片子| 国产精品不卡视频一区二区| 免费看光身美女| 亚洲国产精品一区二区三区在线| 永久免费av网站大全| 午夜av观看不卡| 免费观看性生交大片5| 日韩欧美精品免费久久| 男人爽女人下面视频在线观看| av又黄又爽大尺度在线免费看| 亚洲成人av在线免费| 亚洲成人手机| 丁香六月天网| 人体艺术视频欧美日本| 国产精品秋霞免费鲁丝片| av在线播放精品| 久久久久久久亚洲中文字幕| 久久国产精品大桥未久av| 久久精品久久久久久噜噜老黄| 国产亚洲午夜精品一区二区久久| 最近手机中文字幕大全| 飞空精品影院首页| 少妇的逼水好多| 丝袜人妻中文字幕| 一级a做视频免费观看| 欧美 亚洲 国产 日韩一| 久久久国产欧美日韩av| 夫妻性生交免费视频一级片| 欧美日韩视频精品一区| 亚洲成人手机| 日韩一区二区三区影片| 国产精品久久久av美女十八| 亚洲人成网站在线观看播放| 亚洲av电影在线进入| 欧美人与性动交α欧美精品济南到 | 街头女战士在线观看网站| 欧美少妇被猛烈插入视频| 日韩免费高清中文字幕av| av线在线观看网站| 欧美最新免费一区二区三区| 国产一区二区在线观看日韩| 在线观看美女被高潮喷水网站| 最近的中文字幕免费完整| 久久久国产欧美日韩av| av线在线观看网站| 久久av网站| 欧美日本中文国产一区发布| 亚洲av国产av综合av卡| 777米奇影视久久| 99九九在线精品视频| 飞空精品影院首页| 中国美白少妇内射xxxbb| 母亲3免费完整高清在线观看 | 大片免费播放器 马上看| 美女国产高潮福利片在线看| 久久精品国产亚洲av天美| 久久亚洲国产成人精品v| 人成视频在线观看免费观看| 国产精品国产三级国产av玫瑰| 丝袜美足系列| 午夜精品国产一区二区电影| 免费在线观看完整版高清| 黄片无遮挡物在线观看| 免费av中文字幕在线| 搡老乐熟女国产| 国内精品宾馆在线| 九九爱精品视频在线观看| 国产成人精品在线电影| 搡女人真爽免费视频火全软件| 精品少妇久久久久久888优播| 免费少妇av软件| 天天操日日干夜夜撸| 高清av免费在线| 在线 av 中文字幕| 国产亚洲精品第一综合不卡 | 国产综合精华液| 侵犯人妻中文字幕一二三四区| 久久精品aⅴ一区二区三区四区 | 丁香六月天网| 久久精品人人爽人人爽视色| 最新中文字幕久久久久| 高清av免费在线| 晚上一个人看的免费电影| 看免费成人av毛片| 免费观看a级毛片全部| 91在线精品国自产拍蜜月| 女人被躁到高潮嗷嗷叫费观| 国产精品.久久久| 亚洲欧美成人综合另类久久久| 国产1区2区3区精品| 九九爱精品视频在线观看| 大陆偷拍与自拍| 国产一区二区在线观看av| 嫩草影院入口| 亚洲av在线观看美女高潮| 欧美日韩综合久久久久久| 亚洲欧美一区二区三区黑人 | 亚洲丝袜综合中文字幕| 日本免费在线观看一区| 十分钟在线观看高清视频www| 丝袜美足系列| 韩国精品一区二区三区 | 精品人妻熟女毛片av久久网站| 日产精品乱码卡一卡2卡三| 中文字幕最新亚洲高清| 国产精品久久久久久精品电影小说| 久久久久精品性色| 午夜精品国产一区二区电影| 日韩,欧美,国产一区二区三区| 哪个播放器可以免费观看大片| 色视频在线一区二区三区| 99热这里只有是精品在线观看| 午夜福利,免费看| 久久99蜜桃精品久久| 人人澡人人妻人| 黄色 视频免费看| 国产白丝娇喘喷水9色精品| 中文字幕人妻熟女乱码| 亚洲精品av麻豆狂野| 亚洲精品成人av观看孕妇| 国产不卡av网站在线观看| 香蕉国产在线看| av天堂久久9| 黄网站色视频无遮挡免费观看| 亚洲精华国产精华液的使用体验| 黑丝袜美女国产一区| 国产男女超爽视频在线观看| 搡女人真爽免费视频火全软件| 国产免费一区二区三区四区乱码| 一二三四在线观看免费中文在 | 少妇被粗大的猛进出69影院 | 夫妻性生交免费视频一级片| 婷婷色麻豆天堂久久| av在线播放精品| 91成人精品电影| 中文天堂在线官网| 国产高清国产精品国产三级| 亚洲第一av免费看| av网站免费在线观看视频| 久久热在线av| 国产深夜福利视频在线观看| 精品酒店卫生间| 免费大片18禁| 亚洲图色成人| 天堂中文最新版在线下载| 91久久精品国产一区二区三区| 久久久久精品性色| 午夜视频国产福利| 日本wwww免费看| av不卡在线播放| 美女中出高潮动态图| 极品人妻少妇av视频| 曰老女人黄片| 99国产精品免费福利视频| 国产精品三级大全| 另类精品久久| 国产成人精品在线电影| 久久久久视频综合| 亚洲欧美日韩另类电影网站| 亚洲精品久久成人aⅴ小说| 最近最新中文字幕大全免费视频 | 两个人免费观看高清视频| 久久99一区二区三区| 9191精品国产免费久久| 亚洲精品久久成人aⅴ小说| 只有这里有精品99| 啦啦啦在线观看免费高清www| 成人18禁高潮啪啪吃奶动态图| 免费黄网站久久成人精品| 满18在线观看网站| 国产精品三级大全| 欧美丝袜亚洲另类| 精品一区二区三卡| 色94色欧美一区二区| 插逼视频在线观看| 日韩大片免费观看网站| 91精品国产国语对白视频| av黄色大香蕉| 99国产综合亚洲精品| 午夜久久久在线观看| 成人二区视频| 国产精品.久久久| 人人妻人人澡人人看| 搡老乐熟女国产| 欧美 日韩 精品 国产| 欧美最新免费一区二区三区| 久久人人爽人人爽人人片va| 一边摸一边做爽爽视频免费| 熟女av电影| 亚洲欧美成人综合另类久久久| 婷婷色综合www| 久久av网站| 欧美精品av麻豆av| 热re99久久精品国产66热6| 久久久久国产网址| 亚洲av电影在线观看一区二区三区| 晚上一个人看的免费电影| 日本-黄色视频高清免费观看| 亚洲情色 制服丝袜| 久久精品国产鲁丝片午夜精品| 男的添女的下面高潮视频| 久久ye,这里只有精品| 午夜视频国产福利| 少妇被粗大猛烈的视频| av不卡在线播放| 亚洲成av片中文字幕在线观看 | 午夜精品国产一区二区电影| 亚洲精品aⅴ在线观看| 天美传媒精品一区二区| 国产亚洲精品第一综合不卡 | 国产免费一级a男人的天堂| 狂野欧美激情性xxxx在线观看| 久久久久久久精品精品| 如何舔出高潮| 免费大片黄手机在线观看| 国产激情久久老熟女| 看免费av毛片| 黑人巨大精品欧美一区二区蜜桃 | 最后的刺客免费高清国语| 久久久久国产精品人妻一区二区| 欧美人与善性xxx| 色吧在线观看| 国语对白做爰xxxⅹ性视频网站| 80岁老熟妇乱子伦牲交| 如日韩欧美国产精品一区二区三区| 国产免费一区二区三区四区乱码| 久久毛片免费看一区二区三区| 日本wwww免费看| 丝袜脚勾引网站| 少妇高潮的动态图| 亚洲欧洲日产国产| 久久久久久人妻| av网站免费在线观看视频| 免费大片18禁| 成年av动漫网址| 精品国产露脸久久av麻豆| 成人影院久久| 在线观看人妻少妇| 国产一区二区三区av在线| 9191精品国产免费久久| 国产精品国产av在线观看| 久久这里只有精品19| 99热全是精品| 你懂的网址亚洲精品在线观看| 亚洲综合色惰| 久久精品国产综合久久久 | 女的被弄到高潮叫床怎么办| 免费在线观看完整版高清| 91国产中文字幕| 9191精品国产免费久久| 欧美精品高潮呻吟av久久| 亚洲欧美一区二区三区黑人 | 你懂的网址亚洲精品在线观看| 男女国产视频网站| 考比视频在线观看| av片东京热男人的天堂| 最近2019中文字幕mv第一页| 久久久久久久久久久免费av| 欧美精品国产亚洲| 中文字幕免费在线视频6| 亚洲欧美一区二区三区黑人 | 午夜日本视频在线| 一二三四中文在线观看免费高清| 男女无遮挡免费网站观看| 欧美bdsm另类| 亚洲精品国产av蜜桃| a级毛片在线看网站| 日韩成人av中文字幕在线观看| 久久久久久久国产电影| 免费高清在线观看日韩| 夫妻午夜视频| 久热久热在线精品观看| 高清黄色对白视频在线免费看| 777米奇影视久久| 久久久久久人人人人人| 亚洲精品视频女| 午夜福利影视在线免费观看| 日本黄大片高清| 亚洲精品成人av观看孕妇| 在线观看免费日韩欧美大片| 香蕉精品网在线| 精品久久久久久电影网| 精品人妻一区二区三区麻豆| 成人国语在线视频| 成人综合一区亚洲| 狂野欧美激情性xxxx在线观看| 国产在线视频一区二区| 午夜91福利影院| 天天操日日干夜夜撸| a级毛色黄片| 日韩av在线免费看完整版不卡| 99热6这里只有精品| 亚洲av中文av极速乱| 亚洲精品456在线播放app| 人人妻人人澡人人看| 国产成人精品一,二区| 亚洲精品国产色婷婷电影| 午夜福利,免费看| 美女福利国产在线| 亚洲欧美精品自产自拍| 久久人人爽人人片av| 精品国产乱码久久久久久小说| 日本91视频免费播放| 久久久精品94久久精品| 国产在线视频一区二区| 丰满少妇做爰视频| 欧美精品人与动牲交sv欧美| 妹子高潮喷水视频| 日韩精品有码人妻一区| 亚洲欧美色中文字幕在线| 多毛熟女@视频| 午夜福利,免费看| 久久青草综合色| www.av在线官网国产| 国产av国产精品国产| 精品国产露脸久久av麻豆| 国产精品熟女久久久久浪| 国产毛片在线视频| 黄色怎么调成土黄色| 久久婷婷青草| 午夜激情久久久久久久| 亚洲欧美成人精品一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 高清毛片免费看| 亚洲一区二区三区欧美精品| 亚洲一级一片aⅴ在线观看| 国产成人精品一,二区| 精品久久久久久电影网| 另类精品久久| 亚洲av欧美aⅴ国产| 亚洲精品第二区| 91午夜精品亚洲一区二区三区| 又黄又粗又硬又大视频| 国产一区二区激情短视频 | 黄片无遮挡物在线观看| 中文字幕精品免费在线观看视频 | 丰满迷人的少妇在线观看| 国产高清三级在线| 国产一区二区在线观看av| 一区二区日韩欧美中文字幕 | 男女边吃奶边做爰视频| 老女人水多毛片| 国产国语露脸激情在线看| 免费看av在线观看网站| 欧美日韩精品成人综合77777| 97在线视频观看| 久久久国产欧美日韩av| 亚洲天堂av无毛| 我要看黄色一级片免费的| 美女大奶头黄色视频| 欧美 亚洲 国产 日韩一| 久久女婷五月综合色啪小说| 亚洲 欧美一区二区三区| 日韩精品免费视频一区二区三区 | 黑人猛操日本美女一级片| 91成人精品电影| 一区在线观看完整版| 在线观看人妻少妇| 欧美日韩视频精品一区| 久久久久精品性色| 满18在线观看网站| 日日爽夜夜爽网站| 免费看光身美女| 人人妻人人添人人爽欧美一区卜| 午夜福利乱码中文字幕| 搡老乐熟女国产| 国产xxxxx性猛交| 成人国产av品久久久| 精品熟女少妇av免费看| 嫩草影院入口| 久久久国产精品麻豆| 搡老乐熟女国产| 国产淫语在线视频| 日日撸夜夜添| 亚洲国产精品国产精品| 99热国产这里只有精品6| 久久久久久久大尺度免费视频|