• <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)收獲
    久久久久久久久大av| 国产精品人妻久久久影院| 亚洲欧美日韩东京热| 国产精华一区二区三区| 婷婷六月久久综合丁香| 成人高潮视频无遮挡免费网站| 久99久视频精品免费| 国产成人免费观看mmmm| 两性午夜刺激爽爽歪歪视频在线观看| 在线免费十八禁| 国产私拍福利视频在线观看| 在线观看av片永久免费下载| kizo精华| 蜜桃亚洲精品一区二区三区| 深爱激情五月婷婷| 搡女人真爽免费视频火全软件| 午夜福利网站1000一区二区三区| 九九在线视频观看精品| 午夜日本视频在线| 天天一区二区日本电影三级| 精品99又大又爽又粗少妇毛片| 人体艺术视频欧美日本| 国产高清三级在线| 直男gayav资源| 在线观看美女被高潮喷水网站| 美女大奶头视频| 日本免费在线观看一区| 波多野结衣巨乳人妻| 亚洲一区高清亚洲精品| 亚洲av一区综合| 国产精品国产三级专区第一集| 99热精品在线国产| 熟女人妻精品中文字幕| 舔av片在线| 美女大奶头视频| 国产精品.久久久| 尤物成人国产欧美一区二区三区| 亚洲va在线va天堂va国产| 亚洲av福利一区| 亚洲精品久久久久久婷婷小说 | 国产激情偷乱视频一区二区| 日韩一区二区三区影片| 国产在视频线精品| 精品一区二区三区视频在线| 久久人人爽人人片av| 特大巨黑吊av在线直播| 成人亚洲精品av一区二区| 国产精品蜜桃在线观看| 亚洲欧美精品综合久久99| 特大巨黑吊av在线直播| 国产在视频线在精品| 国产激情偷乱视频一区二区| 国产午夜精品论理片| 成年版毛片免费区| 一个人看的www免费观看视频| 色网站视频免费| 精品酒店卫生间| 搞女人的毛片| 女人久久www免费人成看片 | 色哟哟·www| 综合色丁香网| 中文乱码字字幕精品一区二区三区 | 在线天堂最新版资源| 午夜福利高清视频| 欧美又色又爽又黄视频| 日本五十路高清| 亚洲精品,欧美精品| 人人妻人人澡欧美一区二区| 亚洲精品久久久久久婷婷小说 | 色综合站精品国产| 91精品伊人久久大香线蕉| 99久久精品国产国产毛片| 简卡轻食公司| 热99re8久久精品国产| 久久精品国产亚洲网站| 男人的好看免费观看在线视频| 51国产日韩欧美| 啦啦啦啦在线视频资源| av黄色大香蕉| 国产在视频线在精品| 啦啦啦观看免费观看视频高清| 美女被艹到高潮喷水动态| 黄色欧美视频在线观看| 看片在线看免费视频| 国产色婷婷99| 18+在线观看网站| 国产精品一区二区在线观看99 | 国产乱来视频区| a级毛色黄片| 午夜爱爱视频在线播放| 欧美丝袜亚洲另类| 国产乱来视频区| 国产老妇伦熟女老妇高清| 亚洲欧美中文字幕日韩二区| 亚洲最大成人手机在线| 国模一区二区三区四区视频| 亚洲高清免费不卡视频| 卡戴珊不雅视频在线播放| 九九热线精品视视频播放| 观看美女的网站| 亚洲精品色激情综合| ponron亚洲| 国内精品美女久久久久久| 日韩,欧美,国产一区二区三区 | 国产黄a三级三级三级人| 91久久精品国产一区二区三区| 日韩成人av中文字幕在线观看| 成人无遮挡网站| 麻豆乱淫一区二区| 乱码一卡2卡4卡精品| 中文字幕av在线有码专区| 亚洲美女搞黄在线观看| 久久久久久久久大av| 美女脱内裤让男人舔精品视频| 国产亚洲av片在线观看秒播厂 | 日日啪夜夜撸| 国语自产精品视频在线第100页| 中国美白少妇内射xxxbb| 亚洲天堂国产精品一区在线| 久久精品久久久久久久性| 最近手机中文字幕大全| 成年女人永久免费观看视频| 黄片wwwwww| 国内少妇人妻偷人精品xxx网站| 国产伦一二天堂av在线观看| 美女被艹到高潮喷水动态| 精品久久久久久久久久久久久| 久久精品夜色国产| 丰满少妇做爰视频| 淫秽高清视频在线观看| 亚洲国产精品专区欧美| 人人妻人人澡欧美一区二区| 人体艺术视频欧美日本| 少妇熟女aⅴ在线视频| a级一级毛片免费在线观看| 精品人妻一区二区三区麻豆| 99热这里只有是精品在线观看| 亚洲欧美中文字幕日韩二区| 女人十人毛片免费观看3o分钟| 精品久久久久久久久av| 亚洲成人精品中文字幕电影| 久久精品91蜜桃| 国产亚洲精品av在线| 亚洲av成人精品一区久久| 麻豆av噜噜一区二区三区| 少妇人妻精品综合一区二区| 日韩成人av中文字幕在线观看| 亚洲av福利一区| 美女黄网站色视频| 国产精品无大码| 中文天堂在线官网| 亚洲第一区二区三区不卡| 免费无遮挡裸体视频| 亚洲欧美成人精品一区二区| 一级毛片aaaaaa免费看小| 亚洲不卡免费看| 桃色一区二区三区在线观看| 国产美女午夜福利| 春色校园在线视频观看| 一区二区三区免费毛片| 久久精品久久精品一区二区三区| or卡值多少钱| 国国产精品蜜臀av免费| 免费不卡的大黄色大毛片视频在线观看 | 三级国产精品片| 久久久久网色| 一本久久精品| 丰满少妇做爰视频| 18禁裸乳无遮挡免费网站照片| 女人十人毛片免费观看3o分钟| 嫩草影院新地址| or卡值多少钱| 天堂√8在线中文| av免费在线看不卡| www日本黄色视频网| 日本五十路高清| 国产三级在线视频| 日日撸夜夜添| 观看美女的网站| 亚洲最大成人手机在线| 99热这里只有精品一区| 嘟嘟电影网在线观看| 简卡轻食公司| 少妇猛男粗大的猛烈进出视频 | 成人性生交大片免费视频hd| 亚洲av免费高清在线观看| 国产精品久久视频播放| 少妇裸体淫交视频免费看高清| 国产精品女同一区二区软件| 久久精品综合一区二区三区| 亚洲国产欧美在线一区| 七月丁香在线播放| 国产成人免费观看mmmm| 久久久久久伊人网av| 美女国产视频在线观看| 热99在线观看视频| 亚洲美女视频黄频| 国产精品人妻久久久影院| 国产成人精品婷婷| 国国产精品蜜臀av免费| 午夜老司机福利剧场| 国产视频内射| 亚洲欧美精品综合久久99| 国产老妇伦熟女老妇高清| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久久久| 精品不卡国产一区二区三区| 亚洲欧美日韩无卡精品| 久久久久久国产a免费观看| 性插视频无遮挡在线免费观看| 欧美3d第一页| 中文在线观看免费www的网站| 亚洲国产精品国产精品| 秋霞伦理黄片| 赤兔流量卡办理| 亚洲成av人片在线播放无| 人妻制服诱惑在线中文字幕| 国产极品精品免费视频能看的| 美女国产视频在线观看| 中国国产av一级| 久久草成人影院| 国产伦一二天堂av在线观看| 亚洲国产精品sss在线观看| 狠狠狠狠99中文字幕| 午夜精品国产一区二区电影 | 国产视频首页在线观看| 99热这里只有是精品在线观看| 亚洲,欧美,日韩| 成人毛片a级毛片在线播放| 国产精品国产三级国产av玫瑰| 国产黄片视频在线免费观看| 日本一二三区视频观看| 久久精品国产99精品国产亚洲性色| 久久久久久久久久久丰满| 国产一区有黄有色的免费视频 | 成年女人永久免费观看视频| 欧美性猛交╳xxx乱大交人| 亚洲aⅴ乱码一区二区在线播放| 国产成人aa在线观看| 久久这里有精品视频免费| 成人毛片60女人毛片免费| 丰满人妻一区二区三区视频av| 91久久精品国产一区二区成人| 建设人人有责人人尽责人人享有的 | 国模一区二区三区四区视频| 亚洲三级黄色毛片| 日韩欧美 国产精品| 日本wwww免费看| 性色avwww在线观看| 蜜桃久久精品国产亚洲av| 日本黄色视频三级网站网址| 我的女老师完整版在线观看| 久久久色成人| 久久精品熟女亚洲av麻豆精品 | 一个人看的www免费观看视频| 亚洲欧洲国产日韩| 99久久精品国产国产毛片| 久久久久久久久久久丰满| 亚洲国产高清在线一区二区三| 久久久久久久久久黄片| 午夜精品在线福利| 国产大屁股一区二区在线视频| 中文字幕久久专区| 亚洲激情五月婷婷啪啪| 我要搜黄色片| 免费av毛片视频| 乱人视频在线观看| 欧美丝袜亚洲另类| 看片在线看免费视频| 国产av在哪里看| 女人十人毛片免费观看3o分钟| 国产av一区在线观看免费| 三级国产精品欧美在线观看| 男女国产视频网站| 亚洲久久久久久中文字幕| 夜夜看夜夜爽夜夜摸| 久久久久久久午夜电影| 韩国高清视频一区二区三区| 成年版毛片免费区| 亚洲精品色激情综合| 少妇人妻精品综合一区二区| 亚洲在久久综合| 久久婷婷人人爽人人干人人爱| 男插女下体视频免费在线播放| 国产伦精品一区二区三区四那| 亚洲国产精品sss在线观看| 午夜福利在线观看吧| 一级毛片aaaaaa免费看小| 日本午夜av视频| 99热全是精品| 欧美又色又爽又黄视频| 三级毛片av免费| 日本黄大片高清| 久久午夜福利片| 人妻制服诱惑在线中文字幕| 久久精品久久精品一区二区三区| 国产免费男女视频| 日本免费a在线| 国产精品野战在线观看| 有码 亚洲区| 国产大屁股一区二区在线视频| 村上凉子中文字幕在线| 国模一区二区三区四区视频| 淫秽高清视频在线观看| 国产精品一区www在线观看| 在线播放无遮挡| 草草在线视频免费看| 热99re8久久精品国产| 不卡视频在线观看欧美| 精品一区二区免费观看| 亚洲高清免费不卡视频| 精品国内亚洲2022精品成人| 亚洲精品成人久久久久久| videossex国产| 老女人水多毛片| 亚洲av二区三区四区| 久久精品久久久久久噜噜老黄 | 国产精品久久视频播放| 午夜久久久久精精品| 国语自产精品视频在线第100页| 色噜噜av男人的天堂激情| 观看美女的网站| 国产精品国产三级国产专区5o | 亚洲av.av天堂| 美女内射精品一级片tv| 色尼玛亚洲综合影院| 又粗又爽又猛毛片免费看| 国产国拍精品亚洲av在线观看| 在现免费观看毛片| 美女被艹到高潮喷水动态| ponron亚洲| 午夜亚洲福利在线播放| 国产精品久久久久久久久免| 国产伦一二天堂av在线观看| www日本黄色视频网| 国产片特级美女逼逼视频| 不卡视频在线观看欧美| 热99在线观看视频| 深爱激情五月婷婷| 热99re8久久精品国产| 国产一区有黄有色的免费视频 | 精品久久久久久久久久久久久| 国产午夜福利久久久久久| АⅤ资源中文在线天堂| 美女高潮的动态| 国产毛片a区久久久久| 亚洲欧美成人综合另类久久久 | 三级男女做爰猛烈吃奶摸视频| 国产成人精品婷婷| 亚洲成人久久爱视频| 日韩欧美在线乱码| 97人妻精品一区二区三区麻豆| 性插视频无遮挡在线免费观看| 久久久精品大字幕| 亚洲激情五月婷婷啪啪| 国产黄色小视频在线观看| 观看免费一级毛片| 国产三级中文精品| 91精品伊人久久大香线蕉| 欧美97在线视频| av国产免费在线观看| 国产精品熟女久久久久浪| 久久久欧美国产精品| 精品人妻视频免费看| 少妇裸体淫交视频免费看高清| 又爽又黄a免费视频| 国产乱人视频| 亚洲aⅴ乱码一区二区在线播放| 最近最新中文字幕大全电影3| 婷婷色麻豆天堂久久 | 久久久久久久久久久丰满| 中文亚洲av片在线观看爽| 日韩制服骚丝袜av| 成人无遮挡网站| 91狼人影院| 国产私拍福利视频在线观看| 秋霞在线观看毛片| 久久久久国产网址| 五月玫瑰六月丁香| 国产精品国产三级国产专区5o | 狠狠狠狠99中文字幕| 国产色爽女视频免费观看| 国产熟女欧美一区二区| 亚洲内射少妇av| 免费看av在线观看网站| 国产成人午夜福利电影在线观看| 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| 久久久久久久国产电影| 永久网站在线| 草草在线视频免费看| 99久国产av精品国产电影| 韩国av在线不卡| 日韩视频在线欧美| 三级国产精品欧美在线观看| 69人妻影院| 久久精品国产亚洲av天美| 伦理电影大哥的女人| 亚洲欧美日韩东京热| 国产精品,欧美在线| 国产真实乱freesex| 欧美激情国产日韩精品一区| 中文欧美无线码| 国产精品久久视频播放| 可以在线观看毛片的网站| 国产不卡一卡二| 久久久久网色| 精品少妇黑人巨大在线播放 | 男人的好看免费观看在线视频| 蜜桃久久精品国产亚洲av| 色综合色国产| 成人特级av手机在线观看| 亚洲欧美精品综合久久99| av.在线天堂| 嫩草影院新地址| 高清在线视频一区二区三区 | АⅤ资源中文在线天堂| 晚上一个人看的免费电影| 在线观看av片永久免费下载| 国产亚洲av嫩草精品影院| 免费看光身美女| 国产高清国产精品国产三级 | 久久鲁丝午夜福利片| 欧美性猛交黑人性爽| 蜜桃久久精品国产亚洲av| 亚洲精品久久久久久婷婷小说 | 观看免费一级毛片| 久久99蜜桃精品久久| 一级毛片我不卡| 成人午夜精彩视频在线观看| 精品国产三级普通话版| 亚洲国产色片| 免费av不卡在线播放| 一级毛片久久久久久久久女| 99热精品在线国产| 草草在线视频免费看| 建设人人有责人人尽责人人享有的 | 日韩欧美精品免费久久| 人人妻人人看人人澡| 亚洲乱码一区二区免费版| 1024手机看黄色片| 亚洲,欧美,日韩| 亚洲自偷自拍三级| 亚洲精品乱码久久久v下载方式| 91午夜精品亚洲一区二区三区| videossex国产| 亚洲国产日韩欧美精品在线观看| 中文欧美无线码| 亚洲欧美日韩高清专用| h日本视频在线播放| 久久人人爽人人爽人人片va| 99久国产av精品| 亚洲av日韩在线播放| 少妇裸体淫交视频免费看高清| 精品久久久噜噜| 久久99热这里只有精品18| 99热全是精品| 久久久久久伊人网av| 热99re8久久精品国产| 国语自产精品视频在线第100页| 国产国拍精品亚洲av在线观看| 最近2019中文字幕mv第一页| 亚洲乱码一区二区免费版| 淫秽高清视频在线观看| 精品欧美国产一区二区三| 亚洲电影在线观看av| 最后的刺客免费高清国语| 亚洲色图av天堂| 97热精品久久久久久| 九九爱精品视频在线观看| 日韩欧美 国产精品| 成人无遮挡网站| 少妇熟女aⅴ在线视频| 欧美一区二区亚洲| av在线老鸭窝| eeuss影院久久| 精品人妻熟女av久视频| 国产精品国产高清国产av| 天美传媒精品一区二区| 欧美+日韩+精品| 一个人免费在线观看电影| 高清在线视频一区二区三区 | 国内揄拍国产精品人妻在线| 老司机福利观看| 哪个播放器可以免费观看大片| 久久99热6这里只有精品| 在线天堂最新版资源| 日韩亚洲欧美综合| 亚洲av二区三区四区| 国产黄片视频在线免费观看| 国产又色又爽无遮挡免| 激情 狠狠 欧美| 国产高清三级在线| 国产精品一及| 久久人妻av系列| 丰满人妻一区二区三区视频av| 一级av片app| 日产精品乱码卡一卡2卡三| 免费在线观看成人毛片| 亚洲国产精品sss在线观看| 国产精品蜜桃在线观看| 国产一区二区三区av在线| 欧美xxxx黑人xx丫x性爽| 欧美激情久久久久久爽电影| 岛国在线免费视频观看| 99久久九九国产精品国产免费| 亚洲婷婷狠狠爱综合网| 亚洲成人av在线免费| 亚洲va在线va天堂va国产| 国产极品精品免费视频能看的| 午夜精品在线福利| 夜夜看夜夜爽夜夜摸| 日日摸夜夜添夜夜爱| 欧美激情久久久久久爽电影| 一夜夜www| 久久久久国产网址| 日日摸夜夜添夜夜爱| 成年版毛片免费区| 欧美成人a在线观看| 色网站视频免费| 亚洲精品国产av成人精品| 亚洲电影在线观看av| 国产老妇伦熟女老妇高清| 成人毛片a级毛片在线播放| 国产麻豆成人av免费视频| 精品人妻偷拍中文字幕| 狂野欧美白嫩少妇大欣赏| 神马国产精品三级电影在线观看| 成年免费大片在线观看| 国产 一区精品| 欧美另类亚洲清纯唯美| 久久精品夜夜夜夜夜久久蜜豆| 日本三级黄在线观看| 高清日韩中文字幕在线| 国产成人精品一,二区| 国产成人精品婷婷| 精品久久久久久久人妻蜜臀av| 欧美性猛交黑人性爽| 国产精品不卡视频一区二区| 舔av片在线| 国产伦在线观看视频一区| 最近视频中文字幕2019在线8| 91精品一卡2卡3卡4卡| 听说在线观看完整版免费高清| 国产精品99久久久久久久久| 欧美一区二区亚洲| 亚洲欧美日韩东京热| 久久韩国三级中文字幕| 成人亚洲欧美一区二区av| 久久国内精品自在自线图片| 成人毛片60女人毛片免费| 国产麻豆成人av免费视频| 自拍偷自拍亚洲精品老妇| 国产v大片淫在线免费观看| 黄片wwwwww| 天堂影院成人在线观看| 中文精品一卡2卡3卡4更新| 国产在线男女| 我要看日韩黄色一级片| 国产精品乱码一区二三区的特点| 女人十人毛片免费观看3o分钟| .国产精品久久| 午夜福利视频1000在线观看| 国产极品天堂在线| 十八禁国产超污无遮挡网站| 久久99热这里只频精品6学生 | 超碰av人人做人人爽久久| 国产成人精品久久久久久| 嘟嘟电影网在线观看| 国产黄a三级三级三级人| 一个人看的www免费观看视频| 少妇熟女欧美另类| 亚州av有码| 黄色一级大片看看| 99热6这里只有精品| 国产精品国产三级国产av玫瑰| 一个人看视频在线观看www免费| av女优亚洲男人天堂| 热99在线观看视频| 日日摸夜夜添夜夜添av毛片| 在现免费观看毛片| 亚洲美女视频黄频| 91狼人影院| 亚洲国产精品成人久久小说| 久久婷婷人人爽人人干人人爱| 尾随美女入室| 中文乱码字字幕精品一区二区三区 | 精品久久久久久电影网 | 一级毛片我不卡| 国产真实伦视频高清在线观看| 国产精品电影一区二区三区| 国模一区二区三区四区视频| 亚洲国产精品sss在线观看| 亚洲丝袜综合中文字幕| 一边摸一边抽搐一进一小说| 成人三级黄色视频| 赤兔流量卡办理| 国产亚洲精品久久久com| 最新中文字幕久久久久| 人人妻人人看人人澡| 亚洲四区av| 久久这里只有精品中国| 麻豆成人av视频| 国产精品嫩草影院av在线观看| 亚洲欧美成人综合另类久久久 | 一本久久精品| 国产 一区精品| 中文字幕av在线有码专区| 国产亚洲精品久久久com| 成人国产麻豆网| 高清视频免费观看一区二区 | 国产成人午夜福利电影在线观看| 大香蕉97超碰在线| 草草在线视频免费看| 精品少妇黑人巨大在线播放 |