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

    光譜特征波長的SPA選取和基于SVM的玉米顆粒霉變程度定性判別

    2016-06-15 16:35:24喜明杰
    光譜學(xué)與光譜分析 2016年1期
    關(guān)鍵詞:波長預(yù)處理光譜

    袁 瑩,王 偉,褚 璇,喜明杰

    中國農(nóng)業(yè)大學(xué)工學(xué)院,現(xiàn)代農(nóng)業(yè)裝備優(yōu)化設(shè)計北京市重點實驗室,北京 100083

    光譜特征波長的SPA選取和基于SVM的玉米顆粒霉變程度定性判別

    袁 瑩,王 偉*,褚 璇,喜明杰

    中國農(nóng)業(yè)大學(xué)工學(xué)院,現(xiàn)代農(nóng)業(yè)裝備優(yōu)化設(shè)計北京市重點實驗室,北京 100083

    利用波長范圍在833~2 500 nm的傅里葉變換近紅外光譜(Fourier transform near infrared spectroscopy, FT-NIR)對不同霉變程度的玉米顆粒進行檢測區(qū)分。首先,為避免光譜數(shù)據(jù)首尾噪聲影響,對比四種常見的預(yù)處理方法,最終選擇移動平均平滑法對原始光譜數(shù)據(jù)進行預(yù)處理;然后為選出合適的樣本集劃分方法以提高模型預(yù)測性能,對常見的四種方法進行對比,最終利用SPXY(sample set partitioning based on joint x-y distance)法進行樣本集劃分;進一步為減少數(shù)據(jù)量,降低維度,使用連續(xù)投影算法(successive projections algorithm, SPA)提取出7個特征波長,分別為833,927,1 208,1 337,1 454,1 861和2 280 nm;最后,將七個特征波長數(shù)據(jù)作為輸入,選取徑向基函數(shù)(radial basis function, RBF)作為支持向量機(support vector machine, SVM)核函數(shù),取參數(shù)C=7 760 469,γ=0.017 003建立判別模型。SVM模型對訓(xùn)練集和測試集的預(yù)測準確率分別達到97.78%和93.33%。另取不同品種的玉米顆粒,以同樣的標準挑選樣品組成獨立驗證集,所建立的判別模型對獨立驗證集的預(yù)測準確率達到91.11%。結(jié)果表明基于SPA和SVM能有效地對玉米顆粒霉變程度進行判別,所選取的7個特征波長為實現(xiàn)在線霉變玉米顆粒近紅外檢測提供了理論依據(jù)。

    霉變;玉米顆粒;FT-NIR;SPA;SVM

    引 言

    我國是玉米生產(chǎn)和消費大國,玉米除食用之外,還用作飼料和工業(yè)原料。然而目前因霉變造成的玉米產(chǎn)后損失過于嚴重,這不僅影響著玉米產(chǎn)量,而且一旦霉變玉米被人畜食用,也會嚴重危害人畜健康,因此,及時檢出霉變玉米尤為重要。

    目前,對霉變玉米的檢出大多采用人工感官鑒定方法進行品質(zhì)檢測,不僅效率低下,工作量大,而且對于品質(zhì)的檢測不能標準化,同時雜質(zhì)、霉變、蟲蝕等缺陷有時不易進行肉眼的判斷;而傳統(tǒng)的生物培養(yǎng)方法需要破壞大量樣本,不利于大群體篩選,而且測定程序過于復(fù)雜不能及時分析等。

    近紅外光譜技術(shù)是近些年來發(fā)展較快的分析方法之一,它依據(jù)化學(xué)成分對近紅外光譜的吸收特性而進行測定。其優(yōu)點在于無需復(fù)雜的樣品前處理,是一種快速無損的分析檢測方法,結(jié)合光纖可以被運用到在線檢測中,因此廣泛應(yīng)用于谷物品質(zhì)和營養(yǎng)的定性及定量檢測[1, 2]。許多研究表明利用近紅外結(jié)合適當(dāng)?shù)幕瘜W(xué)計量學(xué)方法可識別谷物的品種[3],并對其內(nèi)部的營養(yǎng)物質(zhì)進行定性定量的檢測,而且谷物顆粒的霉變檢測也已經(jīng)成為西方發(fā)達國家的研究熱點[4-5]。

    本工作在運用移動平均平滑法進行光譜預(yù)處理和運用SPXY(sample set partitioning based on joint x-y distance)法進行樣本集劃分的基礎(chǔ)上,采用連續(xù)投影算法(successive projections algorithm,SPA)[6]提取光譜中的特征波長,降低了光譜信息的維度,剔除了冗余信息,提取出的特征波長作為支持向量機(support vector machine, SVM)的輸入用于建立不同霉變程度玉米顆粒分類模型,探索玉米顆粒霉變程度近紅外檢測方法的可行性。

    1 實驗部分

    1.1 儀器設(shè)備

    實驗使用儀器為德國布魯克公司(BRUKER,德國)的MPA型傅里葉變換近紅外光譜儀,其采集模式為積分球反射模式,波長范圍為833~2 500 nm(12 000~4 000 cm-1),分辨率為4 cm-1。利用光譜儀自帶的OPUS軟件進行光譜數(shù)據(jù)獲取。并利用Matlab2012b平臺結(jié)合臺灣大學(xué)林智仁(C.J Lin)等開發(fā)的LIBSVM工具箱對光譜數(shù)據(jù)進行處理分析。

    1.2 材料和樣品

    選取150粒2012年收獲的自然感染霉菌的豫玉32玉米顆粒作為樣品,其形狀和大小基本一致。參照文獻[7]的方法進行分組,每組挑選50粒玉米顆粒組成樣品集:

    (1)無癥狀組:玉米顆粒表面無明顯的霉菌損害的跡象,除了非常輕微的變色。

    (2)中度霉變組:可見的真菌生長占玉米顆粒表面的30%~70%,且嚴重變色。

    (3)重度霉變組:真菌生長覆蓋整個嚴重變色玉米顆粒的表面。

    1.3 光譜數(shù)據(jù)預(yù)處理

    由于樣品形狀不同及表面凹凸不平,由此產(chǎn)生的散射會對光譜的采集造成一定的影響,并且由光譜儀得到的光譜信號中不僅含有有用信息,同時也疊加著隨機誤差,使得光譜曲線中存在噪聲干擾,因而需要對原始光譜數(shù)據(jù)進行預(yù)處理。常見的預(yù)處理方法有標準化(autoscaling)、移動平均平滑法(moving average smoothing)、局部加權(quán)回歸散點平滑法(locally weighted scatterplot smoothing,LOWESS)、S-G卷積平滑法(Savitzky-Golay smoothing)等。

    標準化,可以給光譜中所有波長變量以相同的權(quán)重,但是不能反映光譜的趨勢及變化。移動平均平滑法,優(yōu)點為計算量少,能較好的反映光譜的趨勢及變化,但是當(dāng)預(yù)處理數(shù)據(jù)較大時,存儲的數(shù)據(jù)量也較大,且平滑點數(shù)對預(yù)處理的影響較大。LOWESS平滑法,其特點是使得線性回歸模型中的參數(shù)可以隨著自變量的不同取值改變,但是也從而加大了計算量。S-G卷積平滑法,可以有效減少光譜數(shù)據(jù)中的隨機噪聲,平滑點數(shù)的影響較小,但計算相對復(fù)雜。

    為選擇適合本研究光譜數(shù)據(jù)的預(yù)處理方法,利用Matlab2012b平臺對四種常見的光譜數(shù)據(jù)預(yù)處理方法進行對比。取各組的前30個樣品,共90個樣品數(shù)據(jù)組成訓(xùn)練集,取各組剩余20個樣品共60個樣品數(shù)據(jù)組成測試集,分別利用標準化、移動平均平滑法、LOWESS平滑法以及S-G卷積平滑法四種方法進行預(yù)處理,并分別將預(yù)處理后的數(shù)據(jù)利用SVM進行初期的預(yù)測,對比預(yù)測結(jié)果,從而選擇出適合的預(yù)處理方法。

    1.4 樣本集劃分

    樣本集劃分即采用某種特定方法將總體樣本劃分為訓(xùn)練集和測試集,用于進行模型建立及對模型的預(yù)測能力進行驗證,它在一定程度上決定著所建模型的預(yù)測性能。在進行樣本集劃分時,數(shù)據(jù)需具有一定的代表性,并且要保證訓(xùn)練集和測試集樣本數(shù)量的合理分配,既要滿足建模的要求,也不會引起數(shù)據(jù)冗余或模型過擬合。常見的樣本集劃分方法主要包括隨機法(random sampling, RS), Kennard-Stone算法(KS)及SPXY算法等。

    隨機法,即按照簡單規(guī)則或完全無規(guī)律劃分樣本集,這種方法簡單、速度快,但是很難得到較為理想的樣本集,且每次劃分都是隨機的,訓(xùn)練集和測試集的樣本差異大,從而會影響模型的穩(wěn)定性。KS法是基于變量之間的歐氏距離,在特征空間中均勻選取樣本。該方法所選的訓(xùn)練集樣本分布均勻,但因需要計算兩兩樣本之間的歐式距離,計算量大。SPXY法是Galvao等[8]在KS法的基礎(chǔ)上提出的,其逐步選擇過程與KS法相同,只是在計算過程中,將變量x和y均考慮在內(nèi),這樣有效兼顧了樣本的光譜和待測量信息。

    為選擇出預(yù)測性能好的樣本集劃分方法,本文利用Matlab2012b平臺對順序劃分、隨機法、KS法以及SPXY法四種方法進行樣本集劃分方法的對比。

    1.5 SPA

    SPA是一種前向循環(huán)波長選取方法,從一個波長開始,每次循環(huán)都計算其在未選入波長上的投影,將投影向量最大的波長引入到波長組合中,每個新選入的波長,都與前一個線性關(guān)系最小。這種方法能夠有效地消除波長變量之間的共線性影響,降低樣品光譜數(shù)據(jù)維數(shù),減少計算量,在多種光譜的多元定量和定性分析中得到了廣泛的應(yīng)用。

    對于樣品數(shù)M和波長數(shù)K組成的光譜矩陣XM×K,N(N

    (1)第一次迭代(p=1)開始前,在光譜矩陣中任選一列向量xj,記為xk(0),即k(0)=j,j∈1, …,m;

    (2)將未被選人的列向量位置的集合記為s,s={j, 1≤j≤m,j?{k(0), …,k(p-1)}};

    (3)計算剩余列向量xj(j∈s)與當(dāng)前所選向量xk(p-1)的投影:

    (4)提取投影向量最大波長的變量序號:

    k(p)=arg[max(‖Pxj‖)],j∈s;

    (5)令xj=Pxj,j∈s;p=p+1,如果p

    最終選取的波長變量組合為{k(p),p=0, …,h-1}。對于每一個初始k(0),循環(huán)一次后利用多元線性回歸(multple linear regresion,MLR)進行交互驗證分析, 最小均方根誤差(root mean square error,RMSE)對應(yīng)的k(p)即為最終的選擇結(jié)果[11]。

    2 結(jié)果與討論

    2.1 原始光譜數(shù)據(jù)獲取

    利用傅里葉變換近紅外光譜儀采集每粒玉米胚芽面的光譜數(shù)據(jù),最終得到的光譜數(shù)據(jù)為掃描64次所得的平均光譜,所有樣品的光譜曲線如圖1所示。從圖中可以看出,無癥狀組和重度霉變組基本可以分開,而中度霉變組與其他兩組略有重疊,因而需要運用化學(xué)計量學(xué)方法進行分類判別,樣品的原始光譜結(jié)果與1.2節(jié)分組基本一致。

    2.2 光譜預(yù)處理結(jié)果

    由于原始光譜曲線在波長1 000 nm之前存在噪聲干擾,同時為保留光譜曲線中的有用信息,分別利用標準化、移動平均平滑法、LOWESS平滑法及S-G卷積平滑法對光譜數(shù)據(jù)進行預(yù)處理,取各組的前30個樣品,共90個樣品數(shù)據(jù)組成訓(xùn)練集,取各組剩余20個樣品共60個樣品數(shù)據(jù)組成測試集,預(yù)處理后得到的數(shù)據(jù)分別作為SVM的輸入進行初步的預(yù)測,對比預(yù)測結(jié)果,從而選取合適的預(yù)處理方法,便于后續(xù)的數(shù)據(jù)分析。四種預(yù)處理方法對應(yīng)的訓(xùn)練集及測試集的預(yù)測準確率如表1所示。

    Fig.1 Original spectra

    Table 1 Prediction accuracy got by different preprocessing methods/%

    由表1可知,采用移動平均平滑法對光譜數(shù)據(jù)進行預(yù)處理,訓(xùn)練集及測試集的預(yù)測結(jié)果都較高。通過試驗,最終采用21個窗口的移動平均平滑法對樣本數(shù)據(jù)進行預(yù)處理。

    2.3 樣本集劃分方法選擇

    由于樣本集的劃分在一定程度上決定著所建模型的預(yù)測性能,為選出合適的樣本集劃分方法,對常見的四種方法進行對比。共采用了150個玉米顆粒樣品,取其中的90個樣品作為訓(xùn)練集,60個樣品作為測試集。使用2.2節(jié)中最終選擇的21個窗口的移動平均平滑法對樣本數(shù)據(jù)進行預(yù)處理,然后分別利用順序劃分、隨機法、KS法以及SPXY法四種方法進行樣本集劃分。其中順序劃分即取各組的前30個樣品數(shù)據(jù)組成訓(xùn)練集,取各組剩余20個樣品數(shù)據(jù)組成測試集。劃分后得到的訓(xùn)練集和測試集輸入到SVM進行劃分效果的對比,四種樣本集劃分方法的預(yù)測準確率如表2所示。

    Table 2 Prediction accuracy got by different grouping method

    利用SPXY法進行樣本集劃分時,訓(xùn)練集及測試集的預(yù)測準確率均高于其他方法,因此,采用SPXY法對樣本集進行劃分。運用SPXY法的樣本集劃分結(jié)果如表3所示。從表中可以看出,各組的訓(xùn)練集測試集的樣品數(shù)差異不大。

    Table 3 Result of grouping by SPXY

    2.4 波長優(yōu)選

    樣本數(shù)據(jù)共含有2 100個光譜數(shù)據(jù)點,由于數(shù)據(jù)量大,運算時間長,且存在多重共線性,在進行分類模型建立時,會降低模型的準確性。因此利用SPA對樣本數(shù)據(jù)進行降維,提取出樣本數(shù)據(jù)中重要的波長點。

    對預(yù)處理后的數(shù)據(jù)用SPA進行波長的優(yōu)選,然后利用優(yōu)選的波長數(shù)據(jù)建立MLR校正模型,最小RMSE值對應(yīng)的波長變量個數(shù)即為最終的選擇結(jié)果。圖2為提取不同個數(shù)波長變量時RMSE的走勢圖,從圖中可知,提取7個波長變量時模型的RMSE最小為0.246 6。

    Fig.2 RMSE in different wavelength variables

    圖3表示使用SPA篩選得到的7個特征波長,在圖中以白色方塊表示。7個特征波長分別為833,927,1 208,1 337,1 454,1 861和2 280 nm。其中,1 208 nm代表C—H二階倍頻伸縮,即CH2;1 454 nm代表O-H一階倍頻伸縮,表明存在淀粉成分;1 861 nm代表C—Cl六階倍頻伸縮;2 280 nm代表C—H伸縮和CH2變形組合,同樣表明存在淀粉成分。

    Fig.3 Optimal characteristis wavelengths selected by SPA

    2.5 SVM模型驗證

    以SPA提取的7個波長數(shù)據(jù)作為輸入,利用SVM進行分類模型的建立和驗證。利用RBF函數(shù)作為SVM的核函數(shù),在懲罰參數(shù)C=7 760 469,核函數(shù)參數(shù)γ=0.017 003的條件下,建立分類模型。利用建立好的模型對訓(xùn)練集和測試集進行預(yù)測,預(yù)測結(jié)果如表4所示。

    可以看出訓(xùn)練集的90個樣品中,有1粒重度霉變的顆粒被錯判為無癥狀組,1粒被錯判為中度霉變,判別準確率為97.78%。測試集的60個樣品中,有1粒無癥狀組的顆粒被錯判為中度霉變,1粒中度霉變組的顆粒被錯判為重度霉變,2粒重度霉變的顆粒被錯判為無癥狀組,判別準確率為93.33%。分析錯判原因,在對顆粒進行分組挑選時,可能會由于部分霉變無法用肉眼識別,從而造成錯分。

    Table 5 Classification result

    Table 5 Predictive result of the independent validation set

    同時,為了驗證模型的通用性,另取2013年收獲的先玉335玉米顆粒,并以相同的分組標準每組挑選30粒共90粒作為獨立驗證集,利用21個窗口的移動平均平滑法進行預(yù)處理,然后將7個特征波長的信息輸入到已建立的判別模型中,判別結(jié)果如表5所示。可以看出,獨立驗證集的預(yù)測準確率達到91.11%,表明本文所建立的判別模型具有通用性。

    3 結(jié) 論

    按照霉變程度將同年收獲的自然感染霉菌的豫玉32玉米顆粒分成三組,運用波長范圍在833~2 500 nm的FT-NIR系統(tǒng)進行光譜數(shù)據(jù)的采集,首先采用21個窗口的移動平均平滑法對樣本數(shù)據(jù)進行預(yù)處理,并利用SPXY算法對樣本集進行劃分;然后為了對光譜數(shù)據(jù)進行降維,利用SPA篩選特征波長,選取7個特征波長數(shù)據(jù)作為輸入,利用SVM進行分類模型的建立和驗證。最終建立的模型對訓(xùn)練集及測試集的預(yù)測準確率分別達到97.78%和93.33%。另取2013年收獲的先玉335玉米顆粒,以同樣的標準挑選90粒作為獨立驗證集,按照相同的方法進行預(yù)處理,所建立的判別模型對獨立驗證集的預(yù)測準確率達到91.11%。結(jié)果表明基于SPA和SVM能有效地對霉變玉米顆粒進行判別,所選取的7個特征波長為在線霉變玉米顆粒近紅外檢測提供了理論依據(jù)。

    [2] LIU Xin-ru, ZHANG Li-ping, WANG Jian-fu, et al(劉心如, 張黎平, 王建福, 等). Spectroscopy and Spectral Analysis(光譜學(xué)與光譜分析), 2013, 33(8): 2092.

    [3] LIANG Jian, LIU Bin-mei, TAO Liang-zhi, et al(梁 劍,劉斌美,陶亮之,等). The Journal of Light Scattering(光散射學(xué)報), 2013, 25(4): 423.

    [4] Singh C B, Jayas D S, Paliwal J, et al. International Journal of Food Properties, 2012, 15(1): 11.

    [5] YAO Hai-bo, Hruska Z, Kincaid R, et al. Biosystems Engineering, 2013, 115: 125.

    [6] Bregman L M. Akademiia Nauk SSSR, Doklady, 1965, 162: 487.

    [7] Tallada J G, Wicklow D T, Pearson T C, et al. Transactions of the ASABE, 2011, 54(3): 1151.

    [8] Galv?o R K H, Araujo M C U, José G E, et al. Talanta, 2005, 67(4): 736.

    [9] CHU Xiao-li(褚小立). Molecular Spectroscopy Analytical Technology Combined with Chemometrics and Its Applications(化學(xué)計量學(xué)方法與分子光譜分析技術(shù)). Beijing: Chemical Industry Press(北京: 化學(xué)工業(yè)出版社), 2011.

    *Corresponding author

    Selection of Characteristic Wavelengths Using SPA and Qualitative Discrimination of Mildew Degree of Corn Kernels Based on SVM

    YUAN Ying, WANG Wei*, CHU Xuan, XI Ming-jie

    Beijing Key Laboratory of Optimization Design for Modern Agricultural Equipment, College of Engineering, China Agricultural University, Beijing 100083, China

    The feasibility of Fourier transform near infrared (FT-NIR) spectroscopy with spectral range between 833 and 2 500 nm to detect the moldy corn kernels with different levels of mildew was verified in this paper. Firstly, to avoid the influence of noise, moving average smoothing was used for spectral data preprocessing after four common pretreatment methods were compared. Then to improve the prediction performance of the model, SPXY (sample set partitioning based on joint x-y distance) was selected and used for sample set partition. Furthermore, in order to reduce the dimensions of the original spectral data, successive projection algorithm (SPA) was adopted and ultimately 7 characteristic wavelengths were extracted, the characteristic wavelengths were 833, 927, 1 208, 1 337, 1 454, 1 861, 2 280 nm. The experimental results showed when the spectrum data of the 7 characteristic wavelengths were taken as the input of SVM, the radial basic function (RBF) used as the kernel function, and kernel parameterC=7 760 469,γ=0.017 003, the classification accuracies of the established SVM model were 97.78% and 93.33% for the training and testing sets respectively. In addition, the independent validation set was selected in the same standard, and used to verify the model. At last, the classification accuracy of 91.11% for the independent validation set was achieved. The result indicated that it is feasible to identify and classify different degree of moldy corn grain kernels using SPA and SVM, and characteristic wavelengths selected by SPA in this paper also lay a foundation for the online NIR detection of mildew corn kernels.

    Mildew grain;Corn kernels;FT-NIR;SPA;SVM

    Aug. 21, 2014; accepted Dec. 5, 2014)

    2014-08-21,

    2014-12-05

    國家科技支撐計劃項目(2012BAK08B04)資助

    袁 瑩,女,1991年生,中國農(nóng)業(yè)大學(xué)機械工程碩士研究生 e-mail:410708790@qq.com *通訊聯(lián)系人 e-mail: playerwxw@cau.edu.cn

    S132

    A

    10.3964/j.issn.1000-0593(2016)01-0226-05

    猜你喜歡
    波長預(yù)處理光譜
    HPLC-PDA雙波長法同時測定四季草片中沒食子酸和槲皮苷的含量
    基于三維Saab變換的高光譜圖像壓縮方法
    基于預(yù)處理MUSIC算法的分布式陣列DOA估計
    雙波長激光治療慢性牙周炎的療效觀察
    日本研發(fā)出可完全覆蓋可見光波長的LED光源
    中國照明(2016年4期)2016-05-17 06:16:15
    淺談PLC在預(yù)處理生產(chǎn)線自動化改造中的應(yīng)用
    星載近紅外高光譜CO2遙感進展
    絡(luò)合萃取法預(yù)處理H酸廢水
    便攜式多用途光波波長測量儀
    物理實驗(2015年9期)2015-02-28 17:36:46
    基于自適應(yīng)預(yù)處理的改進CPF-GMRES算法
    国产精品98久久久久久宅男小说| 看片在线看免费视频| 午夜福利在线观看免费完整高清在 | 高清日韩中文字幕在线| 亚洲五月婷婷丁香| 在线观看免费午夜福利视频| 啦啦啦韩国在线观看视频| 婷婷精品国产亚洲av在线| 免费人成视频x8x8入口观看| 看免费av毛片| av中文乱码字幕在线| 成人午夜高清在线视频| 久久久久久久亚洲中文字幕 | 色吧在线观看| 一进一出抽搐gif免费好疼| 亚洲男人的天堂狠狠| 精品久久久久久久毛片微露脸| 久久精品91蜜桃| 国产精品久久久久久久久免 | 亚洲欧美日韩东京热| 高清日韩中文字幕在线| 免费看日本二区| 日韩中文字幕欧美一区二区| 91麻豆精品激情在线观看国产| 成人av在线播放网站| 午夜影院日韩av| 夜夜躁狠狠躁天天躁| 两个人视频免费观看高清| 久久精品人妻少妇| 国产综合懂色| av女优亚洲男人天堂| 国产成人a区在线观看| 亚洲国产精品999在线| 亚洲av免费高清在线观看| 丰满乱子伦码专区| 国产久久久一区二区三区| 久久久久久九九精品二区国产| 久久亚洲真实| 欧美日韩国产亚洲二区| 亚洲av美国av| 国产精品国产高清国产av| 久久婷婷人人爽人人干人人爱| 一个人观看的视频www高清免费观看| 国产成人啪精品午夜网站| 嫩草影视91久久| 最后的刺客免费高清国语| 欧美日本视频| 99精品欧美一区二区三区四区| 国产淫片久久久久久久久 | 在线观看66精品国产| 国产精品精品国产色婷婷| 亚洲男人的天堂狠狠| 嫩草影视91久久| 日韩欧美免费精品| 国产精品一区二区免费欧美| 日本与韩国留学比较| 一本久久中文字幕| 欧美精品啪啪一区二区三区| 国产黄a三级三级三级人| 亚洲avbb在线观看| 亚洲国产高清在线一区二区三| 90打野战视频偷拍视频| 99久久精品国产亚洲精品| 成年女人看的毛片在线观看| 免费高清视频大片| 国产在视频线在精品| 日韩 欧美 亚洲 中文字幕| 99久久精品国产亚洲精品| 国内精品美女久久久久久| 亚洲精品久久国产高清桃花| 99久久精品一区二区三区| 高清日韩中文字幕在线| 99久久综合精品五月天人人| 欧美+亚洲+日韩+国产| 一边摸一边抽搐一进一小说| 亚洲欧美日韩东京热| av在线天堂中文字幕| 免费av观看视频| 日本精品一区二区三区蜜桃| 国产精品 欧美亚洲| 国产成人系列免费观看| 久久久国产成人精品二区| 日本五十路高清| 美女 人体艺术 gogo| 在线观看av片永久免费下载| 国产淫片久久久久久久久 | 国产男靠女视频免费网站| 亚洲国产色片| 每晚都被弄得嗷嗷叫到高潮| 一本久久中文字幕| 精品国产超薄肉色丝袜足j| 最后的刺客免费高清国语| 国产视频一区二区在线看| 在线观看美女被高潮喷水网站 | 搡女人真爽免费视频火全软件 | 精品久久久久久,| 一本综合久久免费| 97超视频在线观看视频| 夜夜躁狠狠躁天天躁| 90打野战视频偷拍视频| 国产又黄又爽又无遮挡在线| 男女做爰动态图高潮gif福利片| 少妇人妻精品综合一区二区 | 精品一区二区三区视频在线 | 桃红色精品国产亚洲av| 小蜜桃在线观看免费完整版高清| 国产精品影院久久| 亚洲中文日韩欧美视频| 亚洲成人久久爱视频| 亚洲成人中文字幕在线播放| 成人一区二区视频在线观看| 久久久久国内视频| 神马国产精品三级电影在线观看| 麻豆一二三区av精品| 天堂√8在线中文| 欧美日韩黄片免| 午夜免费男女啪啪视频观看 | 性色avwww在线观看| 免费电影在线观看免费观看| 在线观看免费视频日本深夜| 亚洲av电影不卡..在线观看| 欧美激情在线99| 三级男女做爰猛烈吃奶摸视频| 国产日本99.免费观看| 麻豆国产97在线/欧美| 久久精品人妻少妇| 色综合婷婷激情| 日韩亚洲欧美综合| 国产精品嫩草影院av在线观看 | 一进一出抽搐gif免费好疼| 色视频www国产| 18禁在线播放成人免费| 18禁黄网站禁片免费观看直播| 国产伦在线观看视频一区| 最近在线观看免费完整版| 俄罗斯特黄特色一大片| 美女高潮喷水抽搐中文字幕| 精品人妻1区二区| 亚洲无线观看免费| 99久久精品一区二区三区| 99热这里只有是精品50| 国产成人欧美在线观看| 欧美色视频一区免费| 99国产综合亚洲精品| 五月玫瑰六月丁香| 男女做爰动态图高潮gif福利片| 99国产极品粉嫩在线观看| 日本熟妇午夜| 天堂影院成人在线观看| 成年人黄色毛片网站| 在线视频色国产色| xxx96com| 亚洲18禁久久av| 国内少妇人妻偷人精品xxx网站| 国产精品野战在线观看| 成年女人毛片免费观看观看9| 国产成人影院久久av| 99国产综合亚洲精品| 国产高清视频在线观看网站| 精华霜和精华液先用哪个| 99久久久亚洲精品蜜臀av| 99久久精品一区二区三区| 亚洲欧美日韩高清专用| 夜夜爽天天搞| 日本免费a在线| 男女床上黄色一级片免费看| 搡女人真爽免费视频火全软件 | 欧美+日韩+精品| 国产精品 欧美亚洲| 欧美色欧美亚洲另类二区| 性色avwww在线观看| 亚洲黑人精品在线| 欧美日韩中文字幕国产精品一区二区三区| 黄色女人牲交| 色在线成人网| 日本免费一区二区三区高清不卡| 国产伦一二天堂av在线观看| 亚洲在线自拍视频| 在线观看日韩欧美| 久久久久九九精品影院| 舔av片在线| 波多野结衣高清作品| 欧美在线黄色| 国产野战对白在线观看| 国产精品影院久久| 变态另类成人亚洲欧美熟女| 无限看片的www在线观看| 免费观看精品视频网站| xxx96com| 精品欧美国产一区二区三| 日韩 欧美 亚洲 中文字幕| 亚洲专区中文字幕在线| 精品电影一区二区在线| 日韩免费av在线播放| www日本在线高清视频| 亚洲中文字幕日韩| 99久久久亚洲精品蜜臀av| 一本综合久久免费| 国产亚洲精品综合一区在线观看| 最新在线观看一区二区三区| 精品无人区乱码1区二区| 国产伦人伦偷精品视频| 国产欧美日韩精品亚洲av| 国产国拍精品亚洲av在线观看 | 最新美女视频免费是黄的| 久久午夜亚洲精品久久| 午夜亚洲福利在线播放| 国产精品av视频在线免费观看| 国产成人av激情在线播放| 男女视频在线观看网站免费| 欧美成人a在线观看| 日本 欧美在线| 天天添夜夜摸| 88av欧美| 精品电影一区二区在线| 亚洲av二区三区四区| 国产精品一区二区三区四区免费观看 | 老熟妇仑乱视频hdxx| 亚洲色图av天堂| 久久久久免费精品人妻一区二区| 制服人妻中文乱码| 一级黄片播放器| 非洲黑人性xxxx精品又粗又长| 亚洲欧美激情综合另类| 亚洲精品久久国产高清桃花| 亚洲五月婷婷丁香| 日韩亚洲欧美综合| 性色av乱码一区二区三区2| 性色avwww在线观看| 国产成人av教育| 久久精品国产99精品国产亚洲性色| 天堂影院成人在线观看| a级一级毛片免费在线观看| 天堂av国产一区二区熟女人妻| 人妻丰满熟妇av一区二区三区| 两个人看的免费小视频| 成人永久免费在线观看视频| 亚洲电影在线观看av| 村上凉子中文字幕在线| 久久精品影院6| 国产精品久久视频播放| 亚洲男人的天堂狠狠| 亚洲av中文字字幕乱码综合| ponron亚洲| 亚洲精品色激情综合| 欧美成人a在线观看| 99热这里只有精品一区| 日韩国内少妇激情av| 最好的美女福利视频网| 亚洲精品美女久久久久99蜜臀| 亚洲欧美日韩高清专用| 精品99又大又爽又粗少妇毛片 | 丰满人妻一区二区三区视频av | 久久久色成人| 国产高清videossex| 国产伦一二天堂av在线观看| 国产高潮美女av| 男女做爰动态图高潮gif福利片| 少妇裸体淫交视频免费看高清| 桃红色精品国产亚洲av| 国产精品久久久久久久久免 | 熟女人妻精品中文字幕| 99热这里只有精品一区| 日日摸夜夜添夜夜添小说| 国产精品1区2区在线观看.| 嫩草影视91久久| 色吧在线观看| 国内毛片毛片毛片毛片毛片| 99在线视频只有这里精品首页| 99热这里只有精品一区| 9191精品国产免费久久| 亚洲专区国产一区二区| 国产男靠女视频免费网站| 久久香蕉国产精品| 男女视频在线观看网站免费| 国产成人欧美在线观看| 亚洲在线观看片| 搞女人的毛片| 在线观看av片永久免费下载| 色噜噜av男人的天堂激情| 久久香蕉国产精品| 国产精品久久久久久人妻精品电影| 亚洲成人久久爱视频| 国产成人a区在线观看| 欧美不卡视频在线免费观看| 最好的美女福利视频网| 久久久久久国产a免费观看| 最近在线观看免费完整版| 九色国产91popny在线| 一级毛片女人18水好多| 国产黄片美女视频| 最好的美女福利视频网| 国产欧美日韩精品亚洲av| 91在线精品国自产拍蜜月 | 日本熟妇午夜| 日韩欧美国产一区二区入口| 老汉色av国产亚洲站长工具| 国产精华一区二区三区| 国产精品久久电影中文字幕| 欧美+日韩+精品| 亚洲专区中文字幕在线| 亚洲18禁久久av| 免费观看的影片在线观看| 久久国产精品人妻蜜桃| 老熟妇乱子伦视频在线观看| 91av网一区二区| 一个人免费在线观看电影| 欧美成狂野欧美在线观看| 欧美国产日韩亚洲一区| 在线国产一区二区在线| 国产午夜福利久久久久久| 国产伦精品一区二区三区四那| 国产精品国产高清国产av| 91av网一区二区| 黄色视频,在线免费观看| 美女 人体艺术 gogo| 人人妻人人看人人澡| 国产三级中文精品| 欧美一区二区国产精品久久精品| 最近在线观看免费完整版| 看免费av毛片| 亚洲国产高清在线一区二区三| 九九久久精品国产亚洲av麻豆| 69人妻影院| 首页视频小说图片口味搜索| 欧美乱色亚洲激情| 国产野战对白在线观看| 亚洲精品国产精品久久久不卡| 一级作爱视频免费观看| 99热这里只有是精品50| 精品国产超薄肉色丝袜足j| av在线天堂中文字幕| 欧美乱妇无乱码| 亚洲精品乱码久久久v下载方式 | 免费人成视频x8x8入口观看| 欧美大码av| 午夜福利在线观看免费完整高清在 | 国产精品自产拍在线观看55亚洲| 少妇人妻精品综合一区二区 | h日本视频在线播放| 少妇人妻一区二区三区视频| 女人十人毛片免费观看3o分钟| 日韩欧美在线二视频| 成年女人看的毛片在线观看| 亚洲精品国产精品久久久不卡| 老司机午夜十八禁免费视频| 日韩欧美在线乱码| 国产一区二区在线观看日韩 | 熟女少妇亚洲综合色aaa.| 亚洲人成网站在线播放欧美日韩| 国产成年人精品一区二区| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 精品一区二区三区人妻视频| 国产精品免费一区二区三区在线| 日韩欧美精品免费久久 | 特大巨黑吊av在线直播| 黄色日韩在线| 十八禁网站免费在线| 久久性视频一级片| 真人做人爱边吃奶动态| 久久久久久人人人人人| 制服丝袜大香蕉在线| 丰满人妻熟妇乱又伦精品不卡| 少妇熟女aⅴ在线视频| 黄色女人牲交| 欧美乱妇无乱码| 欧美乱码精品一区二区三区| 一进一出抽搐gif免费好疼| 麻豆国产97在线/欧美| 精品午夜福利视频在线观看一区| 99视频精品全部免费 在线| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| 好看av亚洲va欧美ⅴa在| 91麻豆精品激情在线观看国产| 日韩亚洲欧美综合| 国产精品国产高清国产av| 免费电影在线观看免费观看| bbb黄色大片| 特级一级黄色大片| www.色视频.com| 黄色片一级片一级黄色片| 少妇熟女aⅴ在线视频| 国产91精品成人一区二区三区| 少妇人妻精品综合一区二区 | 亚洲,欧美精品.| 一级作爱视频免费观看| 首页视频小说图片口味搜索| 成人高潮视频无遮挡免费网站| 国产真实伦视频高清在线观看 | 午夜日韩欧美国产| 日本五十路高清| 性色avwww在线观看| 成人永久免费在线观看视频| 亚洲精品色激情综合| 色视频www国产| 小说图片视频综合网站| 99在线人妻在线中文字幕| 亚洲欧美日韩高清专用| 中文字幕久久专区| 亚洲成人久久性| 性色av乱码一区二区三区2| 欧美极品一区二区三区四区| 欧美日韩综合久久久久久 | 亚洲av熟女| 亚洲男人的天堂狠狠| 国语自产精品视频在线第100页| 亚洲精品国产精品久久久不卡| 国产精品永久免费网站| 人人妻人人澡欧美一区二区| 亚洲欧美日韩高清专用| 亚洲国产欧美网| av中文乱码字幕在线| 日日干狠狠操夜夜爽| 高清在线国产一区| 精品人妻偷拍中文字幕| 国产伦人伦偷精品视频| 法律面前人人平等表现在哪些方面| 精品福利观看| 国产三级黄色录像| 级片在线观看| 亚洲美女黄片视频| 啪啪无遮挡十八禁网站| 脱女人内裤的视频| 香蕉av资源在线| 日本黄色片子视频| 精品一区二区三区人妻视频| 欧美日韩乱码在线| 怎么达到女性高潮| 国产av在哪里看| 9191精品国产免费久久| 欧美区成人在线视频| 久久精品综合一区二区三区| 日本撒尿小便嘘嘘汇集6| 欧美绝顶高潮抽搐喷水| 黄片小视频在线播放| 国产中年淑女户外野战色| 中文资源天堂在线| 99久久精品热视频| 男人和女人高潮做爰伦理| 变态另类成人亚洲欧美熟女| 一区二区三区免费毛片| 高清毛片免费观看视频网站| 国产欧美日韩精品一区二区| avwww免费| 欧美3d第一页| 久久久久亚洲av毛片大全| 日韩欧美精品免费久久 | 蜜桃亚洲精品一区二区三区| 婷婷精品国产亚洲av| 美女大奶头视频| 亚洲欧美一区二区三区黑人| 国产精品一区二区免费欧美| 女同久久另类99精品国产91| 精品国产美女av久久久久小说| 99久久精品国产亚洲精品| 日韩中文字幕欧美一区二区| 亚洲精品影视一区二区三区av| 国产精品国产高清国产av| 狠狠狠狠99中文字幕| 中文字幕av在线有码专区| 久久精品综合一区二区三区| 国产av一区在线观看免费| 国产真人三级小视频在线观看| 男插女下体视频免费在线播放| 国产精品1区2区在线观看.| 欧美一区二区亚洲| 久久久精品大字幕| 在线视频色国产色| 欧美色视频一区免费| 久久久久国产精品人妻aⅴ院| 国产成人福利小说| 欧美国产日韩亚洲一区| 人人妻人人澡欧美一区二区| a在线观看视频网站| 无限看片的www在线观看| 亚洲国产欧洲综合997久久,| www国产在线视频色| 欧美日韩黄片免| 日本a在线网址| 欧美高清成人免费视频www| 亚洲,欧美精品.| www.熟女人妻精品国产| 久久精品国产自在天天线| 首页视频小说图片口味搜索| 中文字幕av成人在线电影| 一本一本综合久久| 老司机在亚洲福利影院| 免费大片18禁| 天堂网av新在线| 国产高清激情床上av| 亚洲成a人片在线一区二区| 精品一区二区三区人妻视频| 夜夜爽天天搞| 国产伦一二天堂av在线观看| 变态另类丝袜制服| 最近最新中文字幕大全电影3| 国产成人系列免费观看| 国产中年淑女户外野战色| 在线天堂最新版资源| 欧美成人性av电影在线观看| 国产欧美日韩精品一区二区| 中文字幕人妻丝袜一区二区| 成熟少妇高潮喷水视频| 午夜福利18| 亚洲成av人片免费观看| 香蕉久久夜色| 欧美日韩中文字幕国产精品一区二区三区| 亚洲国产精品久久男人天堂| 午夜福利18| 精品无人区乱码1区二区| 国产高清视频在线播放一区| 国产亚洲av嫩草精品影院| 黄色日韩在线| 69av精品久久久久久| 又黄又粗又硬又大视频| 90打野战视频偷拍视频| 精品国内亚洲2022精品成人| 岛国视频午夜一区免费看| 日韩欧美国产在线观看| 色综合亚洲欧美另类图片| 亚洲不卡免费看| 日韩免费av在线播放| 欧美色视频一区免费| 亚洲国产高清在线一区二区三| 亚洲五月天丁香| 成年版毛片免费区| 天天添夜夜摸| 国产欧美日韩精品一区二区| 天堂√8在线中文| 午夜精品一区二区三区免费看| 国产欧美日韩一区二区三| 成熟少妇高潮喷水视频| 亚洲欧美日韩无卡精品| 美女高潮的动态| 搡老岳熟女国产| 性欧美人与动物交配| 亚洲精品成人久久久久久| 最新在线观看一区二区三区| 嫩草影视91久久| 国产日本99.免费观看| 欧美日韩瑟瑟在线播放| 狠狠狠狠99中文字幕| 午夜免费男女啪啪视频观看 | 久久久久久久精品吃奶| 欧美日本视频| 亚洲无线在线观看| 免费人成在线观看视频色| 免费高清视频大片| 免费在线观看日本一区| 久久久成人免费电影| 老司机福利观看| 五月伊人婷婷丁香| 黄色片一级片一级黄色片| 日韩人妻高清精品专区| 操出白浆在线播放| 一进一出抽搐gif免费好疼| 国产淫片久久久久久久久 | 欧美大码av| 色综合亚洲欧美另类图片| 欧洲精品卡2卡3卡4卡5卡区| 最新中文字幕久久久久| 老司机福利观看| www日本黄色视频网| 精品日产1卡2卡| www日本黄色视频网| 欧美成人免费av一区二区三区| 亚洲精品456在线播放app | 美女高潮的动态| 无遮挡黄片免费观看| 国产单亲对白刺激| 色在线成人网| 久久精品国产99精品国产亚洲性色| 午夜福利18| 国产精品久久久久久久久免 | 亚洲va日本ⅴa欧美va伊人久久| 免费观看人在逋| 国产亚洲精品一区二区www| 在线观看美女被高潮喷水网站 | 欧洲精品卡2卡3卡4卡5卡区| 欧美午夜高清在线| 观看免费一级毛片| 日本精品一区二区三区蜜桃| 级片在线观看| 丰满人妻一区二区三区视频av | 亚洲一区二区三区不卡视频| 亚洲狠狠婷婷综合久久图片| 免费av观看视频| 欧美一区二区亚洲| 久久久久免费精品人妻一区二区| 综合色av麻豆| 国产亚洲精品综合一区在线观看| 观看美女的网站| 免费av不卡在线播放| 久久久久九九精品影院| 91久久精品电影网| 成人鲁丝片一二三区免费| 午夜福利欧美成人| 天堂影院成人在线观看| 亚洲av成人av| 一本综合久久免费| 免费搜索国产男女视频| 亚洲国产高清在线一区二区三| 十八禁人妻一区二区| h日本视频在线播放| 桃红色精品国产亚洲av| 精品人妻偷拍中文字幕| 欧美中文日本在线观看视频| 精品国产三级普通话版| 一本精品99久久精品77| 国产精品三级大全| 天天一区二区日本电影三级| 亚洲aⅴ乱码一区二区在线播放| 少妇高潮的动态图|