• <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)收獲
    九九爱精品视频在线观看| 国产 一区 欧美 日韩| 高清毛片免费看| 亚洲精华国产精华液的使用体验| 久久久亚洲精品成人影院| 精品一区二区免费观看| 精品久久久久久久末码| 91精品国产九色| 人人妻人人看人人澡| 欧美日韩精品成人综合77777| 久久久久久久亚洲中文字幕| 亚洲av免费在线观看| 精品久久久久久成人av| 亚洲国产精品国产精品| 亚洲欧美清纯卡通| 69av精品久久久久久| 99久国产av精品国产电影| 老司机影院成人| 一个人看的www免费观看视频| 好男人在线观看高清免费视频| 黄色日韩在线| 欧美成人免费av一区二区三区| 亚洲一级一片aⅴ在线观看| 精品久久久久久久人妻蜜臀av| 国产成人福利小说| 欧美不卡视频在线免费观看| 色尼玛亚洲综合影院| 久久综合国产亚洲精品| 身体一侧抽搐| 国产91av在线免费观看| 热99在线观看视频| 极品教师在线视频| 国产欧美另类精品又又久久亚洲欧美| 欧美丝袜亚洲另类| 国产午夜精品久久久久久一区二区三区| 日韩一区二区视频免费看| 久久韩国三级中文字幕| 久久99精品国语久久久| 亚洲在线自拍视频| 亚洲高清免费不卡视频| 国产一区二区在线观看日韩| 欧美潮喷喷水| 亚洲国产精品成人综合色| 亚洲精品,欧美精品| 韩国高清视频一区二区三区| 亚洲真实伦在线观看| av线在线观看网站| 日韩av在线大香蕉| 国产精品无大码| 亚洲人成网站在线观看播放| 一二三四中文在线观看免费高清| 午夜福利高清视频| 日韩欧美在线乱码| 联通29元200g的流量卡| 国语对白做爰xxxⅹ性视频网站| 成人国产麻豆网| 国产精品1区2区在线观看.| 亚洲欧美精品综合久久99| 国产视频首页在线观看| 国产免费福利视频在线观看| 精品午夜福利在线看| 久久亚洲国产成人精品v| 久久久国产成人免费| 婷婷色麻豆天堂久久 | 午夜亚洲福利在线播放| 久久精品国产鲁丝片午夜精品| 亚洲中文字幕一区二区三区有码在线看| 69人妻影院| 麻豆成人午夜福利视频| 国产极品天堂在线| 欧美日本视频| 一本一本综合久久| 亚洲熟妇中文字幕五十中出| 久久久亚洲精品成人影院| 看免费成人av毛片| 直男gayav资源| 春色校园在线视频观看| 变态另类丝袜制服| 亚洲欧洲国产日韩| 国产极品精品免费视频能看的| 可以在线观看毛片的网站| 国产真实乱freesex| 免费大片18禁| 亚洲av成人精品一二三区| 日日啪夜夜撸| 欧美不卡视频在线免费观看| 亚洲天堂国产精品一区在线| 一本一本综合久久| 九九热线精品视视频播放| 成人毛片60女人毛片免费| 午夜免费男女啪啪视频观看| 亚洲av电影在线观看一区二区三区 | 伊人久久精品亚洲午夜| 亚洲av免费在线观看| 少妇熟女aⅴ在线视频| 美女高潮的动态| av国产久精品久网站免费入址| 有码 亚洲区| 久久精品国产自在天天线| 久久久亚洲精品成人影院| 高清视频免费观看一区二区 | a级毛片免费高清观看在线播放| 日韩中字成人| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 乱人视频在线观看| 床上黄色一级片| 欧美日韩精品成人综合77777| 麻豆久久精品国产亚洲av| 黄色配什么色好看| 欧美一级a爱片免费观看看| 亚洲真实伦在线观看| 国产成人午夜福利电影在线观看| 99久国产av精品| 18禁在线无遮挡免费观看视频| 成年女人永久免费观看视频| 1000部很黄的大片| 搞女人的毛片| 久久久精品大字幕| 青春草国产在线视频| 亚洲精品国产av成人精品| 久久久久久国产a免费观看| av在线观看视频网站免费| 亚洲av二区三区四区| 国产综合懂色| 亚洲一区高清亚洲精品| 一区二区三区乱码不卡18| 淫秽高清视频在线观看| 国内精品美女久久久久久| 99热全是精品| 欧美日本亚洲视频在线播放| 免费黄网站久久成人精品| 国产精品一及| 成人美女网站在线观看视频| 日韩欧美精品v在线| 国产精品伦人一区二区| 长腿黑丝高跟| 成人综合一区亚洲| 国产成人91sexporn| 日本三级黄在线观看| 日本wwww免费看| 精品久久久久久成人av| 国内揄拍国产精品人妻在线| 晚上一个人看的免费电影| 精品久久国产蜜桃| 国产av在哪里看| 国内精品一区二区在线观看| 免费人成在线观看视频色| 女人被狂操c到高潮| 久久久国产成人免费| 六月丁香七月| 搡老妇女老女人老熟妇| 免费不卡的大黄色大毛片视频在线观看 | 男人和女人高潮做爰伦理| 久久精品久久精品一区二区三区| av免费在线看不卡| 国产不卡一卡二| 亚洲精品久久久久久婷婷小说 | 日本与韩国留学比较| 九九在线视频观看精品| 国产 一区 欧美 日韩| 欧美高清性xxxxhd video| 热99re8久久精品国产| 嫩草影院新地址| www.色视频.com| 22中文网久久字幕| av天堂中文字幕网| 亚洲欧美成人精品一区二区| 亚洲在线观看片| 国产日韩欧美在线精品| 国产91av在线免费观看| 免费在线观看成人毛片| 人妻少妇偷人精品九色| 2021少妇久久久久久久久久久| 性色avwww在线观看| 一级毛片我不卡| 免费观看性生交大片5| 一个人免费在线观看电影| 一边亲一边摸免费视频| www日本黄色视频网| 婷婷色麻豆天堂久久 | 岛国在线免费视频观看| 久久久精品欧美日韩精品| 久久精品综合一区二区三区| 日本爱情动作片www.在线观看| 欧美激情在线99| 久久精品综合一区二区三区| 亚洲av成人av| 日日撸夜夜添| 国产一区二区在线观看日韩| 又粗又爽又猛毛片免费看| 中文天堂在线官网| 亚洲国产高清在线一区二区三| 18+在线观看网站| 别揉我奶头 嗯啊视频| 欧美高清性xxxxhd video| 中文字幕免费在线视频6| 国产亚洲一区二区精品| 亚洲国产最新在线播放| 一个人免费在线观看电影| 午夜亚洲福利在线播放| 久久综合国产亚洲精品| 日日摸夜夜添夜夜爱| 成人无遮挡网站| 性色avwww在线观看| av在线播放精品| 性插视频无遮挡在线免费观看| 青春草亚洲视频在线观看| 久久久久久久午夜电影| 中文字幕熟女人妻在线| 七月丁香在线播放| 久久久久久伊人网av| 久久久亚洲精品成人影院| 身体一侧抽搐| 午夜福利网站1000一区二区三区| 国产老妇女一区| 欧美日本视频| 夜夜爽夜夜爽视频| 久久久久久久亚洲中文字幕| 如何舔出高潮| 午夜激情福利司机影院| 非洲黑人性xxxx精品又粗又长| 热99re8久久精品国产| 国产真实乱freesex| 99久国产av精品国产电影| 色播亚洲综合网| 一级毛片电影观看 | 国产高清视频在线观看网站| 国产黄a三级三级三级人| 久久精品影院6| 国产精品爽爽va在线观看网站| 超碰av人人做人人爽久久| 免费观看性生交大片5| 国产成人a∨麻豆精品| 久久国内精品自在自线图片| 精品一区二区三区视频在线| 免费在线观看成人毛片| 亚洲国产成人一精品久久久| 国产精品久久电影中文字幕| 国产免费又黄又爽又色| 可以在线观看毛片的网站| 免费av观看视频| 亚洲在久久综合| 国产精品一区二区在线观看99 | 久久久久九九精品影院| 一区二区三区乱码不卡18| 午夜免费激情av| 亚洲电影在线观看av| 国产精品不卡视频一区二区| 岛国毛片在线播放| 你懂的网址亚洲精品在线观看 | 国产欧美另类精品又又久久亚洲欧美| 亚洲av免费在线观看| 久久久久久久久久黄片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久精品久久久久久噜噜老黄 | 最近手机中文字幕大全| 免费观看性生交大片5| 亚洲综合精品二区| 乱系列少妇在线播放| 乱系列少妇在线播放| 久久亚洲国产成人精品v| www日本黄色视频网| 人人妻人人看人人澡| 亚洲av日韩在线播放| 少妇熟女aⅴ在线视频| 不卡视频在线观看欧美| 国产白丝娇喘喷水9色精品| 亚洲最大成人中文| 久久久久久久久久久免费av| 男人舔女人下体高潮全视频| 最近中文字幕高清免费大全6| 99久久无色码亚洲精品果冻| 一级毛片久久久久久久久女| 综合色av麻豆| 亚洲五月天丁香| 春色校园在线视频观看| 女人被狂操c到高潮| av在线老鸭窝| 国产又色又爽无遮挡免| 国产大屁股一区二区在线视频| 成年免费大片在线观看| 亚洲av不卡在线观看| 中文在线观看免费www的网站| 国产精品野战在线观看| 岛国毛片在线播放| 日本黄色视频三级网站网址| 天天躁夜夜躁狠狠久久av| 97人妻精品一区二区三区麻豆| 草草在线视频免费看| 亚洲精品色激情综合| 久久久久久久午夜电影| 69av精品久久久久久| 国语对白做爰xxxⅹ性视频网站| 国产v大片淫在线免费观看| 成人综合一区亚洲| 免费av不卡在线播放| 亚洲乱码一区二区免费版| 又粗又硬又长又爽又黄的视频| 毛片一级片免费看久久久久| 一级黄片播放器| 国产成人精品一,二区| 日韩,欧美,国产一区二区三区 | av免费在线看不卡| 久久婷婷人人爽人人干人人爱| 日韩精品有码人妻一区| 免费观看a级毛片全部| 在线观看av片永久免费下载| 亚洲欧美精品自产自拍| 欧美区成人在线视频| 一个人看视频在线观看www免费| 男女那种视频在线观看| 精品久久久久久久久久久久久| 精品99又大又爽又粗少妇毛片| 欧美3d第一页| 国产视频内射| 最近2019中文字幕mv第一页| av女优亚洲男人天堂| 美女cb高潮喷水在线观看| 一区二区三区乱码不卡18| 亚洲综合精品二区| 日本一二三区视频观看| 午夜福利在线观看吧| 日韩av在线大香蕉| 人妻系列 视频| 亚洲丝袜综合中文字幕| 91久久精品电影网| 婷婷六月久久综合丁香| 91久久精品国产一区二区三区| 男人和女人高潮做爰伦理| 别揉我奶头 嗯啊视频| 国产精品一区二区三区四区免费观看| 国产三级中文精品| 又粗又硬又长又爽又黄的视频| 成人无遮挡网站| 午夜精品国产一区二区电影 | 成人欧美大片| 久久鲁丝午夜福利片| 一区二区三区高清视频在线| 亚洲国产精品专区欧美| 日日摸夜夜添夜夜爱| 蜜桃亚洲精品一区二区三区| 看黄色毛片网站| 一级黄色大片毛片| 免费观看a级毛片全部| 18禁在线播放成人免费| 国内精品宾馆在线| 色视频www国产| 国产欧美日韩精品一区二区| 成人亚洲欧美一区二区av| 神马国产精品三级电影在线观看| 亚洲av免费在线观看| 亚洲va在线va天堂va国产| 欧美成人免费av一区二区三区| 国产淫语在线视频| 亚洲第一区二区三区不卡| 亚洲欧洲国产日韩| 久久6这里有精品| 国产精品一二三区在线看| 日韩一区二区视频免费看| 成人毛片60女人毛片免费| 亚洲av男天堂| 日本av手机在线免费观看| 麻豆乱淫一区二区| 欧美97在线视频| 国语自产精品视频在线第100页| www日本黄色视频网| 亚洲av一区综合| 久久久国产成人精品二区| 日韩在线高清观看一区二区三区| 男的添女的下面高潮视频| 久久人人爽人人爽人人片va| 亚洲高清免费不卡视频| 男人狂女人下面高潮的视频| 搞女人的毛片| 国产精品无大码| 欧美性感艳星| 建设人人有责人人尽责人人享有的 | 久久精品91蜜桃| 国产精品av视频在线免费观看| 乱人视频在线观看| 99热6这里只有精品| 国产精品人妻久久久久久| 国产亚洲5aaaaa淫片| 噜噜噜噜噜久久久久久91| 国产中年淑女户外野战色| 久久精品熟女亚洲av麻豆精品 | 国产69精品久久久久777片| 久久人妻av系列| 亚洲经典国产精华液单| 日韩中字成人| 国产黄片视频在线免费观看| 少妇熟女欧美另类| 国产又色又爽无遮挡免| 久久久午夜欧美精品| 少妇的逼好多水| 欧美色视频一区免费| 天堂网av新在线| 亚洲国产高清在线一区二区三| 草草在线视频免费看| 少妇猛男粗大的猛烈进出视频 | 欧美激情久久久久久爽电影| 中文字幕av成人在线电影| 变态另类丝袜制服| 亚洲国产色片| 欧美97在线视频| 永久网站在线| 亚洲av成人av| 在现免费观看毛片| 联通29元200g的流量卡| ponron亚洲| 欧美性猛交黑人性爽| 久久精品国产亚洲网站| 91精品一卡2卡3卡4卡| 男人舔奶头视频| 成人特级av手机在线观看| 精品一区二区三区视频在线| videos熟女内射| 成年女人看的毛片在线观看| 亚洲国产精品专区欧美| 日本欧美国产在线视频| 国产在线男女| 日本午夜av视频| 一区二区三区免费毛片| 精品久久久噜噜| 国产爱豆传媒在线观看| 亚洲精品国产成人久久av| 51国产日韩欧美| 亚洲人成网站在线观看播放| 亚洲真实伦在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品永久免费网站| 国产免费一级a男人的天堂| 亚洲最大成人av| 成年av动漫网址| 你懂的网址亚洲精品在线观看 | 亚洲人成网站在线播| 亚洲电影在线观看av| 丰满乱子伦码专区| 又爽又黄无遮挡网站| 欧美最新免费一区二区三区| 国产男人的电影天堂91| 一级黄片播放器| 欧美激情国产日韩精品一区| 精品国产三级普通话版| 亚洲欧美清纯卡通| 亚洲欧美成人精品一区二区| 男女下面进入的视频免费午夜| 亚洲综合精品二区| 日韩在线高清观看一区二区三区| 国产一区二区在线观看日韩| 性色avwww在线观看| 久久久久久久亚洲中文字幕| 22中文网久久字幕| 成人午夜高清在线视频| 日本一本二区三区精品| 久久99蜜桃精品久久| 国产高清不卡午夜福利| 人人妻人人澡欧美一区二区| 黑人高潮一二区| 狠狠狠狠99中文字幕| 亚洲国产色片| 亚洲精品一区蜜桃| 特级一级黄色大片| 国内精品美女久久久久久| 亚洲美女视频黄频| 欧美性猛交╳xxx乱大交人| 国产精品嫩草影院av在线观看| 亚洲国产精品久久男人天堂| 高清视频免费观看一区二区 | 免费黄色在线免费观看| 乱码一卡2卡4卡精品| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧美日韩无卡精品| 国内精品宾馆在线| 国产精品不卡视频一区二区| 久久婷婷人人爽人人干人人爱| 午夜爱爱视频在线播放| 欧美+日韩+精品| 中文精品一卡2卡3卡4更新| 亚洲精品456在线播放app| 日本免费在线观看一区| 激情 狠狠 欧美| 级片在线观看| 免费观看a级毛片全部| 99热6这里只有精品| .国产精品久久| 七月丁香在线播放| 国产中年淑女户外野战色| 女人久久www免费人成看片 | 18禁动态无遮挡网站| 欧美三级亚洲精品| 国产黄色小视频在线观看| 黄色欧美视频在线观看| 成人三级黄色视频| 亚洲最大成人手机在线| 欧美日本视频| 亚洲精品乱久久久久久| 97超碰精品成人国产| 少妇被粗大猛烈的视频| 精品一区二区三区视频在线| 国产黄片美女视频| 变态另类丝袜制服| 中文字幕亚洲精品专区| 22中文网久久字幕| 日韩高清综合在线| 秋霞伦理黄片| 久久久成人免费电影| 精品无人区乱码1区二区| 国语自产精品视频在线第100页| 国产av不卡久久| 亚洲欧美精品综合久久99| 欧美精品一区二区大全| 91在线精品国自产拍蜜月| 看片在线看免费视频| 亚洲怡红院男人天堂| 麻豆精品久久久久久蜜桃| 日本免费在线观看一区| 在线观看美女被高潮喷水网站| 最近最新中文字幕大全电影3| 久久99蜜桃精品久久| 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| av免费在线看不卡| 久久人人爽人人片av| 亚洲成av人片在线播放无| 精品一区二区三区人妻视频| 国产精品无大码| 久久亚洲精品不卡| 少妇的逼好多水| 好男人视频免费观看在线| 欧美丝袜亚洲另类| 一级爰片在线观看| 亚洲国产高清在线一区二区三| 午夜亚洲福利在线播放| 尤物成人国产欧美一区二区三区| .国产精品久久| 免费看美女性在线毛片视频| 青春草视频在线免费观看| 99热网站在线观看| 超碰97精品在线观看| 不卡视频在线观看欧美| 久久婷婷人人爽人人干人人爱| 亚洲精品456在线播放app| 精品国产一区二区三区久久久樱花 | 久久久久性生活片| 亚洲av日韩在线播放| 一本久久精品| 夫妻性生交免费视频一级片| 狠狠狠狠99中文字幕| 亚洲无线观看免费| 精品久久国产蜜桃| 黄色欧美视频在线观看| 国产成人免费观看mmmm| 国产精品久久久久久久久免| 亚洲欧美精品专区久久| 亚洲人与动物交配视频| 亚洲av日韩在线播放| 国产精华一区二区三区| 亚洲经典国产精华液单| 欧美激情在线99| 亚洲无线观看免费| 色综合站精品国产| 午夜精品一区二区三区免费看| 国产精品,欧美在线| av线在线观看网站| 最新中文字幕久久久久| 国产精品久久久久久久久免| 91aial.com中文字幕在线观看| 亚洲国产精品成人久久小说| 熟女电影av网| 97热精品久久久久久| 精品酒店卫生间| 岛国在线免费视频观看| 欧美激情国产日韩精品一区| 欧美性感艳星| 国产伦理片在线播放av一区| or卡值多少钱| kizo精华| 我要看日韩黄色一级片| 免费一级毛片在线播放高清视频| 国产高清视频在线观看网站| 成人美女网站在线观看视频| 亚洲国产欧美在线一区| 99热精品在线国产| 亚洲av.av天堂| 免费黄网站久久成人精品| 2021天堂中文幕一二区在线观| 久久久久性生活片| 只有这里有精品99| 日本黄色片子视频| 亚洲欧美精品专区久久| 欧美丝袜亚洲另类| 国产精品国产三级专区第一集| 51国产日韩欧美| 神马国产精品三级电影在线观看| 热99re8久久精品国产| 精品人妻偷拍中文字幕| 99在线视频只有这里精品首页| 天堂中文最新版在线下载 | 老司机影院成人| 成人毛片60女人毛片免费| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 国产老妇伦熟女老妇高清| 国产午夜福利久久久久久| 亚洲精品国产av成人精品| 亚洲成色77777| 好男人视频免费观看在线| 色综合色国产| 日韩av在线大香蕉| 国产精品.久久久| 男女视频在线观看网站免费| 少妇高潮的动态图| 岛国在线免费视频观看| 高清视频免费观看一区二区 |