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

    基于優(yōu)化SVM和BP神經(jīng)網(wǎng)絡(luò)檢測雙孢蘑菇早期病害

    2022-11-29 03:09:52陳子涵溫志強(qiáng)
    關(guān)鍵詞:雙孢染病蘑菇

    黃 亮, 魏 萱,2,3, 陳子涵, 溫志強(qiáng)

    (1.福建農(nóng)林大學(xué)機(jī)電工程學(xué)院,福建 福州 350002;2.福建農(nóng)林大學(xué)食品學(xué)院,福建 福州 350002;3.福建省農(nóng)業(yè)信息感知技術(shù)重點(diǎn)實(shí)驗(yàn)室,福建 福州 350002;4.福建水利電力職業(yè)技術(shù)學(xué)院,福建 永安 366000;5.福建農(nóng)林大學(xué)生命科學(xué)學(xué)院,福建 福州 350002)

    雙孢蘑菇(Agaricusbisporus)又名白蘑菇、洋蘑菇,是傘菌目(Agaricales)蘑菇屬(Agaricus)食用菌.雙孢蘑菇含有豐富的氨基酸、核苷酸、超氧化物歧化酶、胰蛋白酶等,蛋白質(zhì)含量高達(dá)42%,具有較高的食用價(jià)值和醫(yī)藥價(jià)值[1-2].有害疣孢霉菌(MycogoneperniciosaMagn.)屬于子囊菌門(Ascomytota)盤菌亞門(Pezizomycotina),在20~28 ℃時(shí)生長較快[3].雙孢蘑菇被疣孢霉菌侵染引起的疣孢霉病,又稱褐腐病、白腐病,輕度侵染后其表面會(huì)形成褐色斑點(diǎn)和白色菌絲,重度侵染可能會(huì)形成畸形菇,且不同生長發(fā)育階段感染癥狀的差異明顯[4].

    多酚氧化酶(polyphenol oxidase,PPO)又稱兒茶酚氧化酶,是一種存在于植物、食用菌中的金屬蛋白酶.它能夠引起果蔬酶促褐變,從而影響果蔬的外觀品質(zhì)和口感風(fēng)味[5].PPO也能作為一種參與黑色素細(xì)胞形成的氧化還原酶,它能抑制黑色素的沉積,從而避免色斑、褐斑等疾病的產(chǎn)生[6].相反,黑色素可以治療白化病等疾病,故防御活性酶可能與食用菌病害存在一定聯(lián)系.已有許多專家研究了果蔬防御活性酶與真菌病害的相關(guān)性,指出防御酶對(duì)大麥抗病起到一定作用[7].傳統(tǒng)的疣孢霉菌檢測方法有柯赫氏法則致病性檢驗(yàn)、核糖體轉(zhuǎn)錄間隔區(qū)(rDNA-ITS)序列分析[8],這些分子生物學(xué)技術(shù)方法往往會(huì)破壞樣本的外部品質(zhì)和內(nèi)部結(jié)構(gòu).因此,采用無損、快速的高光譜成像技術(shù)檢測樣本的品質(zhì)尤為重要.

    近年來光譜成像技術(shù)已經(jīng)大量運(yùn)用到水果[9]、蔬菜[10]、食用菌[11]等檢測中,具有檢測范圍廣、精度高等優(yōu)點(diǎn).馬淏等[12]基于光譜技術(shù)對(duì)雙孢蘑菇的新鮮度進(jìn)行分級(jí)檢測,結(jié)果表明,結(jié)合粒子群算法(particle swarm optimization,PSO)和海姆優(yōu)化算法(seagull optimization algorithm,SOA)極限學(xué)習(xí)機(jī)(extreme learning machine,ELM)檢測模型精度較高,預(yù)測精度分別為92.5%、94%.Yang et al[13]采用反向傳播神經(jīng)網(wǎng)絡(luò)(back propagation, BP),構(gòu)建了不同生育期杏鮑菇MSC-CARS-BP菌絲體檢測模型,不同生育期菌絲體檢測結(jié)果準(zhǔn)確率均為99.67%.Xiao et al[14]基于高光譜技術(shù)預(yù)測滲透脫水中雙孢蘑菇的可溶性固形物(Soluble Solids Content,SSC)含量,結(jié)果表明經(jīng)過正交信號(hào)校正(orthogonal signal correction,OSC)的支持向量機(jī)(support vector machine,SVM)模型預(yù)測集R2p為0.883,且OSC-CARS-SVM預(yù)測模型的精度最佳.Mollazade[15]基于高光譜成像技術(shù)對(duì)4種褐變雙孢蘑菇進(jìn)行了分類,其中經(jīng)過競爭自適應(yīng)加權(quán)算法(competitive adaptive reweighted sampling,CARS)后的偏最小二乘回歸(partial least squares regression,PLSR)模型的校正集和預(yù)測集的分類準(zhǔn)確率分別為80.6%、80.3%,并將全光譜點(diǎn)的主成分分析(principal components analysis,PCA)圖與常規(guī)紅、綠、藍(lán)(red,green,blue, RGB)成像波段分類圖進(jìn)行了比較.

    早期雙孢蘑菇(原基期至小菇期)菌蓋上的菌絲細(xì)小,很難提取到疣孢霉菌進(jìn)行病害判別.又由于疣孢霉菌入侵雙孢蘑菇后,會(huì)合成大量PPO來抵制病菌的侵害[16].故PPO與疣孢霉病可能存在某種相關(guān)性.目前利用高光譜成像技術(shù)預(yù)測雙孢蘑菇早期PPO的研究較少,對(duì)雙孢蘑菇疣孢霉病判別鮮有研究.因此,本研究先提取與PPO相關(guān)的特征光譜,然后進(jìn)一步篩選光譜特征判別雙孢蘑菇的早期病害,旨在為食用菌早期病害判別提供一種新方法,為進(jìn)一步開發(fā)早期雙孢蘑菇病害的快速無損鑒別設(shè)備提供參考.

    1 材料與方法

    1.1 樣本的培養(yǎng)與接菌

    試驗(yàn)材料為福建省農(nóng)業(yè)科學(xué)院食用菌研究所提供的W192型雙孢蘑菇菌菌種、My.p0012型疣孢霉菌菌種.以小麥播種的方式將菌種播撒在栽培料中混合均勻,采用無菌食用菌栽培袋以400 g·袋-1的標(biāo)準(zhǔn)分裝;將雙孢蘑菇放入22 ℃、90%相對(duì)濕度、無光照條件的人工氣候培養(yǎng)箱中培養(yǎng);當(dāng)菌絲大約生長到栽培料的2/3時(shí)覆土,覆土高度約3 cm.覆土后適當(dāng)調(diào)高環(huán)境溫度并進(jìn)行輕噴,以利于菌絲爬土.待覆土后10 d左右,溫度調(diào)低至20 ℃,觀察雙孢蘑菇出菇情況.培養(yǎng)期間每天通風(fēng)2次,每次1 h[17].

    采用孢子懸浮液接種雙孢蘑菇.先將疣孢霉菌放入馬鈴薯葡萄糖瓊脂(potato dextrose agar,PDA)中培養(yǎng)7 d,再用無菌水沖洗PDA活化疣孢霉菌孢子,并用血球計(jì)數(shù)板檢測疣孢霉菌孢子懸液濃度,其值為(1.0±0.4)×105個(gè)·mL-1.然后在2/3覆土層的表面均勻噴灑5 mL疣孢霉菌孢子懸液接菌[18].

    1.2 防御酶活性值的測定

    選用蘇州科銘生物技術(shù)有限公司生產(chǎn)的試劑盒測定指標(biāo)值,測定方法參考文獻(xiàn)[19].PPO的特征波段光密度用酶標(biāo)儀測定,計(jì)算PPO活性值.酶標(biāo)儀(SpectraMax?Plus 384)配備高精度光柵型單色器和溫度控制器,波長可調(diào)節(jié)至190~1 000 nm,步進(jìn)1 nm,帶寬2 nm.

    1.3 高光譜圖像的獲取

    高光譜成像系統(tǒng)主要包括ImSpector V10E光譜儀、OLE 23型成像鏡頭、2個(gè)150 W光纖鹵素?zé)?、電控升降臺(tái)和計(jì)算機(jī)等[20].采用HSI Analyzer軟件進(jìn)行白板校正、圖像存儲(chǔ)等.設(shè)置合理的曝光時(shí)間等參數(shù)以避免信息飽和以及圖像失真,調(diào)整后物鏡高度為350 mm,CCD相機(jī)曝光時(shí)間為0.13 s,相機(jī)推掃移動(dòng)速度為2.20 m·s-1.為了消除傳感器暗電流、光線不均勻等產(chǎn)生的影響,對(duì)高光譜圖像進(jìn)行黑白校正:

    (1)

    式中:I0表示校正后的高光譜圖像;A表示白板圖像;B表示黑板圖像;IC表示原始高光譜圖像.

    采集原基期、菇蕾期、幼菇期和小菇期的染病雙孢蘑菇高光譜圖像以實(shí)現(xiàn)早期檢測.接種后0~7 d為原基期,8~9 d為菇蕾期,10~11 d為幼菇期,11~12 d為小菇期,各個(gè)生長周期為2 d左右.原基期普遍呈白色米粒狀;菇蕾期則為黃豆大小,且菌柄、菌蓋已成雛形;幼菇期、小菇期的菌蓋分別呈蓋半球狀、扁蓋半球狀.小菇期后的病害特征較為明顯,菌蓋表面會(huì)出現(xiàn)褐變等患病癥狀,故可將接種后7~11 d的樣本認(rèn)定為早期.每個(gè)周期健康和染病雙孢蘑菇的數(shù)量各50個(gè).以每個(gè)圖像菌蓋周圍部分區(qū)域?yàn)楦信d趣區(qū)域(region of interest,ROI),采用新區(qū)域疊加所有樣本的數(shù)據(jù)提取區(qū)域,提取形狀有橢圓形、矩形、八邊形等.最終得到4×2×50個(gè)高光譜圖像的ROI光譜數(shù)據(jù).

    1.4 特征波長的優(yōu)選方法

    CARS采用自適應(yīng)重加權(quán)采樣(adaptive reweighted sampling,ARS)和指數(shù)衰減函數(shù)(exponential decay function,EDF)來建立多個(gè)PLSR模型,從而篩選出極值回歸系數(shù),然后利用十字交叉算法優(yōu)選出均方誤差較小的PLS模型[21].連續(xù)投影算法(successive projections algorithm,SPA)通過引入最小化變量解決矢量空間的共線性來選擇最優(yōu)變量,從而提高建模效率.SPA先對(duì)所有光譜波段進(jìn)行投影分析,然后將相對(duì)較大的投影向量值作為特征波長,從全波段范圍內(nèi)篩選出無關(guān)信息量最少的一組光譜信息來建立多元回歸模型.所篩選的變量不僅共線性較小,而且避免了光譜波段的共線性和信息重疊,提高了建模精度[22].

    1.5 預(yù)測及判別模型

    SVM是一種典型的二分類有監(jiān)督學(xué)習(xí)算法.SVM通過采用序列最小優(yōu)化法(sequential minimal optimization, SMO)來不斷迭代出2個(gè)最合適的拉格朗日算子α,從而確認(rèn)最終模型.在本研究中采用Sigmoid核函數(shù)來線性劃分?jǐn)?shù)據(jù)集.

    ELM與單隱含層前饋神經(jīng)網(wǎng)絡(luò)(single hidden layer feed-forward neural network,SLFN)相似.與反向傳播的神經(jīng)網(wǎng)絡(luò)不同,ELM隨機(jī)選取輸入層權(quán)重和隱含層偏置,輸出層權(quán)重由最小化訓(xùn)練誤差項(xiàng)和輸出層權(quán)重范數(shù)的正則項(xiàng)構(gòu)成,具有學(xué)習(xí)速度快、訓(xùn)練參數(shù)少、泛化能力強(qiáng)等特點(diǎn)[23].

    BP神經(jīng)網(wǎng)絡(luò)具有正、反向傳播的特點(diǎn),是目前應(yīng)用最廣的多層前饋神經(jīng)網(wǎng)絡(luò)之一.其結(jié)構(gòu)主要由輸入層、隱含層和輸出層構(gòu)成,每層網(wǎng)絡(luò)的神經(jīng)元節(jié)點(diǎn)數(shù)由實(shí)際問題確定.BP神經(jīng)網(wǎng)絡(luò)迭代訓(xùn)練時(shí)首先進(jìn)行正向傳播,輸入層的每個(gè)神經(jīng)元通過權(quán)值連接隱含層的每個(gè)節(jié)點(diǎn),將隱含層每個(gè)節(jié)點(diǎn)的輸入值帶入激活函數(shù)(sigmoid)中可得到每個(gè)節(jié)點(diǎn)的輸出值,sigmoid函數(shù)的導(dǎo)數(shù)能很好地逼近非線性關(guān)系,并將值控制在0~1.每個(gè)隱含層節(jié)點(diǎn)輸出值又通過權(quán)值連接到每個(gè)輸出層的神經(jīng)元.然后計(jì)算輸出值與實(shí)際值的誤差,將誤差帶入反向傳播中,不斷更新網(wǎng)絡(luò)層之間的連接權(quán)值和隱含層、輸出層的閾值,直到輸出誤差達(dá)到最小[24].本研究的關(guān)鍵算法步驟表示如下:

    Oy=sigmoid(x×u+hy)

    (2)

    Oz=sigmoid(y×v+hz)

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    式中,Oy、OZ分別為隱含層、輸出層節(jié)點(diǎn)的輸出值,x、y分別為輸入層、隱含層神經(jīng)元節(jié)點(diǎn)值,u、v分別為輸入層與隱含層之間的權(quán)值、隱含層與輸出層之間的權(quán)值,hy、hz分別為隱含層、輸出層的閾值.Ix、Iy、Iz分別為輸入層、隱含層、輸出層的輸入節(jié)點(diǎn)值.δ為輸出值與實(shí)際值的誤差.σ′為sigmoid函數(shù)的導(dǎo)數(shù).

    式(4)、(6)、(8)中的LRA、LRB、LRC為學(xué)習(xí)率.式(2)、(3)分別是正向傳播中隱含層、輸出層的網(wǎng)絡(luò)節(jié)點(diǎn)輸出公式,式(5)為損失函數(shù)的計(jì)算公式,式(4)、(6)、(7)、(8)分別是反向傳播中權(quán)值、閾值的更新公式.

    1.6 粒子群優(yōu)化算法

    PSO的基本思想是共享群體中個(gè)體的信息來尋求適應(yīng)度函數(shù)的最優(yōu)解.群體中每個(gè)粒子有且僅有速度和位置2個(gè)屬性[25].若計(jì)算適應(yīng)度函數(shù)最小值,將個(gè)體極小值賦予全局最優(yōu)解;反之,將個(gè)體極大值賦予全局最優(yōu)解.然后將當(dāng)前全局最優(yōu)位置和個(gè)體歷史最優(yōu)位置帶入式(9)~(10)中迭代速度、位置,最后根據(jù)目標(biāo)函數(shù)問題性質(zhì)更新全局最優(yōu)位置解.

    (9)

    (10)

    式中,ω為慣性因子,取值范圍為0.8~0.2.其值越大,全局尋優(yōu)能力越強(qiáng),局部尋優(yōu)能力越弱.Vj為第j個(gè)粒子的速度,Xj為第j個(gè)粒子的位置.C1和C2為學(xué)習(xí)率.Pbest為個(gè)體極值,Gbest為全局最優(yōu)解.將個(gè)體的位置帶入目標(biāo)函數(shù)求得個(gè)體極值,然后根據(jù)目標(biāo)函數(shù)最優(yōu)解性質(zhì)將個(gè)體極值(Pbest)賦予全局最優(yōu)解(Gbest).

    采用PSO優(yōu)化支持向量機(jī)算法的步驟為:(1)歸一化數(shù)據(jù)集,初始化粒子群的速度、位置;(2)將核寬度(g)和懲罰參數(shù)(c)分別作為粒子的位置(X1,X2),并帶入式(9)、(10)迭代更新速度、位置,計(jì)算粒子的個(gè)體最優(yōu)值和群體極值;(3)若滿足終止條件,將獲得的最優(yōu)參數(shù)c、g帶入SVM模型中迭代訓(xùn)練;(4)若不滿足終止條件,重復(fù)更新粒子的速度、位置,尋找當(dāng)前個(gè)體最優(yōu)極值,直至滿足條件,輸出c、g.

    表1 400個(gè)樣本的PPO活性真實(shí)值統(tǒng)計(jì)結(jié)果Table 1 Statistics of PPO values of 400 A. bisporus samples

    圖1 健康、染病雙孢蘑菇高光譜平均光譜曲線Fig.1 Average hyperspectral curve of healthy and infected A. bisporus

    2 結(jié)果與分析

    2.1 光譜數(shù)據(jù)與樣本

    采用柯爾莫哥洛夫-斯摩洛夫(kolmogorov-smirnov,KS)算法[26]將400個(gè)樣本PPO活性值按2∶1的比例分為訓(xùn)練集和測試集,統(tǒng)計(jì)結(jié)果如表1所示.訓(xùn)練集、測試集的PPO值范圍分別為109.560~474.360 U·g-1、115.247~433.861 U·g-1.測試集指標(biāo)范圍在訓(xùn)練集內(nèi),數(shù)據(jù)劃分合理.健康、染病雙孢蘑菇所有時(shí)期的平均光譜曲線如圖1所示.健康、染病雙孢蘑菇的平均光譜曲線變化相似.相對(duì)于表面有褐斑、霉菌的染病雙孢蘑菇,健康雙孢蘑菇表面光滑、白凈,總體光譜反射率較高.去除首位噪聲,選取450~1 000 nm高光譜波長范圍(307個(gè)波段)進(jìn)行研究.

    2.2 基于PPO特征波段的活性值預(yù)測

    2.2.1 競爭自適應(yīng)加權(quán)算法優(yōu)選波段 采用CARS法對(duì)所有樣本全波段(307個(gè)波段)進(jìn)行PPO活性值特征波段的選取.如圖2a所示,設(shè)定采樣次數(shù)N=50,隨著采樣次數(shù)的增加,曲線斜率不斷減小,表明從“粗選”到“細(xì)選”過程中高光譜波段數(shù)由原來的307個(gè)減少至24個(gè).如圖2b所示,當(dāng)采樣次數(shù)為26次時(shí)均方誤差最小,其值為39.采樣次數(shù)從26次增加到50次的過程中,均方誤差不斷地增大或減小,表明與PPO敏感的特征波段被剔除或保留.為了建立穩(wěn)定的預(yù)測模型,選取均方誤差最小時(shí)采樣次數(shù)所對(duì)應(yīng)的特征波段數(shù).即當(dāng)采樣數(shù)為26次時(shí),提取到了與雙孢蘑菇PPO活性值敏感的24個(gè)特征波段.特征波段依次為540.8、542.6、549.5、565.2、566.9、570.4、579.1、616.0、617.7、619.5、621.3、644.2、656.7、658.4、764.5、801.0、802.9、810.2、823.1、908.4、925.3、942.2、949.7、980.0 nm.

    a:光譜波段數(shù)變化過程;b:均方誤差變化過程.圖2 CARS算法優(yōu)選波段過程圖Fig.2 Optimization of spectral bands by CARS algorithm

    2.2.2 連續(xù)投影算法優(yōu)選波段 采用SPA對(duì)所有樣本全波段(307個(gè)波段)進(jìn)行PPO活性值特征波段的選取,選取過程如圖3所示.圖3a表示隨著模型內(nèi)變量的引入,均方誤差的變化過程.當(dāng)引入10個(gè)變量時(shí),均方誤差最小.圖3b表示各個(gè)波段的反射率.最終提取到了10個(gè)特征波段,依次為499.4、566.9、653.1、726.5、823.1、910.3、929.0、962.9、995.2、1 000.9 nm.

    a:均方誤差變化過程;b:波段篩選過程.圖3 SPA算法優(yōu)選波段過程圖Fig.3 Optimization of spectral bands by SPA algorithm

    分別將經(jīng)過SPA、CARS提取的特征波段帶入PLSR模型中預(yù)測PPO活性值.SPA-PLSR模型的訓(xùn)練集的相關(guān)系數(shù)(Rc)和均方誤差(RMSEC)分別為0.921、0.303,測試集的Rp、RMSEP分別為0.917、0.349.CARS-PLSR模型訓(xùn)練集的Rc、RMSEC分別為0.923、0.306,測試集的Rp、RMSEP分別為0.904、0.340.

    2.3 光譜反射率與疣孢霉病的相關(guān)性

    本研究所提取的高光譜信息位于發(fā)病區(qū)域,可采用所提取特征信息對(duì)雙孢蘑菇的染病情況進(jìn)行判別.為了在對(duì)PPO敏感的波段中選取對(duì)病害敏感的特征波段,擬采用逐步回歸法和相關(guān)性分析在經(jīng)過CARS所提取的24個(gè)特征波段中篩選出與染病相關(guān)的光譜波段.

    首先篩選出對(duì)染病情況貢獻(xiàn)最大的特征波段建立回歸模型,然后逐步加入剩余特征波段,最后在加入波段的過程中及時(shí)剔除不顯著和共線性的特征波段,共建立8個(gè)回歸模型.對(duì)染病情況貢獻(xiàn)最大的566.9 nm波段最先被引入第1個(gè)回歸模型,其次566.9、656.7 nm波段被引入第2個(gè)回歸模型.以此類推,最后有8個(gè)特征波段按貢獻(xiàn)率依次被引入第8個(gè)回歸模型中,引入波段的順序?yàn)?66.9、656.7、617.7、658.4、801.0、980.0、908.4和942.2 nm波段.引入標(biāo)準(zhǔn)為F檢驗(yàn)概率小于0.05.各回歸模型的統(tǒng)計(jì)參數(shù)如表2所示.

    隨著蘑菇染病日的延長,貢獻(xiàn)較低的特征波段樣本內(nèi)水分的流失會(huì)導(dǎo)致其吸收峰發(fā)生偏移,616.0、619.5、621.3 nm波段可能是水分特征波段617 nm發(fā)生偏移所致.此外雙孢蘑菇在培養(yǎng)、接菌后,其內(nèi)部物質(zhì)含量會(huì)發(fā)生改變,802.9、810.2、823.1 nm可能是維C特征波段801、810、821 nm發(fā)生位置偏移所致.540.8、542.6和549.5 nm附近的吸收峰可能與雙孢蘑菇內(nèi)PPO分子中的O-H二級(jí)倍頻(540 nm)、C-H二級(jí)倍頻(543 nm)有關(guān),565.2 nm附近的吸收峰可能與其C-H四級(jí)倍頻(566 nm)有關(guān)[27].從表3可知,隨著特征波段的引入,各回歸模型的復(fù)測定系數(shù)(R)逐漸變大.回歸模型8的R為0.891,表明所篩選的8個(gè)特征波段與染病判別的擬合度最大,R上標(biāo)的a~h分別表示各個(gè)模型被引入特征波段的光密度,a表示第1個(gè)引入的波段光密度,h表示被引入的8個(gè)波段光密度.其回歸模型8的決定系數(shù)(R2)為0.793,表示8個(gè)波段的光密度能解釋雙孢蘑菇染病情況的79%,其染病情況的21%由其他變量來解釋.

    表2 染病判別回歸模型的統(tǒng)計(jì)參數(shù)Table 2 Statistical parameters for regression analysis of infection recognition model

    表3 特征波段與染病情況的顯著性分析結(jié)果1)Table 3 Significance analysis results of characteristic bands and infection

    將經(jīng)過逐步回歸篩選出的特征波段與經(jīng)過SPA篩選出的特征波段結(jié)合,采用皮爾曼相關(guān)系數(shù)對(duì)疣孢霉病進(jìn)行顯著性分析(表3).從表3可知,929.0 nm波段與染病情況在0.05水平上有顯著相關(guān)性,其他光譜波段都在0.01水平上具有極顯著相關(guān)性,且449.4 nm波段與染病情況呈負(fù)相關(guān).

    2.4 PSO算法優(yōu)化及BP判別模型的建立

    將篩選出的12個(gè)與染病相關(guān)的高光譜波段分別帶入ELM、SVM、PSO-SVM和BP神經(jīng)網(wǎng)絡(luò)中進(jìn)行建模比較分析.

    從圖4可知,設(shè)置迭代次數(shù)為100,種群規(guī)模為30,采用8折交叉驗(yàn)證來優(yōu)化參數(shù).為了平衡尋優(yōu)能力,慣性權(quán)重w設(shè)為1.以核寬度和懲罰參數(shù)為目標(biāo)函數(shù)的輸入變量進(jìn)行迭代,優(yōu)化后的懲罰參數(shù)、核寬度分別為78.149、0.01,然后將最優(yōu)解帶入SVM計(jì)算雙孢蘑菇的正確判別率.

    圖4 PSO優(yōu)化后的適應(yīng)度函數(shù)值變化圖Fig.4 Changes on fitting values after PSO optimization

    在圖5的3組學(xué)習(xí)率下的誤差走勢圖中,迭代次數(shù)設(shè)為500,反向傳播更新權(quán)值、閾值后,均方誤差不斷降低.學(xué)習(xí)率A=0.1,學(xué)習(xí)率B=0.1時(shí),學(xué)習(xí)率較小,曲線很平滑,收斂速度較慢.A=0.5,B=0.6時(shí),學(xué)習(xí)率適當(dāng)增加后,收斂速度較第1組有所降低,曲線較平滑.A=0.8,B=0.9時(shí),學(xué)習(xí)率快到上限時(shí),收斂速度未發(fā)生明顯改變,但曲線發(fā)散,迭代過程不穩(wěn)定.輸入層節(jié)點(diǎn)數(shù)為光譜波段數(shù)(12個(gè)),輸出層節(jié)點(diǎn)數(shù)為輸出類型數(shù)(2個(gè)).隱含層節(jié)點(diǎn)數(shù)可由下列3個(gè)經(jīng)驗(yàn)公式計(jì)算:

    (11)

    m=log2n

    (12)

    (13)

    式中,m、n、l分別是隱含層、輸入層、輸出層的節(jié)點(diǎn)數(shù),α為1~10的常數(shù).

    從式(11)可得到隱含層節(jié)點(diǎn)數(shù)范圍為4.7~13.7,從式(12)可得到隱含層節(jié)點(diǎn)數(shù)為3.6,從式(13)可得到隱含層節(jié)點(diǎn)數(shù)為4.9.隱含層節(jié)點(diǎn)數(shù)分別取式(12)、(13)的值以及式(11)的最大值.迭代后均方誤差的變化趨勢如圖6所示,當(dāng)m=4時(shí),在迭代次數(shù)200次以內(nèi)誤差曲線較為發(fā)散,大于200次迭代次數(shù)誤差趨向平穩(wěn).當(dāng)m=5時(shí),迭代次數(shù)在250次以內(nèi)誤差曲線較平滑,迭代次數(shù)大于250次曲線較為發(fā)散;且與m=4相比較,收斂后的均方誤差較大.當(dāng)m=14時(shí),迭代過程中誤差曲線較為平滑,且與m=4,m=5相比較,收斂后的均方誤差有所減小.當(dāng)m=40時(shí),迭代2次達(dá)到過擬合的狀態(tài),訓(xùn)練模型停止.

    圖5 不同學(xué)習(xí)率下的BP網(wǎng)絡(luò)誤差圖Fig.5 Error diagram of BP network under different learning rates

    表4為6種模型的建模效果,PSO-SVM、BP1、BP2測試集健康識(shí)別率最高,其值為93.333%,BP2測試集染病識(shí)別率、總體識(shí)別率最高,其值分別為95.890%、94.737%.

    表4 6種模型的建模效果Table 4 Modeling effects of 6 models

    3 小結(jié)

    本文采用粒子群算法優(yōu)化支持向量機(jī)和BP神經(jīng)網(wǎng)絡(luò)檢測雙孢蘑菇的早期品質(zhì).所建立的PPO預(yù)測模型SPA-PLSR精度最高,其RP、RMSEP分別為0.917、0.348.在所提取的特征波段中有8個(gè)被引入到回歸模型中,決定系數(shù)可達(dá)78.6%.所建立的BP模型對(duì)雙孢蘑菇病害的判別率最高,總體判別率可達(dá)94.74%.結(jié)果表明,優(yōu)化后的機(jī)器學(xué)習(xí)模型可較好地實(shí)現(xiàn)雙孢蘑菇外部品質(zhì)的判別.

    猜你喜歡
    雙孢染病蘑菇
    偶感
    均勻網(wǎng)絡(luò)上SIR模型三種不同逼近方法比較
    雙孢蘑菇菌種質(zhì)量檢驗(yàn)規(guī)程探討
    雙孢蘑菇低溫保鮮試驗(yàn)初報(bào)
    食用菌(2017年3期)2017-05-24 06:52:18
    愛 情
    詩選刊(2016年9期)2016-11-26 13:47:43
    爆圖團(tuán)
    婦女之友(2016年9期)2016-11-07 19:39:58
    西藏野生雙孢蘑菇栽培種營養(yǎng)成分研究
    西藏科技(2016年8期)2016-09-26 09:00:59
    蘑菇
    蘑菇傘
    食餌染病生態(tài)流行病系統(tǒng)的穩(wěn)定性與最優(yōu)收獲
    五月天丁香电影| 亚洲欧美中文字幕日韩二区| 亚洲精品,欧美精品| 午夜日本视频在线| 人人妻人人澡人人看| 亚洲伊人久久精品综合| 色哟哟·www| 亚洲情色 制服丝袜| 免费久久久久久久精品成人欧美视频 | 亚洲精品视频女| 久久久久国产网址| 国产精品熟女久久久久浪| 午夜老司机福利剧场| 男男h啪啪无遮挡| 亚洲av免费高清在线观看| a级毛色黄片| 亚洲高清免费不卡视频| 亚洲综合色惰| 免费不卡的大黄色大毛片视频在线观看| 制服诱惑二区| 高清午夜精品一区二区三区| av专区在线播放| 亚洲欧美日韩另类电影网站| videosex国产| 久久精品国产a三级三级三级| 一区二区三区免费毛片| 亚洲av成人精品一二三区| 在线天堂最新版资源| 大又大粗又爽又黄少妇毛片口| 国内精品宾馆在线| 一边亲一边摸免费视频| 日韩人妻高清精品专区| 欧美人与善性xxx| 在线亚洲精品国产二区图片欧美 | 男人爽女人下面视频在线观看| 国产男人的电影天堂91| 国产精品熟女久久久久浪| 久久人人爽人人爽人人片va| 男女边摸边吃奶| 人妻一区二区av| 亚洲成人一二三区av| 王馨瑶露胸无遮挡在线观看| 午夜激情久久久久久久| 男女边摸边吃奶| 成年人免费黄色播放视频| 乱人伦中国视频| 在现免费观看毛片| 国产片内射在线| 日本午夜av视频| 国产av码专区亚洲av| 久热这里只有精品99| 美女视频免费永久观看网站| 嫩草影院入口| av又黄又爽大尺度在线免费看| 美女大奶头黄色视频| 伦精品一区二区三区| 精品99又大又爽又粗少妇毛片| 亚洲婷婷狠狠爱综合网| 又粗又硬又长又爽又黄的视频| 久久热精品热| 精品人妻熟女毛片av久久网站| 99久久人妻综合| 国产成人aa在线观看| 久久 成人 亚洲| 亚洲精品乱久久久久久| 免费观看av网站的网址| 久久久精品免费免费高清| 热re99久久精品国产66热6| 日日爽夜夜爽网站| 国产一区亚洲一区在线观看| 国产女主播在线喷水免费视频网站| 老司机影院毛片| 久久精品国产亚洲av天美| 国内精品宾馆在线| 18禁在线播放成人免费| 日本wwww免费看| 亚洲精品亚洲一区二区| 伊人久久精品亚洲午夜| 尾随美女入室| 高清午夜精品一区二区三区| 国产午夜精品久久久久久一区二区三区| 男人操女人黄网站| 三级国产精品欧美在线观看| 亚洲久久久国产精品| 久久女婷五月综合色啪小说| 国产精品国产三级专区第一集| 国产午夜精品久久久久久一区二区三区| 欧美xxⅹ黑人| 少妇高潮的动态图| 久久毛片免费看一区二区三区| 久久精品国产亚洲av天美| 建设人人有责人人尽责人人享有的| 国模一区二区三区四区视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产成人一精品久久久| 日韩 亚洲 欧美在线| 嫩草影院入口| 国产视频首页在线观看| av免费在线看不卡| 一区在线观看完整版| 波野结衣二区三区在线| 尾随美女入室| 久久狼人影院| 99热这里只有精品一区| 人妻系列 视频| 日本欧美国产在线视频| 国产精品久久久久久av不卡| 久久人人爽av亚洲精品天堂| 多毛熟女@视频| 一级片'在线观看视频| 少妇人妻久久综合中文| 精品久久久久久久久亚洲| 亚洲一区二区三区欧美精品| 日本黄大片高清| 国产白丝娇喘喷水9色精品| 国产精品不卡视频一区二区| 成年女人在线观看亚洲视频| 亚洲图色成人| videos熟女内射| 国产成人aa在线观看| 熟女人妻精品中文字幕| 王馨瑶露胸无遮挡在线观看| 国产精品一区二区在线不卡| 高清av免费在线| 中文字幕最新亚洲高清| 国产精品国产三级国产av玫瑰| 最新中文字幕久久久久| 在线观看免费日韩欧美大片 | 成人漫画全彩无遮挡| 日韩大片免费观看网站| 国产一区二区在线观看日韩| 久久99精品国语久久久| av.在线天堂| 国产成人精品福利久久| 日韩在线高清观看一区二区三区| 精品少妇黑人巨大在线播放| 看十八女毛片水多多多| 亚洲欧美成人精品一区二区| av线在线观看网站| 建设人人有责人人尽责人人享有的| 婷婷色综合www| 亚洲精品国产色婷婷电影| 七月丁香在线播放| 日本av免费视频播放| 国产av码专区亚洲av| 黄色一级大片看看| 精品人妻在线不人妻| 婷婷成人精品国产| 日韩亚洲欧美综合| av天堂久久9| 搡女人真爽免费视频火全软件| 你懂的网址亚洲精品在线观看| 欧美人与善性xxx| 亚洲人成网站在线观看播放| 一个人免费看片子| 欧美成人午夜免费资源| 秋霞在线观看毛片| 少妇猛男粗大的猛烈进出视频| 国产黄片视频在线免费观看| 青春草亚洲视频在线观看| 亚洲综合色网址| 免费播放大片免费观看视频在线观看| 97超视频在线观看视频| av一本久久久久| 好男人视频免费观看在线| 九色亚洲精品在线播放| 飞空精品影院首页| 亚洲少妇的诱惑av| 欧美+日韩+精品| 中国美白少妇内射xxxbb| 女性生殖器流出的白浆| 亚洲第一区二区三区不卡| 免费高清在线观看视频在线观看| 少妇被粗大猛烈的视频| 99久久中文字幕三级久久日本| 中文精品一卡2卡3卡4更新| 夫妻午夜视频| 久久99热这里只频精品6学生| 日韩欧美精品免费久久| 水蜜桃什么品种好| 欧美丝袜亚洲另类| 成人综合一区亚洲| 久久久国产精品麻豆| 亚洲色图 男人天堂 中文字幕 | av专区在线播放| 免费少妇av软件| 你懂的网址亚洲精品在线观看| 又大又黄又爽视频免费| 大码成人一级视频| 我的女老师完整版在线观看| 国产精品人妻久久久影院| 亚洲内射少妇av| 亚洲欧洲精品一区二区精品久久久 | 80岁老熟妇乱子伦牲交| 在线播放无遮挡| 不卡视频在线观看欧美| 高清午夜精品一区二区三区| 我的女老师完整版在线观看| 男男h啪啪无遮挡| 欧美成人午夜免费资源| 大香蕉久久网| 少妇丰满av| 观看美女的网站| 高清欧美精品videossex| 免费黄网站久久成人精品| 欧美激情极品国产一区二区三区 | 久久精品久久久久久噜噜老黄| 中文字幕最新亚洲高清| 国产 精品1| av不卡在线播放| 最近的中文字幕免费完整| 人妻人人澡人人爽人人| 日本91视频免费播放| 免费高清在线观看日韩| 欧美最新免费一区二区三区| 99国产综合亚洲精品| 亚洲精品久久久久久婷婷小说| 日韩中文字幕视频在线看片| 国产成人a∨麻豆精品| 麻豆精品久久久久久蜜桃| 香蕉精品网在线| 麻豆乱淫一区二区| 18禁动态无遮挡网站| 亚洲丝袜综合中文字幕| 男的添女的下面高潮视频| 亚洲国产欧美在线一区| 亚洲成人手机| 啦啦啦啦在线视频资源| 久久 成人 亚洲| 国产亚洲精品久久久com| 国产无遮挡羞羞视频在线观看| 在线观看人妻少妇| av视频免费观看在线观看| 亚洲第一区二区三区不卡| 777米奇影视久久| 久久 成人 亚洲| 国产精品久久久久久精品古装| 亚洲精品一二三| 欧美精品人与动牲交sv欧美| a级毛片在线看网站| 欧美人与性动交α欧美精品济南到 | 又粗又硬又长又爽又黄的视频| 我要看黄色一级片免费的| 乱人伦中国视频| 免费观看av网站的网址| 久久久久久久久大av| 亚洲欧美清纯卡通| 天堂俺去俺来也www色官网| 国产乱人偷精品视频| 国产欧美亚洲国产| 亚洲成色77777| 视频区图区小说| 亚洲综合色网址| av天堂久久9| 久久 成人 亚洲| 亚洲精品乱码久久久v下载方式| 最近中文字幕2019免费版| 纵有疾风起免费观看全集完整版| 国产 一区精品| 少妇被粗大猛烈的视频| 夜夜骑夜夜射夜夜干| 蜜桃久久精品国产亚洲av| 99视频精品全部免费 在线| 18禁裸乳无遮挡动漫免费视频| 有码 亚洲区| 国产成人免费无遮挡视频| 久久精品国产亚洲av天美| videosex国产| 精品午夜福利在线看| 精品人妻在线不人妻| 免费av不卡在线播放| 欧美人与性动交α欧美精品济南到 | 亚洲av综合色区一区| 欧美激情国产日韩精品一区| 满18在线观看网站| 亚洲内射少妇av| 中文精品一卡2卡3卡4更新| 久久久久精品久久久久真实原创| 涩涩av久久男人的天堂| 满18在线观看网站| av天堂久久9| 亚洲精品国产av蜜桃| 老司机影院毛片| 九色成人免费人妻av| 热re99久久精品国产66热6| 色94色欧美一区二区| 亚洲国产欧美在线一区| 建设人人有责人人尽责人人享有的| 看免费成人av毛片| 国产精品国产三级国产av玫瑰| 亚洲人与动物交配视频| www.av在线官网国产| 日韩中文字幕视频在线看片| 国产精品国产三级国产专区5o| 国产白丝娇喘喷水9色精品| 亚洲精品乱久久久久久| 韩国av在线不卡| 女人久久www免费人成看片| 七月丁香在线播放| 熟女人妻精品中文字幕| 欧美激情极品国产一区二区三区 | 国产一区二区在线观看av| 国产精品久久久久久精品电影小说| 久久精品国产鲁丝片午夜精品| 成人毛片a级毛片在线播放| 国产片内射在线| 国产无遮挡羞羞视频在线观看| 伦理电影免费视频| 欧美日韩av久久| 久久av网站| 男人操女人黄网站| 精品少妇内射三级| 97在线人人人人妻| 狠狠婷婷综合久久久久久88av| 亚洲伊人久久精品综合| 亚洲国产最新在线播放| 国产精品不卡视频一区二区| 成人国语在线视频| 简卡轻食公司| 欧美+日韩+精品| 自拍欧美九色日韩亚洲蝌蚪91| 一边摸一边做爽爽视频免费| 人妻制服诱惑在线中文字幕| 国产国语露脸激情在线看| 热99国产精品久久久久久7| 精品人妻在线不人妻| 哪个播放器可以免费观看大片| 亚洲四区av| 考比视频在线观看| 亚洲经典国产精华液单| 亚洲人成网站在线观看播放| 飞空精品影院首页| 你懂的网址亚洲精品在线观看| 丰满乱子伦码专区| 在线观看www视频免费| 成人亚洲精品一区在线观看| 多毛熟女@视频| 乱码一卡2卡4卡精品| 亚洲精品日韩在线中文字幕| 国产高清有码在线观看视频| 日本欧美视频一区| 一区在线观看完整版| 亚洲欧洲精品一区二区精品久久久 | 草草在线视频免费看| 国语对白做爰xxxⅹ性视频网站| 在线观看人妻少妇| 亚洲av中文av极速乱| 国产精品女同一区二区软件| 亚洲国产毛片av蜜桃av| 精品国产露脸久久av麻豆| 菩萨蛮人人尽说江南好唐韦庄| 日韩亚洲欧美综合| 久久久精品免费免费高清| 亚洲精品色激情综合| 欧美xxxx性猛交bbbb| 伊人久久精品亚洲午夜| av免费在线看不卡| 性色av一级| 久久免费观看电影| 久久97久久精品| 一级毛片aaaaaa免费看小| 日韩不卡一区二区三区视频在线| 精品人妻熟女av久视频| 久久久国产精品麻豆| 日韩av免费高清视频| 超碰97精品在线观看| 欧美精品人与动牲交sv欧美| 热re99久久国产66热| 午夜免费男女啪啪视频观看| 亚洲一级一片aⅴ在线观看| 久久久久国产精品人妻一区二区| 女人久久www免费人成看片| 国产一区亚洲一区在线观看| 久久久久久久久大av| 桃花免费在线播放| 99久国产av精品国产电影| 伊人久久精品亚洲午夜| 精品国产一区二区久久| 国内精品宾馆在线| 国产一区二区三区综合在线观看 | 国产国拍精品亚洲av在线观看| 亚洲av综合色区一区| 欧美三级亚洲精品| 国产成人av激情在线播放 | a级毛片在线看网站| 亚洲四区av| .国产精品久久| 丰满乱子伦码专区| 亚洲av在线观看美女高潮| 日本色播在线视频| 曰老女人黄片| 免费高清在线观看日韩| 久久久久人妻精品一区果冻| 少妇精品久久久久久久| 亚洲欧洲国产日韩| 99热这里只有是精品在线观看| 中文天堂在线官网| 91久久精品国产一区二区成人| 国产色爽女视频免费观看| 亚洲美女黄色视频免费看| 国产精品久久久久久精品电影小说| 免费大片黄手机在线观看| 少妇的逼好多水| 在线观看人妻少妇| 蜜臀久久99精品久久宅男| 一二三四中文在线观看免费高清| 日韩电影二区| 欧美日韩视频精品一区| 日韩 亚洲 欧美在线| 自线自在国产av| 欧美一级a爱片免费观看看| 亚洲精品久久久久久婷婷小说| 国产免费视频播放在线视频| 婷婷色av中文字幕| 大话2 男鬼变身卡| 精品卡一卡二卡四卡免费| 精品视频人人做人人爽| 日韩中字成人| 国产免费视频播放在线视频| 国产精品不卡视频一区二区| 亚洲av.av天堂| 精品亚洲乱码少妇综合久久| 国产精品无大码| videossex国产| 97超碰精品成人国产| 99久久中文字幕三级久久日本| 热re99久久精品国产66热6| 你懂的网址亚洲精品在线观看| 午夜久久久在线观看| 十八禁高潮呻吟视频| 亚洲av国产av综合av卡| 91在线精品国自产拍蜜月| 大又大粗又爽又黄少妇毛片口| 26uuu在线亚洲综合色| 亚洲五月色婷婷综合| 老司机影院毛片| 国产免费一区二区三区四区乱码| 成人午夜精彩视频在线观看| 日韩不卡一区二区三区视频在线| 国产精品久久久久久av不卡| 夜夜爽夜夜爽视频| 26uuu在线亚洲综合色| 日韩一区二区三区影片| 女的被弄到高潮叫床怎么办| 美女cb高潮喷水在线观看| 精品少妇黑人巨大在线播放| 一级毛片电影观看| 黑人巨大精品欧美一区二区蜜桃 | 最近最新中文字幕免费大全7| 久久午夜综合久久蜜桃| 黄色怎么调成土黄色| 乱人伦中国视频| 美女主播在线视频| a 毛片基地| 免费观看av网站的网址| 一二三四中文在线观看免费高清| 色5月婷婷丁香| 一级毛片电影观看| 精品午夜福利在线看| 一区二区av电影网| videossex国产| 在线播放无遮挡| 国产一区亚洲一区在线观看| 国产黄频视频在线观看| 午夜福利网站1000一区二区三区| 我的老师免费观看完整版| 狠狠婷婷综合久久久久久88av| 纵有疾风起免费观看全集完整版| 9色porny在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精华国产精华液的使用体验| 国产熟女午夜一区二区三区 | a级毛片在线看网站| 天堂中文最新版在线下载| 久久免费观看电影| 一区二区三区四区激情视频| av黄色大香蕉| av在线观看视频网站免费| 91久久精品国产一区二区成人| 中文字幕久久专区| 日本av免费视频播放| 一本一本综合久久| 日日摸夜夜添夜夜添av毛片| 最新中文字幕久久久久| 午夜福利视频精品| 久久久国产欧美日韩av| freevideosex欧美| 国产免费一区二区三区四区乱码| 久久亚洲国产成人精品v| 一级毛片 在线播放| 亚洲不卡免费看| 男女免费视频国产| 嫩草影院入口| 免费大片黄手机在线观看| 青青草视频在线视频观看| 熟女av电影| 亚洲一级一片aⅴ在线观看| 精品人妻一区二区三区麻豆| 91国产中文字幕| 久久精品夜色国产| 亚洲av欧美aⅴ国产| 大片免费播放器 马上看| 岛国毛片在线播放| 久久午夜综合久久蜜桃| 黄色欧美视频在线观看| 久久久久视频综合| 久久久a久久爽久久v久久| av免费在线看不卡| 婷婷色av中文字幕| 街头女战士在线观看网站| 久久精品久久久久久久性| 蜜臀久久99精品久久宅男| 秋霞在线观看毛片| 国产精品国产三级专区第一集| a 毛片基地| 蜜臀久久99精品久久宅男| 一边摸一边做爽爽视频免费| 日产精品乱码卡一卡2卡三| 国产精品久久久久成人av| 国产女主播在线喷水免费视频网站| 免费日韩欧美在线观看| 国产精品久久久久久久电影| 我要看黄色一级片免费的| 美女脱内裤让男人舔精品视频| 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 精品卡一卡二卡四卡免费| 美女脱内裤让男人舔精品视频| 春色校园在线视频观看| 制服丝袜香蕉在线| 国产成人aa在线观看| 日产精品乱码卡一卡2卡三| 亚洲av综合色区一区| 国产欧美日韩综合在线一区二区| 中文字幕免费在线视频6| 日本wwww免费看| tube8黄色片| 黑人高潮一二区| 国产精品麻豆人妻色哟哟久久| 亚洲av.av天堂| 99热这里只有精品一区| av免费观看日本| 青春草国产在线视频| av网站免费在线观看视频| 国模一区二区三区四区视频| 蜜桃在线观看..| 日韩中文字幕视频在线看片| 亚洲人成网站在线播| 亚洲av国产av综合av卡| 国产精品熟女久久久久浪| 久久99一区二区三区| 亚洲高清免费不卡视频| 五月开心婷婷网| 国产亚洲精品久久久com| a级毛片在线看网站| 久久精品久久久久久久性| 哪个播放器可以免费观看大片| 99热全是精品| 亚洲av日韩在线播放| 大香蕉97超碰在线| 亚洲人成77777在线视频| 一本大道久久a久久精品| 制服丝袜香蕉在线| 曰老女人黄片| 街头女战士在线观看网站| 乱人伦中国视频| 啦啦啦在线观看免费高清www| 少妇的逼水好多| 日韩成人av中文字幕在线观看| 国产精品.久久久| 久久午夜福利片| 亚洲人成77777在线视频| 色婷婷久久久亚洲欧美| 91久久精品电影网| 水蜜桃什么品种好| 亚洲美女黄色视频免费看| 亚洲国产精品国产精品| 日韩不卡一区二区三区视频在线| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产三级国产av玫瑰| 免费人成在线观看视频色| a级毛片免费高清观看在线播放| 赤兔流量卡办理| 亚洲情色 制服丝袜| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av一本久久久久| 日韩精品免费视频一区二区三区 | 男女边摸边吃奶| 亚洲欧美中文字幕日韩二区| 日本-黄色视频高清免费观看| 大片免费播放器 马上看| 亚洲人成77777在线视频| 99久久精品国产国产毛片| av有码第一页| 亚洲成人手机| 永久网站在线| 欧美 亚洲 国产 日韩一| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 黑丝袜美女国产一区| av天堂久久9| 另类亚洲欧美激情| 男女国产视频网站| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡 | 亚洲精品日韩在线中文字幕| 美女国产高潮福利片在线看| 国产免费福利视频在线观看| 黄色视频在线播放观看不卡| 日韩一本色道免费dvd| 水蜜桃什么品种好| 国产黄频视频在线观看| 亚洲av不卡在线观看| 搡女人真爽免费视频火全软件| 久久精品国产亚洲av天美|