陸玉玲，謝錢(qián)姣，朱家明，李德政

(安徽財(cái)經(jīng)大學(xué) 統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院，安徽 蚌埠 233030)

基于多元回歸對(duì)糧食產(chǎn)量的研究

陸玉玲，謝錢(qián)姣，朱家明，李德政

(安徽財(cái)經(jīng)大學(xué) 統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院，安徽 蚌埠 233030)

為建立關(guān)于糧食產(chǎn)量的模型，并運(yùn)用模型計(jì)算安徽省2020年化肥、農(nóng)藥使用量零增長(zhǎng)時(shí)化肥利用率提高的百分點(diǎn)，利用時(shí)間序列分析預(yù)測(cè)得到了2015-2020年糧食產(chǎn)量、農(nóng)藥和化肥使用量，建立糧食產(chǎn)量與有效化肥使用量和農(nóng)藥使用量的回歸模型，并利用其解決化肥、農(nóng)藥使用量零增長(zhǎng)時(shí)化肥利用率提高的百分點(diǎn)，運(yùn)用MATLAB進(jìn)行求解，得出化肥利用率提高的百分點(diǎn)為7.3%。

糧食產(chǎn)量；多元回歸；時(shí)間序列；MATLAB

糧食產(chǎn)量是我國(guó)第一產(chǎn)業(yè)的重要組成部分，是國(guó)家經(jīng)濟(jì)發(fā)展的基礎(chǔ)，受到多重因素的影響，比如農(nóng)藥使用量以及化肥使用量。有效的糧食產(chǎn)量預(yù)測(cè)可以為政府決策提供科學(xué)依據(jù)，同時(shí)，在減少農(nóng)藥和化肥使用量下，如何實(shí)現(xiàn)糧食的高產(chǎn)出也是我們一直追求的目標(biāo)。

近年來(lái)，隨著科技的不斷發(fā)展，糧食產(chǎn)量的預(yù)測(cè)方法也得到了不斷的發(fā)展。文獻(xiàn)[1]225-228+267采用了灰色預(yù)測(cè)模型預(yù)測(cè)糧食產(chǎn)量的變化趨勢(shì)，然后運(yùn)用BP神經(jīng)網(wǎng)絡(luò)對(duì)灰色預(yù)測(cè)模型所得到的結(jié)果進(jìn)行了修正，從而提高了預(yù)測(cè)精度。但是，此方法僅能提高單一的糧食產(chǎn)量預(yù)測(cè)精度。文獻(xiàn)[2]6-10利用支持向量機(jī)模型對(duì)糧食產(chǎn)量進(jìn)行預(yù)測(cè)，提高了小樣本糧食產(chǎn)量預(yù)測(cè)精度。但此方法僅考慮糧食歷史產(chǎn)量數(shù)據(jù)，沒(méi)有考慮到其他影響因素。文獻(xiàn)[3]429-433利用灰色預(yù)測(cè)和馬爾可夫預(yù)測(cè)模型相結(jié)合來(lái)對(duì)糧食產(chǎn)量進(jìn)行預(yù)測(cè)，從而大大提高了糧食產(chǎn)量的預(yù)測(cè)精度。

1 模型假設(shè)

為了便于解決和研究問(wèn)題，提出以下幾條假設(shè)：(1)假設(shè)所有數(shù)據(jù)的來(lái)源真實(shí)可靠，具有科學(xué)性；(2)假設(shè)糧食產(chǎn)量只與化肥、農(nóng)藥的使用量及其利用率有關(guān)，不受其他因素的影響。

2 預(yù)測(cè)2015-2020年安徽省的糧食產(chǎn)量、化肥和農(nóng)藥使用量

2.1 研究思路

首先，收集安徽省2005—2014年糧食產(chǎn)量、化肥施用量、農(nóng)藥使用量數(shù)據(jù)，其次根據(jù)所收集到的數(shù)據(jù)建立時(shí)間序列的預(yù)測(cè)模型，并利用MATLAB軟件進(jìn)行求解，從而得到安徽省在2015—2020年糧食產(chǎn)量、化肥施用量、農(nóng)藥使用量的具體數(shù)值。

2.2 研究方法

(1)數(shù)據(jù)收集

①根據(jù)國(guó)家統(tǒng)計(jì)局公布的官方數(shù)據(jù)，將其倒入Excel中得到2005—2014年安徽省的糧食產(chǎn)量、化肥使用量、農(nóng)藥使用量指標(biāo)的具體數(shù)據(jù)見(jiàn)表1。

表1 2005-2014年農(nóng)業(yè)相關(guān)數(shù)據(jù) (單位：萬(wàn)噸)

②根據(jù)安徽省農(nóng)業(yè)委員會(huì)公布的數(shù)據(jù)，我們可以得到2005—2014年安徽省的化肥利用率和農(nóng)藥利用率數(shù)據(jù)，見(jiàn)表2。

表2 2005-2014年化肥、農(nóng)藥利用率

為了反映化肥與農(nóng)藥利用率變動(dòng)的情形及其好壞趨向，根據(jù)表2中的數(shù)據(jù)，做出化肥和農(nóng)藥利用率的雷達(dá)圖，如圖1所示。

圖1 化肥與農(nóng)藥利用率的雷達(dá)圖

(2)數(shù)據(jù)預(yù)處理

有效化肥使用量：指化肥使用量與化肥利用率的乘積。假設(shè)用H表示有效化肥使用量，則H=π1·π2,π1表示化肥使用量,π2表示化肥利用率。

根據(jù)化肥使用量和化肥利用率的數(shù)據(jù)，運(yùn)用EXCEL軟件計(jì)算得到2005—2014年安徽省各年的有效化肥使用量，見(jiàn)表3。

(3)模型的建立與求解

根據(jù)安徽省2005—2014年的糧食產(chǎn)量及化肥、農(nóng)藥的使用量，建立時(shí)間序列預(yù)測(cè)模型，[4]并結(jié)合MATLAB軟件，得到了2005—2014年的糧食產(chǎn)量及化肥、農(nóng)藥使用量的預(yù)測(cè)值，具體數(shù)據(jù)見(jiàn)表4。

表3 2005-2014年有效化肥使用量 (單位：萬(wàn)噸)

表4 糧食產(chǎn)量、化肥和農(nóng)藥使用量預(yù)測(cè)值(單位：萬(wàn)噸)

表5 2005-2014年糧食產(chǎn)量、化肥使用量和農(nóng)藥使用量實(shí)際與預(yù)測(cè)值(單位：萬(wàn)噸)

2.3 結(jié)果分析

一般來(lái)說(shuō)，相對(duì)誤差的絕對(duì)值在0.2以?xún)?nèi)說(shuō)明模型有較高的合理性，且相對(duì)誤差越小，模型的精度越高由表5中的數(shù)據(jù)可知，相對(duì)誤差的絕對(duì)值均在0.2以?xún)?nèi)，這說(shuō)明建立的模型是相對(duì)合理的，可靠性較高。

3 糧食產(chǎn)量模型

3.1 研究思路

首先，根據(jù)表1中2005-2014年間各年的糧食產(chǎn)量、農(nóng)藥使用量和表3中2005-2014年間各年的有效化肥使用量的具體數(shù)據(jù)，以糧食產(chǎn)量為因變量，有效化肥使用量和農(nóng)藥使用量為自變量，建立多元線(xiàn)性回歸模型，[5]利用MATLAB軟件進(jìn)行求解得到回歸方程，并對(duì)模型進(jìn)行檢驗(yàn)。其次，根據(jù)方程，以2015年的化肥使用量為基準(zhǔn)，在二者用量保持不變的情況下，求解2020年化肥利用率所提高的百分點(diǎn)。

3.2 具體過(guò)程和結(jié)果

(1)以2005年為第1年，分別用xt,yt,gt(t=1,2,…16)表示第t年的安徽省糧食產(chǎn)量、農(nóng)藥使用量，有效化肥使用量，它們之間的關(guān)系用公式表為：

xt=a0+a1yt+a2gt+εt

(2)矩陣形式

X=(x1,x2,…x10),

B=(β1,β2,…β10),

N=(ε1,ε2,…ε10)

則糧食產(chǎn)量與各變量的關(guān)系可用矩陣表示為：

(3)未知參數(shù)的估計(jì)值可以由最小二乘來(lái)估計(jì)，表示為：

所以可得回歸方程：

根據(jù)相關(guān)數(shù)據(jù)，并結(jié)合MATLAB軟件得到第t年的糧食產(chǎn)量xt與第t年的農(nóng)藥使用量yt和有效化肥使用量gt之間的二元線(xiàn)性回歸方程為：

xt=1748.6985+14.1856gt-0.1075yt

可決系數(shù)為0.9665。

表6 多元線(xiàn)性回歸預(yù)測(cè)值相對(duì)誤差

由表6可知，相對(duì)誤差的絕對(duì)值均小于0.03，故模型的可靠度較高。

(5)對(duì)殘差進(jìn)行正態(tài)性檢驗(yàn)和t檢驗(yàn)：由MATLAB運(yùn)行結(jié)果可以知道，正態(tài)性檢驗(yàn)的測(cè)試結(jié)果為0，說(shuō)明殘差通過(guò)正態(tài)性檢驗(yàn)，且t檢驗(yàn)的測(cè)試結(jié)果為0，說(shuō)明殘差通過(guò)t檢驗(yàn)。因此模型的建立是比較合理的，利用此模型得到的數(shù)據(jù)和結(jié)論比較真實(shí)可靠。

(6)由表4可知2020年的糧食產(chǎn)量的預(yù)測(cè)值為3864.8萬(wàn)噸，2015年的農(nóng)藥使用量的預(yù)測(cè)值為12.4107萬(wàn)噸，2015年的化肥使用量預(yù)測(cè)值為351.0127萬(wàn)噸，以2015年的糧食產(chǎn)量，農(nóng)藥、化肥使用量為基準(zhǔn)，來(lái)求解2020年的化肥利用率。[6]

由于農(nóng)藥使用量零增長(zhǎng)，將2020糧食產(chǎn)量和2015年農(nóng)藥使用量具體數(shù)值代入方程：

x16=1748.6985+14.1856g16-0.1075y11

求解得到2020年的有效化肥使用量為149.2666萬(wàn)噸。

根據(jù)第t年化肥利用率ft、化肥使用量yt、有效化肥使用量gt之間的關(guān)系，即第t年化肥利用率為：

通過(guò)查閱中國(guó)農(nóng)業(yè)委員會(huì)公布的數(shù)據(jù)，得到2015年化肥利用率為35.2%，所以2020年化肥利用率提高的百分點(diǎn)為42.5%-35.2%=7.3%。

結(jié)束語(yǔ)

在安徽省2020年實(shí)現(xiàn)化肥、農(nóng)藥使用量零增長(zhǎng)的前提下，為了預(yù)測(cè)化肥利用率提高的百分點(diǎn)，本文建立了時(shí)間序列預(yù)測(cè)模型，根據(jù)相關(guān)數(shù)據(jù)預(yù)測(cè)得到了安徽省2015-2020年間各年的糧食產(chǎn)量、農(nóng)藥使用量以及化肥使用量的預(yù)測(cè)值，同時(shí)引入了有效化肥使用量這個(gè)概念，建立了糧食產(chǎn)量關(guān)于農(nóng)藥使用量和有效化肥使用量的二元線(xiàn)性回歸方程，并對(duì)模型進(jìn)行了相對(duì)誤差檢驗(yàn)，對(duì)殘差進(jìn)行了正態(tài)性檢驗(yàn)和t檢驗(yàn)，從而可知模型具有較高的可靠性，但是本文只考慮了農(nóng)藥和化肥使用量對(duì)糧食產(chǎn)量的影響，沒(méi)有考慮其它因素的影響，故結(jié)果具有一定的誤差。此外，本文使用的時(shí)間序列預(yù)測(cè)模型不僅可以預(yù)測(cè)糧食產(chǎn)量，還可以與當(dāng)前先進(jìn)的科學(xué)技術(shù)相聯(lián)系應(yīng)用于測(cè)量領(lǐng)域，比如多時(shí)相遙感圖像數(shù)據(jù)分析，即遙感時(shí)間序列分析，用于監(jiān)測(cè)各種資源動(dòng)態(tài)，通用性較強(qiáng)。

[1]林芳.灰色神經(jīng)網(wǎng)絡(luò)在糧食產(chǎn)量預(yù)測(cè)中的應(yīng)用[J].計(jì)算機(jī)仿真，2012(4):225-228+267.

[2]向昌盛,周子英,武麗娜.糧食產(chǎn)量預(yù)測(cè)的支持向量機(jī)模型研究[J].湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版),2010,(1):6-10.

[3]陳煥珍.基于灰色馬爾科夫模型的青島市糧食產(chǎn)量預(yù)測(cè)[J].計(jì)算機(jī)仿真,2013(5):429-433.

[4]羅鳳曼.時(shí)間序列預(yù)測(cè)模型及其算法研究[D].四川大學(xué),2006.

[5]李瑩.基于時(shí)間序列與多元線(xiàn)性回歸綜合模型的農(nóng)村卷煙銷(xiāo)量預(yù)測(cè)[D].云南大學(xué),2015.

[6]李亮科.生產(chǎn)要素利用對(duì)糧食增產(chǎn)和環(huán)境影響研究[D].中國(guó)農(nóng)業(yè)大學(xué),2015.

Class No.:F323 Document Mark：A

(責(zé)任編輯：蔡雪嵐)

Study on Grain Yield Based on Multiple Regressions

Lu Yuling, Xie Qianjiao, Zhu Jiaming, Li Dezheng

(School of Statistics and Applied Mathematics, Anhui Finance and Economics University, Bengbu, Anhui 233030，China)

In order to establish the food production model and calculate the percentage increase of fertilizer utilization when the use of chemical fertilizers and pesticides achieve zero growth in Anhui Province in 2020, we used the time series analysis approach to predict the grain yield and the use of pesticide and fertilizer in 2015-2020. We established a regression model of grain yield and use of fertilizers and pesticides. Based on the model, we concluded that the utilization ratio of fertilizer is 7.3% solved by MATLAB.

grain yield；multiple regressions；time series；MATLAB

陸玉玲，本科學(xué)生，安徽財(cái)經(jīng)大學(xué)統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院。研究方向:數(shù)學(xué)與應(yīng)用數(shù)學(xué)。 朱家明，碩士，副教授，安徽財(cái)經(jīng)大學(xué)統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院。研究方向：應(yīng)用數(shù)學(xué)與數(shù)學(xué)建模。

國(guó)家自然科學(xué)基金(編號(hào)：11301001)； 國(guó)家級(jí)大學(xué)生創(chuàng)新訓(xùn)練計(jì)劃(編號(hào)：201510378470)。

1672-6758(2017)06-0090-4

F323

A