• <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光源
    日韩免费av在线播放| 欧美黑人精品巨大| 欧美日韩成人在线一区二区| 中文字幕高清在线视频| av视频免费观看在线观看| 汤姆久久久久久久影院中文字幕| 黄片播放在线免费| 免费在线观看影片大全网站| 久久精品国产亚洲av高清一级| 法律面前人人平等表现在哪些方面| 久久久久久久大尺度免费视频| 51午夜福利影视在线观看| 一本色道久久久久久精品综合| 国产视频一区二区在线看| 天天添夜夜摸| 亚洲 欧美一区二区三区| 少妇粗大呻吟视频| 一级片免费观看大全| 黑人操中国人逼视频| 亚洲成人免费av在线播放| 成在线人永久免费视频| 老司机午夜福利在线观看视频 | 国产在线免费精品| 成人18禁在线播放| 丁香六月天网| 日本av免费视频播放| 国产视频一区二区在线看| 男女免费视频国产| 精品久久久精品久久久| 亚洲中文av在线| 怎么达到女性高潮| 亚洲精品一二三| 夫妻午夜视频| 精品一区二区三卡| 久久久久视频综合| 日本五十路高清| av福利片在线| 好男人电影高清在线观看| 自线自在国产av| 男人舔女人的私密视频| 亚洲色图 男人天堂 中文字幕| 久久九九热精品免费| 性色av乱码一区二区三区2| 午夜福利免费观看在线| 啦啦啦 在线观看视频| 国产av精品麻豆| 久久人妻熟女aⅴ| 狠狠婷婷综合久久久久久88av| 狠狠婷婷综合久久久久久88av| 最近最新免费中文字幕在线| 在线看a的网站| 国产亚洲精品久久久久5区| 久久久久久久久久久久大奶| 夜夜爽天天搞| 亚洲精品美女久久av网站| a级片在线免费高清观看视频| 人妻 亚洲 视频| 亚洲第一欧美日韩一区二区三区 | 又黄又粗又硬又大视频| 亚洲人成77777在线视频| 久久 成人 亚洲| 国产亚洲欧美在线一区二区| 99国产精品99久久久久| 制服人妻中文乱码| 99久久精品国产亚洲精品| 亚洲精华国产精华精| 国产在线观看jvid| 五月天丁香电影| 国产精品一区二区在线不卡| 中亚洲国语对白在线视频| www日本在线高清视频| 国产精品久久久久久人妻精品电影 | 亚洲国产毛片av蜜桃av| 老司机福利观看| 午夜福利一区二区在线看| 久久午夜综合久久蜜桃| 久久久久视频综合| 欧美成人免费av一区二区三区 | 啪啪无遮挡十八禁网站| 99精品久久久久人妻精品| 欧美黄色淫秽网站| 欧美国产精品一级二级三级| www.熟女人妻精品国产| 日日夜夜操网爽| 国产精品熟女久久久久浪| 美国免费a级毛片| 亚洲成人手机| 一本久久精品| 国产精品99久久99久久久不卡| 久久久精品区二区三区| 国产一区二区三区在线臀色熟女 | 国产区一区二久久| 视频在线观看一区二区三区| 丝袜人妻中文字幕| 两性夫妻黄色片| 久久中文字幕人妻熟女| 91成年电影在线观看| 99精国产麻豆久久婷婷| 国产免费福利视频在线观看| 国产主播在线观看一区二区| 一本一本久久a久久精品综合妖精| 国产精品一区二区在线观看99| 男女免费视频国产| 99香蕉大伊视频| 精品国产一区二区久久| 老司机午夜十八禁免费视频| 亚洲人成伊人成综合网2020| 亚洲精品自拍成人| 欧美人与性动交α欧美精品济南到| 搡老熟女国产l中国老女人| 久久人妻熟女aⅴ| 天堂中文最新版在线下载| 国产伦人伦偷精品视频| 丝袜人妻中文字幕| 免费女性裸体啪啪无遮挡网站| 狠狠狠狠99中文字幕| 久久九九热精品免费| 婷婷丁香在线五月| 亚洲精品久久午夜乱码| xxxhd国产人妻xxx| 国产成人免费无遮挡视频| 亚洲精品av麻豆狂野| 久热爱精品视频在线9| 国产亚洲一区二区精品| 国产又爽黄色视频| 亚洲成人国产一区在线观看| 精品久久久久久久毛片微露脸| 视频在线观看一区二区三区| 99精品久久久久人妻精品| 欧美另类亚洲清纯唯美| 免费在线观看日本一区| 亚洲精品久久午夜乱码| av网站在线播放免费| 高潮久久久久久久久久久不卡| 国产在视频线精品| 水蜜桃什么品种好| 丝袜在线中文字幕| 国产成人欧美| 日韩视频在线欧美| 精品一区二区三卡| 久久天堂一区二区三区四区| 精品乱码久久久久久99久播| 十八禁网站网址无遮挡| 亚洲九九香蕉| 中文字幕人妻丝袜制服| 肉色欧美久久久久久久蜜桃| 这个男人来自地球电影免费观看| 免费观看a级毛片全部| 欧美精品av麻豆av| 国产精品98久久久久久宅男小说| 一区福利在线观看| 久久精品成人免费网站| 男女高潮啪啪啪动态图| 午夜成年电影在线免费观看| av有码第一页| av福利片在线| 我的亚洲天堂| 纯流量卡能插随身wifi吗| 国产午夜精品久久久久久| 女警被强在线播放| 亚洲中文日韩欧美视频| 老熟女久久久| 亚洲性夜色夜夜综合| 性少妇av在线| 欧美性长视频在线观看| 成年动漫av网址| 欧美黑人欧美精品刺激| 国产av又大| av一本久久久久| 一级黄色大片毛片| 亚洲久久久国产精品| 最新的欧美精品一区二区| 久久久久精品国产欧美久久久| 国产亚洲av高清不卡| 黄片小视频在线播放| 久久九九热精品免费| 一级毛片女人18水好多| 91字幕亚洲| 亚洲精品粉嫩美女一区| 汤姆久久久久久久影院中文字幕| 亚洲成人手机| 少妇精品久久久久久久| 交换朋友夫妻互换小说| 两个人看的免费小视频| 激情视频va一区二区三区| 国产精品二区激情视频| 99热国产这里只有精品6| 99精品在免费线老司机午夜| 久久久国产欧美日韩av| 男女边摸边吃奶| 欧美中文综合在线视频| 在线观看免费视频日本深夜| 十八禁网站网址无遮挡| 69精品国产乱码久久久| 成人精品一区二区免费| tube8黄色片| 天天操日日干夜夜撸| 亚洲自偷自拍图片 自拍| 狠狠精品人妻久久久久久综合| 91成年电影在线观看| 首页视频小说图片口味搜索| 水蜜桃什么品种好| 亚洲欧美色中文字幕在线| 欧美日韩成人在线一区二区| 黄色毛片三级朝国网站| 91成人精品电影| 久久99一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 成人永久免费在线观看视频 | 国产一区二区三区综合在线观看| 亚洲人成伊人成综合网2020| av不卡在线播放| 国产真人三级小视频在线观看| 99国产精品一区二区蜜桃av | 国产av国产精品国产| 悠悠久久av| 久久精品人人爽人人爽视色| 亚洲精品乱久久久久久| 黄色成人免费大全| 捣出白浆h1v1| 精品国产乱码久久久久久小说| 中文字幕高清在线视频| 国产成人系列免费观看| 嫁个100分男人电影在线观看| 日日爽夜夜爽网站| 日韩制服丝袜自拍偷拍| 亚洲性夜色夜夜综合| 欧美精品高潮呻吟av久久| 美女高潮喷水抽搐中文字幕| 亚洲国产av影院在线观看| 高清欧美精品videossex| 亚洲精品成人av观看孕妇| 国产不卡一卡二| 又大又爽又粗| 亚洲中文日韩欧美视频| 成年动漫av网址| 免费人妻精品一区二区三区视频| e午夜精品久久久久久久| 色播在线永久视频| 欧美 日韩 精品 国产| 国产精品免费一区二区三区在线 | 久久久久网色| 天堂动漫精品| 精品午夜福利视频在线观看一区 | 亚洲全国av大片| 国内毛片毛片毛片毛片毛片| 国产片内射在线| 亚洲熟女毛片儿| 捣出白浆h1v1| 777米奇影视久久| 精品国产超薄肉色丝袜足j| 女人久久www免费人成看片| 麻豆av在线久日| 精品一区二区三区av网在线观看 | 国产av一区二区精品久久| 一区二区三区精品91| 精品少妇久久久久久888优播| 国产xxxxx性猛交| 女警被强在线播放| 另类精品久久| 免费看十八禁软件| 老司机亚洲免费影院| 成年人免费黄色播放视频| 婷婷丁香在线五月| 免费看a级黄色片| 男人操女人黄网站| 欧美成人午夜精品| 久久99一区二区三区| 精品少妇内射三级| 国产一区二区三区综合在线观看| 国产高清国产精品国产三级| 久久国产精品男人的天堂亚洲| 精品国产乱码久久久久久小说| 亚洲精品一卡2卡三卡4卡5卡| 岛国毛片在线播放| 美女扒开内裤让男人捅视频| 日本av手机在线免费观看| 久久人人97超碰香蕉20202| 国产一区二区三区综合在线观看| 亚洲黑人精品在线| 在线观看免费午夜福利视频| 欧美性长视频在线观看| 少妇被粗大的猛进出69影院| 一区二区av电影网| 欧美精品亚洲一区二区| 欧美黑人欧美精品刺激| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久精品古装| 免费观看av网站的网址| 一级毛片精品| 国产一区有黄有色的免费视频| 成年人免费黄色播放视频| 99久久国产精品久久久| 午夜激情av网站| av天堂在线播放| 精品国产一区二区三区四区第35| 亚洲第一av免费看| 国产成人精品在线电影| 免费不卡黄色视频| 亚洲欧美日韩另类电影网站| 黑丝袜美女国产一区| 成人18禁在线播放| 97在线人人人人妻| netflix在线观看网站| 91字幕亚洲| 国产福利在线免费观看视频| 国产精品电影一区二区三区 | 考比视频在线观看| 日本vs欧美在线观看视频| 变态另类成人亚洲欧美熟女 | 大码成人一级视频| 超碰97精品在线观看| 一夜夜www| 精品一区二区三区av网在线观看 | 欧美日韩成人在线一区二区| 一本久久精品| 精品少妇黑人巨大在线播放| 久久久久精品国产欧美久久久| 欧美老熟妇乱子伦牲交| 久久久久久亚洲精品国产蜜桃av| 国产一区二区三区视频了| 亚洲国产精品一区二区三区在线| 精品少妇内射三级| 国产精品久久久人人做人人爽| 男人操女人黄网站| 亚洲欧美精品综合一区二区三区| 久久国产精品大桥未久av| 亚洲欧洲精品一区二区精品久久久| 男女午夜视频在线观看| 99精国产麻豆久久婷婷| 男女免费视频国产| avwww免费| 久久久国产精品麻豆| 欧美在线一区亚洲| av一本久久久久| 成人特级黄色片久久久久久久 | 欧美黑人欧美精品刺激| 欧美日韩视频精品一区| 国产99久久九九免费精品| 女人被躁到高潮嗷嗷叫费观| 精品人妻熟女毛片av久久网站| 免费观看人在逋| 久久精品国产99精品国产亚洲性色 | 免费在线观看影片大全网站| 成人亚洲精品一区在线观看| 老司机午夜福利在线观看视频 | 悠悠久久av| 十八禁网站网址无遮挡| 久久精品91无色码中文字幕| 日韩一卡2卡3卡4卡2021年| 亚洲中文日韩欧美视频| 亚洲第一av免费看| 电影成人av| 国产人伦9x9x在线观看| 在线亚洲精品国产二区图片欧美| 五月开心婷婷网| 两性夫妻黄色片| 亚洲av片天天在线观看| 蜜桃在线观看..| av又黄又爽大尺度在线免费看| 亚洲国产成人一精品久久久| 久久精品成人免费网站| a级毛片黄视频| 国产成人精品在线电影| 精品国内亚洲2022精品成人 | 午夜老司机福利片| 久久国产精品影院| 91av网站免费观看| 日日摸夜夜添夜夜添小说| 日韩欧美一区视频在线观看| 国产单亲对白刺激| 国产日韩欧美视频二区| 91国产中文字幕| 在线 av 中文字幕| 精品少妇黑人巨大在线播放| 久久久精品94久久精品| 丁香欧美五月| 日韩有码中文字幕| 久热爱精品视频在线9| 性高湖久久久久久久久免费观看| 午夜成年电影在线免费观看| 国产精品麻豆人妻色哟哟久久| 麻豆国产av国片精品| 水蜜桃什么品种好| 日本a在线网址| 日日夜夜操网爽| 国产不卡av网站在线观看| 一本—道久久a久久精品蜜桃钙片| 91九色精品人成在线观看| 宅男免费午夜| 中文字幕色久视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲精华国产精华精| 在线观看舔阴道视频| 亚洲熟妇熟女久久| 国产高清国产精品国产三级| 日韩精品免费视频一区二区三区| 亚洲成国产人片在线观看| 一区二区三区激情视频| 一级,二级,三级黄色视频| 一本综合久久免费| 国产97色在线日韩免费| 国产日韩欧美视频二区| 天天操日日干夜夜撸| 午夜福利影视在线免费观看| 免费高清在线观看日韩| 丁香六月欧美| 中文字幕人妻丝袜一区二区| 亚洲国产中文字幕在线视频| 多毛熟女@视频| 国产一区二区三区视频了| 视频在线观看一区二区三区| 久久精品亚洲熟妇少妇任你| 国产一区二区三区在线臀色熟女 | 香蕉丝袜av| 国产精品电影一区二区三区 | 国产日韩欧美亚洲二区| 人人妻人人澡人人看| 精品国产超薄肉色丝袜足j| 亚洲精品自拍成人| 午夜激情av网站| 狠狠狠狠99中文字幕| 香蕉久久夜色| 叶爱在线成人免费视频播放| www.熟女人妻精品国产| 美女高潮喷水抽搐中文字幕| 国产1区2区3区精品| 亚洲精品久久成人aⅴ小说| 十分钟在线观看高清视频www| 欧美av亚洲av综合av国产av| 国产精品久久电影中文字幕 | 国产欧美日韩一区二区三| av天堂久久9| 国产高清激情床上av| 黑人猛操日本美女一级片| 80岁老熟妇乱子伦牲交| 淫妇啪啪啪对白视频| 日本黄色日本黄色录像| 王馨瑶露胸无遮挡在线观看| 亚洲av第一区精品v没综合| 国产老妇伦熟女老妇高清| 亚洲黑人精品在线| 动漫黄色视频在线观看| 热99re8久久精品国产| 国产深夜福利视频在线观看| 少妇的丰满在线观看| 精品国产一区二区三区四区第35| 91大片在线观看| 人人妻人人爽人人添夜夜欢视频| 18禁国产床啪视频网站| 精品少妇久久久久久888优播| 性高湖久久久久久久久免费观看| 国产精品久久久久成人av| 久久久久久久久久久久大奶| 不卡一级毛片| 日本wwww免费看| 在线观看人妻少妇| 免费久久久久久久精品成人欧美视频| 成人特级黄色片久久久久久久 | 丰满人妻熟妇乱又伦精品不卡| 久久久久久久大尺度免费视频| 女同久久另类99精品国产91| 高清视频免费观看一区二区| 免费av中文字幕在线| 国产精品香港三级国产av潘金莲| 成人国语在线视频| 午夜视频精品福利| 一区二区三区精品91| 亚洲欧美激情在线| 18在线观看网站| 亚洲专区国产一区二区| cao死你这个sao货| 国产精品国产高清国产av | 少妇被粗大的猛进出69影院| 亚洲成国产人片在线观看| 成人国产一区最新在线观看| 亚洲人成电影观看| 国产精品久久久久久精品电影小说| 90打野战视频偷拍视频| 国精品久久久久久国模美| 宅男免费午夜| 婷婷丁香在线五月| 亚洲专区中文字幕在线| 老鸭窝网址在线观看| 亚洲欧美一区二区三区久久| www.999成人在线观看| 成人影院久久| 十八禁高潮呻吟视频| 亚洲精品自拍成人| 欧美日韩av久久| 亚洲人成77777在线视频| 亚洲欧美日韩高清在线视频 | 国产精品亚洲av一区麻豆| 日本黄色视频三级网站网址 | 亚洲av日韩在线播放| 免费在线观看日本一区| 国产精品国产av在线观看| 久久久久精品国产欧美久久久| 成年人免费黄色播放视频| 无限看片的www在线观看| 纵有疾风起免费观看全集完整版| 在线观看www视频免费| 老司机午夜十八禁免费视频| 考比视频在线观看| 99国产精品99久久久久| 亚洲黑人精品在线| 中文字幕制服av| 91麻豆av在线| 黄片小视频在线播放| 精品卡一卡二卡四卡免费| 久久亚洲真实| 老熟妇仑乱视频hdxx| e午夜精品久久久久久久| 午夜福利一区二区在线看| 黄片播放在线免费| 国产成人精品在线电影| 美女国产高潮福利片在线看| 一区二区三区精品91| 97在线人人人人妻| 国产高清国产精品国产三级| 国产一区二区三区在线臀色熟女 | 搡老熟女国产l中国老女人| av网站免费在线观看视频| 蜜桃国产av成人99| 精品欧美一区二区三区在线| 一二三四社区在线视频社区8| 99九九在线精品视频| 午夜福利在线观看吧| 三上悠亚av全集在线观看| 久久久国产一区二区| 亚洲精品美女久久av网站| 欧美性长视频在线观看| 波多野结衣一区麻豆| 搡老岳熟女国产| 人人妻人人添人人爽欧美一区卜| 天堂中文最新版在线下载| 国产日韩欧美视频二区| e午夜精品久久久久久久| 1024香蕉在线观看| 搡老熟女国产l中国老女人| 亚洲国产欧美一区二区综合| 自线自在国产av| 欧美国产精品一级二级三级| 亚洲精品国产精品久久久不卡| www日本在线高清视频| 满18在线观看网站| 亚洲欧美精品综合一区二区三区| 香蕉国产在线看| 一个人免费看片子| 亚洲欧美日韩高清在线视频 | 激情视频va一区二区三区| 操美女的视频在线观看| 热99久久久久精品小说推荐| 丰满少妇做爰视频| 亚洲五月婷婷丁香| 亚洲国产毛片av蜜桃av| 91老司机精品| 免费观看av网站的网址| 美女视频免费永久观看网站| 久久久国产一区二区| 美女高潮到喷水免费观看| 下体分泌物呈黄色| 操出白浆在线播放| 在线永久观看黄色视频| 一区二区三区乱码不卡18| 91av网站免费观看| 黄色 视频免费看| 亚洲avbb在线观看| 国产三级黄色录像| 亚洲欧美日韩另类电影网站| 桃花免费在线播放| 十八禁网站免费在线| 美女国产高潮福利片在线看| 人人妻人人澡人人看| 精品视频人人做人人爽| 啦啦啦视频在线资源免费观看| 最近最新免费中文字幕在线| 十八禁高潮呻吟视频| 国产在线精品亚洲第一网站| 最近最新中文字幕大全免费视频| 亚洲国产成人一精品久久久| 高清在线国产一区| 日韩欧美免费精品| 欧美日韩黄片免| 777米奇影视久久| 亚洲中文av在线| 亚洲精品成人av观看孕妇| 一本色道久久久久久精品综合| 男女床上黄色一级片免费看| 亚洲午夜精品一区,二区,三区| 国产精品秋霞免费鲁丝片| 国产精品免费一区二区三区在线 | 女同久久另类99精品国产91| 人人妻人人澡人人看| 欧美老熟妇乱子伦牲交| 国产精品熟女久久久久浪| 精品国产乱码久久久久久小说| 两个人免费观看高清视频| 亚洲av成人一区二区三| 成年人黄色毛片网站| av欧美777| 午夜日韩欧美国产| 成人手机av| 国产日韩欧美亚洲二区| 亚洲精品国产一区二区精华液| 欧美激情极品国产一区二区三区| 侵犯人妻中文字幕一二三四区| 波多野结衣一区麻豆| 免费黄频网站在线观看国产| 一级毛片女人18水好多| 麻豆av在线久日| 国产成人系列免费观看| 99国产精品99久久久久| 成人手机av| 精品人妻在线不人妻|