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

    基于大數(shù)據(jù)的玉米田四代棉鈴蟲發(fā)生量的預(yù)測(cè)模型

    2016-04-07 01:34:53趙雷楊波劉勇牟少敏溫孚江
    大數(shù)據(jù) 2016年1期
    關(guān)鍵詞:發(fā)生量玉米田棉鈴蟲

    趙雷,楊波,劉勇,牟少敏,溫孚江

    山東農(nóng)業(yè)大學(xué)農(nóng)業(yè)大數(shù)據(jù)研究中心,山東 泰安 271018

    基于大數(shù)據(jù)的玉米田四代棉鈴蟲發(fā)生量的預(yù)測(cè)模型

    趙雷,楊波,劉勇,牟少敏,溫孚江

    山東農(nóng)業(yè)大學(xué)農(nóng)業(yè)大數(shù)據(jù)研究中心,山東 泰安 271018

    提出了一種基于支持向量機(jī)的預(yù)測(cè)模型。根據(jù)山東省1999-2013年玉米田第四代棉鈴蟲發(fā)生程度采集的數(shù)據(jù),采用支持向量回歸(SVR)算法,構(gòu)建了玉米田第四代棉鈴蟲發(fā)生程度與其關(guān)聯(lián)因子間的非線性關(guān)系模型,并對(duì)該方法進(jìn)行了測(cè)試與分析。結(jié)果表明,由SVR預(yù)測(cè)模型得到的預(yù)測(cè)發(fā)生量與實(shí)際發(fā)生量基本一致,預(yù)測(cè)的平均絕對(duì)百分比誤差為4.36%,預(yù)測(cè)值與實(shí)際值的相關(guān)系數(shù)為0.960 6,為玉米田第四代棉鈴蟲的有效防控提供了科學(xué)指導(dǎo)。

    農(nóng)業(yè)大數(shù)據(jù);棉鈴蟲;支持向量回歸;監(jiān)測(cè)預(yù)警;玉米

    1 引言

    隨著云計(jì)算、物聯(lián)網(wǎng)等技術(shù)的迅猛發(fā)展,數(shù)據(jù)正以前所未有的速度不斷地增長(zhǎng)和累積,大數(shù)據(jù)時(shí)代已經(jīng)來臨。依靠數(shù)據(jù)分析獲得洞察力,做出更好的業(yè)務(wù)決策是數(shù)據(jù)分析挖掘的主要目的[1,2]。

    預(yù)測(cè)是大數(shù)據(jù)的核心。跳出傳統(tǒng)的因果關(guān)系的思維局限,通過對(duì)大量數(shù)據(jù)的搜集、挖掘和分析,發(fā)現(xiàn)數(shù)據(jù)間隱藏的相關(guān)關(guān)系,從而揭示事物發(fā)生和發(fā)展的內(nèi)在規(guī)律,能做出更快、更符合實(shí)際的預(yù)測(cè)[3,4]。如英國(guó)和美國(guó)科學(xué)家在《Science》發(fā)文指出,全球變暖會(huì)導(dǎo)致非洲和南美洲高原地帶瘧疾的流行[5];基于多年的數(shù)據(jù)挖掘和分析,荷蘭科學(xué)家指明,新煙堿類殺蟲劑吡蟲啉的應(yīng)用是本地食蟲鳥類種群數(shù)量減少的主要原因[6]。因此,借助于大數(shù)據(jù)的研究手段和方法,能夠使已有的農(nóng)業(yè)數(shù)據(jù)“活起來”,認(rèn)識(shí)其內(nèi)在的關(guān)聯(lián)性,預(yù)測(cè)發(fā)展趨勢(shì),使其在指導(dǎo)生產(chǎn)活動(dòng)中產(chǎn)生價(jià)值[7-10]。

    棉鈴蟲屬鱗翅目、夜蛾科,是一種重要的農(nóng)業(yè)害蟲,寄主范圍廣。近年來,隨Bt棉的大范圍種植,第四代棉鈴蟲對(duì)玉米葉片和果穗,特別是果穗的危害逐漸加重,造成玉米產(chǎn)量下降,品質(zhì)降低[11]。截至目前,國(guó)內(nèi)外對(duì)棉鈴蟲的監(jiān)測(cè)預(yù)警研究主要是依靠有限的氣象因素,如溫度、降雨和光周期等[12-14],采用線性回歸分析,建立相關(guān)模型,開展對(duì)棉花田棉鈴蟲的預(yù)測(cè)研究。而對(duì)玉米棉鈴蟲的監(jiān)測(cè)預(yù)警未見報(bào)道。本文基于大數(shù)據(jù)理念,依據(jù)已有的數(shù)據(jù)積累,采用支持向量回歸(SVR)算法,建立了預(yù)測(cè)玉米田第四代棉鈴蟲發(fā)生量的支持向量機(jī)模型,為指導(dǎo)玉米田第四代棉鈴蟲的發(fā)生預(yù)測(cè)及科學(xué)防治打下了基礎(chǔ)。

    2 資料與方法

    2.1 數(shù)據(jù)來源

    本文涉及的變量主要有1999-2013年山東省濱州地區(qū)玉米田第四代棉鈴蟲的發(fā)生量及氣象數(shù)據(jù),主要包括7月中旬到8月上旬的平均氣溫、降水量、最高氣溫(≥35℃)的日數(shù)、降水(≥10 mm)的日數(shù)、極大風(fēng)速、平均本站氣壓、平均風(fēng)速、平均水汽壓、平均相對(duì)濕度、日最低本站氣壓、日最低氣溫、日最高本站氣壓、日最高氣溫、最大風(fēng)速、最大風(fēng)速的風(fēng)向和最小相對(duì)濕度,分別計(jì)算出每年月氣象因子的平均數(shù)。1999-2013年第四代棉鈴蟲的發(fā)生量資料來自山東省濱州地區(qū)植保部門;該時(shí)期的逐日氣象觀測(cè)資料來自國(guó)家氣象信息中心。

    2.2 支持向量機(jī)

    2.2.1 支持向量機(jī)的基本思想

    支持向量機(jī)(support vector machines,SVM)是Vapnik等人根據(jù)統(tǒng)計(jì)學(xué)習(xí)理論提出的機(jī)器學(xué)習(xí)方法[15]?;舅枷胧峭ㄟ^一個(gè)非線性映射把樣本空間映射到一個(gè)高維特征空間中,將尋找最優(yōu)線性回歸超平面的算法歸結(jié)為求解一個(gè)凸約束特性下的凸規(guī)劃問題,并得到全局最優(yōu)解。同時(shí)支持向量機(jī)通過定義核函數(shù)(kernel function),將高維空間中的內(nèi)積運(yùn)算轉(zhuǎn)化為原空間中的核函數(shù)運(yùn)算(如圖1所示)。由于棉鈴蟲的發(fā)生具有非線性、不穩(wěn)定、多變量的特點(diǎn),對(duì)于其蟲害發(fā)生系統(tǒng),很難用確切的公式和解析方法將棉鈴蟲發(fā)生的規(guī)律表達(dá)出來[16]。而這種信息處理方式正是支持向量機(jī)所具備的,因此,將支持向量回歸用于棉鈴蟲發(fā)生量的建模和預(yù)測(cè)分析。

    圖1 支持向量機(jī)結(jié)構(gòu)示意

    2.2.2 支持向量回歸預(yù)測(cè)方法

    支持向量機(jī)的回歸[17]函數(shù)擬合分為線性擬合和非線性擬合,首先考慮線性擬合函數(shù)。假設(shè)有一個(gè)樣本集為:(y1,x1),(y2,x2),…,(yi,xi),y∈R,回歸函數(shù)線性方程表示如下:

    跟基礎(chǔ)定義的線性可分原理一樣,通過函數(shù)的最小值找到最佳的回歸函數(shù),得出:

    其中,w代表W的維數(shù),ξ和ξ*為松弛變量,ξ為上限,ξ*為下限。Vapnik定義不敏感耗損函數(shù)為:

    得到對(duì)偶優(yōu)化方程為:

    于是得到支持向量機(jī)的回歸函數(shù):

    其中,ai、ai*將只有小部分不為0,它們對(duì)應(yīng)的樣本就是支持向量。

    對(duì)于非線性的支持向量機(jī)的回歸,通過一個(gè)非線性映射把數(shù)據(jù)x映射到高維特征空間,然后可以在這個(gè)空間進(jìn)行線性回歸,也就是類似于分類問題。跟支持向量機(jī)定義的線性不可分的原理類似,在它的基礎(chǔ)上進(jìn)行回歸,需用到一個(gè)非敏感性損耗函數(shù),且目標(biāo)函數(shù)為:

    其約束條件是跟式(4)的約束條件一樣。可通過求得的ai和ai*,得到回歸函數(shù):

    2.2.3 核函數(shù)

    核函數(shù)在支持向量回歸中起著重要的作用,它不僅可以解決非線性問題,克服維數(shù)災(zāi)難問題,而且還可以代替高維特征空間中的內(nèi)積運(yùn)算,避免高維度運(yùn)算的復(fù)雜性。

    支持向量機(jī)核函數(shù)的回歸預(yù)測(cè)在諸如疾病預(yù)測(cè)、天氣預(yù)測(cè)、市場(chǎng)預(yù)測(cè)、股價(jià)預(yù)測(cè)等很多方面都得到了廣泛應(yīng)用。核函數(shù)主要有4種。本文將徑向基核函數(shù)(RBF)應(yīng)用于支持向量回歸模型中,支持向量回歸對(duì)玉米田第四代棉鈴蟲的發(fā)生量預(yù)測(cè)從線性到非線性轉(zhuǎn)換,是通過徑向基核函數(shù)。徑向基核函數(shù)能夠針對(duì)棉鈴蟲的發(fā)生具有非線性、不穩(wěn)定性和多變量的特點(diǎn)進(jìn)行處理。

    徑向基核函數(shù)為:

    3 棉鈴蟲發(fā)生量模型實(shí)例研究

    3.1 數(shù)據(jù)的預(yù)處理

    本文是以濱州地區(qū)玉米地1999-2010年玉米田第四代棉鈴蟲的實(shí)際發(fā)生量作為訓(xùn)練樣本(見表1),基于支持向量機(jī)理論,建立支持向量回歸第四代棉鈴蟲發(fā)生量的預(yù)測(cè)模型,以2011-2013年第四代棉鈴蟲發(fā)生量數(shù)據(jù)進(jìn)行預(yù)測(cè)(見表2)。為了提高數(shù)據(jù)之間的可比性和收斂速度、縮短訓(xùn)練時(shí)間,本文在對(duì)數(shù)據(jù)處理時(shí),先對(duì)原始數(shù)據(jù)進(jìn)行歸一化處理。

    表1 支持向量回歸擬合結(jié)果與實(shí)際值對(duì)比

    表2 支持向量回歸預(yù)測(cè)結(jié)果與實(shí)際值對(duì)比

    其中,xi是原始數(shù)據(jù);xi′為歸一化后的數(shù)據(jù);xmin和xmax分別代表數(shù)據(jù)的最小值和最大值。將數(shù)據(jù)壓縮到[0,1]。

    3.2 模型參數(shù)的選取

    本研究中,采用LIBSVM-3.20軟件包來建立SVR模型。通過對(duì)各種核函數(shù)進(jìn)行測(cè)試,最終確定預(yù)測(cè)模型的核函數(shù)為徑向基核函數(shù)。模型的參數(shù)按照相應(yīng)的標(biāo)準(zhǔn)來選取,本文選定的滑動(dòng)窗口寬度(嵌入維數(shù))為16,用gridregression.py自動(dòng)搜索最佳懲罰參數(shù)、靈敏度及徑向?qū)挾鹊群撕瘮?shù)參數(shù)。

    3.3 基于SVR的第四代棉鈴蟲發(fā)生量預(yù)測(cè)模型

    根據(jù)式(10)得到實(shí)際值與擬合值之間的相關(guān)系數(shù)[18],得到的數(shù)據(jù)表明實(shí)際值和擬合值之間有高度的相關(guān)性,相關(guān)系數(shù)接近1,擬合效果好。數(shù)據(jù)表明實(shí)際值與預(yù)測(cè)值之間的相關(guān)系數(shù)為0.96,有高度的相關(guān)性,預(yù)測(cè)結(jié)果與實(shí)際值相符合。

    通過回歸模型得到的訓(xùn)練集樣本的擬合值與實(shí)際值相符合(如圖2所示),而測(cè)試集樣本的預(yù)測(cè)值與實(shí)際值相匹配(如圖3所示)。

    圖2 訓(xùn)練集擬合結(jié)果

    圖3 樣本的預(yù)測(cè)結(jié)果

    其中,yi為樣本的實(shí)際值,為樣本的預(yù)測(cè)值,n為預(yù)測(cè)樣本數(shù)。

    支持向量機(jī)的預(yù)測(cè)準(zhǔn)確率用偏差來表示,主要包括平均絕對(duì)誤差(MAE)、平均相對(duì)誤差(MRE)、均方根誤差(RMSE),表達(dá)式如下:

    其中,yi為樣本的實(shí)際值,為樣本的預(yù)測(cè)值,n為預(yù)測(cè)樣本數(shù)。

    用式(11)~式(13)來計(jì)算偏差,計(jì)算結(jié)果見表3。結(jié)果表明,擬合誤差水平相當(dāng)?shù)停A(yù)測(cè)誤差水平高于擬合誤差??傮w而言,誤差水平總體較低,MAE低于2,MRE低于5%,RMSE低于4。而由表2可知,除了2013年的相對(duì)誤差達(dá)到8.01%之外,其他樣本的相對(duì)誤差均較小。由此可知,基于SVR的玉米田第四代棉鈴蟲發(fā)生量預(yù)測(cè)模型具有可行性。

    表3 預(yù)測(cè)樣本的MAE、MRE、RMSE值

    4 應(yīng)用分析與討論

    玉米田棉鈴蟲的發(fā)生具有復(fù)雜的非線性變化規(guī)律,要建立精確的數(shù)學(xué)模型相當(dāng)困難,因此,根據(jù)大數(shù)據(jù)的理念,將支持向量機(jī)引入第四代棉鈴蟲發(fā)生量的預(yù)測(cè)中。支持向量機(jī)不需要對(duì)數(shù)據(jù)分布性質(zhì)做任何假設(shè),通用性較強(qiáng),實(shí)驗(yàn)結(jié)果表明,將SVR用于玉米田第四代棉鈴蟲發(fā)生量建模與預(yù)測(cè),能較好地揭示玉米田第四代棉鈴蟲的發(fā)生規(guī)律。傳統(tǒng)的模型主要集中在對(duì)棉田棉鈴蟲的發(fā)生預(yù)測(cè)上。傳統(tǒng)的方法主要有多元線性回歸分析[12]、二級(jí)分辨率模型、Fisher模型、期距法等數(shù)十種。近年來,隨著物聯(lián)網(wǎng)、機(jī)器學(xué)習(xí)等技術(shù)的迅速發(fā)展,出現(xiàn)了更為方便、準(zhǔn)確率更高的預(yù)測(cè)方法。如黃健等人[14]利用人工神經(jīng)網(wǎng)絡(luò)(ANN)模型對(duì)新疆地區(qū)1990-2007年第二代棉鈴蟲的發(fā)生等級(jí)進(jìn)行預(yù)測(cè)的結(jié)果表明,BP神經(jīng)網(wǎng)絡(luò)預(yù)報(bào)模型的擬合精度和預(yù)報(bào)精度高于逐步回歸模型。朱軍生等人[19]利用1966-1995年山東惠民縣棉鈴蟲的監(jiān)測(cè)數(shù)據(jù)建立了基于徑向基小波神經(jīng)網(wǎng)絡(luò)的第二代棉鈴蟲卵峰日期預(yù)測(cè)模型,結(jié)果表明,在5年的預(yù)測(cè)數(shù)據(jù)中,4年的預(yù)測(cè)數(shù)據(jù)偏差在3天以內(nèi),另外一年的預(yù)測(cè)數(shù)據(jù)偏差為4天。這對(duì)于卵峰日的預(yù)測(cè)偏差較大。以前對(duì)棉鈴蟲的預(yù)測(cè),大都是對(duì)棉田第二、三代棉鈴蟲發(fā)生期、發(fā)生等級(jí)的預(yù)測(cè),而未見對(duì)玉米田主害代棉鈴蟲進(jìn)行預(yù)測(cè)。雖然基于人工神經(jīng)網(wǎng)絡(luò)的模型比傳統(tǒng)的預(yù)測(cè)模型效果好,但是人工神經(jīng)網(wǎng)絡(luò)容易出現(xiàn)過擬合、維數(shù)災(zāi)難等問題;徑向基小波神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)復(fù)雜。而支持向量機(jī)模型可以很好地解決小樣本、非線性、過擬合、維數(shù)災(zāi)難和局部最優(yōu)的問題,結(jié)構(gòu)簡(jiǎn)單,便于應(yīng)用。因此,首次將SVR應(yīng)用到玉米田第四代棉鈴蟲發(fā)生量的預(yù)測(cè)上。

    從理論上講,SVR算法得到全局最優(yōu),可解決其他神經(jīng)網(wǎng)絡(luò)算法無法避免的局部最優(yōu)問題;SVR預(yù)測(cè)是把線性回歸轉(zhuǎn)為非線性,需要將內(nèi)機(jī)核函數(shù)轉(zhuǎn)換為高維空間的非線性映射,計(jì)算的復(fù)雜性取決于支持向量機(jī)的數(shù)目,不是樣本空間的維數(shù),因此,從一定層面上避免了維數(shù)災(zāi)難的問題。

    但是,SVR算法是半監(jiān)督式學(xué)習(xí)算法,其模型具有一定的局限性,由表1可知,2011-2013年的預(yù)測(cè)值的相對(duì)誤差逐漸增大,與實(shí)際值相差逐漸偏大,說明預(yù)測(cè)的年份距建立模型的年份越遠(yuǎn),預(yù)測(cè)結(jié)果偏差越大。這需要連續(xù)跟蹤現(xiàn)實(shí)數(shù)據(jù)的采集,不斷優(yōu)化預(yù)測(cè)模型。

    氣象因子是影響玉米第四代棉鈴蟲種群發(fā)生和發(fā)展的重要因子。本文通過對(duì)氣象因子進(jìn)行分析,建立基于SVR的玉米第四代棉鈴蟲的預(yù)測(cè)模型,得到較為準(zhǔn)確的結(jié)果。但在2006年,預(yù)測(cè)相對(duì)誤差達(dá)到4.07%,推測(cè)原因主要有以下兩點(diǎn)。

    ● 歷史數(shù)據(jù)的限制。過去十幾年采集的數(shù)據(jù)量小、數(shù)據(jù)范圍比較窄。2014年以后,擴(kuò)展了數(shù)據(jù)采集面,將農(nóng)田中諸如生物因子(包括自然天敵、周邊生物環(huán)境等)、非生物因子(包括土壤性質(zhì)等以及地塊類型、施肥、澆水等)都列于采集范圍之內(nèi)。隨著數(shù)據(jù)采集的不斷完善以及數(shù)據(jù)分析方法的逐步改進(jìn),將會(huì)使預(yù)測(cè)結(jié)果與實(shí)際情況更加接近。

    ● SVR本身的不足,如參數(shù)優(yōu)化等。

    5 結(jié)束語

    本研究首次將SVR應(yīng)用到玉米田第四代棉鈴蟲發(fā)生量的預(yù)測(cè)上,根據(jù)1999-2010年第四代棉鈴蟲采集的數(shù)據(jù)構(gòu)建了玉米田第四代棉鈴蟲發(fā)生量的SVR模型,并對(duì)2011-2013年進(jìn)行了測(cè)試,得到的預(yù)測(cè)發(fā)生量與實(shí)際發(fā)生量基本一致,呈現(xiàn)高度相關(guān)性。這一模型的應(yīng)用,能及時(shí)和準(zhǔn)確地發(fā)布第四代棉鈴蟲監(jiān)測(cè)預(yù)警信息,有效地指導(dǎo)玉米田棉鈴蟲的科學(xué)防控。

    本文首次將SVR用于第四代棉鈴蟲發(fā)生量建模與預(yù)測(cè)中,回歸精度與泛化能力都較高。以上研究表明,SVR應(yīng)用于玉米田棉鈴蟲發(fā)生量預(yù)測(cè)是可行的。

    [1] GANDOMI A, HAIDER M. Beyond the hype: big data concepts, methods, and analytics[J]. International Journal of Information Management, 2015, 35(2): 137-144.

    [2] KAMBATLA K, KOLLIAS G, KUMAR V, et al. Trends in big data analytics[J]. Journal ofParallel and Distributed Computing, 2014, 74(7): 2561-2573.

    [3] LAZER D, KENNEDY R, KING G, et al. The parable of Google flu: traps in big data analysis[J]. Science, 2014, 343(6176): 1203-1205.

    [4] SUNZF, DUKM, ZHENG F X. Perspectives of research and application of big data on smart agriculture[J]. Journal of Agricultural Science and Technology, 2013, 15(6): 63-71.

    [5] SIRAJ A S, SANTOS-VEGA M, BOUMA M J, et al. Altitudinal changes in malaria incidence in highlands of ethiopia and colombia[J]. Science, 2014, 343(6175): 1154-1158.

    [6] HALLM A N NCA, FOPPENRPB, TURNHOUT C AMV, et al. Declines in insectivorous birds are associated with high neonicotinoid concentrations[J]. Nature, 2014, 511(7509): 341-343.

    [7] 宋長(zhǎng)青, 牟少敏, 周虎, 等. 農(nóng)業(yè)大數(shù)據(jù)研究中心的建設(shè)、研究與實(shí)踐[J]. 中國(guó)現(xiàn)代教育裝備, 2015(3): 8-11. SONG C Q, MU S M, ZHOU H, et al. Construction, research and practice of agricultural big data research center[J]. China Modern Educational Equipment, 2015(3): 8-11.

    [8] 辛妍. 大數(shù)據(jù)在農(nóng)業(yè)中的應(yīng)用[J]. 新經(jīng)濟(jì)導(dǎo)刊, 2015(4): 46-51. XIN Y. Bigdataapp lications in agriculture[J]. New Economy Weekly, 2015(4): 46-51.

    [9] 溫孚江. 農(nóng)業(yè)大數(shù)據(jù)與發(fā)展新機(jī)遇[J]. 中國(guó)農(nóng)村科技, 2013(3): 4-7. WEN F J. Agricultural big data and development of new opportunities[J]. Agricultural Sciences in China, 2013(3): 4-7.

    [10] 楊波, 劉勇, 牟少敏, 等. 大數(shù)據(jù)背景下山東省二代玉米螟發(fā)生程度預(yù)測(cè)模型的構(gòu)建[J]. 計(jì)算機(jī)研究與發(fā)展, 2014(S2): 160-165. YA NGB, LIUY, MUSM, etal. Based on big data: the establishment of meteorological forecast model for the occurrence degree of the second generation of corn borer in Shandong[J]. Journal of Computer Research and Development, 2014(S2): 160-165.

    [11] 陳廣泉. 河西走廊玉米田棉鈴蟲發(fā)生規(guī)律與藥劑防治技術(shù)研究[D]. 咸陽:西北農(nóng)林科技大學(xué), 2004: 12-16. CHEN G Q. The occurrence regularity and cottrol technology research of cotton bollworm in the hexi corridor cornfield[D]. Xianyang: Northwest A&F University, 2004: 12-16.

    [12] CHAN GNOND, SAND STRO MM, ASTOLFI J, et al. Using climatology to predict the first majorsummer corn earworm (lepidoptera: noctuidae) catch in north central Illinois[J]. Meteorological Applications, 2010, 17(3): 321-328.

    [13] FENG H, GOULD F, HUANG Y, et al. Modeling the population dynamics of cotton bollworm Helicoverpa armigera (Hübner) (lepidoptera: noctuidae) over a wide area in northern China[J]. Ecological Modelling, 2010, 221(15): 1819-1830.

    [14] RAJ K R, KARDAM A, ARORA J K, et al. Artificial neural network (ANN) design for Hg-Se interactions and their effect on reduction of Hg uptake by radish plant[J]. Journal of Radioanalytical and Nuclear Chemistry, 2010, 283(3): 797-801.

    [15] VAPNIK V N. The nature of statistical learning theory[J]. IEEE Transactions on Neural Networks, 1995, 10(5): 988-999.

    [16] 趙仲華, 沈佐銳. 昆蟲種群動(dòng)態(tài)非線性建模理論與應(yīng)用[J]. 生物數(shù)學(xué)學(xué)報(bào), 2001, 16(4):439-447. ZHAO Z H, SHEN Z R. Insect population dynamics of nonlinear modeling theory and application[J]. Journal of Biomathematics, 2001, 16(4): 439-447.

    [17] 李永娜. 基于支持向量機(jī)的回歸預(yù)測(cè)綜述[J].信息通信, 2014(11):32-33. LI Y N. Regression forecast review based on support vector machine[J]. Message Communication, 2014(11): 32-33.

    [18] ERDAL H I. Two-level and hybrid ensembles of decision trees for high performance concrete compressive strength prediction[J]. Engineering Applications of Artificial Intelligence, 2013, 26(7): 1689-1697.

    [19] 朱軍生, 翟保平, 董保信. 基于徑向基小波網(wǎng)絡(luò)的二代棉鈴蟲卵峰日預(yù)測(cè)模型[J]. 昆蟲學(xué)報(bào),2010, 53(12): 1429-1435. ZHU J S, ZHAI B P, DONG B X. Forecasting model for the oviposition peak day in thesecond generation of helicoverpa armigera (lepidopeter: noctuidae) based on radial basis wavelet network[J]. Acta Entomologica Sinica, 2010, 53(12): 1429-1435.

    Forecasting model for the fourth generation of cotton bollworm in corn fi elds based on big data

    ZHAO Lei, YANG Bo, LIU Yong, MU Shaomin, WEN Fujiang
    Agricultural Big Data Research Center, Shandong Agricultural University, Taian 271018, China

    The monitoring and forecasting model was put forward based on support vector machine program. According to the data collection of the fourth generation occurrence degree of the corn bollworm in Shandong province from 1999 to 2013, the support vector regression (SVR) method was adopted to build the nonlinear correlation model between the occurrence degree of the fourth generation bollworm and the associated factors. The method and the model were tested and analyzed. The results showed that the SVR forecasting model for prediction was almost in accord with the actual insect occurrence situation. The mean absolute percentage error was 4.36%, and the actual and estimated value of the correlation coefficient was 0.960 6. It could provide effective and accurate guidance to the cotton bollworm control in corn fields.

    agricultural big data, cotton bollworm, support vector regression, monitoring and forecasting, corn

    S431.9

    A

    10.11959/j.issn.2096-0271.2016008

    趙雷(1990-),女,山東農(nóng)業(yè)大學(xué)碩士生,主要研究方向?yàn)檗r(nóng)業(yè)科研與大數(shù)據(jù)。

    楊波(1988-),女,山東農(nóng)業(yè)大學(xué)碩士生,主要研究方向?yàn)檗r(nóng)業(yè)大數(shù)據(jù)。

    劉勇(1968-),男,山東農(nóng)業(yè)大學(xué)教授、博士生導(dǎo)師,主要研究方向?yàn)楹οx綠色防控和農(nóng)業(yè)大數(shù)據(jù)。

    牟少敏(1964-),男,博士,山東農(nóng)業(yè)大學(xué)教授,主要研究方向?yàn)榇髷?shù)據(jù)、機(jī)器學(xué)習(xí)和模式識(shí)別。

    溫孚江(1955-),男,現(xiàn)任山東農(nóng)業(yè)大學(xué)校長(zhǎng)、教授,農(nóng)業(yè)大數(shù)據(jù)創(chuàng)新戰(zhàn)略聯(lián)盟理事長(zhǎng),全國(guó)人民代表大會(huì)常務(wù)委員會(huì)委員。早年留學(xué)美國(guó),并獲博士學(xué)位。主要從事植物保護(hù)研究和宏觀農(nóng)業(yè)研究工作。發(fā)表論文210余篇,專著5部。最近一部專著《大數(shù)據(jù)農(nóng)業(yè)》由中國(guó)農(nóng)業(yè)出版社于2015年9月出版。目前主要從事農(nóng)業(yè)大數(shù)據(jù)應(yīng)用研究工作,是我國(guó)農(nóng)業(yè)大數(shù)據(jù)研究主要發(fā)起人之一。

    2015-10-28

    劉勇,liuyong@sdau.edu.cn;溫孚江,fjw@sdau.edu.cn

    山東省農(nóng)業(yè)重大應(yīng)用技術(shù)創(chuàng)新課題基金資助項(xiàng)目

    Foundation Item: Major Innovation of Applied Technology in Agriculture of Shandong Province

    猜你喜歡
    發(fā)生量玉米田棉鈴蟲
    基于眼高刻槽數(shù)預(yù)測(cè)松褐天牛發(fā)生量分析
    六盤山地區(qū)落葉松紅腹葉蜂發(fā)生量對(duì)溫度與降水的響應(yīng)研究
    不同誘芯、誘捕器誘集棉鈴蟲效果試驗(yàn)
    高效低毒新型藥劑防治棉鈴蟲試驗(yàn)
    棉鈴蟲幾種常見藥劑大田試驗(yàn)篩選
    玉米田除草劑的那些事
    高溫干旱下的夏播玉米田害蟲防治?
    1%苦皮藤素防治棉鈴蟲試驗(yàn)
    玉米田不同控釋肥料效用研討
    玉米田化學(xué)除草技術(shù)
    免费av不卡在线播放| 国产乱人偷精品视频| 曰老女人黄片| 精品国产国语对白av| 久久人人爽av亚洲精品天堂| 日韩大片免费观看网站| 免费人妻精品一区二区三区视频| 日产精品乱码卡一卡2卡三| 亚洲精品,欧美精品| 免费高清在线观看视频在线观看| 亚洲精品乱码久久久久久按摩| 人成视频在线观看免费观看| 国产有黄有色有爽视频| 在线天堂最新版资源| 伦理电影大哥的女人| 国产白丝娇喘喷水9色精品| 国产成人一区二区在线| 国产又色又爽无遮挡免| 国产淫语在线视频| 考比视频在线观看| 乱码一卡2卡4卡精品| 新久久久久国产一级毛片| 日韩制服丝袜自拍偷拍| 中文字幕亚洲精品专区| 国产精品蜜桃在线观看| 建设人人有责人人尽责人人享有的| 在线观看www视频免费| 青春草国产在线视频| 国产亚洲一区二区精品| 久久久久久人妻| 韩国精品一区二区三区 | 天天躁夜夜躁狠狠躁躁| 内地一区二区视频在线| 国产亚洲av片在线观看秒播厂| tube8黄色片| 国产高清三级在线| 男女午夜视频在线观看 | 草草在线视频免费看| 一区二区av电影网| 亚洲av成人精品一二三区| 免费不卡的大黄色大毛片视频在线观看| 在线观看三级黄色| 黑人猛操日本美女一级片| 欧美精品av麻豆av| 人人澡人人妻人| 精品少妇久久久久久888优播| av国产精品久久久久影院| 婷婷色综合www| 久久精品夜色国产| 丰满饥渴人妻一区二区三| 欧美 日韩 精品 国产| 亚洲国产成人一精品久久久| 成年av动漫网址| 国产免费又黄又爽又色| 欧美激情 高清一区二区三区| 亚洲国产日韩一区二区| 伦理电影大哥的女人| 亚洲国产精品一区二区三区在线| 亚洲av欧美aⅴ国产| 亚洲精品自拍成人| 亚洲一区二区三区欧美精品| 欧美日韩精品成人综合77777| 日韩精品有码人妻一区| av福利片在线| 国产熟女午夜一区二区三区| 成人国产av品久久久| 国产高清国产精品国产三级| 黄色毛片三级朝国网站| 日日摸夜夜添夜夜爱| 亚洲欧美一区二区三区国产| 国产精品久久久久久精品古装| 亚洲少妇的诱惑av| 国产xxxxx性猛交| 人人妻人人爽人人添夜夜欢视频| 成人亚洲欧美一区二区av| 日本黄大片高清| 成人黄色视频免费在线看| xxx大片免费视频| 青青草视频在线视频观看| 99久久综合免费| 成年人午夜在线观看视频| 丝袜在线中文字幕| 又黄又爽又刺激的免费视频.| 宅男免费午夜| 成人免费观看视频高清| 日韩人妻精品一区2区三区| 免费在线观看黄色视频的| 久久久国产一区二区| 国产精品女同一区二区软件| 久久韩国三级中文字幕| 久久精品夜色国产| 十八禁高潮呻吟视频| 亚洲,欧美精品.| 欧美日韩一区二区视频在线观看视频在线| 国产极品天堂在线| av免费在线看不卡| 亚洲av电影在线观看一区二区三区| 精品卡一卡二卡四卡免费| 美女视频免费永久观看网站| 精品国产一区二区三区四区第35| 18禁国产床啪视频网站| 中文字幕另类日韩欧美亚洲嫩草| 视频在线观看一区二区三区| 精品久久久久久电影网| av不卡在线播放| 超碰97精品在线观看| 水蜜桃什么品种好| 纯流量卡能插随身wifi吗| 午夜免费观看性视频| 人妻一区二区av| 免费观看av网站的网址| 91aial.com中文字幕在线观看| 黑人高潮一二区| av在线播放精品| 国产成人aa在线观看| 亚洲国产精品国产精品| 色5月婷婷丁香| 欧美精品一区二区大全| 免费在线观看黄色视频的| 亚洲内射少妇av| 大香蕉97超碰在线| 一二三四在线观看免费中文在 | 乱人伦中国视频| 在线观看一区二区三区激情| 少妇的逼水好多| 看免费成人av毛片| 日韩人妻精品一区2区三区| 一本色道久久久久久精品综合| 亚洲国产精品一区三区| 中文字幕最新亚洲高清| 在线观看免费高清a一片| 亚洲精品第二区| 国产在线一区二区三区精| 久久久精品免费免费高清| 波野结衣二区三区在线| 人成视频在线观看免费观看| 最黄视频免费看| 大码成人一级视频| 国产成人av激情在线播放| 一级毛片我不卡| 三级国产精品片| 免费黄频网站在线观看国产| 免费看不卡的av| 久久国产精品大桥未久av| 日本色播在线视频| 精品久久国产蜜桃| 亚洲成人手机| 在线观看www视频免费| 日日啪夜夜爽| 国产成人一区二区在线| 男女免费视频国产| av福利片在线| 国产片内射在线| 欧美国产精品va在线观看不卡| 啦啦啦啦在线视频资源| 丁香六月天网| 亚洲精品美女久久久久99蜜臀 | 亚洲国产精品一区二区三区在线| av福利片在线| 在线观看国产h片| 国产免费一区二区三区四区乱码| 亚洲精品成人av观看孕妇| 国产日韩欧美亚洲二区| 91精品三级在线观看| 伊人亚洲综合成人网| 成人手机av| 中文精品一卡2卡3卡4更新| 久久韩国三级中文字幕| 麻豆精品久久久久久蜜桃| 国产色爽女视频免费观看| 亚洲在久久综合| 午夜福利网站1000一区二区三区| 不卡视频在线观看欧美| 大香蕉久久网| 亚洲性久久影院| 免费av不卡在线播放| 只有这里有精品99| 18禁裸乳无遮挡动漫免费视频| 日日爽夜夜爽网站| 麻豆精品久久久久久蜜桃| 国产色婷婷99| 性高湖久久久久久久久免费观看| 精品福利永久在线观看| 中文欧美无线码| av又黄又爽大尺度在线免费看| 国产免费视频播放在线视频| 国产乱来视频区| 国产日韩欧美在线精品| 看免费成人av毛片| 亚洲精品自拍成人| av又黄又爽大尺度在线免费看| 午夜老司机福利剧场| 中文字幕制服av| 另类亚洲欧美激情| 高清黄色对白视频在线免费看| 少妇人妻久久综合中文| 久久人人爽av亚洲精品天堂| 51国产日韩欧美| 午夜影院在线不卡| 成人亚洲精品一区在线观看| 亚洲av中文av极速乱| 久久久久视频综合| 高清不卡的av网站| 狂野欧美激情性bbbbbb| 国产日韩欧美亚洲二区| 黄色一级大片看看| 日韩av在线免费看完整版不卡| 亚洲伊人久久精品综合| 人人澡人人妻人| 国产精品一二三区在线看| 亚洲精品久久成人aⅴ小说| 人人妻人人澡人人爽人人夜夜| 免费人妻精品一区二区三区视频| 在线观看三级黄色| videossex国产| 中文字幕免费在线视频6| av视频免费观看在线观看| 日本免费在线观看一区| 精品一区二区三区四区五区乱码 | 一本—道久久a久久精品蜜桃钙片| 国产1区2区3区精品| 中文字幕亚洲精品专区| 色5月婷婷丁香| 国产精品 国内视频| 欧美精品国产亚洲| 亚洲四区av| 午夜影院在线不卡| 国产精品久久久久久久电影| 国产欧美日韩综合在线一区二区| xxxhd国产人妻xxx| 天堂俺去俺来也www色官网| 国产精品人妻久久久久久| 亚洲精品,欧美精品| av天堂久久9| a级毛色黄片| 欧美激情 高清一区二区三区| 日韩av不卡免费在线播放| 久久人人爽人人爽人人片va| 久久久久精品性色| 欧美97在线视频| 亚洲精品aⅴ在线观看| 国产精品人妻久久久影院| videossex国产| 高清毛片免费看| 女人久久www免费人成看片| 免费播放大片免费观看视频在线观看| 亚洲精品久久午夜乱码| 麻豆精品久久久久久蜜桃| 国产又色又爽无遮挡免| 视频在线观看一区二区三区| 国产精品三级大全| 国产精品久久久久久精品电影小说| 欧美97在线视频| 免费观看av网站的网址| 国产男人的电影天堂91| 免费日韩欧美在线观看| 久久久a久久爽久久v久久| 狠狠婷婷综合久久久久久88av| av片东京热男人的天堂| 一二三四中文在线观看免费高清| 香蕉精品网在线| 捣出白浆h1v1| 男人操女人黄网站| 黄色怎么调成土黄色| 亚洲国产欧美日韩在线播放| 国产探花极品一区二区| 免费观看a级毛片全部| 国精品久久久久久国模美| 久久久久久人妻| 久久女婷五月综合色啪小说| 免费播放大片免费观看视频在线观看| 最近中文字幕高清免费大全6| 99热这里只有是精品在线观看| 日韩欧美一区视频在线观看| 男女免费视频国产| 啦啦啦在线观看免费高清www| 老司机影院成人| 成年女人在线观看亚洲视频| 国产成人精品无人区| 啦啦啦在线观看免费高清www| 美女国产高潮福利片在线看| 亚洲精品久久成人aⅴ小说| 一级,二级,三级黄色视频| 欧美另类一区| 免费日韩欧美在线观看| 午夜免费男女啪啪视频观看| freevideosex欧美| 免费大片18禁| 国产爽快片一区二区三区| 国产av精品麻豆| 中文精品一卡2卡3卡4更新| 曰老女人黄片| 男人添女人高潮全过程视频| 永久免费av网站大全| 亚洲欧美日韩另类电影网站| 成年动漫av网址| 男女边摸边吃奶| 精品少妇久久久久久888优播| 国产精品蜜桃在线观看| 成年人免费黄色播放视频| 精品一区二区免费观看| 久久精品国产综合久久久 | av免费观看日本| 少妇猛男粗大的猛烈进出视频| 自线自在国产av| 亚洲欧美精品自产自拍| 亚洲av福利一区| 久久99热6这里只有精品| 久久鲁丝午夜福利片| 亚洲综合色惰| 2021少妇久久久久久久久久久| 青春草视频在线免费观看| 丝袜人妻中文字幕| 色视频在线一区二区三区| 天天躁夜夜躁狠狠躁躁| 免费高清在线观看日韩| 精品久久蜜臀av无| 国产 精品1| 精品视频人人做人人爽| 91国产中文字幕| 天美传媒精品一区二区| 中文字幕亚洲精品专区| 精品国产露脸久久av麻豆| 国产精品人妻久久久影院| 人妻系列 视频| 性高湖久久久久久久久免费观看| 各种免费的搞黄视频| 狂野欧美激情性bbbbbb| 久久久久久久久久成人| 免费播放大片免费观看视频在线观看| 日韩免费高清中文字幕av| 精品一区二区三卡| 欧美亚洲 丝袜 人妻 在线| 看免费av毛片| 日本91视频免费播放| 母亲3免费完整高清在线观看 | 亚洲国产av影院在线观看| 在线 av 中文字幕| 亚洲精品久久成人aⅴ小说| 日本与韩国留学比较| av卡一久久| 午夜影院在线不卡| 欧美日韩成人在线一区二区| 男男h啪啪无遮挡| 亚洲少妇的诱惑av| av网站免费在线观看视频| 欧美性感艳星| 中文字幕制服av| 国产免费福利视频在线观看| 欧美另类一区| 成人二区视频| 中文字幕另类日韩欧美亚洲嫩草| 成人18禁高潮啪啪吃奶动态图| 熟妇人妻不卡中文字幕| 午夜福利在线观看免费完整高清在| 亚洲欧洲日产国产| 国产69精品久久久久777片| 国产av精品麻豆| 97在线视频观看| 欧美 日韩 精品 国产| 国产av码专区亚洲av| 国国产精品蜜臀av免费| 飞空精品影院首页| 欧美xxⅹ黑人| 国产精品女同一区二区软件| 久久精品久久久久久久性| 99视频精品全部免费 在线| 夫妻午夜视频| 国产爽快片一区二区三区| 纵有疾风起免费观看全集完整版| 男女国产视频网站| 欧美亚洲日本最大视频资源| 国产成人精品在线电影| 伦理电影免费视频| 狠狠精品人妻久久久久久综合| 高清黄色对白视频在线免费看| 综合色丁香网| 午夜福利视频在线观看免费| 亚洲av国产av综合av卡| 亚洲色图综合在线观看| 久久这里有精品视频免费| 国产一区二区三区av在线| 高清不卡的av网站| 国产色婷婷99| 观看av在线不卡| 国产精品嫩草影院av在线观看| 久久av网站| 成人国语在线视频| 在线免费观看不下载黄p国产| 久久鲁丝午夜福利片| 九色成人免费人妻av| 国产欧美日韩一区二区三区在线| 你懂的网址亚洲精品在线观看| 丰满迷人的少妇在线观看| 老女人水多毛片| 在线观看免费日韩欧美大片| 亚洲婷婷狠狠爱综合网| 久久毛片免费看一区二区三区| 日本wwww免费看| 国产日韩欧美视频二区| 人体艺术视频欧美日本| 精品第一国产精品| 极品少妇高潮喷水抽搐| 欧美亚洲 丝袜 人妻 在线| 黄色毛片三级朝国网站| 成年女人在线观看亚洲视频| 免费在线观看黄色视频的| 午夜激情av网站| 国产免费现黄频在线看| 91久久精品国产一区二区三区| 九九爱精品视频在线观看| 亚洲国产精品一区二区三区在线| 最近中文字幕高清免费大全6| 久热久热在线精品观看| 日本黄大片高清| 久久97久久精品| 精品一区二区三卡| av免费观看日本| 免费黄色在线免费观看| 久久人妻熟女aⅴ| 一二三四中文在线观看免费高清| 久久精品aⅴ一区二区三区四区 | 99精国产麻豆久久婷婷| 成人影院久久| videos熟女内射| 日本av免费视频播放| 男人爽女人下面视频在线观看| 精品人妻熟女毛片av久久网站| 大片电影免费在线观看免费| 男人舔女人的私密视频| 欧美日韩成人在线一区二区| 大香蕉97超碰在线| 欧美激情极品国产一区二区三区 | 免费人成在线观看视频色| 免费播放大片免费观看视频在线观看| 一二三四在线观看免费中文在 | 色婷婷av一区二区三区视频| 在线免费观看不下载黄p国产| 国产av码专区亚洲av| 久久午夜综合久久蜜桃| 一个人免费看片子| 狠狠婷婷综合久久久久久88av| 久久影院123| 成年av动漫网址| 国产男女超爽视频在线观看| 国产熟女午夜一区二区三区| 国产精品一区二区在线不卡| 国产精品久久久久久av不卡| 春色校园在线视频观看| 久久国内精品自在自线图片| 国内精品宾馆在线| 黄色 视频免费看| 亚洲精品美女久久av网站| 久久这里只有精品19| 热99国产精品久久久久久7| 国产精品.久久久| 美女xxoo啪啪120秒动态图| 黄片无遮挡物在线观看| 精品少妇久久久久久888优播| 人妻 亚洲 视频| 色婷婷av一区二区三区视频| 99久久精品国产国产毛片| 色94色欧美一区二区| 免费观看性生交大片5| 国产成人午夜福利电影在线观看| 大片电影免费在线观看免费| 国产高清三级在线| 日本wwww免费看| 18禁观看日本| 一区二区三区精品91| 美女国产高潮福利片在线看| 国产av精品麻豆| 午夜福利视频精品| 欧美日本中文国产一区发布| 黄片播放在线免费| 在线观看免费视频网站a站| 一本色道久久久久久精品综合| 亚洲成人一二三区av| 欧美日韩成人在线一区二区| 国产精品久久久久久精品电影小说| 免费大片黄手机在线观看| 久久影院123| 国产免费又黄又爽又色| 我的女老师完整版在线观看| 韩国精品一区二区三区 | 性色av一级| 亚洲激情五月婷婷啪啪| 日本欧美国产在线视频| 波野结衣二区三区在线| 亚洲精品av麻豆狂野| 亚洲成av片中文字幕在线观看 | 亚洲精品久久久久久婷婷小说| 在现免费观看毛片| 久久久久国产网址| 亚洲国产精品一区三区| 人妻 亚洲 视频| 2018国产大陆天天弄谢| 男女午夜视频在线观看 | 男人舔女人的私密视频| 一区在线观看完整版| 99热全是精品| 男女啪啪激烈高潮av片| a级毛片黄视频| 中文字幕人妻丝袜制服| 久久精品国产综合久久久 | 亚洲,一卡二卡三卡| 成年人免费黄色播放视频| 精品少妇久久久久久888优播| 国产不卡av网站在线观看| 欧美少妇被猛烈插入视频| 欧美97在线视频| 国产男人的电影天堂91| 日本91视频免费播放| 久久久久久伊人网av| 熟妇人妻不卡中文字幕| 看非洲黑人一级黄片| 女人被躁到高潮嗷嗷叫费观| 亚洲熟女精品中文字幕| 欧美成人精品欧美一级黄| 只有这里有精品99| 亚洲国产色片| 欧美精品一区二区大全| 亚洲av免费高清在线观看| 亚洲第一区二区三区不卡| 最后的刺客免费高清国语| www.av在线官网国产| 高清av免费在线| 国产淫语在线视频| 亚洲成国产人片在线观看| 午夜福利视频在线观看免费| 国产成人免费观看mmmm| 丁香六月天网| 免费在线观看黄色视频的| 韩国高清视频一区二区三区| 亚洲一区二区三区欧美精品| 欧美xxxx性猛交bbbb| 亚洲精品自拍成人| 欧美xxxx性猛交bbbb| 伊人亚洲综合成人网| 午夜激情av网站| 母亲3免费完整高清在线观看 | 青青草视频在线视频观看| 午夜福利,免费看| 校园人妻丝袜中文字幕| 欧美+日韩+精品| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放| 视频中文字幕在线观看| 99re6热这里在线精品视频| 亚洲av电影在线观看一区二区三区| 午夜免费男女啪啪视频观看| 亚洲精品国产色婷婷电影| 亚洲精品久久成人aⅴ小说| 国产精品三级大全| 久久久精品免费免费高清| 精品久久蜜臀av无| 男人添女人高潮全过程视频| 97在线视频观看| 欧美老熟妇乱子伦牲交| 中国国产av一级| 深夜精品福利| 免费观看无遮挡的男女| 人人妻人人澡人人爽人人夜夜| 夜夜骑夜夜射夜夜干| 久久精品人人爽人人爽视色| 亚洲 欧美一区二区三区| 五月开心婷婷网| 9191精品国产免费久久| 欧美日韩综合久久久久久| 卡戴珊不雅视频在线播放| 久久女婷五月综合色啪小说| 国语对白做爰xxxⅹ性视频网站| 国产熟女欧美一区二区| 午夜福利影视在线免费观看| 高清黄色对白视频在线免费看| 国产精品一区www在线观看| 亚洲av国产av综合av卡| 在线观看www视频免费| 少妇被粗大猛烈的视频| 蜜桃国产av成人99| 最黄视频免费看| 伊人久久国产一区二区| 色吧在线观看| 成年女人在线观看亚洲视频| 欧美精品av麻豆av| 久久久久久久久久成人| 啦啦啦视频在线资源免费观看| 成人免费观看视频高清| 一二三四中文在线观看免费高清| 国产av精品麻豆| 中文乱码字字幕精品一区二区三区| 亚洲精品乱码久久久久久按摩| 久久精品久久久久久久性| 日韩中字成人| 丁香六月天网| 少妇的逼好多水| 亚洲欧美成人精品一区二区| 一区二区三区四区激情视频| 黑丝袜美女国产一区| 2021少妇久久久久久久久久久| 国产精品 国内视频| 热re99久久精品国产66热6| 免费在线观看黄色视频的| 国产日韩欧美亚洲二区| 五月玫瑰六月丁香| 黄网站色视频无遮挡免费观看| 亚洲熟女精品中文字幕| 国产乱人偷精品视频| xxx大片免费视频| 国产精品人妻久久久久久| 日本av手机在线免费观看| 国产高清不卡午夜福利| 免费日韩欧美在线观看| 久久久亚洲精品成人影院| 黄色 视频免费看|