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

    光譜特征變量和BP神經(jīng)網(wǎng)絡(luò)構(gòu)建油用牡丹種子含水率估算模型

    2021-01-19 05:00:34劉秀英余俊茹王世華
    關(guān)鍵詞:油用牡丹波長(zhǎng)

    劉秀英,余俊茹,王世華

    (1.河南科技大學(xué)農(nóng)學(xué)院,洛陽 471023;2.洛陽市共微生物與綠色發(fā)展重點(diǎn)實(shí)驗(yàn)室,洛陽 471023;3.洛陽市植物營(yíng)養(yǎng)與環(huán)境生態(tài)重點(diǎn)實(shí)驗(yàn)室,洛陽 471023)

    0 引 言

    種子水分是評(píng)定種子質(zhì)量的重要指標(biāo)之一,通過種子水分測(cè)定,能夠及時(shí)采取措施,防止種子發(fā)熱、霉變、生蟲,避免造成種子變質(zhì),是保證種子質(zhì)量的重要手段[1]。傳統(tǒng)測(cè)定種子含水率的方法是室內(nèi)烘干稱重法,該方法測(cè)定結(jié)果準(zhǔn)確,但比較繁瑣,測(cè)定時(shí)間較長(zhǎng)。其他的直接或間接含水率測(cè)定方法雖然也較多,但都存在某些方面的缺點(diǎn),具有不同的適用范圍[2-3]。已有研究表明水分速測(cè)儀具有使用簡(jiǎn)單、測(cè)定迅速的特點(diǎn),但其測(cè)試結(jié)果的準(zhǔn)確性與適宜的溫度條件、種子的品種及質(zhì)量范圍有密切的關(guān)系[4]。目前,市場(chǎng)上種子水分快速檢測(cè)設(shè)備較多,但多數(shù)儀器適合測(cè)定特定的品種,需要適宜的工作環(huán)境,測(cè)定水分的范圍因品種而異。且油脂類種子含有不飽和脂肪酸,在研磨、剪切、加熱過程中容易氧化,因而較多方法和儀器并不適合進(jìn)行油脂類種子含水率測(cè)定[2]。

    近紅外光譜分析技術(shù)(Near Infrared Reflectance Spectroscopy,NIRS)是20世紀(jì)80年代發(fā)展起來的一項(xiàng)物理檢測(cè)技術(shù)。與傳統(tǒng)分析方法相比,該技術(shù)具有檢測(cè)速度快、非破壞性、無損樣品等優(yōu)點(diǎn),尤其適用于農(nóng)產(chǎn)品品質(zhì)分析[5-7]和食品品質(zhì)檢測(cè)[8-9]領(lǐng)域。由于在近紅外光譜區(qū),水的O-H基團(tuán)泛音和組合帶非常明顯,因此該技術(shù)的發(fā)展為種子含水率的測(cè)定提供了新的方法[10]。國(guó)內(nèi)外學(xué)者采用近紅外光譜技術(shù)結(jié)合化學(xué)計(jì)量學(xué)方法特別是偏最小二乘方法對(duì)種子含水率的快速檢測(cè)進(jìn)行了研究,其估算精度差異較大[10-15]。

    人工神經(jīng)網(wǎng)絡(luò)具有學(xué)習(xí)性、容錯(cuò)性、實(shí)時(shí)性等優(yōu)點(diǎn),能夠?qū)W習(xí)和模擬任意復(fù)雜的非線性函數(shù),因而被越來越多的應(yīng)用于可見/近紅外光譜分析[13,15-16]。牡丹籽油富含人體需要的氨基酸、維生素、不飽和脂肪酸(含量高達(dá)92%以上)等多種成分,其中“植物腦黃金”α-亞麻酸占42%左右[17],并且油用牡丹種子出油率較高,因此油用牡丹成為了一種新興的高端木本油料作物。含水率是衡量油用牡丹種子質(zhì)量、影響牡丹籽油品質(zhì)的重要指標(biāo),因而在銷售、儲(chǔ)存及加工過程中需要進(jìn)行快速、準(zhǔn)確測(cè)量。然而以油用牡丹種子為研究對(duì)象,使用BP神經(jīng)網(wǎng)絡(luò)模型反演種子含水率的研究目前尚不多見。為此,本研究以油用牡丹種子為研究對(duì)象,通過相關(guān)分析選擇特征變量作為自變量,建立一元線性回歸模型,然后將從中優(yōu)選的參數(shù)作為BP神經(jīng)網(wǎng)絡(luò)的輸入,油用牡丹種子含水率作為輸出,建立含水率BP神經(jīng)網(wǎng)絡(luò)估算模型,并對(duì)比分析一元線性回歸(SLR,Simple Linear Regression)模型、逐步多元線性回歸(SMLR,Stepwise Multiple Linear Regression)模型、偏最小二乘回歸(PLSR,Partial Least-squares Regression)模型和BP神經(jīng)網(wǎng)絡(luò)(BPNN,BP Neural Network)模型的估算精度,旨在探索油用牡丹種子含水率高光譜估算的最優(yōu)模型,提高油用牡丹種子含水率反演精度,為促進(jìn)中國(guó)牡丹產(chǎn)業(yè)發(fā)展及糧油生產(chǎn)、保障國(guó)家糧油安全提供理論依據(jù)和技術(shù)支持。

    1 材料與方法

    1.1 供試材料

    本研究收集了2個(gè)年份(即:2016年和2017年)、不同產(chǎn)地(即:安徽安慶、毫州、銅陵;河南洛陽;山東曹縣、吏口、平陰;四川德陽、眉山、資陽)、不同品種(“鳳丹”和“紫斑”)的油用牡丹種子 156份。分年份收集種子樣品后,首先進(jìn)行人工目視篩選,選擇顆粒比較飽滿、大小均勻一致的種子充分進(jìn)行自然干燥,盡可能的去除種子中的含水率。然后分裝在密封袋內(nèi),其中,山東曹縣、平陰、吏口的種子根據(jù)年份、品種不同分成6個(gè)平行樣,共獲得72個(gè)樣本,設(shè)置3個(gè)梯度;而其他省份不同地區(qū)的樣本根據(jù)年份、品種不同分成 3個(gè)平行樣,共獲得84個(gè)樣本,同樣設(shè)置3個(gè)梯度,每份種子樣本約30 g,( 135±10)顆,正好裝滿75 mm培養(yǎng)皿(見圖1),并根據(jù)年份、產(chǎn)地、品種不同進(jìn)行編號(hào)。然后將不同梯度樣本放入烘箱干燥,通過調(diào)節(jié)烘干的時(shí)間,形成不同的水分梯度,樣本含水率范圍為 0.22%~5.77%,詳細(xì)參數(shù)見表1。

    圖1 單個(gè)油用牡丹種子樣本Fig.1 Single sample of oil tree peony seeds

    表1 油用牡丹種子樣本含水率的統(tǒng)計(jì)特征Table 1 Statistical characteristic of moisture content measured oil tree peony seed samples

    1.2 測(cè)定方法

    含水率測(cè)定:種子含水率測(cè)定采用常壓烘箱干燥法,依照GB/T 5009.3—2003進(jìn)行。

    光譜數(shù)據(jù)采集:采用美國(guó)Spectra Vista公司的SVC HR-1024i便攜式光譜儀測(cè)定置于黑色背景上裝滿75 mm培養(yǎng)皿的油用牡丹種子光譜,該儀器的光譜范圍 350~2 500 nm。在波長(zhǎng)350~1 000 nm,光譜分辨率≤3.5 nm;1 000~1 850 nm,光譜分辨率≤9.5 nm;1 850~2 500 nm,光譜分辨率≤6.5 nm。光源為功率50 W鹵素?zé)?,光源入射角?45°;探頭視場(chǎng)角為 8°,垂直放置于目標(biāo)物正上方。儀器開機(jī)預(yù)熱30 min,然后將油用牡丹種子樣本裝盤,從4個(gè)方向?qū)悠愤M(jìn)行掃描,每個(gè)方向掃描3次,一個(gè)樣品總共掃描12次,去掉異常線后求平均值作為該樣品的光譜值。

    1.3 數(shù)據(jù)處理、模型構(gòu)建及驗(yàn)證

    通過對(duì)種子近紅外吸收光譜、一階微分值與含水率的相關(guān)分析,選擇相關(guān)系數(shù)絕對(duì)值大于0.85的波段即特征波長(zhǎng)吸收光譜及其特征波長(zhǎng)一階微分光譜作為自變量構(gòu)建含水率的一元線性回歸估算模型;將水分吸收特征參數(shù)與含水率進(jìn)行相關(guān)分析,選擇相關(guān)系數(shù)絕對(duì)值大于0.85的參數(shù)作為自變量構(gòu)建含水率的一元線性回歸估算模型。為了進(jìn)一步簡(jiǎn)化及優(yōu)化模型,構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型、多元線性回歸模型和偏最小二乘回歸模型時(shí),選用的輸入變量為構(gòu)建一元線性回歸模型時(shí),建模及驗(yàn)?zāi)?shù)R2大于0.8,驗(yàn)?zāi)MSE小于0.5%,而RPD大于2.5的變量。

    為了驗(yàn)證模型的預(yù)測(cè)精度,采用獨(dú)立驗(yàn)證樣本將模型預(yù)測(cè)值和實(shí)測(cè)值進(jìn)行回歸擬合,以預(yù)測(cè)決定系數(shù)(R2)、均方根誤差(RMSE)及相對(duì)預(yù)測(cè)偏差(RPD)作為指標(biāo)評(píng)價(jià)模型的學(xué)習(xí)能力和預(yù)測(cè)能力。RPD值可以解釋模型的預(yù)測(cè)能力,其評(píng)價(jià)采用 Chang等[18]提出的閾值劃分標(biāo)準(zhǔn)。R2及RPD值越大,RMSE越小,表明模型的精度越高。綜合建模和驗(yàn)?zāi)>葯z驗(yàn)結(jié)果,選擇出最優(yōu)估算模型。

    2 結(jié)果與分析

    2.1 油用牡丹種子的近紅外吸收光譜特征與水分吸收特征

    圖2和圖3為不同含水率油用牡丹種子樣本的近紅外吸收光譜及吸收特征曲線。由圖2可以看出,油用牡丹種子的近紅外吸收光譜隨含水率增加而增大,在整個(gè)波段范圍內(nèi),光譜曲線總體變化比較明顯,具有多個(gè)波峰、波谷;此外,光譜曲線具有相似的變化規(guī)律,在1 200、1 440、1 930、2 140 nm波長(zhǎng)附近有強(qiáng)烈的吸收峰,隨含水率增加吸收峰強(qiáng)度具有增大的趨勢(shì),而吸收峰位置具有朝長(zhǎng)波方向偏移的趨勢(shì)。依據(jù)前人研究可知,1 400、1 900和2 100 nm波長(zhǎng)附近的吸收峰是由于水分的強(qiáng)吸收引起的[10,19-20],因而采用包絡(luò)線法制作了吸收深度曲線圖(圖3)。從圖3中可以看出,吸收深度最大的波長(zhǎng)位于3個(gè)水分吸收帶,具體位置為1 440、1 930、2 140 nm波長(zhǎng)附近。當(dāng)油用牡丹種子含水率發(fā)生變化時(shí),1 930 nm波長(zhǎng)處的波段寬度最窄;2 140 nm波長(zhǎng)處的吸收深度及吸收面積變化最明顯,1 440 nm波長(zhǎng)處的吸收深度參數(shù)變化最小。因而可以提取1 440、1 930、2 140 nm波長(zhǎng)附近的水分吸收特征參數(shù)進(jìn)行油用牡丹種子含水率反演。

    圖2 油用牡丹種子的近紅外吸收光譜Fig.2 Near infrared spectrogram of oil tree peony seeds

    圖3 油用牡丹種子的光譜吸收深度曲線Fig.3 Spectral absorption depth curves of oil tree peony seeds

    2.2 光譜特征變量選擇

    對(duì)油用牡丹種子近紅外吸收光譜及一階微分光譜與含水率進(jìn)行相關(guān)分析,其相關(guān)系數(shù)見圖4。依據(jù)相關(guān)系數(shù)較大且位于 3個(gè)水分吸收帶附近的原則,分別選取了 3個(gè)特征波長(zhǎng)光譜及 3個(gè)特征波長(zhǎng)一階微分光譜變量(見圖4)用于后續(xù)分析建模,其相關(guān)系數(shù)大小見表2。分析表2可知,特征波長(zhǎng)光譜及特征波長(zhǎng)一階微分光譜與含水率的相關(guān)性較好,相關(guān)系數(shù)較高,均通過了0.01極顯著性檢驗(yàn)。特征波長(zhǎng)光譜變量與含水率的相關(guān)系數(shù)介于0.83~0.90之間,而特征波長(zhǎng)一階微分光譜變量與含水率的相關(guān)系數(shù)均大于 0.91,介于 0.91~0.96之間,后者的相關(guān)性明顯高于前者。選擇這3個(gè)特征波長(zhǎng)光譜及3個(gè)特征波長(zhǎng)一階微分光譜變量參與一元線性含水率估算模型的構(gòu)建。

    對(duì)光譜進(jìn)行包絡(luò)線去除后,提取水分吸收特征波段1 440、1 930、2 140 nm處的水分吸收特征參數(shù),包括:最大吸收深度(D)、吸收總面積(A)、吸收峰右面積(RA)、吸收峰左面積(LA)、面積歸一化最大吸收深度(AD),及吸收波段波長(zhǎng)位置(P)。將提取的吸收特征參數(shù)與含水率進(jìn)行相關(guān)分析,其相關(guān)系數(shù)的大小見表3。從表3可知,除P1930無法計(jì)算,LA1440、D2140、RA2140與含水率的相關(guān)系數(shù)在0.05水平上相關(guān)外,其余各參數(shù)與含水率的相關(guān)性均達(dá)到了 0.01水平極顯著相關(guān)。其中,面積歸一化最大吸收深度(AD)與含水率的相關(guān)性最好,兩者間的相關(guān)系數(shù)均大于0.9。比較3個(gè)波段位置的吸收特征參數(shù)與含水率的相關(guān)系數(shù)可知,2 140 nm波長(zhǎng)處的吸收特征參數(shù)與含水率的相關(guān)性明顯優(yōu)于另外兩個(gè)波長(zhǎng)位置的參數(shù)??傮w來看,吸收特征參數(shù)與含水率的相關(guān)系數(shù)高于0.85的有7個(gè),說明這些吸收特征參數(shù)對(duì)含水率變化更敏感,因而選擇這 7個(gè)吸收特征參數(shù)作為變量參與一元線性含水率估算模型的構(gòu)建。

    圖4 近紅外吸收光譜、一階微分光譜與含水率的相關(guān)性Fig.4 Correlation between near infrared absorption spectra and first derivative of spectra and moisture content

    表2 油用牡丹種子特征波長(zhǎng)光譜與含水率的相關(guān)系數(shù)Table 2 Correlation coefficients between characteristic wavelength spectrum of oil tree peony seed and moisture content

    2.3 基于近紅外光譜的油用牡丹種子含水率估算

    2.3.1 含水率的一元線性回歸估算模型

    建立 13個(gè)光譜特征變量與含水率的一元線性估算模型,然后采用獨(dú)立樣本進(jìn)行驗(yàn)證,其建模及驗(yàn)?zāi)?shù)見表4。分析表4可知,除R1410和S1930這2個(gè)變量建立的線性估算模型的驗(yàn)證RPD值小于2.0外,其余變量建立的一元線性估算模型的驗(yàn)證RPD值均大于2.0,校正及驗(yàn)證R2均大于0.75,說明這些變量建立的模型均能用于油用牡丹種子含水率快速、無損檢測(cè)。3個(gè)特征波長(zhǎng)一階微分光譜建立的估算模型效果均較好,建模及驗(yàn)?zāi)2和RPD值均較大,而RMSE值均較小,尤其以DF1950、DF2080這2個(gè)參數(shù)建立的模型估算效果更佳,其建模及驗(yàn)?zāi)2均高于0.92,RPD值高于4.0,而RMSE值小于0.4%。

    表3 油用牡丹種子水分吸收特征參數(shù)與含水率的相關(guān)系數(shù)Table 3 Correlation coefficients between absorptioncharacteristics parameters of oil tree peony seed and moisture content

    吸收特征參數(shù)中以面積歸一化吸收深度參數(shù)建立的估算模型效果最好,尤其以AD2140建立的估算模型效果最佳,其建模及驗(yàn)?zāi)2高于0.94,RPD值高于4.0,而RMSE值小于0.4%。綜合分析比較一元線性模型的建模及驗(yàn)?zāi)5母黜?xiàng)參數(shù)值可知,以DF2080和AD2140為自變量建立的模型優(yōu)于其他光譜特征變量建立的模型,都能對(duì)含水率進(jìn)行很好的估算;而且這2個(gè)模型建模及驗(yàn)?zāi)8黜?xiàng)參數(shù)幾乎相同,相比較而言,以DF2080為自變量建立的模型均方根誤差RMSEv較小,而RPDv值較大所以該模型為一元線性模型中的最優(yōu)估算模型。由于一元線性模型只是應(yīng)用了一個(gè)光譜特征變量建模,為了進(jìn)一步提高含水率預(yù)測(cè)精度,依據(jù)建模及驗(yàn)?zāi)?shù)R2大于0.8,RMSEv小于0.5%,而RPDv大于2.5的原則,從13個(gè)光譜特征變量中進(jìn)一步優(yōu)選出 5個(gè)光譜特征變量(即 DF1150、DF1950、DF2080、AD1930、AD2140)用于逐步多元線性回歸、偏最小二乘回歸及人工神經(jīng)網(wǎng)絡(luò)建模。

    表4 基于光譜特征變量的油用牡丹種子含水率估算模型參數(shù)Table 4 The model parameters of moisture content in oil peony seed based on spectral characteristic variable

    2.3.2 含水率的BP神經(jīng)網(wǎng)絡(luò)模型

    以優(yōu)選出來的5個(gè)光譜特征變量(DF1150、DF1950、DF2080、AD1930、AD2140)作為輸入層,油用牡丹種子實(shí)測(cè)含水率作為輸出層,隱含層節(jié)點(diǎn)數(shù)q根據(jù)經(jīng)驗(yàn)公式(1)[21]給定的范圍,通過試錯(cuò)法多次訓(xùn)練后選出最佳節(jié)點(diǎn)數(shù)為5,建立一個(gè)網(wǎng)絡(luò)結(jié)構(gòu)為5-5-1的BP神經(jīng)網(wǎng)絡(luò)模型。

    式中k為輸入層單元數(shù);m為輸出層單元數(shù);α為[1,10]之間的常數(shù)。

    利用Matlab編程進(jìn)行BP神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練,為了使網(wǎng)絡(luò)在訓(xùn)練階段容易收斂,首先對(duì)輸入輸出變量進(jìn)行歸一化處理[22];網(wǎng)絡(luò)中間層神經(jīng)元傳遞函數(shù)采用 S型正切函數(shù)Tansig,輸出層神經(jīng)元傳遞函數(shù)采用線性函數(shù) Purelin,訓(xùn)練函數(shù)為 trainlm[23]。通過 BP神經(jīng)網(wǎng)絡(luò)分析得到含水率的預(yù)測(cè)值,然后對(duì)獨(dú)立驗(yàn)證樣本的實(shí)測(cè)值與預(yù)測(cè)值進(jìn)行擬合及精度檢驗(yàn),結(jié)果見表5。模型的訓(xùn)練值和目標(biāo)值均方根誤差RMSE為0.220%,相關(guān)系數(shù)為0.989,相對(duì)分析誤差RPD值為6.478,模擬結(jié)果很好。預(yù)測(cè)值和實(shí)測(cè)值的均方根誤差RMSE為0.242%,預(yù)測(cè)值和實(shí)測(cè)值相關(guān)系數(shù)達(dá)到 0.986,相對(duì)分析誤差RPD值為5.889,模型預(yù)測(cè)效果也非常好。綜上所述,BP神經(jīng)網(wǎng)絡(luò)模型的建模及驗(yàn)證精度都非常高,可進(jìn)行油用牡丹種子含水率精確估算,幾乎可以與實(shí)驗(yàn)室測(cè)量結(jié)果相媲美。

    2.3.3 BP神經(jīng)網(wǎng)絡(luò)模型與回歸模型的比較

    將BP神經(jīng)網(wǎng)絡(luò)模型與一元線性回歸模型、逐步多元線性回歸模型以及目前在含水率預(yù)測(cè)中應(yīng)用非常廣泛的偏最小二乘回歸模型進(jìn)行比較,其結(jié)果見表5。為了保證模型具有可比性,建立SMLR模型和PLSR模型時(shí),選取與構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型相同的自變量(即優(yōu)選的 5個(gè)光譜特征變量 DF1150、DF1950、DF2080、AD1930、AD2140)作為輸入變量,并采用相同的獨(dú)立樣本進(jìn)行驗(yàn)證。表5列出了不同估算模型的建模及驗(yàn)證結(jié)果。結(jié)果表明,BPNN模型的建模效果最好,建模決定系數(shù)R2為0.978,RMSE為0.220%;其次為多元線性模型,建模決定系數(shù)R2為0.961,RMSE為0.289%。就模型檢驗(yàn)效果而言,BP神經(jīng)網(wǎng)絡(luò)模型驗(yàn)證效果仍然最佳,其次為多元線性模型,均優(yōu)于PLSR模型和一元線性模型。PLSR模型驗(yàn)證決定系數(shù)Rv2為0.955,RMSEv和 RPDv分別為 0.242%和 5.889,預(yù)測(cè)精度略低于 BP神經(jīng)網(wǎng)絡(luò)模型和逐步多元線性模型;而一元線性模型無論是建模還是驗(yàn)?zāi)PЧ畈睢5?,一元線性模型僅用了一個(gè)光譜特征變量建模,從建模及驗(yàn)?zāi)8黜?xiàng)參數(shù)來看估算效果也比較好,精度比較高。由于該方法最簡(jiǎn)單且容易操作,在能夠滿足油用牡丹種子含水率估算的精度情況下,不失為優(yōu)先選擇的估算方法。綜合比較4類模型的建模及驗(yàn)?zāi)8黜?xiàng)參數(shù)可知,4類模型的預(yù)測(cè)精度都比較高,而BPNN模型的預(yù)測(cè)效果最好,其次是多元線性估算模型。

    表5 不同估算模型的建模及驗(yàn)證結(jié)果Table 5 Modeling and testing results under different estimation

    3 討 論

    種子的近紅外光譜是其各種組成成分及其內(nèi)、外結(jié)構(gòu)的綜合表現(xiàn),兩者之間存在著密切關(guān)系。油用牡丹種子樣品中含有大量的C-H、O-H、N-H等含氫基團(tuán),因而在近紅外光譜圖中出現(xiàn)很強(qiáng)烈的吸收,形成吸收峰。但是高光譜數(shù)據(jù)具有波段多、數(shù)據(jù)量大、冗余性強(qiáng)等特點(diǎn),若直接用全波段數(shù)據(jù)進(jìn)行數(shù)據(jù)建模就會(huì)導(dǎo)致建模效率低、模型的性能差[24],因此本研究通過相關(guān)分析篩選出了不同油用牡丹種子含水率的敏感波長(zhǎng),有利于構(gòu)建較高精度的、比較簡(jiǎn)單的估測(cè)模型。將近紅外吸收光譜、一階微分光譜與含水率進(jìn)行相關(guān)分析,結(jié)果表明一階微分光譜與含水率的相關(guān)性明顯高于近紅外吸收光譜。通過相關(guān)分析選擇了相關(guān)系數(shù)較大且位于水分吸收波長(zhǎng)附近的3個(gè)特征波長(zhǎng)1 410、1 900、1 990 nm處的吸收光譜,和3個(gè)特征波長(zhǎng)1 150、1 950、2 080 nm處的一階微分光譜作為變量,參與一元線性回歸建模,結(jié)果表明特征波長(zhǎng)一階微分光譜建立的模型預(yù)測(cè)效果更好,究其原因可能是一階微處分光譜可以部分消除線性或接近線性噪聲的影響[23]。已有研究表明,含水率不同時(shí)吸收光譜的吸收峰的位置會(huì)產(chǎn)生偏移,通過微分變換可以進(jìn)一步增強(qiáng)高頻信息,因而微分變換可以調(diào)整波長(zhǎng)光譜,使吸收峰的位置更加明顯[10],同樣,通過相關(guān)分析可以使不同含水率的吸收光譜的吸收峰的位置更加明顯。本研究分析吸收光譜曲線探究水分的吸收峰位置與通過相關(guān)分析和微分變換選擇的特征吸收峰的具體波段位置雖然稍有差異,但大多都位于1 400、1 900、2 100 nm附近,可見不同方法的分析結(jié)果具有一致性。

    通過包絡(luò)線消除法可以去除那些不感興趣的吸收特征,孤立單個(gè)感興趣的吸收特征,并將其歸一到一個(gè)一致的光譜背景上,從而具有可以擴(kuò)大弱吸收特征信息的優(yōu)勢(shì)[25]。故本研究采用包絡(luò)線消除法提取了 1 440、1 930和2 140 nm 3個(gè)水分吸收波段附近的18個(gè)吸收特征參數(shù),通過相關(guān)分析選擇了7個(gè)相關(guān)系數(shù)大于0.85的水分吸收特征參數(shù)參與一元線性回歸分析建模。通過比較分析建模及預(yù)測(cè)結(jié)果可知,3個(gè)吸收深度參數(shù)為變量構(gòu)建的模型優(yōu)于其他吸收特征參數(shù)為變量構(gòu)建的模型;并且2 100 nm波長(zhǎng)位置附近的參數(shù)(比如:AD2140)建立的模型優(yōu)于1 400和1 900 nm波長(zhǎng)位置附近參數(shù)建立的模型。分析原因可能是:從吸收深度圖可以發(fā)現(xiàn),不同含水率的種子樣本在 2 100 nm波長(zhǎng)位置附近吸收峰的吸收深度和吸收面積的變化最明顯,說明含水率的變化對(duì)這個(gè)波長(zhǎng)位置的光譜峰值變化影響更大,從而這個(gè)位置的參數(shù)建立的模型效果更優(yōu),這一結(jié)論有待進(jìn)一步研究。

    建模方法是影響可見/近紅外光譜分析法預(yù)測(cè)精度的一個(gè)主要因素,因此采用多變量校正算法來構(gòu)建光譜數(shù)據(jù)與分析目標(biāo)之間的關(guān)系并進(jìn)行預(yù)測(cè)精度比較是非常重要的[24]。本研究比較了SLR、SMLR、PLSR和BPNN方法建立的含水率估算模型的預(yù)測(cè)精度,結(jié)果表明 BPNN模型的預(yù)測(cè)精度最高,其次為SMLR模型,這與其他研究者的結(jié)論一致[23,26]。多元回歸模型可以實(shí)現(xiàn)降維和剔除冗余信息[27],顯著提升光譜分析的精度和可靠性[28],因而模型的預(yù)測(cè)精度相對(duì)較高。而BPNN模型預(yù)測(cè)結(jié)果優(yōu)于PLSR模型,原因可能是BPNN模型能夠解釋光譜變量與含水率間存在的非線性關(guān)系,而偏最小二乘回歸模型是一種線性算法,沒有考慮光譜變量中某些潛在的非線性信息。BP神經(jīng)網(wǎng)絡(luò)是一種采用誤差反向傳播算法對(duì)網(wǎng)絡(luò)權(quán)值進(jìn)行訓(xùn)練的多層前向網(wǎng)絡(luò),通過反向傳播途徑能夠不斷的調(diào)整網(wǎng)絡(luò)的閾值與權(quán)值,直至滿足誤差最小精度條件,輸出最優(yōu)結(jié)果[29]。在數(shù)據(jù)分析時(shí),理論上BP神經(jīng)網(wǎng)絡(luò)模型能夠無限逼近任意復(fù)雜的非線性函數(shù),因而對(duì)非線性問題有很好的解釋性。但是,BP神經(jīng)網(wǎng)絡(luò)算法的網(wǎng)絡(luò)權(quán)值初始化是隨機(jī)的,神經(jīng)網(wǎng)絡(luò)程序在每次運(yùn)行后得到的訓(xùn)練結(jié)果是不同的[23,30],因此BP 神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)過程很難取得全局最佳效果。而且BP神經(jīng)網(wǎng)絡(luò)的隱含層的層數(shù)及節(jié)點(diǎn)數(shù)對(duì)網(wǎng)絡(luò)的模擬及預(yù)測(cè)精度影響較大[24],所以與其他回歸模型相比,其在實(shí)用性方面還有待于進(jìn)一步的研究。此外,為了驗(yàn)證所建模型的實(shí)用性,下一步將進(jìn)行另外批次自然含水率的種子驗(yàn)證試驗(yàn);其次,本研究油用牡丹種子含水率比較低,范圍比較窄,為了拓寬模型的適用范圍,進(jìn)一步研究將擴(kuò)大種子含水率的范圍。

    4 結(jié) 論

    本研究針對(duì)油用牡丹種子的含水率,采用包絡(luò)線消除法提取水分吸收特征參數(shù),分析了含水率與近紅外吸收光譜、一階微分光譜及水分吸收特征參數(shù)的相關(guān)關(guān)系,建立了油用牡丹種子含水率估算的一元線性回歸模型、逐步多元線性回歸模型、偏最小二乘回歸模型和BP神經(jīng)網(wǎng)絡(luò)模型,并比較了不同模型的預(yù)測(cè)精度,得出如下結(jié)論:

    1)通過相關(guān)分析確定了含水率的吸收光譜特征波長(zhǎng)為 14 10、1 900、1 990 nm,一階微分光譜特征波長(zhǎng)為1 150、1 950、2 080 nm;

    2)水分吸收特征參數(shù)與含水率的相關(guān)性較好,相關(guān)系數(shù)大于 0.9的水分吸收特征參數(shù)為吸收深度參數(shù),即AD1930、AD2140、AD1440;

    3)以DF2080和AD2140為自變量建立的一元線性回歸模型預(yù)測(cè)效果明顯優(yōu)于其他模型;

    4)將構(gòu)建一元線性回歸模型時(shí)篩選出的光譜特征變量作為輸入,實(shí)測(cè)含水率值作為輸出,構(gòu)建含水率估算的逐步多元線性回歸模型、偏最小二乘回歸模型及BP神經(jīng)網(wǎng)絡(luò)模型。其中BP神經(jīng)網(wǎng)絡(luò)模型的建模與驗(yàn)?zāi)2分別為0.978和0.973,RMSE分別為0.220%和0.242%,而RPD值分別為6.478和5.889,與其他回歸模型相比,BP神經(jīng)網(wǎng)絡(luò)模型的建模及驗(yàn)?zāi)>染罡?,且比較穩(wěn)定,是估算油用牡丹種子含水率的最優(yōu)模型,其次為逐步多元線性回歸模型。在能夠滿足油用牡丹種子含水率估算的精度情況下,一元線性回歸模型估算方法可以作為優(yōu)先選擇。

    猜你喜歡
    油用牡丹波長(zhǎng)
    HPLC-PDA雙波長(zhǎng)法同時(shí)測(cè)定四季草片中沒食子酸和槲皮苷的含量
    “三不夠”牡丹節(jié)
    油用牡丹助農(nóng)增收
    油用牡丹在園林中的運(yùn)用
    牡丹的整形修剪
    雙波長(zhǎng)激光治療慢性牙周炎的療效觀察
    綠牡丹
    河北省油用牡丹發(fā)展前景分析及對(duì)策
    “牡丹”情
    北方音樂(2016年12期)2016-08-23 03:20:03
    日本研發(fā)出可完全覆蓋可見光波長(zhǎng)的LED光源
    少妇被粗大猛烈的视频| 中文乱码字字幕精品一区二区三区| 中文字幕亚洲精品专区| 80岁老熟妇乱子伦牲交| 国产午夜福利久久久久久| 草草在线视频免费看| 噜噜噜噜噜久久久久久91| 99久久精品国产国产毛片| 别揉我奶头 嗯啊视频| 伊人久久精品亚洲午夜| 只有这里有精品99| 国产精品国产三级国产av玫瑰| 丝袜美腿在线中文| 简卡轻食公司| 秋霞在线观看毛片| 丝瓜视频免费看黄片| 国产精品人妻久久久久久| 欧美精品一区二区大全| 国产成年人精品一区二区| 丰满少妇做爰视频| 久久精品国产亚洲av天美| 十八禁网站网址无遮挡 | 亚洲人成网站在线观看播放| 亚洲av欧美aⅴ国产| 亚洲自偷自拍三级| 春色校园在线视频观看| 又爽又黄a免费视频| 成人免费观看视频高清| 3wmmmm亚洲av在线观看| 黄色怎么调成土黄色| 国产精品国产三级国产专区5o| 亚洲丝袜综合中文字幕| 国产免费一级a男人的天堂| 免费不卡的大黄色大毛片视频在线观看| 美女内射精品一级片tv| 久热久热在线精品观看| 亚洲精品一二三| 22中文网久久字幕| 高清av免费在线| 男人爽女人下面视频在线观看| 亚洲人成网站高清观看| 免费观看在线日韩| 色吧在线观看| 国产淫片久久久久久久久| 国产免费又黄又爽又色| 精品熟女少妇av免费看| 亚洲国产精品专区欧美| 91久久精品国产一区二区成人| 免费看日本二区| 国产午夜精品一二区理论片| 国产午夜精品久久久久久一区二区三区| 噜噜噜噜噜久久久久久91| 99热6这里只有精品| 国产伦精品一区二区三区四那| 午夜福利高清视频| 亚洲自偷自拍三级| 日日摸夜夜添夜夜爱| 亚洲高清免费不卡视频| 免费播放大片免费观看视频在线观看| 我要看日韩黄色一级片| 精品一区二区三区视频在线| 中文在线观看免费www的网站| 伊人久久精品亚洲午夜| 婷婷色麻豆天堂久久| 国产成人免费观看mmmm| 久久97久久精品| 97精品久久久久久久久久精品| 免费黄频网站在线观看国产| www.av在线官网国产| 久久国产乱子免费精品| 精品一区在线观看国产| 下体分泌物呈黄色| 老女人水多毛片| 国产精品.久久久| 国产亚洲5aaaaa淫片| 少妇人妻久久综合中文| a级毛片免费高清观看在线播放| 国产精品精品国产色婷婷| 视频中文字幕在线观看| av国产久精品久网站免费入址| 国产精品久久久久久精品电影小说 | 亚洲精品,欧美精品| 久久久久久久久大av| 国产视频首页在线观看| 六月丁香七月| 水蜜桃什么品种好| 久久精品国产亚洲av天美| 91精品一卡2卡3卡4卡| 99精国产麻豆久久婷婷| 久久精品久久久久久噜噜老黄| 亚洲,一卡二卡三卡| 丝瓜视频免费看黄片| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看| 男人和女人高潮做爰伦理| 联通29元200g的流量卡| 天天躁夜夜躁狠狠久久av| 亚洲精品国产色婷婷电影| 性色av一级| 日韩国内少妇激情av| 日韩伦理黄色片| 97热精品久久久久久| 国产 一区 欧美 日韩| 天堂中文最新版在线下载 | 九九爱精品视频在线观看| 老司机影院毛片| 国产精品国产三级国产av玫瑰| 亚洲国产色片| 亚洲在久久综合| 99热国产这里只有精品6| 一区二区三区免费毛片| 国产成人a区在线观看| 欧美精品一区二区大全| 听说在线观看完整版免费高清| 久热久热在线精品观看| 三级男女做爰猛烈吃奶摸视频| 日韩大片免费观看网站| 丝袜脚勾引网站| 国产精品不卡视频一区二区| 国产精品久久久久久精品电影| 禁无遮挡网站| 尾随美女入室| 日日撸夜夜添| 国产美女午夜福利| 特级一级黄色大片| 久久99热这里只有精品18| 亚洲精品国产av成人精品| 一级毛片aaaaaa免费看小| 亚洲av中文av极速乱| 嫩草影院新地址| 直男gayav资源| 国产午夜福利久久久久久| 男女国产视频网站| 久热久热在线精品观看| 欧美区成人在线视频| 亚洲av福利一区| 少妇人妻精品综合一区二区| 亚洲国产日韩一区二区| 26uuu在线亚洲综合色| 最后的刺客免费高清国语| 免费人成在线观看视频色| 国产亚洲5aaaaa淫片| 亚洲国产精品成人综合色| av线在线观看网站| 我的老师免费观看完整版| 亚洲精品国产成人久久av| 男人狂女人下面高潮的视频| 国产av码专区亚洲av| 一边亲一边摸免费视频| 国产精品一二三区在线看| 岛国毛片在线播放| 一本色道久久久久久精品综合| 乱系列少妇在线播放| 亚洲在久久综合| 日日撸夜夜添| 亚洲成人一二三区av| av福利片在线观看| 国产黄a三级三级三级人| 日本黄色片子视频| 亚洲av二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91 | 大香蕉97超碰在线| 18+在线观看网站| 免费观看a级毛片全部| 午夜老司机福利剧场| 午夜免费观看性视频| 亚洲欧美日韩无卡精品| 丰满少妇做爰视频| 国产又色又爽无遮挡免| 国内揄拍国产精品人妻在线| 少妇 在线观看| 白带黄色成豆腐渣| 国产男女内射视频| 嘟嘟电影网在线观看| 精品国产露脸久久av麻豆| 亚洲欧洲日产国产| 亚洲国产成人一精品久久久| 中文天堂在线官网| 久久6这里有精品| 一级爰片在线观看| 中文欧美无线码| 亚洲综合精品二区| 在线播放无遮挡| 亚洲国产av新网站| 狂野欧美白嫩少妇大欣赏| 国产毛片a区久久久久| 大码成人一级视频| 亚洲精品456在线播放app| 插逼视频在线观看| 狂野欧美激情性xxxx在线观看| 成年av动漫网址| 91久久精品电影网| 啦啦啦中文免费视频观看日本| 观看免费一级毛片| 免费高清在线观看视频在线观看| 欧美成人a在线观看| 黄色欧美视频在线观看| 国产一区二区三区av在线| 性色avwww在线观看| 另类亚洲欧美激情| 亚洲av成人精品一二三区| 亚洲成人久久爱视频| 女人久久www免费人成看片| 爱豆传媒免费全集在线观看| 国产亚洲一区二区精品| 国产av码专区亚洲av| av在线天堂中文字幕| 久久久久久久久久久丰满| 日韩av在线免费看完整版不卡| 国产免费视频播放在线视频| 国产伦理片在线播放av一区| 高清av免费在线| 亚洲人成网站在线播| 精品人妻偷拍中文字幕| 色播亚洲综合网| 可以在线观看毛片的网站| 少妇被粗大猛烈的视频| 国产成人aa在线观看| 亚洲经典国产精华液单| 国产视频首页在线观看| 国产乱人偷精品视频| 亚洲精品久久久久久婷婷小说| 国产免费一区二区三区四区乱码| 91狼人影院| 亚州av有码| 免费看日本二区| 永久免费av网站大全| 日日摸夜夜添夜夜爱| 亚洲av中文字字幕乱码综合| 精品午夜福利在线看| 国产精品久久久久久久久免| 午夜福利在线在线| 国产亚洲av嫩草精品影院| 高清视频免费观看一区二区| 欧美另类一区| 欧美一级a爱片免费观看看| 免费观看a级毛片全部| 国产精品久久久久久精品电影小说 | 国产老妇女一区| 嫩草影院精品99| 免费黄频网站在线观看国产| 草草在线视频免费看| 国产乱来视频区| 99精国产麻豆久久婷婷| 亚洲精品久久午夜乱码| 成人欧美大片| 国产男人的电影天堂91| 欧美日韩亚洲高清精品| 国产乱人视频| 好男人在线观看高清免费视频| 欧美激情久久久久久爽电影| 七月丁香在线播放| 一本一本综合久久| 大香蕉久久网| 自拍欧美九色日韩亚洲蝌蚪91 | 2021少妇久久久久久久久久久| 日韩,欧美,国产一区二区三区| 大片电影免费在线观看免费| 又爽又黄a免费视频| 少妇人妻一区二区三区视频| 亚洲欧美日韩另类电影网站 | 亚洲精品日韩在线中文字幕| 2021少妇久久久久久久久久久| 国产成人精品福利久久| 青春草亚洲视频在线观看| 一级爰片在线观看| 亚洲怡红院男人天堂| 热re99久久精品国产66热6| 午夜福利网站1000一区二区三区| 美女cb高潮喷水在线观看| 久久鲁丝午夜福利片| 97在线视频观看| 久久久a久久爽久久v久久| av免费观看日本| 亚洲欧美一区二区三区国产| av在线蜜桃| 一区二区三区免费毛片| 亚洲最大成人中文| 中国国产av一级| 在线观看一区二区三区激情| 日韩不卡一区二区三区视频在线| 免费看av在线观看网站| 亚洲精品自拍成人| 精品一区二区三区视频在线| 99视频精品全部免费 在线| 久久综合国产亚洲精品| 美女被艹到高潮喷水动态| 久久久久久久久大av| 女的被弄到高潮叫床怎么办| 日本午夜av视频| 一级黄片播放器| 两个人的视频大全免费| 国语对白做爰xxxⅹ性视频网站| 国产欧美亚洲国产| 亚洲精品成人av观看孕妇| 亚洲欧美一区二区三区黑人 | 亚洲av一区综合| 在线观看一区二区三区激情| 97超视频在线观看视频| 内地一区二区视频在线| 久久国产乱子免费精品| 国产精品一二三区在线看| 午夜福利高清视频| 五月天丁香电影| 99视频精品全部免费 在线| 日日摸夜夜添夜夜爱| 高清视频免费观看一区二区| 国产精品女同一区二区软件| 91精品一卡2卡3卡4卡| 麻豆国产97在线/欧美| 老女人水多毛片| 一区二区av电影网| 蜜桃亚洲精品一区二区三区| 亚洲真实伦在线观看| 国产成人aa在线观看| 99精国产麻豆久久婷婷| 国产精品久久久久久av不卡| 亚洲成色77777| av在线播放精品| 日韩成人伦理影院| 激情五月婷婷亚洲| 成人一区二区视频在线观看| 高清av免费在线| 亚洲人成网站高清观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲真实伦在线观看| 水蜜桃什么品种好| 中文天堂在线官网| 亚洲国产色片| 国产永久视频网站| 成年人午夜在线观看视频| av.在线天堂| 国模一区二区三区四区视频| 人人妻人人澡人人爽人人夜夜| 久久97久久精品| 亚洲av中文字字幕乱码综合| 九九久久精品国产亚洲av麻豆| 国产乱来视频区| 成人毛片a级毛片在线播放| 观看美女的网站| 亚洲在线观看片| 丰满人妻一区二区三区视频av| 别揉我奶头 嗯啊视频| 久久久久精品久久久久真实原创| 亚洲不卡免费看| 街头女战士在线观看网站| 久久久久国产精品人妻一区二区| 亚洲精品日韩在线中文字幕| 建设人人有责人人尽责人人享有的 | 肉色欧美久久久久久久蜜桃 | 亚洲精品视频女| 国产精品一区www在线观看| 身体一侧抽搐| 人妻少妇偷人精品九色| av网站免费在线观看视频| 亚洲四区av| 夜夜看夜夜爽夜夜摸| 国产在视频线精品| 插逼视频在线观看| 高清欧美精品videossex| 美女脱内裤让男人舔精品视频| 天天躁日日操中文字幕| 三级男女做爰猛烈吃奶摸视频| 亚洲精品影视一区二区三区av| 伦理电影大哥的女人| 免费观看性生交大片5| 内射极品少妇av片p| 午夜免费鲁丝| 免费大片黄手机在线观看| 亚洲,一卡二卡三卡| 久久久久国产网址| 国产成人a∨麻豆精品| 国产毛片a区久久久久| 日日摸夜夜添夜夜爱| 五月开心婷婷网| 日日摸夜夜添夜夜爱| 看十八女毛片水多多多| 美女内射精品一级片tv| 九九在线视频观看精品| 2021少妇久久久久久久久久久| 成年女人看的毛片在线观看| 午夜爱爱视频在线播放| 国产乱人视频| 欧美丝袜亚洲另类| 免费高清在线观看视频在线观看| 男人狂女人下面高潮的视频| 最新中文字幕久久久久| 亚洲精品一区蜜桃| 日韩大片免费观看网站| 中文字幕久久专区| 又粗又硬又长又爽又黄的视频| 亚洲在线观看片| 少妇的逼好多水| 80岁老熟妇乱子伦牲交| av在线天堂中文字幕| 91狼人影院| 久久99热这里只频精品6学生| 伊人久久精品亚洲午夜| 欧美成人a在线观看| 在线a可以看的网站| 国模一区二区三区四区视频| 日本黄大片高清| 麻豆成人午夜福利视频| 中文字幕av成人在线电影| 免费黄色在线免费观看| 久久99精品国语久久久| 久久精品久久久久久久性| 国产真实伦视频高清在线观看| 蜜桃亚洲精品一区二区三区| 99久久精品热视频| 99热这里只有是精品50| 精品久久国产蜜桃| 欧美日本视频| 国产av国产精品国产| h日本视频在线播放| 欧美日韩国产mv在线观看视频 | 国产成人精品婷婷| 小蜜桃在线观看免费完整版高清| 99九九线精品视频在线观看视频| 国产毛片a区久久久久| 嫩草影院入口| 永久网站在线| 少妇人妻精品综合一区二区| 街头女战士在线观看网站| 女的被弄到高潮叫床怎么办| av免费在线看不卡| 国产欧美日韩精品一区二区| 一级二级三级毛片免费看| 麻豆成人av视频| 一二三四中文在线观看免费高清| 插逼视频在线观看| 熟女人妻精品中文字幕| 久久99蜜桃精品久久| 欧美成人精品欧美一级黄| 少妇人妻 视频| 国产色爽女视频免费观看| 国产成人福利小说| 麻豆成人av视频| 卡戴珊不雅视频在线播放| 精品国产露脸久久av麻豆| 好男人在线观看高清免费视频| 亚洲精品自拍成人| 国产精品蜜桃在线观看| 少妇人妻久久综合中文| 美女被艹到高潮喷水动态| 青春草亚洲视频在线观看| 日韩,欧美,国产一区二区三区| 六月丁香七月| 亚洲精品aⅴ在线观看| 日韩一区二区三区影片| 国内揄拍国产精品人妻在线| 天堂网av新在线| 国产视频首页在线观看| 大香蕉久久网| 男人和女人高潮做爰伦理| 男人舔奶头视频| 晚上一个人看的免费电影| 美女国产视频在线观看| 少妇人妻一区二区三区视频| 国产久久久一区二区三区| 国产精品人妻久久久影院| 看黄色毛片网站| av在线播放精品| 在线天堂最新版资源| 亚洲国产色片| 亚洲天堂av无毛| 直男gayav资源| 大片电影免费在线观看免费| 国产伦在线观看视频一区| 在线观看免费高清a一片| 天堂中文最新版在线下载 | 三级国产精品欧美在线观看| 性色av一级| 国产精品av视频在线免费观看| 又黄又爽又刺激的免费视频.| 色哟哟·www| 国产探花极品一区二区| 一区二区三区乱码不卡18| 永久免费av网站大全| 久久精品久久久久久久性| 久久99热这里只有精品18| 日韩电影二区| 综合色av麻豆| 国产毛片a区久久久久| 大香蕉97超碰在线| 国产精品福利在线免费观看| 又爽又黄a免费视频| 国产精品三级大全| 欧美激情久久久久久爽电影| 老司机影院成人| 亚洲自拍偷在线| 中文字幕制服av| 日韩大片免费观看网站| 91精品一卡2卡3卡4卡| 毛片一级片免费看久久久久| 肉色欧美久久久久久久蜜桃 | 精品久久久久久电影网| 联通29元200g的流量卡| 啦啦啦中文免费视频观看日本| av.在线天堂| a级毛色黄片| 少妇人妻久久综合中文| 另类亚洲欧美激情| 大片免费播放器 马上看| 亚洲av欧美aⅴ国产| 欧美成人精品欧美一级黄| 国产成年人精品一区二区| 熟女av电影| 18禁在线无遮挡免费观看视频| 两个人的视频大全免费| 女的被弄到高潮叫床怎么办| 午夜激情久久久久久久| 久久鲁丝午夜福利片| 男人和女人高潮做爰伦理| 国产午夜精品久久久久久一区二区三区| 新久久久久国产一级毛片| 97在线视频观看| 伦理电影大哥的女人| av女优亚洲男人天堂| 日韩欧美精品v在线| av在线app专区| 国产亚洲最大av| 亚洲经典国产精华液单| 最近2019中文字幕mv第一页| 成人国产麻豆网| 久久精品熟女亚洲av麻豆精品| 高清av免费在线| 亚洲av中文字字幕乱码综合| 国产黄片视频在线免费观看| 制服丝袜香蕉在线| 少妇熟女欧美另类| 深夜a级毛片| 干丝袜人妻中文字幕| 黑人高潮一二区| 精品少妇黑人巨大在线播放| 国产亚洲精品久久久com| 午夜日本视频在线| 直男gayav资源| 国产亚洲最大av| 99精国产麻豆久久婷婷| 国产精品国产三级国产专区5o| 国产成人免费观看mmmm| 色5月婷婷丁香| 免费观看的影片在线观看| 国产精品福利在线免费观看| 国产一区有黄有色的免费视频| 深爱激情五月婷婷| 一级毛片黄色毛片免费观看视频| 国产亚洲一区二区精品| 国产欧美另类精品又又久久亚洲欧美| 一个人看视频在线观看www免费| 午夜福利视频1000在线观看| 欧美日韩视频精品一区| 亚洲av成人精品一区久久| 日本爱情动作片www.在线观看| 人妻夜夜爽99麻豆av| 国产欧美日韩一区二区三区在线 | 人妻制服诱惑在线中文字幕| 国产人妻一区二区三区在| 亚洲欧美成人综合另类久久久| 成人午夜精彩视频在线观看| 成人无遮挡网站| 国产精品福利在线免费观看| 国产日韩欧美在线精品| 亚洲欧美日韩卡通动漫| 男女啪啪激烈高潮av片| 国产日韩欧美亚洲二区| 国产精品女同一区二区软件| av在线蜜桃| 成年av动漫网址| 亚洲最大成人av| 国产有黄有色有爽视频| 婷婷色av中文字幕| 大片电影免费在线观看免费| 国产黄片美女视频| 99热这里只有是精品在线观看| 亚洲自偷自拍三级| 精品酒店卫生间| 亚洲欧美日韩另类电影网站 | 18禁在线播放成人免费| 亚洲av日韩在线播放| 精品99又大又爽又粗少妇毛片| 91久久精品国产一区二区三区| 全区人妻精品视频| 晚上一个人看的免费电影| 久久久久久久午夜电影| 人体艺术视频欧美日本| 一级av片app| 夫妻午夜视频| 国产男女超爽视频在线观看| 久久精品久久久久久久性| 一级片'在线观看视频| 少妇高潮的动态图| 少妇 在线观看| 各种免费的搞黄视频| 国产精品一区二区三区四区免费观看| 国产精品99久久99久久久不卡 | 久久久久久久久久成人| 丝袜脚勾引网站| 美女国产视频在线观看| 十八禁网站网址无遮挡 | 热99国产精品久久久久久7| 精品久久久久久久久亚洲| 人妻 亚洲 视频| 能在线免费看毛片的网站| 免费观看a级毛片全部| 80岁老熟妇乱子伦牲交| 日韩免费高清中文字幕av| 久久精品国产自在天天线| kizo精华| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品电影| 国产探花极品一区二区| 97精品久久久久久久久久精品| 亚洲精品色激情综合|