朱新明，梁玉清

（甘肅省酒泉市農(nóng)業(yè)科學(xué)研究院，甘肅 酒泉 735000）

酒泉是全國(guó)重要的農(nóng)作物種子生產(chǎn)基地，2013年全市制種面積達(dá)2.86萬(wàn)hm2。種子生產(chǎn)面積的大小，直接影響著市場(chǎng)種子的供求，也涉及到地方的產(chǎn)業(yè)規(guī)劃和企業(yè)的銷售策略，因此，建立制種面積的時(shí)間序列分析預(yù)測(cè)模型，對(duì)科學(xué)合理規(guī)劃種子生產(chǎn)面積，制定地方種業(yè)發(fā)展規(guī)劃，避免種子市場(chǎng)的大幅波動(dòng)具有重要的意義。

1 數(shù)據(jù)來(lái)源及分析方法

2000—2013年全市各類農(nóng)作物種子生產(chǎn)面積數(shù)據(jù)，來(lái)源于酒泉市種子管理站上報(bào)政府和上級(jí)主管部門的統(tǒng)計(jì)數(shù)據(jù)。用統(tǒng)計(jì)分析軟件SPSS19.0對(duì)數(shù)據(jù)進(jìn)行匯總，對(duì)制種作物面積變化規(guī)律、組織成分進(jìn)行初步分析，尋找數(shù)據(jù)的變化趨勢(shì)；用時(shí)間序列模塊對(duì)制種面積的變化趨勢(shì)進(jìn)行模型擬合選出最優(yōu)模型，預(yù)測(cè)未來(lái)全市制種業(yè)發(fā)展?fàn)顩r。

2 結(jié)果與分析

2.1 制種作物種類及總面積變化趨勢(shì)分析

酒泉市2000—2013年的制種作物有小麥、玉米、棉花、瓜類、蔬菜、花卉、向日葵、牧草、甜菜九大類。各類作物不同年份制種面積見表1。表1數(shù)據(jù)顯示，酒泉市制種面積總體呈現(xiàn)上升趨勢(shì)，玉米、蔬菜、棉花和瓜類2000—2013年的年平均制種面積分別是0.99萬(wàn)、0.39萬(wàn)、0.26萬(wàn)、0.20萬(wàn)hm2，分別占作物制種面積的47.37%、18.40%、12.21%、9.36%。

表1 2000—2013年酒泉市制種面積 萬(wàn)hm2

圖1 2000—2013年酒泉市制種作物面積分布箱

用統(tǒng)計(jì)分析軟件SPSS19.0對(duì)表1各項(xiàng)數(shù)據(jù)分析，制作箱圖1。圖1表明，在9類作物中，影響制種面積波動(dòng)變化較大的制種作物是玉米、蔬菜、棉花和瓜類。波動(dòng)范圍為制種玉米0.48萬(wàn)～1.46萬(wàn)hm2，制種蔬菜0.14萬(wàn)～0.84萬(wàn)hm2，繁種棉花0.04萬(wàn)～0.46萬(wàn)hm2，制種瓜類0.08萬(wàn)～0.33萬(wàn)hm2。其中，玉米制種面積因子是影響2000—2013年全市制種面積的波動(dòng)的決定因子，其因子成分占到47.37%。

2.2 最優(yōu)模型的選擇

通過時(shí)間序列分析選擇建立數(shù)據(jù)模型，其原理就是利用時(shí)間順序的一組數(shù)字序列，應(yīng)用數(shù)理統(tǒng)計(jì)方法加以處理，以預(yù)測(cè)未來(lái)事物的發(fā)展。時(shí)間序列預(yù)測(cè)一般反映趨勢(shì)變化、周期性變化、隨機(jī)性變化3種實(shí)際變化規(guī)律。時(shí)間序列分析模型有指數(shù)平滑模型和ARIMA模型。我們用統(tǒng)計(jì)分析軟件SPSS19.0時(shí)間序列模塊，對(duì)2000—2013年全市制種面積數(shù)據(jù)進(jìn)行時(shí)間序列分析，建模擬合優(yōu)化。

從殘差序列（圖2）的自相關(guān)（ACF）和偏相關(guān)（PACF）序列圖看出，兩個(gè)圖形都沒有顯著的趨勢(shì)特征（拖尾和截尾現(xiàn)象），故選用指數(shù)平滑模型擬合較為適合。

圖2 殘差相關(guān)函數(shù)序列

從表2看出，平穩(wěn)R和R2均大于0（0.790，0.420），在指數(shù)平滑模型中，Holt線性模型要優(yōu)于其它的均值模型。

2.3 Holt線性趨勢(shì)模型

式中t為當(dāng)前日期；τ為預(yù)測(cè)超前日期數(shù)；Tt、Tt-1為利用t期或前t-1期數(shù)據(jù)，對(duì)t期或t-1期趨勢(shì)的估計(jì)；bt、bt-1為利用t期或前t-1期數(shù)據(jù)，對(duì)趨勢(shì)增加量b的估計(jì)；xt為第t期的實(shí)際觀察值；x^t+τ為利用t期數(shù)據(jù)，對(duì)第t+τ期的預(yù)測(cè)值；α、β為平滑常數(shù)，滿足0＜α、β＜1。

表2 模型擬合

圖3 Holt預(yù)測(cè)結(jié)果

2.4 Holt模型擬合分析預(yù)測(cè)

通過SPSS19.0軟件，用Holt模型計(jì)算2014－2017年酒泉市制種面積的預(yù)測(cè)值及置信區(qū)間，結(jié)果見表3。到2014年全市制種面積預(yù)測(cè)值為2.56萬(wàn)hm2，置信區(qū)間 1.95萬(wàn)～3.18萬(wàn)hm2；2017年為2.74萬(wàn)hm2，置信區(qū)間2.11萬(wàn)～3.37萬(wàn)hm2。

表3 2014—2017年酒泉市制種面積Holt模型預(yù)測(cè)萬(wàn)hm2

圖3為實(shí)際面積序列和Holt模型擬合的變化趨勢(shì)線形圖。從圖中看出，酒泉市制種面積變化波動(dòng)較大，但總體趨勢(shì)呈上升的態(tài)勢(shì)。

3 小結(jié)與分析

1）酒泉市農(nóng)作物制種組成主要是雜交玉米、蔬菜和瓜類三大作物，占全部農(nóng)作物制種面積的75.13%。用SPSS19.0統(tǒng)計(jì)分析時(shí)間序列模塊對(duì)2000—2013年的制種面積數(shù)據(jù)趨勢(shì)擬合，指數(shù)平滑模型Holt模型擬合良好，用該模型對(duì)2014—2017年的制種面積分析預(yù)測(cè)，制種面積呈上升趨勢(shì)。模型預(yù)測(cè)，到2017年，酒泉市制種面積預(yù)計(jì)將達(dá)到2.74萬(wàn)hm2。

2） Holt模型預(yù)測(cè)表明，酒泉市種子生產(chǎn)面積依然呈逐年增長(zhǎng)的趨勢(shì)。2013年，酒泉市肅州區(qū)已被農(nóng)業(yè)部確定為國(guó)家級(jí)玉米種子生產(chǎn)基地，隨著國(guó)家政策對(duì)種業(yè)支持力度的加大，酒泉制種產(chǎn)業(yè)也將由單一的繁種向育繁推銷一體化轉(zhuǎn)變，面積將穩(wěn)定增加。

3） 區(qū)間大波動(dòng)是酒泉制種業(yè)上升趨勢(shì)的次要特征。引起較大波動(dòng)的主要原因是雜交玉米制種和蔬菜制種兩頭在外，即生產(chǎn)訂單在外，銷售市場(chǎng)在外。從全國(guó)看，目前約有種子公司6 200家，公司多，自主知識(shí)產(chǎn)權(quán)品種和品牌少，企業(yè)間同質(zhì)性強(qiáng)，種子市場(chǎng)競(jìng)爭(zhēng)激烈，各企業(yè)安排和從事種子生產(chǎn)任務(wù)隨意性和盲從性大，導(dǎo)致種子市場(chǎng)特別是玉米種子市場(chǎng)周期性波動(dòng)較大。