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

    可見(jiàn)/近紅外技術(shù)與感官評(píng)價(jià)信息融合預(yù)測(cè)冷藏條件下羅非魚(yú)TVB-N的研究

    2017-11-21 10:50:50,,,,,,,*
    食品工業(yè)科技 2017年21期
    關(guān)鍵詞:新鮮度羅非魚(yú)尾部

    ,,,, ,,,*

    (1.北京農(nóng)業(yè)信息技術(shù)研究中心,北京 100097; 2.天津科技大學(xué)電子信息與自動(dòng)化學(xué)院,天津 300222; 3.國(guó)家農(nóng)業(yè)信息化工程技術(shù)研究中心,北京 100097; 4.農(nóng)產(chǎn)品質(zhì)量安全追溯技術(shù)及應(yīng)用國(guó)家工程實(shí)驗(yàn)室,北京 100097)

    可見(jiàn)/近紅外技術(shù)與感官評(píng)價(jià)信息融合預(yù)測(cè)冷藏條件下羅非魚(yú)TVB-N的研究

    史策1,3,4,孫立濤1,2,錢(qián)建平1,3,4,韓帥1,3,4,范蓓蕾1,3,4,楊信廷1,3,4,*

    (1.北京農(nóng)業(yè)信息技術(shù)研究中心,北京 100097; 2.天津科技大學(xué)電子信息與自動(dòng)化學(xué)院,天津 300222; 3.國(guó)家農(nóng)業(yè)信息化工程技術(shù)研究中心,北京 100097; 4.農(nóng)產(chǎn)品質(zhì)量安全追溯技術(shù)及應(yīng)用國(guó)家工程實(shí)驗(yàn)室,北京 100097)

    應(yīng)用可見(jiàn)/近紅外光譜技術(shù)與感官評(píng)價(jià)信息融合實(shí)現(xiàn)冷藏條件下有鱗羅非魚(yú)不同部位魚(yú)肉揮發(fā)性鹽基氮(TVB-N)的預(yù)測(cè)。通過(guò)便攜式近紅外光譜儀采集有鱗羅非魚(yú)胸部、中部和尾部魚(yú)肉在340~1063 nm的光譜數(shù)據(jù),分別采用卷積平滑法、變量標(biāo)準(zhǔn)化、一階(1st Der)和二階(2nd Der)導(dǎo)數(shù)進(jìn)行光譜預(yù)處理,利用連續(xù)投影算法(SPA)提取羅非魚(yú)不同部位魚(yú)肉的特征波長(zhǎng),建立魚(yú)肉光譜與TVB-N偏最小二乘回歸(PLSR)模型,結(jié)果表明尾部魚(yú)肉的1st Der-SPA-PLSR模型預(yù)測(cè)均方根誤差(RMSEP)=1.1295 mg/100 g,預(yù)測(cè)相關(guān)系數(shù)(Rp2)=0.8998,預(yù)測(cè)結(jié)果高于其胸部和中部魚(yú)肉模型,并稍高于尾部魚(yú)肉全波段模型。因此,選擇尾部魚(yú)肉作為羅非魚(yú)光譜采樣區(qū)域。為進(jìn)一步提高模型預(yù)測(cè)準(zhǔn)確性,將尾部魚(yú)肉特征光譜數(shù)據(jù)與感官評(píng)價(jià)進(jìn)行信息融合。通過(guò)對(duì)比尾部魚(yú)肉光譜、尾部魚(yú)肉光譜與感官評(píng)價(jià)融合的SPA-PLSR、SPA-BP神經(jīng)網(wǎng)絡(luò)和SPA-偏最小二乘支持向量機(jī)(LS-SVM)模型,結(jié)果表明,尾部光譜和感官評(píng)價(jià)信息融合的SPA-LS-SVM模型預(yù)測(cè)結(jié)果為RMSEP=0.9701 mg/100 g,Rp2=0.9255,能更準(zhǔn)確預(yù)測(cè)羅非魚(yú)冷藏條件下TVB-N變化。為冷藏過(guò)程中羅非魚(yú)新鮮度預(yù)測(cè)方法提供了新的思路。

    可見(jiàn)/近紅外光譜技術(shù),羅非魚(yú),新鮮度,預(yù)測(cè)模型,貯藏,信息融合

    羅非魚(yú)(Oreochromis),又名非洲鯽魚(yú),自2000年以來(lái)我國(guó)羅非魚(yú)養(yǎng)殖產(chǎn)量發(fā)展迅猛,現(xiàn)已成為世界上最大的羅非魚(yú)養(yǎng)殖國(guó)家[1]。羅非魚(yú)含有豐富的蛋白質(zhì)和多種不飽和脂肪酸,有望成為未來(lái)動(dòng)物性蛋白質(zhì)的主要來(lái)源之一。然而由于富含蛋白質(zhì)和水分,導(dǎo)致貯藏過(guò)程中極易腐敗變質(zhì)。目前主要采用理化分析、微生物檢測(cè)來(lái)鑒別魚(yú)類(lèi)的新鮮度,而傳統(tǒng)的理化、微生物等分析方法費(fèi)時(shí)費(fèi)力、具有破壞性,因此,迫切需要一種無(wú)損的新鮮度檢測(cè)技術(shù)對(duì)魚(yú)類(lèi)進(jìn)行品質(zhì)評(píng)價(jià)。

    表1 羅非魚(yú)感官評(píng)價(jià)標(biāo)準(zhǔn)Table 1 The criterion of sensory index for tilapia

    可見(jiàn)/近紅外光譜技術(shù)是近年來(lái)發(fā)展起來(lái)的一種高效、快速、無(wú)損的檢測(cè)技術(shù)[2],目前廣泛的應(yīng)用于食品領(lǐng)域。陳偉華等[3]采集絞碎前后羅非魚(yú)片背肉及腹肉的近紅外光譜,并將其與揮發(fā)性鹽基氮(total volatile basic nitrogen,TVB-N)含量進(jìn)行擬合,能夠較好的預(yù)測(cè)TVB-N的變化;Nilsen等[4]利用近紅外光譜研究冰藏鱈魚(yú)和三文魚(yú)的冷藏時(shí)間,應(yīng)用多元統(tǒng)計(jì)分析光譜數(shù)據(jù)與冷藏時(shí)間的相關(guān)性,相關(guān)系數(shù)分別為0.97和0.98;Sivertsen等[5]采用可見(jiàn)/近紅外光譜檢測(cè)了解凍鱈魚(yú)在冰藏條件下的存儲(chǔ)時(shí)間,通過(guò)偏最小二乘法對(duì)光譜數(shù)據(jù)和存儲(chǔ)時(shí)間建模,預(yù)測(cè)鱈魚(yú)的貯藏時(shí)間精度在1.6 d左右。目前大部分的研究集中在近紅外光譜與冷藏時(shí)間或TVB-N等之間的相關(guān)性。然而僅憑光譜信息單一指標(biāo)判定魚(yú)類(lèi)冷藏條件下的新鮮度或貨架期,其準(zhǔn)確性有限。因此,將光譜信息與感官評(píng)價(jià)等內(nèi)外品質(zhì)信息融合用于魚(yú)新鮮度無(wú)損檢測(cè),有助于綜合評(píng)價(jià)魚(yú)類(lèi)的新鮮度。

    本研究提出了基于可見(jiàn)/近紅外光譜和感官評(píng)價(jià)信息融合預(yù)測(cè)羅非魚(yú)新鮮度的研究思路。主要采用商業(yè)化便攜式手持近紅外光譜儀,分別測(cè)定并對(duì)比了冷藏條件下有鱗羅非魚(yú)胸部、中部和尾部魚(yú)肉的光譜數(shù)據(jù)與TVB-N的相關(guān)性,選擇較優(yōu)的羅非魚(yú)新鮮度檢測(cè)部位;其次采用光譜數(shù)據(jù)和感官評(píng)價(jià)進(jìn)行信息融合,建立PLSR,BP-ANN和LS-SVM預(yù)測(cè)模型,以期獲得更好的羅非魚(yú)TVB-N預(yù)測(cè)結(jié)果。為進(jìn)一步提高冷藏過(guò)程中羅非魚(yú)新鮮度預(yù)測(cè)準(zhǔn)確率提供依據(jù)。

    1 材料與方法

    1.1材料與儀器

    60條羅非魚(yú)(657±67) g 北京農(nóng)貿(mào)市場(chǎng);氧化鎂,硼酸,甲基紅,次甲基藍(lán)(均為分析純)北京化學(xué)試劑公司。

    D2010W電動(dòng)攪拌器 上海梅穎浦儀器儀表制造有限公司;YP30002越平電子天平 上海右一儀器有限公司;KDY-9820凱氏定氮儀 北京通潤(rùn)源機(jī)電技術(shù)公司;便攜式可見(jiàn)近紅外光譜儀 美國(guó)ASD公司,采集光譜波段范圍是325~1075 nm,光譜采樣間隔為0.81 nm,掃描次數(shù)20次。

    1.2實(shí)驗(yàn)方法

    1.2.1 樣品前處理 鮮活羅非魚(yú)擊斃后,裝入聚乙烯袋中并密封,置于4 ℃培養(yǎng)箱中貯藏。在貯藏0、2、4、6、8、10 d時(shí)取樣,每次取10條羅非魚(yú)用于實(shí)驗(yàn)。每條魚(yú)左右兩側(cè)作為2個(gè)樣本,共120個(gè)樣本。其中80個(gè)樣本作為建模集,40個(gè)樣本作為預(yù)測(cè)集。

    1.2.2 感官評(píng)價(jià) 參照李杉等[6]和史策等[7]關(guān)于魚(yú)體感官指標(biāo)的測(cè)定方法,并做部分修改。經(jīng)過(guò)感官培訓(xùn)的7位食品專(zhuān)業(yè)同學(xué)作為評(píng)審成員,在光線(xiàn)充足無(wú)氣味的環(huán)境中,從魚(yú)體的眼球、體表色澤、肌肉彈性和氣味4個(gè)方面分別進(jìn)行感官評(píng)分,每個(gè)部分滿(mǎn)分為25分。以60分為食用標(biāo)準(zhǔn),100分為新鮮魚(yú),低于60分為不可食用。

    1.2.3 揮發(fā)性鹽基氮(TVB-N) 參照董偉坤[8]的方法,按照半微量蒸餾法進(jìn)行測(cè)定。

    1.2.4 光譜的采集 將羅非魚(yú)體表面水分擦干后置于黑色背景布上,以標(biāo)準(zhǔn)白板(聚四氟乙烯)為參比,以漫反射的方式采集光譜,采集區(qū)域如圖1標(biāo)識(shí)所示。分別采集羅非魚(yú)胸部、中部和尾部各20次(圖1),每個(gè)部位采集區(qū)域不重疊。取每個(gè)樣品得到的20條光譜曲線(xiàn)的平均值。

    圖1 羅非魚(yú)近紅外光譜采集部位Fig.1 Spectra extraction regions of tilapia

    1.3數(shù)據(jù)分析

    光譜預(yù)處理:分別采用Savitzky-Golay(SG)卷積平滑法、變量標(biāo)準(zhǔn)化(Standard Normal Variate,SNV)、一階(1stDer)和二階(2ndDer)求導(dǎo)算法對(duì)光譜進(jìn)行預(yù)處理分析;

    特征波長(zhǎng)選擇:采用連續(xù)投影算法(Successive projections algorithm,SPA)進(jìn)行特征波長(zhǎng)選擇;

    建模方法:采用BP人工神經(jīng)網(wǎng)絡(luò)(Back-propagation Artificial Neural Network,BP-ANN)、偏最小二乘回歸法(Partial Least Squares Regression,PLSR)和最小二乘支持向量機(jī)(Least Squares Support Vector Machine,LS-SVM)進(jìn)行建模。

    實(shí)驗(yàn)數(shù)據(jù)(均方根誤差(RMSE)、平方相關(guān)系數(shù)(R2)等)使用Excel 2007進(jìn)行數(shù)據(jù)處理,采用SAS軟件進(jìn)行方差分析和差異顯著性分析,算法主要采用Unscrambler(CAMO,Process,AS,OSLO,Norway)、ASD View Spec Pro、MATLAB 2014a(The Math Works,Natick,USA)實(shí)現(xiàn)。

    2 結(jié)果與分析

    2.1冷藏過(guò)程中羅非魚(yú)感官值和TVB-N的變化

    圖2對(duì)羅非魚(yú)的眼球、體表色澤、肌肉彈性和氣味4個(gè)方面分別進(jìn)行了評(píng)分,隨著貯藏時(shí)間的延長(zhǎng)感官值逐漸下降。0~6 d時(shí)羅非魚(yú)感官值下降速率較慢,6 d后,羅非魚(yú)感官值下降速率顯著增加。當(dāng)貯藏至第6 d時(shí),對(duì)應(yīng)的感官值為68,6 d后,魚(yú)肉逐漸失去使用價(jià)值。李莎等[9]研究發(fā)現(xiàn)4 ℃貯藏5 d后羅非魚(yú)片品質(zhì)明顯下降;Hong等[10]研究4 ℃冷藏下鳙魚(yú)的感官值發(fā)現(xiàn),8 d后鳙魚(yú)不可食用。以上報(bào)道都與本研究結(jié)果類(lèi)似。

    圖2 羅非魚(yú)4 ℃貯藏過(guò)程中感官值和TVB-N的變化Fig.2 The change of sensory evaluation and TVB-N value of tilapia during storage at 4 ℃

    TVB-N是酶和細(xì)菌分解蛋白質(zhì)產(chǎn)生的氨及胺類(lèi)等堿性含氮物質(zhì)[11],目前被廣泛用于判斷水產(chǎn)品及肉品的新鮮度。其中,一級(jí)鮮度淡水魚(yú)為T(mén)VB-N≤13 mg/100 g,二級(jí)鮮度淡水魚(yú)為T(mén)VB-N≤20 mg/100 g。如圖2所示,新鮮羅非魚(yú)的TVB-N含量為7.6 mg/100 g,隨著貯藏時(shí)間的延長(zhǎng)TVB-N含量逐漸增加,在第6 d達(dá)到感官拒絕值時(shí)的TVB-N含量為17.2 mg/100 g。

    本實(shí)驗(yàn)所采集的羅非魚(yú)感官評(píng)價(jià)和TVB-N的范圍分別為12.54~100和7.6~23.5 mg/100 g,兩種指標(biāo)數(shù)據(jù)范圍均較大,可滿(mǎn)足光譜技術(shù)定量建模的需求。

    2.2冷藏過(guò)程中羅非魚(yú)胸部、中部和尾部的可見(jiàn)/近紅外原始光譜曲線(xiàn)的變化

    由于前后波段噪聲較大,最終選取340~1063 nm波長(zhǎng)范圍的吸光值作為羅非魚(yú)原始光譜曲線(xiàn)。將原始光譜平均后,得到羅非魚(yú)胸部、中部和尾部平均后的光譜曲線(xiàn)(如圖3)。三個(gè)部位羅非魚(yú)肉的光譜峰谷變化趨勢(shì)基本一致,與zhu等[12]測(cè)定的比目魚(yú)片近紅外光譜谷峰趨勢(shì)類(lèi)似。413 nm處光譜吸光度都達(dá)到最高值;3個(gè)部位的魚(yú)肉光譜在452~628 nm波段都出現(xiàn)峰谷波動(dòng),552 nm附近的波峰可能是由于血紅蛋白(Hb)和肌紅蛋白(Mb)分子上血紅素的吸收,512 nm和620 nm處的波谷可能是由于高鐵血紅蛋白(metHb)和高鐵肌紅蛋白(metMb)分子上氧化血紅素的吸收[5];3個(gè)部位魚(yú)肉的光譜曲線(xiàn)在977 nm處的波峰可能是由于水分子中O-H鍵的二級(jí)倍頻吸收帶在970 nm附近有特征吸收[13]。

    圖3 羅非魚(yú)胸部、中部和尾部的平均光譜曲線(xiàn)Fig.3 Average spectra of three regions for tilapia

    2.3不同預(yù)處理方法對(duì)羅非魚(yú)胸部、中部和尾部建模的影響

    光譜除樣品本身的信息外,還會(huì)受到魚(yú)體表面不均勻、紋理等物理因素,甚至背景噪聲等無(wú)關(guān)信息的干擾,因此需要對(duì)光譜進(jìn)行預(yù)處理[2]。SG平滑不僅能夠減少噪聲對(duì)該數(shù)據(jù)點(diǎn)數(shù)值的干擾,還能保持?jǐn)?shù)據(jù)信號(hào)的真實(shí)性[14],SNV處理通常用于校正樣本間由散射引起的光譜誤差[14],1st Der和2nd Der常用于近紅外光譜的預(yù)處理,能夠一定程度上解決光譜信號(hào)中由于倍頻和合頻引起的光譜信號(hào)重疊問(wèn)題[15]。本實(shí)驗(yàn)主要比較了原始光譜(Raw)、SG、SNV、1st Der和2nd Der以及不同方法組合對(duì)建模結(jié)果的影響。

    表2 4 ℃冷藏過(guò)程中經(jīng)過(guò)不同預(yù)處理的羅非魚(yú)胸部、中部和尾部魚(yú)肉光譜與TVB-N的PLS建模結(jié)果Table 2 The result of PLSR models for correlating between spectral measurement of different positions for tilapia by different preprocessing methods and TVB-N

    經(jīng)過(guò)不同預(yù)處理的胸部、中部和尾部魚(yú)肉光譜與TVB-N建立PLSR的模型,如表2所示。對(duì)比胸部、中部和尾部魚(yú)肉的校正集、交互驗(yàn)證集和預(yù)測(cè)集結(jié)果,其中胸部魚(yú)肉經(jīng)過(guò)SNV和1st Der處理的建模結(jié)果總體較優(yōu),中部魚(yú)肉經(jīng)過(guò)2nd Der處理的光譜建模結(jié)果總體較優(yōu),尾部魚(yú)肉經(jīng)過(guò)1st Der和2nd Der處理的光譜建模結(jié)果總體較優(yōu),這可能由于導(dǎo)數(shù)能夠消除光譜基線(xiàn)漂移及魚(yú)表面不平等背景干擾。胸部魚(yú)肉經(jīng)過(guò)1st Der處理的模型校正集均方根誤差(RMSEP=1.1694 mg/100 g)大于其SNV結(jié)果(RMSEP=0.8804 mg/100 g),但交互驗(yàn)證集和預(yù)測(cè)集均方根誤差(RMSEP=1.6437 mg/100 g和1.5092 mg/100 g)小于其SNV處理結(jié)果(RMSEP=1.6837 mg/100 g和1.6788 mg/100 g),由于SNV處理后胸部魚(yú)肉的校正集與交互驗(yàn)證集和預(yù)測(cè)集的均方根誤差及相關(guān)系數(shù)差距較大,可能出現(xiàn)過(guò)擬合情況,因此選擇1st Der處理的胸部魚(yú)肉數(shù)據(jù)進(jìn)行下一步建模。尾部魚(yú)肉光譜經(jīng)過(guò)1st Der處理的模型校正集和交互驗(yàn)證集均方根誤差分別與其2nd Der結(jié)果相近,但尾部魚(yú)肉光譜經(jīng)過(guò)1st Der預(yù)處理的模型預(yù)測(cè)集均方根誤差(1.2293 mg/100 g)小于其2nd Der處理的結(jié)果(1.3998 mg/100 g),因此尾部魚(yú)肉光譜經(jīng)過(guò)1st Der處理的建模結(jié)果總體較優(yōu);而羅非魚(yú)中部建模結(jié)果中,經(jīng)過(guò)2nd Der處理的光譜預(yù)測(cè)集效果較優(yōu),RMSEP=1.0345 mg/100 g,Rp2=0.9088。通過(guò)對(duì)比預(yù)處理后的胸部、中部和尾部魚(yú)肉建模效果發(fā)現(xiàn),胸部和尾部魚(yú)肉經(jīng)過(guò)1st Der處理、中部經(jīng)2nd Der處理的全波段光譜建模結(jié)果總體較優(yōu)。

    2.4基于連續(xù)投影法(SPA)的胸部、中部和尾部魚(yú)肉特征波長(zhǎng)選取及TVB-N預(yù)測(cè)

    為了減少不必要和多余的信息,并優(yōu)化檢測(cè)速度,本文選擇SPA進(jìn)行最優(yōu)波長(zhǎng)的選擇。SPA被廣泛的用于變量的選擇,能最大程度的減少光譜信息的重疊[16],簡(jiǎn)化校正模型和縮短校正時(shí)間。圖4為羅非魚(yú)胸部1st Der、中部2nd Der和尾部1st Der預(yù)處理光譜經(jīng)過(guò)SPA選擇的過(guò)程,胸部、中部和尾部魚(yú)肉分別優(yōu)選出15、15和19個(gè)光譜變量。羅非魚(yú)胸部魚(yú)肉光譜變量主要集中在347~366 nm和976~1062 nm,中部魚(yú)肉光譜變量主要集中在502~652 nm和968~1039 nm,尾部魚(yú)肉光譜變量主要集中在347~369 nm和935~1061 nm。3個(gè)部位選擇的特征波段都集中在900~1060 nm之間,主要由于這部分波段是水分吸收帶,而水分是魚(yú)肉最主要的成分[16]。

    圖4 羅非魚(yú)胸部、中部和尾部魚(yú)肉光譜SPA法選擇變量過(guò)程Fig.4 Selection of optimum spectral variables of different tilapia positions by SPA

    基于SPA算法選擇的特征波長(zhǎng)分別建立了胸部、中部和尾部魚(yú)肉SPA-PLSR模型(表3)。通過(guò)對(duì)比不同部位建模結(jié)果發(fā)現(xiàn),尾部魚(yú)肉SPA-PLSR模型的校正集、交互驗(yàn)證集和預(yù)測(cè)集的相關(guān)系數(shù)(R2)顯著高于胸部和中部。與尾部魚(yú)肉全波段建模結(jié)果相比,經(jīng)過(guò)SPA算法預(yù)處理后的尾部魚(yú)肉PLSR模型RMSEP=1.1295 mg/100 g,Rp2=0.8998,較全波段建模結(jié)果有所提高(RMSEP=1.2293 mg/100 g,Rp2=0.8698),且尾部魚(yú)肉SPA-PLSR模型僅使用19個(gè)變量,建模的速度和效率得到提升。因此,羅非魚(yú)尾部魚(yú)肉光譜經(jīng)過(guò)1st Der、SPA算法預(yù)處理后能夠良好的預(yù)測(cè)4 ℃貯藏條件下羅非魚(yú)TVB-N的變化。

    2.5羅非魚(yú)尾部魚(yú)肉光譜與感官值信息融合的預(yù)測(cè)模型

    表3 經(jīng)過(guò)連續(xù)投影法選擇特征波長(zhǎng)后的PLSR建模結(jié)果Table 3 The result of PLSR models of different positions for tilapia by SPA

    表4 PLSR,BP-ANN和LS-SVM模型結(jié)合多種變量對(duì)TVB-N建模與預(yù)測(cè)結(jié)果的參數(shù)統(tǒng)計(jì)Table 4 Comparison of storage time with different data variables based on PLSR,BP-ANN and LS-SVM

    信息融合根據(jù)融合系統(tǒng)所處理的信息層面主要分為原始數(shù)據(jù)層融合、特征層融合和決策層融合。在本文中,由于光譜和感官評(píng)價(jià)特性不一致,較難直接關(guān)聯(lián);而決策層融合需要與光譜及感官相關(guān)的先驗(yàn)知識(shí)進(jìn)行決策,所需數(shù)據(jù)庫(kù)龐大,因此,本文選擇特征層融合。通過(guò)將光譜數(shù)據(jù)進(jìn)行特征處理后,與感官數(shù)據(jù)聯(lián)合進(jìn)行特征層融合完成對(duì)羅非魚(yú)新鮮度的綜合評(píng)價(jià)。

    為了提高羅非魚(yú)貯藏過(guò)程TVB-N的預(yù)測(cè)精度,將上文中經(jīng)過(guò)1st Der、SPA預(yù)處理的尾部魚(yú)肉光譜特征波長(zhǎng)與感官評(píng)價(jià)聯(lián)合進(jìn)行信息融合,采用PLSR,BP-ANN和LS-SVM建立TVB-N預(yù)測(cè)模型(表4)。

    BP-ANN是一種多層前饋神經(jīng)網(wǎng)絡(luò),從輸出層反向經(jīng)過(guò)隱含層回到輸入層,逐步修正各連接權(quán)值。本文中BP-ANN模型選取的是單隱含層模型,隱含層和輸出層的傳遞函數(shù)分別選取logsig和purelin,經(jīng)過(guò)10000次訓(xùn)練后,光譜、光譜和感官值融合的BP-ANN模型的隱含層的神經(jīng)元分別為10、11時(shí),逼近效果最好。LS-SVM是一種遵循結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則的核函數(shù)學(xué)習(xí)機(jī)器。利用非線(xiàn)性映射函數(shù)將輸入向量映射到高維特征空間,采用拉格朗日乘子計(jì)算向量的偏微分,從而以較快的方法解決多元校正問(wèn)題[17-18],能夠降低計(jì)算機(jī)的復(fù)雜性,提高計(jì)算效率[14]。本文中LS-SVM模型采用的徑向基核函數(shù),并通過(guò)交叉驗(yàn)證優(yōu)化參數(shù)建立LS-SVM模型。

    對(duì)比3種建模方法,光譜和感官信息融合的LS-SVM預(yù)測(cè)模型的結(jié)果總體略?xún)?yōu)。這可能是由于變量較多并且變量之間線(xiàn)性關(guān)系并不明顯,而PLSR模型只能處理光譜數(shù)據(jù)中的線(xiàn)性關(guān)系;在2組BP-ANN模型中RMSEP顯著大于RMSECV,在運(yùn)行過(guò)程中可能出現(xiàn)過(guò)擬合情況;而LS-SVM模型可以利用光譜數(shù)據(jù)的非線(xiàn)性信息,提高了模型預(yù)測(cè)的準(zhǔn)確性[19]。2組LS-SVM模型的交互驗(yàn)證集結(jié)果接近,而光譜和感官值融合模型的預(yù)測(cè)集RMSEP=0.9701 mg/100 g優(yōu)于光譜模型的預(yù)測(cè)集RMSEP=1.2364 mg/100 g,光譜和感官值融合的LS-SVM模型能較好的預(yù)測(cè)羅非魚(yú)的TVB-N含量。羅非魚(yú)的品質(zhì)劣變是個(gè)復(fù)雜的過(guò)程,伴隨著內(nèi)部指標(biāo)和外部指標(biāo)變化,可見(jiàn)/近紅外光譜數(shù)據(jù)反映了羅非魚(yú)內(nèi)部化學(xué)變化,感官評(píng)價(jià)反映了羅非魚(yú)外部感官變化。將光譜數(shù)據(jù)和感官評(píng)價(jià)信息融合的LS-SVM模型能更全面的反映羅非魚(yú)冷藏條件下新鮮度變化。

    3 結(jié)論

    通過(guò)商業(yè)化便攜式手持近紅外光譜儀采集不同部位(胸部、中部和尾部)羅非魚(yú)光譜與TVB-N進(jìn)行擬合,構(gòu)建預(yù)測(cè)模型以評(píng)價(jià)冷藏過(guò)程中新鮮度。通過(guò)SG、SNV、1st Der和2nd Der以及不同方法組合的光譜預(yù)處理及SPA算法的優(yōu)化,結(jié)果表明羅非魚(yú)尾部1st Der-SPA-PLSR模型預(yù)測(cè)效果好于其胸部和中部預(yù)測(cè)模型,并較尾部全波段建模結(jié)果有所提高。因此,選取尾部魚(yú)肉作為羅非魚(yú)光譜采樣區(qū)域。

    為提高羅非魚(yú)冷藏過(guò)程中TVB-N的預(yù)測(cè)準(zhǔn)確性,將經(jīng)過(guò)1st Der、SPA預(yù)處理的尾部魚(yú)肉光譜數(shù)據(jù)與感官評(píng)價(jià)特征層信息融合,對(duì)比光譜、光譜和感官值信息融合的PLSR,BP-ANN和LS-SVM模型效果。尾部魚(yú)肉光譜和感官評(píng)價(jià)信息融合的SPA-LS-SVM模型略?xún)?yōu),預(yù)測(cè)集的RMSEP=0.9701 mg/100 g,Rp2=0.9255。通過(guò)尾部魚(yú)肉光譜和感官評(píng)價(jià)信息融合能更全面的反映羅非魚(yú)冷藏條件下新鮮度變化,為提高冷藏過(guò)程中羅非魚(yú)新鮮度檢測(cè)準(zhǔn)確率提供理論依據(jù)。

    [1]楊弘. 我國(guó)羅非魚(yú)產(chǎn)業(yè)現(xiàn)狀及產(chǎn)業(yè)技術(shù)體系建設(shè)[J]. 中國(guó)水產(chǎn),2010(9):6-10.

    [2]黃濤,李小昱,彭毅,等. 基于近紅外光譜的淡水魚(yú)新鮮度在線(xiàn)檢測(cè)方法研究[J]. 光譜學(xué)與光譜分析,2014(10):2732-2736.

    [3]陳偉華,許長(zhǎng)華,樊玉霞,等. 近紅外光譜技術(shù)快速無(wú)損評(píng)價(jià)羅非魚(yú)片新鮮度[J]. 食品科學(xué),2014,35(24):164-168.

    [4]Nilsen H,Esaiassen M,Heia K,et al. Visible/Near-Infrared Spectroscopy:A New Tool for the Evaluation of Fish Freshness?[J]. Journal of Food Science,2002,67(5):1821-1826.

    [5]Sivertsen A H,Kimiya T,Heia K. Automatic freshness assessment of cod(Gadus morhua)fillets by Vis/Nir spectroscopy[J]. Journal of Food Engineering,2011,103(3):317-323.

    [6]李杉,馬海霞,李來(lái)好,等. 減菌化預(yù)處理對(duì)鮮羅非魚(yú)片質(zhì)量的影響[J]. 食品科學(xué),2009(18):379-384.

    [7]史策,崔建云,王航,等. 反復(fù)冷凍-解凍對(duì)鰱品質(zhì)的影響[J]. 中國(guó)水產(chǎn)科學(xué),2012,19(1):167-173.

    [8]董偉坤. 食品檢驗(yàn)與分析[Z]. 北京:中國(guó)輕工業(yè)出版社,1989.

    [9]李莎,李來(lái)好,楊賢慶,等. 羅非魚(yú)片在冷藏過(guò)程中的品質(zhì)變化研究[J]. 食品科學(xué),2010,31(20):444-447.

    [10]Hong H,Luo Y,Zhou Z,et al. Effects of low concentration of salt and sucrose on the quality of bighead carp(Aristichthysnobilis)fillets stored at 4 C[J]. Food Chemistry,2012,133(1):102-107.

    [11]Olafsdóttir G,Martinsdóttir E,Oehlenschl?ger J,et al. Methods to evaluate fish freshness in research and industry[J]. Trends in Food Science & Technology,1997,8(8):258-265.

    [12]Zhu F,Zhang D,He Y,et al. Application of visible and near infrared hyperspectral imaging to differentiate between fresh and frozen-thawed fish fillets[J]. Food and Bioprocess Technology,2013,6(10):2931-2937.

    [13]B. G. O,T. F. Near Infrared Spectroscopy in Food Analysis[M].New York:Longman Scientific & Technical,1986.

    [14]朱逢樂(lè).基于光譜和高光譜成像技術(shù)的海水魚(yú)品質(zhì)快速無(wú)損檢測(cè)[D]. 浙江大學(xué)生物系統(tǒng)工程,2014.

    [15]Lin M,Cavinato A G,Huang Y,et al. Predicting sodium chloride content in commercial king(Oncorhynchustshawytscha)and chum(O.keta)hot smoked salmon fillet portions by short-wavelength near-infrared(SW-NIR)spectroscopy[J]. Food Research International,2003,36(8):761-766.

    [16]Cheng J,Sun D,Zeng X,et al. Non-destructive and rapid determination of TVB-N content for freshness evaluation of grass carp(Ctenopharyngodonidella)by hyperspectral imaging[J]. Innovative Food Science & Emerging Technologies,2014,21:179-187.

    [17]顧燕萍,趙文杰,吳占松. 基于最小二乘支持向量機(jī)的電站鍋爐燃燒優(yōu)化[J]. 中國(guó)電機(jī)工程學(xué)報(bào),2010,30(17):91-97.

    [18]Suykens J A,Vandewalle J. Least squares support vector machine classifiers[J]. Neural processing letters,1999,9(3):293-300.

    [19]朱逢樂(lè),何勇,邵詠妮. 應(yīng)用近紅外高光譜成像預(yù)測(cè)三文魚(yú)肉的水分含量[J]. 光譜學(xué)與光譜分析,2015(1):113-117.

    TVB-Npredictionoftilapiawithscalesbyinformationfusionofnearinfraredspectrumtechnologyandsensoryevaluationduringchilledstorage

    SHICe1,3,4,SUNLi-tao1,2,QIANJian-ping1,3,4,HANShuai1,3,4,FANBei-lei1,3,4,YANGXin-ting1,3,4,*

    (1.Beijing Research Center for Information Technology in Agriculture,Beijing 100097,China; 2.College of Electronic Information and Automation,Tianjin University of Science & Technology,Tianjin 300222,China; 3.National Engineering Research Center for Information Technology in Agriculture,Beijing 100097,China; 4.National Engineering Laboratory for Agri-product Quality Traceability,Beijing 100097,China)

    The information fusion of near infrared spectrum technology and sensory evaluation was applied to predict the freshness of different parts for tilapia with scales during chilled storage.Spectral signatures ofbreast,middle and tail region in the range of 340~1063 nm were extracted.Smoothing Savitzky-Golay(SG),standard normal variate(SNV),polynomial derivative filters(1st Der and 2nd Der)were used for spectral pre-processing.Partial least square regression(PLSR)was used to correlate the whole wavelengths spectra with total volatile basic nitrogen(TVB-N).Optimal wavelengths of different tilapia positions were selected by successive projections algorithm(SPA)to develop new SPA-PLSR models,and the SPA-PLSR predictive performances of tails position(root mean square error of prediction(RMSEP)=1.1295 mg/100 g,determination coefficient(Rp2)=0.8998)was better than that of breast and middle region,and also better than whole wavelengths model of tails region.Therefore,tail region was selected as spectrum sampling area. In order to evaluate the comprehensively fish freshness and improve the accuracy of model,spectral data and sensory evaluation were integrated for nondestructive measurement of freshness for tail region of tilapia based on PLSR,back-propagation artificial neural network(BP-ANN)and least squares support vector machines(LS-SVM). Compared with single characteristic,information fusion of spectral data and sensory evaluation for LS-SVM had its superiority,which achieved accurate results with Rp2of 0.9255,RMSEP of 0.9701 mg/100 g. This result indicated that information fusion by integrating spectral data and sensory evaluation could significantly improve the TVB-N prediction performance,and it has tremendous potential in prediction of freshness in fish during chilled storage.

    near infrared spectrum technology;tilapia;freshness;predictive model;storage;information fusion

    2017-03-03

    史策(1989-),女,博士,助理研究員,研究方向:農(nóng)產(chǎn)品質(zhì)量安全控制技術(shù),E-mail:shic@nercita.org.cn。

    *

    楊信廷(1974-),男,博士,研究員,主要從事農(nóng)業(yè)信息化關(guān)鍵技術(shù)及農(nóng)產(chǎn)品質(zhì)量安全控制研究,E-mail:yangxt@nercita.org.cn。

    北京市自然科學(xué)基金項(xiàng)目(6174040);北京市農(nóng)林科學(xué)院青年科研基金(QNJJ201720);國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFD0401205)。

    TS254.1

    A

    1002-0306(2017)21-0268-06

    10.13386/j.issn1002-0306.2017.21.053

    猜你喜歡
    新鮮度羅非魚(yú)尾部
    船舶尾部響應(yīng)特性試驗(yàn)與計(jì)算
    超聲及磁共振診斷骶尾部藏毛竇1例
    羅非魚(yú)養(yǎng)殖模式
    貿(mào)易戰(zhàn),羅非魚(yú)首當(dāng)其沖!面臨眾多不利因素,昔日的王者羅非魚(yú)還能打一場(chǎng)翻身戰(zhàn)嗎?
    羅非魚(yú) 年總產(chǎn)量全國(guó)第三位
    彎式尾部接頭注塑模具設(shè)計(jì)
    基于TTI的冷鮮羊肉新鮮度研究
    不同新鮮度金槍魚(yú)肉蒸煮品質(zhì)的研究
    幾種魚(yú)體新鮮度快速檢測(cè)方法的研究進(jìn)展
    羅非魚(yú)肥胖癥的防治
    国产淫片久久久久久久久| 日韩人妻高清精品专区| xxx大片免费视频| 国产精品久久久久久av不卡| 晚上一个人看的免费电影| 日韩成人伦理影院| 身体一侧抽搐| 国产成人精品婷婷| 亚洲四区av| 内地一区二区视频在线| 五月开心婷婷网| 国产黄频视频在线观看| 久久久久人妻精品一区果冻| 欧美丝袜亚洲另类| 精品少妇黑人巨大在线播放| 观看av在线不卡| 久久综合国产亚洲精品| 欧美3d第一页| www.av在线官网国产| 国产精品福利在线免费观看| 亚洲精品久久午夜乱码| 中文字幕免费在线视频6| 欧美+日韩+精品| 久久久亚洲精品成人影院| 中国三级夫妇交换| 高清欧美精品videossex| 涩涩av久久男人的天堂| 男女免费视频国产| 麻豆国产97在线/欧美| 夜夜骑夜夜射夜夜干| 久久久久久久久久久丰满| 日本黄色日本黄色录像| 最新中文字幕久久久久| 欧美成人一区二区免费高清观看| 丝袜喷水一区| av网站免费在线观看视频| 黑人高潮一二区| 97热精品久久久久久| 久久精品国产亚洲av天美| 国产69精品久久久久777片| 直男gayav资源| 欧美激情极品国产一区二区三区 | 内地一区二区视频在线| 成人亚洲欧美一区二区av| 91精品一卡2卡3卡4卡| 日韩欧美 国产精品| 成人国产av品久久久| a级毛片免费高清观看在线播放| av女优亚洲男人天堂| 精品一区二区三卡| 午夜日本视频在线| 美女xxoo啪啪120秒动态图| 国产中年淑女户外野战色| 成人午夜精彩视频在线观看| 国产精品国产三级专区第一集| 男女国产视频网站| 亚洲成人一二三区av| 高清视频免费观看一区二区| 国产一区二区在线观看日韩| 久久午夜福利片| av线在线观看网站| 在线天堂最新版资源| 99久久中文字幕三级久久日本| 亚洲综合精品二区| 国产在线男女| 日本av免费视频播放| 综合色丁香网| 亚洲精品乱码久久久v下载方式| 女人久久www免费人成看片| 又大又黄又爽视频免费| 美女cb高潮喷水在线观看| 丝袜喷水一区| 人妻夜夜爽99麻豆av| 午夜激情久久久久久久| 国产乱来视频区| 97超碰精品成人国产| 日产精品乱码卡一卡2卡三| 一级二级三级毛片免费看| 美女高潮的动态| 18+在线观看网站| 久久99精品国语久久久| 日本av手机在线免费观看| 国产v大片淫在线免费观看| 成人高潮视频无遮挡免费网站| av免费观看日本| 国产高清有码在线观看视频| 丝袜喷水一区| 日产精品乱码卡一卡2卡三| 久久人人爽av亚洲精品天堂 | 色吧在线观看| 久久青草综合色| 一级黄片播放器| 亚洲第一区二区三区不卡| 亚洲精品乱码久久久久久按摩| 国产 一区精品| 久久久久人妻精品一区果冻| 国产精品熟女久久久久浪| 香蕉精品网在线| 国产日韩欧美亚洲二区| 国产精品99久久久久久久久| 亚洲精品自拍成人| 国产一区二区在线观看日韩| 一边亲一边摸免费视频| 国产免费视频播放在线视频| 91久久精品国产一区二区成人| 男人爽女人下面视频在线观看| 嘟嘟电影网在线观看| 又黄又爽又刺激的免费视频.| 中国美白少妇内射xxxbb| 亚洲第一区二区三区不卡| 黄片wwwwww| 国产精品爽爽va在线观看网站| 毛片女人毛片| 亚洲人成网站在线播| 国产成人午夜福利电影在线观看| 国产av码专区亚洲av| 一区在线观看完整版| 午夜激情久久久久久久| 97热精品久久久久久| 国产精品蜜桃在线观看| 亚洲在久久综合| 免费人妻精品一区二区三区视频| 中国国产av一级| 香蕉精品网在线| 大陆偷拍与自拍| 777米奇影视久久| 少妇人妻精品综合一区二区| 看十八女毛片水多多多| xxx大片免费视频| 最近的中文字幕免费完整| 亚洲国产精品999| 搡女人真爽免费视频火全软件| 国产在视频线精品| 色5月婷婷丁香| 又大又黄又爽视频免费| 国产淫片久久久久久久久| 欧美老熟妇乱子伦牲交| 草草在线视频免费看| 我要看黄色一级片免费的| 午夜免费观看性视频| freevideosex欧美| 丰满少妇做爰视频| 国产毛片在线视频| 亚洲一级一片aⅴ在线观看| 国产精品99久久99久久久不卡 | 国产精品一区www在线观看| 91精品一卡2卡3卡4卡| 26uuu在线亚洲综合色| 亚洲va在线va天堂va国产| 99久久精品一区二区三区| 亚洲欧美日韩无卡精品| 日韩免费高清中文字幕av| 秋霞伦理黄片| 成人黄色视频免费在线看| videossex国产| av专区在线播放| 国产一区二区三区av在线| videossex国产| 日韩一区二区视频免费看| 91aial.com中文字幕在线观看| 大片免费播放器 马上看| 男人添女人高潮全过程视频| 在线观看一区二区三区| 中文乱码字字幕精品一区二区三区| 色5月婷婷丁香| 久久久成人免费电影| 久久久久网色| 成人亚洲精品一区在线观看 | 久久精品国产鲁丝片午夜精品| 交换朋友夫妻互换小说| 免费观看性生交大片5| 新久久久久国产一级毛片| 欧美xxxx性猛交bbbb| 日韩,欧美,国产一区二区三区| 国产成人91sexporn| 精品人妻熟女av久视频| 国产黄片美女视频| 久久久国产一区二区| 麻豆乱淫一区二区| 久久国产精品大桥未久av | 亚洲精品成人av观看孕妇| 免费大片黄手机在线观看| 日韩三级伦理在线观看| 国产爽快片一区二区三区| av在线蜜桃| 在线观看美女被高潮喷水网站| 成年av动漫网址| 男的添女的下面高潮视频| 久久99热这里只频精品6学生| 久久久a久久爽久久v久久| 免费人妻精品一区二区三区视频| 亚洲精品日本国产第一区| 欧美+日韩+精品| 97超视频在线观看视频| 亚洲国产成人一精品久久久| 男女免费视频国产| 久久午夜福利片| 国产午夜精品一二区理论片| 我的女老师完整版在线观看| 少妇熟女欧美另类| 制服丝袜香蕉在线| 亚洲高清免费不卡视频| 日日摸夜夜添夜夜添av毛片| 大片电影免费在线观看免费| 香蕉精品网在线| 国产免费一级a男人的天堂| 久久青草综合色| 男人添女人高潮全过程视频| 色婷婷av一区二区三区视频| 久久久久性生活片| 天堂中文最新版在线下载| 偷拍熟女少妇极品色| tube8黄色片| 深夜a级毛片| 一本一本综合久久| 成人免费观看视频高清| 成人18禁高潮啪啪吃奶动态图 | 一本一本综合久久| 波野结衣二区三区在线| 日本爱情动作片www.在线观看| 国产色爽女视频免费观看| 深夜a级毛片| 色婷婷久久久亚洲欧美| 精品国产一区二区三区久久久樱花 | 成人18禁高潮啪啪吃奶动态图 | 亚洲国产精品成人久久小说| 国产熟女欧美一区二区| 美女xxoo啪啪120秒动态图| 丝袜喷水一区| 国产精品秋霞免费鲁丝片| 国产美女午夜福利| 三级国产精品欧美在线观看| 久久午夜福利片| 亚洲无线观看免费| 精品一区在线观看国产| 中文字幕久久专区| 一本—道久久a久久精品蜜桃钙片| 国产亚洲一区二区精品| 久久精品人妻少妇| 成年人午夜在线观看视频| 国产在线免费精品| 国产高清国产精品国产三级 | 久久影院123| 涩涩av久久男人的天堂| 久久精品国产自在天天线| 久久久久国产精品人妻一区二区| 男女免费视频国产| 卡戴珊不雅视频在线播放| 久热这里只有精品99| 少妇人妻 视频| 国产乱来视频区| 三级国产精品欧美在线观看| 婷婷色综合www| 人妻制服诱惑在线中文字幕| 欧美国产精品一级二级三级 | 制服丝袜香蕉在线| 亚洲高清免费不卡视频| 免费大片黄手机在线观看| 久久久精品94久久精品| 青春草视频在线免费观看| 亚洲欧美成人综合另类久久久| 麻豆乱淫一区二区| 亚洲精品视频女| 中文字幕精品免费在线观看视频 | 插阴视频在线观看视频| 国产精品一区二区在线观看99| 美女xxoo啪啪120秒动态图| 精品一区二区三卡| av又黄又爽大尺度在线免费看| 久久久久久久国产电影| 汤姆久久久久久久影院中文字幕| 最近2019中文字幕mv第一页| 两个人的视频大全免费| 男人和女人高潮做爰伦理| 看免费成人av毛片| 日本爱情动作片www.在线观看| a级一级毛片免费在线观看| av免费在线看不卡| 久久久久久久大尺度免费视频| 国产乱人视频| 久久久久精品性色| av在线蜜桃| 能在线免费看毛片的网站| 一个人看视频在线观看www免费| 一级片'在线观看视频| 伊人久久国产一区二区| 国产精品秋霞免费鲁丝片| 成人二区视频| av播播在线观看一区| 国产美女午夜福利| 22中文网久久字幕| 男人和女人高潮做爰伦理| av福利片在线观看| 大片免费播放器 马上看| 国产成人freesex在线| 国产高清国产精品国产三级 | 精品少妇久久久久久888优播| 国产欧美日韩精品一区二区| 大片电影免费在线观看免费| 国产大屁股一区二区在线视频| 日本午夜av视频| 亚洲精品成人av观看孕妇| 纵有疾风起免费观看全集完整版| 成年人午夜在线观看视频| 插逼视频在线观看| 国产在视频线精品| 久久久久久久久久人人人人人人| 久久韩国三级中文字幕| 国产成人精品福利久久| 久久国产亚洲av麻豆专区| 亚洲久久久国产精品| 噜噜噜噜噜久久久久久91| 精品视频人人做人人爽| 国产男女内射视频| 国产伦精品一区二区三区四那| 日韩制服骚丝袜av| 又爽又黄a免费视频| 99热国产这里只有精品6| 久久精品久久久久久久性| 亚洲精品视频女| 亚洲欧美中文字幕日韩二区| 国产爱豆传媒在线观看| 黄色视频在线播放观看不卡| 国产免费福利视频在线观看| 在线 av 中文字幕| 欧美精品一区二区大全| 熟女人妻精品中文字幕| 日韩人妻高清精品专区| 少妇人妻精品综合一区二区| 日产精品乱码卡一卡2卡三| 欧美精品亚洲一区二区| 亚洲精品乱码久久久久久按摩| 亚洲国产色片| 国产爱豆传媒在线观看| 51国产日韩欧美| 永久免费av网站大全| 综合色丁香网| 精品少妇黑人巨大在线播放| 91在线精品国自产拍蜜月| 伦精品一区二区三区| 久久久久人妻精品一区果冻| 日日摸夜夜添夜夜添av毛片| 在线观看美女被高潮喷水网站| 国产精品不卡视频一区二区| 激情五月婷婷亚洲| 熟女av电影| 黑人猛操日本美女一级片| 中文字幕久久专区| 国国产精品蜜臀av免费| 亚洲精品一二三| 亚洲内射少妇av| 欧美成人精品欧美一级黄| 日韩欧美 国产精品| 国产精品熟女久久久久浪| 91狼人影院| 18禁裸乳无遮挡动漫免费视频| 97在线人人人人妻| 亚洲一级一片aⅴ在线观看| 精品亚洲成国产av| 在线观看三级黄色| 久久久精品免费免费高清| 又粗又硬又长又爽又黄的视频| 日韩视频在线欧美| 国国产精品蜜臀av免费| 欧美日本视频| 视频中文字幕在线观看| 免费看光身美女| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩东京热| 啦啦啦中文免费视频观看日本| 亚洲最大成人中文| 国产免费视频播放在线视频| 大陆偷拍与自拍| 一本久久精品| 高清黄色对白视频在线免费看 | 日韩欧美一区视频在线观看 | 久久久久性生活片| 久久国产乱子免费精品| 精品99又大又爽又粗少妇毛片| 日韩精品有码人妻一区| 日日啪夜夜爽| 国产一区亚洲一区在线观看| 日韩一区二区视频免费看| 97在线视频观看| 国产伦精品一区二区三区四那| 久久久精品94久久精品| 国产伦在线观看视频一区| 国产精品女同一区二区软件| 一个人看视频在线观看www免费| 免费观看性生交大片5| 欧美激情极品国产一区二区三区 | 天天躁夜夜躁狠狠久久av| 久久精品夜色国产| 深夜a级毛片| 在线亚洲精品国产二区图片欧美 | 亚洲天堂av无毛| 如何舔出高潮| www.色视频.com| 一个人免费看片子| 有码 亚洲区| 蜜桃在线观看..| 99九九线精品视频在线观看视频| 国产欧美亚洲国产| 欧美成人一区二区免费高清观看| 国产成人精品久久久久久| 国产人妻一区二区三区在| 亚洲精品视频女| 性色av一级| av国产免费在线观看| 狠狠精品人妻久久久久久综合| 免费观看a级毛片全部| 熟妇人妻不卡中文字幕| 最近最新中文字幕免费大全7| 久久国产精品男人的天堂亚洲 | 国产成人精品福利久久| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美一区二区三区黑人 | 蜜桃在线观看..| 视频中文字幕在线观看| 国产精品一区二区在线观看99| 亚洲va在线va天堂va国产| 在线看a的网站| 亚洲人与动物交配视频| 亚洲精品日韩av片在线观看| 精品亚洲乱码少妇综合久久| 中文精品一卡2卡3卡4更新| 成年女人在线观看亚洲视频| 成人漫画全彩无遮挡| 亚洲国产日韩一区二区| 一本久久精品| 日韩一本色道免费dvd| a级毛色黄片| 日本午夜av视频| 久久精品国产鲁丝片午夜精品| 久久久成人免费电影| 纯流量卡能插随身wifi吗| 婷婷色麻豆天堂久久| 男人爽女人下面视频在线观看| 熟女电影av网| 毛片女人毛片| 日产精品乱码卡一卡2卡三| 免费看日本二区| 99久久精品热视频| 男人添女人高潮全过程视频| 欧美一区二区亚洲| 亚洲色图av天堂| 我要看黄色一级片免费的| 国产乱人视频| 女的被弄到高潮叫床怎么办| 日韩大片免费观看网站| 久久久精品免费免费高清| 色视频在线一区二区三区| 久久女婷五月综合色啪小说| 国产老妇伦熟女老妇高清| 亚洲av中文av极速乱| 干丝袜人妻中文字幕| 久久久久精品久久久久真实原创| 国产日韩欧美在线精品| 深爱激情五月婷婷| 三级经典国产精品| 在线 av 中文字幕| 免费人妻精品一区二区三区视频| 国产久久久一区二区三区| av又黄又爽大尺度在线免费看| 久久毛片免费看一区二区三区| 天堂8中文在线网| 亚洲成色77777| 欧美高清性xxxxhd video| 美女主播在线视频| 又黄又爽又刺激的免费视频.| av一本久久久久| 黄色配什么色好看| 97热精品久久久久久| 欧美精品一区二区大全| 91精品国产九色| 99久久中文字幕三级久久日本| 久热这里只有精品99| 亚洲精品久久午夜乱码| 免费观看av网站的网址| 麻豆精品久久久久久蜜桃| 男女下面进入的视频免费午夜| 国产 一区 欧美 日韩| 亚洲欧洲日产国产| 岛国毛片在线播放| 国产淫片久久久久久久久| 欧美精品一区二区免费开放| 国产大屁股一区二区在线视频| 一级毛片电影观看| 女人久久www免费人成看片| 一本一本综合久久| 日本av免费视频播放| 中文在线观看免费www的网站| 久久久精品94久久精品| 看十八女毛片水多多多| 国产精品一区www在线观看| 日本午夜av视频| 99热这里只有是精品50| 青青草视频在线视频观看| 老司机影院毛片| 亚洲人与动物交配视频| 成人国产麻豆网| 成年女人在线观看亚洲视频| 国精品久久久久久国模美| 寂寞人妻少妇视频99o| 99热这里只有是精品在线观看| 亚洲av成人精品一二三区| 一级毛片aaaaaa免费看小| 最近中文字幕高清免费大全6| 欧美三级亚洲精品| 亚洲成人一二三区av| 一本—道久久a久久精品蜜桃钙片| 97超碰精品成人国产| 美女内射精品一级片tv| 日韩av在线免费看完整版不卡| 成年av动漫网址| 有码 亚洲区| 免费大片黄手机在线观看| 在线 av 中文字幕| 看非洲黑人一级黄片| 青春草视频在线免费观看| 久久99热这里只频精品6学生| 亚洲第一区二区三区不卡| 国产69精品久久久久777片| 日韩大片免费观看网站| 亚洲,欧美,日韩| 日韩免费高清中文字幕av| 自拍偷自拍亚洲精品老妇| 亚洲av在线观看美女高潮| 免费av中文字幕在线| 大香蕉久久网| 最近2019中文字幕mv第一页| 草草在线视频免费看| av视频免费观看在线观看| 99久久综合免费| 久久久午夜欧美精品| 精品人妻一区二区三区麻豆| 我要看日韩黄色一级片| 最近手机中文字幕大全| 国产一区二区三区av在线| 亚洲av成人精品一区久久| 校园人妻丝袜中文字幕| 国产女主播在线喷水免费视频网站| 特大巨黑吊av在线直播| 一级a做视频免费观看| 看十八女毛片水多多多| 在线免费观看不下载黄p国产| 精品一区二区免费观看| 中文资源天堂在线| 亚洲婷婷狠狠爱综合网| av在线蜜桃| 久久久久久久国产电影| 一级毛片电影观看| 大陆偷拍与自拍| 色视频www国产| 五月玫瑰六月丁香| 大香蕉久久网| 国产中年淑女户外野战色| 亚洲av.av天堂| 亚洲,欧美,日韩| 女性被躁到高潮视频| 在线观看一区二区三区| 久久毛片免费看一区二区三区| 成人无遮挡网站| 亚洲欧美精品专区久久| 成人综合一区亚洲| 国产成人免费无遮挡视频| 国产精品一区二区性色av| 亚洲av不卡在线观看| 黑人猛操日本美女一级片| 成人亚洲欧美一区二区av| 美女高潮的动态| 欧美精品国产亚洲| 看十八女毛片水多多多| 三级国产精品片| 五月玫瑰六月丁香| 女性生殖器流出的白浆| 最近的中文字幕免费完整| 色吧在线观看| 国产欧美亚洲国产| 亚洲精品第二区| 97超视频在线观看视频| 黄色配什么色好看| 99热6这里只有精品| 一本—道久久a久久精品蜜桃钙片| 国产在线视频一区二区| 亚洲人与动物交配视频| 丰满人妻一区二区三区视频av| 内射极品少妇av片p| 亚洲欧美中文字幕日韩二区| 九九久久精品国产亚洲av麻豆| 91狼人影院| 亚洲国产精品成人久久小说| 久久久成人免费电影| 亚洲久久久国产精品| 乱系列少妇在线播放| 国产成人免费观看mmmm| 久久久久久久久久成人| 国产精品一区www在线观看| 亚洲av男天堂| 亚洲在久久综合| 51国产日韩欧美| 亚洲精品色激情综合| 女人久久www免费人成看片| 国产视频内射| 国产精品免费大片| 九九在线视频观看精品| 久久久久久久大尺度免费视频| 狂野欧美激情性xxxx在线观看| 亚洲欧美一区二区三区黑人 | 久久国产精品大桥未久av | 久久久久久久久久久免费av| 三级国产精品片| 九色成人免费人妻av| 插逼视频在线观看| 午夜精品国产一区二区电影|