（華中師范大學經(jīng)濟與工商管理學院，湖北武漢 430070）

糧食是事關國計民生的戰(zhàn)略物質，糧食問題不僅是治國安邦的大事，而且事關國民經(jīng)濟和社會的發(fā)展。一個國家只有牢牢把握糧食生產的主動權，立足糧食自給自足，才能掌控經(jīng)濟和社會發(fā)展的大局。首先，糧食是人類生存和發(fā)展的最基本物質，如果糧食出了問題，不僅經(jīng)濟建設等自亂陣腳，還會引起社會動蕩，影響政府的執(zhí)政能力。其次，糧食是工業(yè)的重要原料，如果糧食短缺，許多工業(yè)部門就會因為原料不足被迫減產停產，并波及到相關行業(yè)，進而影響到整個國民經(jīng)濟的健康發(fā)展。第三，糧食是農業(yè)發(fā)展、農民致富的依托?？梢娂Z食安全問題不僅僅只是一個經(jīng)濟問題，更是一個社會問題、政治問題。影響糧食產量的因素主要有以下幾個方面：（1）糧食播種面積。糧食播種面積直接和糧食產量掛鉤。種植面積有了保障，糧食就不會出現(xiàn)大問題，即使發(fā)生一些自然災害，也不會導致糧食產量大幅下跌；相反，如果沒有種植面積為保證，即使管理再好，糧食產量也難以有保障。（2）種植技術。科學種田，提高耕地的綜合生產能力，是糧食豐收的重要保障。大力推廣良種和良法，提高糧食單產和品質，是未來轉變糧食生產方式的方向。（3）農田基礎設施。完善的農田水利設施是糧食高產穩(wěn)產的基礎。（4）氣候條件。在當前我國農業(yè)底子薄弱的情況下，糧食生產易受氣候影響。風調雨順，氣候正常，糧食有望獲得豐收；如果出現(xiàn)嚴重的旱澇、病蟲害、極端天氣，糧食就有可能減產。（5）種植結構。不同的糧食品種，單產存在差異，種植結構的調整會影響糧食的產量。雖然影響糧食產量的因素比較多，大部分因素帶有不可知性、不確定性，但只要不掉以輕心，鼓勵糧食生產的政策、資金、技術，措施到位，糧食產量就不會有太大波動。因此，可以從我國糧食本身的運行規(guī)律出發(fā)，把握其變化發(fā)展的趨勢。

重慶是我國的直轄市，轄區(qū)面積大，農業(yè)人口多，城鎮(zhèn)化率相對不高，農業(yè)在國民經(jīng)濟中仍然占有一定的地位。特別是糧食生產，關系到城市的可持續(xù)發(fā)展和社會的持續(xù)穩(wěn)定。重慶農村貧困人口也有一定比例，發(fā)展糧食生產，也是農民增收致富、擺脫貧困的重要途徑之一?？茖W預測重慶市糧食產量，是制定地區(qū)農業(yè)發(fā)展和糧食生產戰(zhàn)略規(guī)劃的基礎，對推進糧食生產，確保糧食穩(wěn)定，維護糧食安全具有重要意義。關于糧食產量預測問題，國內學者進行了大量的研究，主要方法有：自回歸滑動平均模型、指數(shù)平滑法、多元線性回歸、主成分分析、灰色模型、神經(jīng)網(wǎng)絡等[1-6]。但這些方法的預測精度都不盡如人意?？柭鼮V波算法是一種最優(yōu)自回歸估計技術，它以最小均方誤差為估計準則，構造一套遞推估計算法，實現(xiàn)對狀態(tài)變量的發(fā)展趨勢的估計，即利用狀態(tài)變量前一時刻的估計值與現(xiàn)時刻的觀測值來更新對狀態(tài)變量的估計，以獲得狀態(tài)變量現(xiàn)時刻的最佳估計值。該算法預測精度高，計算簡便，在許多工程領域得到了應用[7-14]。文中運用卡爾曼濾波算法預測重慶市的糧食產量，提高了預測的準確性。

一、卡爾曼濾波算法基本原理

卡爾曼濾波算法本質上屬于一種遞推反饋算法?；驹硎且詋-1時刻變量的最優(yōu)估計為基準，預測k時刻變量的狀況，同時對該狀態(tài)進行觀測，通過分析觀測值和預測量的差異，采用觀測值和預測量相結合的辦法對預測值進行修正，從而得到變量k時刻的最優(yōu)狀態(tài)估計?？柭鼮V波算法可分為兩個部分：即時間更新方程與測量狀態(tài)更新方程，前者負責遞推，后者負責反饋，反饋在整個算法中較為關鍵，它將先前的狀態(tài)和新的測量值相結合，經(jīng)過比較分析后，以二者為基礎構造改進規(guī)則，得到其最后的最佳估計。即根據(jù)測量值和估計值這兩個量之前的表現(xiàn)，各自給它們分配一個權重，之前表現(xiàn)越好的量其權重就越高，否則，權重就越低，分配的權重也就是所謂的卡爾曼增益。表現(xiàn)好意思是測量值或估計值方差很小，表現(xiàn)不好就是測量值或估計值方差較大。通過一步一步地修正估計值，漸漸使其達到最佳狀況?？柭鼮V波算法的時間更新方程即為預測方程，測量更新方程即為校正方程，分別可表示為：

預估方程：

校正方程：

運用卡爾曼濾波算法預測時間序列時，就是將其視為時間狀況變量，借鑒卡爾曼濾波算法的特點對其發(fā)展趨勢進行估計。在工程應用上狀態(tài)變量有3個值：實際值、觀察值（測量值）、估計值，預測時間序列時，可將實際值、觀察值合二為一，時間序列的實際值即為觀測值，公式中的參數(shù)A = 1，uk= 0，H = 1，B也就無意義?？梢?，運用卡爾曼濾波算法預測時間序列容易實現(xiàn)，所要確定的參數(shù)僅為P0，Q，R三個；P0一般可任意選取，因為算法能自動更正不合理的取值，但注意P0不能為0，否則，濾波器會認為已經(jīng)沒有誤差了，認定初始值就是系統(tǒng)最優(yōu)估計而停止尋優(yōu)過程，故只需選擇合適的Q，R參數(shù)就行。卡爾曼濾波算法流程如圖1所示。

圖1 卡爾曼濾波算法原理圖

二、重慶市糧食產量預測

重慶是我國的重要工業(yè)基地，農村人口比例大，農業(yè)生產在國民經(jīng)濟中占有一定的位置，特別是糧食收入是農民經(jīng)濟收入的主要來源，也是農民脫貧致富的重要途徑之一。重慶市在堅持走工業(yè)強市道路的同時，也不忘農業(yè)生產發(fā)展，積極推進糧食生產，近幾年糧食產量年年豐收，一步一個新臺階。圖2為2000—2017年重慶市糧食產量統(tǒng)計數(shù)據(jù)（數(shù)據(jù)來源于重慶市歷年國民經(jīng)濟和社會發(fā)展公報，其中2017年數(shù)據(jù)來自國家統(tǒng)計局），從圖2可見，重慶市近些年糧食產量除2006年外，基本穩(wěn)定在1 100萬噸左右，近幾年節(jié)節(jié)攀升，為保證全市乃至全國糧食安全、國民經(jīng)濟發(fā)展及社會穩(wěn)定作出了重要貢獻。

圖2 重慶市糧食產量統(tǒng)計數(shù)據(jù)

以 2000—2017年重慶市糧食產量統(tǒng)計數(shù)據(jù)為觀測值，建立卡爾曼濾波，模型參數(shù)設定P0= 1，經(jīng)嘗試，Q = 0.08，R = 0.01，通過MATLAB2014軟件運算，得到了2000—2017年重慶市糧食產量預測值，結果如表2所示，從表2知，模型具有很高的預測精度，平均預測誤差僅為0.400 87%。說明運用卡爾曼濾波算法預測重慶市糧食產量切實可行，效果顯著。

為了檢驗模型的性能，以上述數(shù)據(jù)建立ARMA模型對重慶市糧食產量進行預測，設重慶市歷年糧食產量為時間序列xt，經(jīng)檢驗xt為非平穩(wěn)時間序列，不符合建模要求，對xt進行一次差分后變d(xt)為平穩(wěn)序列，滿足建模要求，經(jīng)分析比對，模型的最佳階數(shù)為ARMA(2，1)，運用EVIEWS8.0軟件，得到模型的參數(shù)估計，結果如表1所示。

表1 ARMA(2，1)模型參數(shù)估計結果

根據(jù)估計參數(shù)的得到的預測方程為：

式中，εt為白噪聲序列。

根據(jù)方程(3)得到歷年重慶市糧食產量預測值（由于原始數(shù)據(jù)經(jīng)過了差分處理，且與前期預測誤差有關，故前幾期的值不能預測），結果如表2所示。ARMA模型的平均預測誤差為1.460 81%，卡爾曼濾波模型比ARMA模型的平均誤差減小了72.558 4%。可見卡爾曼濾波的優(yōu)勢，對于這樣龐大的數(shù)據(jù)序列，預測誤差很小，可謂高效、實用、可靠。兩種模型的預測曲線如圖3所示，卡爾曼濾波預測曲線與實際曲線幾乎重合。

表2 重慶市糧食產量預測結果及比較

2012 1 138.50 1 135.924 -0.226 24 1 137.608 0.078 32 -65.381 9 2013 1 148.13 1 146.248 -0.163 96 1 147.067 0.092 58 -43.536 8 2014 1 145.50 1 147.186 0.147 15 1 145.658 -0.013 82 -90.608 0 2015 1 154.89 1 155.993 0.095 51 1 153.957 0.080 75 -15.451 9 2016 1 166.00 1 164.971 -0.088 22 1 164.783 0.104 335 18.266 83 2017 1 167.20 1 167.433 0.019 987 1 166.956 0.020 915 4.643 018

圖3 預測曲線及比較

至于對 2018年重慶市糧食產量的預測，由于沒有觀測值，卡爾曼濾波無法獨立實現(xiàn)，此時可以借助ARMA模型，ARMA模型可以預測延后一年的值，因此，可將ARMA模型預測得到的2018年的值作為觀測值，再利用卡爾曼濾波預測2018年的值。從表2來看，以ARMA模型預測2018年的值作為觀測值是可行的，該模型2017年的預測誤差小于0.5%，2018年預測值應該有一定的可信度。由ARMA(2，1)預測得到2018年重慶市糧食產量為1 173.704萬噸，將這一數(shù)值加入到原時間序列，通過卡爾曼濾波得到2018年重慶市糧食產量預測值為1 171.875 2萬噸。

2017年我國糧食雖然也獲得了豐收，但與歷史最高產量比仍有所下降。而重慶逆勢而上，今年糧食產量接近歷史最高值，但決不能沾沾自喜，產生驕傲自滿的情緒，認為我國糧食安全可以高枕無憂。應要繼續(xù)抓好糧食生產，再接再厲，再創(chuàng)佳績，為全國糧食安全作出應有的貢獻。同時，在糧食生產的新形式下，發(fā)展糧食生產也要盡量貼近市場，優(yōu)化種植結構，實現(xiàn)糧食生產產量、質量和效益的統(tǒng)一，經(jīng)積極探索發(fā)展優(yōu)質高效的現(xiàn)代糧食生產模式，提高種植的科學技術水平，提高糧食的質量和品質，更好地滿足人們吃飽、吃好、吃安全的要求，提高糧農的經(jīng)濟效益，促進農民增收致富，早日邁入小康社會。

三、結 語

隨著我國從農業(yè)大國到工業(yè)大國的轉型，農業(yè)在國民經(jīng)濟的比重逐步下降，糧食消費在整個商品消費中的比例也呈逐年下降的態(tài)勢，導致一些人產生了農業(yè)不重要、忽視農業(yè)、特別是放松糧食生產的錯誤觀念。農業(yè)始終是支撐整個國民經(jīng)濟發(fā)展與進步的保障，糧食永遠是國家第一重要的戰(zhàn)略資源。不管在什么時候、在什么形勢下糧食安全都是一切工作之根基。要實現(xiàn)我國糧食高產穩(wěn)產，維護糧食長久安全，首先，要切實保護好基本農田，最低保障性耕地這條底線絕不能逾越。隨著我國工業(yè)化、城鎮(zhèn)化進程的推進，各地建設用地都在快速飆升，對基本農田構成了極大的威脅，有些地方政府片面追求經(jīng)濟的一時發(fā)展，認為農業(yè)用地沒有工業(yè)、商業(yè)等用地的效益高，擠占毀壞農田，這是十分危險的想法，它不僅動搖了糧食安全的基石，也阻礙農業(yè)農村的發(fā)展，并破壞了國民經(jīng)濟持續(xù)健康發(fā)展的基礎。其次，要調動農民種糧的積極性。傳統(tǒng)的種糧模式和現(xiàn)行的糧食購銷體制的確使種糧農民的經(jīng)濟效益下滑，導致許多地方農民不愿意種糧，特別是農村青年，嫌棄種田累、種田苦、種田窮，不愿意與土地打交道，大多選擇外出務工另求生計，導致一些地方的耕地荒廢、農田退化，令人惋惜心痛。要改變這種現(xiàn)狀，必須深化農業(yè)改革，積極推進土地確權工作，加快推動農村土地流轉制度，加大對種糧農民的直接補貼力度，完善糧食等農產品價格形成機制，確保種糧農民的利益，不讓農民吃虧，使農民愿意種糧，樂意種糧。第三，要加強對農田基礎設施建設。長期以來，我國農業(yè)基礎比較薄弱，與經(jīng)濟社會發(fā)展不同步、不相稱，特別是農村實行聯(lián)產責任承包制以來，糧食生產以分散經(jīng)營為主，農村集體經(jīng)濟削弱，農田基礎設施建設無法依靠單門獨戶、單打獨斗來完成，因此，許多地方大型基礎水利設施建設基本空白，既往的一些農田水利設施也年久失修，有的已經(jīng)廢棄，有的不能完全發(fā)揮作用。糧食生產抵抗自然災害能力差，不少地方依然靠天吃飯，糧食收成得不到保障。要大力發(fā)展和壯大農村集體經(jīng)濟，創(chuàng)新農田水利投資方式，形成多元化投資體制機制，按照市場經(jīng)濟規(guī)律，建立政府、農村集體、農民個人、社會各界共同參與的多層次、多渠道、多元化的農田水利投入機制，加強和完善農田基礎設施建設，為糧食生產提供可靠的水利灌溉保證，確保農業(yè)增產增效、農民增收。第四，要加強農業(yè)科技服務。目前，我國農業(yè)科技人才隊伍相對比較薄弱，農業(yè)科技成果轉化率低，僅有的一些農業(yè)技術人員，由于農村條件差，農業(yè)工作待遇低，不肯下到基層從事農業(yè)技術推廣服務，導致我國多數(shù)地區(qū)糧食生產技術還比較落后，農民種田的科學化水平不高，一些糧食生產新技術難以得到推廣，影響了糧食生產的產量和品質，導致糧農賣糧難、賣糧苦。針對這種情況，必須加快農業(yè)人才隊伍建設，積極創(chuàng)造條件，鼓勵農業(yè)科技人員到一線工作，推進農業(yè)科技成果轉化率，吸引大學生到農村就業(yè)創(chuàng)業(yè)，積極落實大學生村官到農村任職的政策措施，改善農村干部隊伍結構，提高農村管理水平和農民耕作技術，不斷提高糧食產品的質量，更好地滿足市場需求。

卡爾曼濾波算法是一種挖掘數(shù)據(jù)內部隱含的信息，并通過不斷改進和不斷完善算法而處理不確定信息的先進技術，它以線性最小方差估計方法給出狀態(tài)的最優(yōu)估計值，能夠得到最接近狀態(tài)真值的估計值，計算方便，預測精度很高，在工程領域得到了廣泛的應用。文中利用卡爾曼濾波算法對重慶市糧食產量進行了預測，取得了很好的效果，平均誤差僅為0.400 87 %，比ARMA模型的1.262 96%減小了68.259 5%。預測得到2018年重慶市糧食產量為1 171.875 2萬噸。

參考文獻：

[1] 舒服華．基于ARMA模型的我國小麥產量預測[J]．現(xiàn)代面粉工業(yè)，2017（4）：26-30．

[2] 曹飛飛．灰色系統(tǒng)理論在糧食產量預測中的應用[J]．數(shù)學的實踐與認識，2017（13）：310-312．

[3] 田秀芹．基于多元線性回歸的糧食產量預測[J]．科技創(chuàng)新與應用，2017（16）：3-4．

[4] 樊超，郭亞菲，曹培格，等．基于主成分分析的糧食產量極限學習機預測模型研究[J]．糧食加工，2017（2）：1-5．

[5] 鄭建安．主成分和BP神經(jīng)網(wǎng)絡在糧食產量預測中的組合應用[J]．計算機系統(tǒng)應用，2016（11）：274-278．

[6] 郭亞菲．基于小波-BP神經(jīng)網(wǎng)絡的糧食產量預測模型[J]．糧食加工，2017（5）：1-4．

[7] 修春波，任曉，李艷晴，等．基于卡爾曼濾波的風速序列短期預測方法[J]．電工技術學報，2014（20）：253-259．

[8] 趙攀，戴義平，夏俊榮，等．卡爾曼濾波修正的風電場短期功率預測模型[J]．西安交通大學學報，2011（5）：47-51．

[9] 謝合亮，張砣．卡爾曼濾波在高頻金融時間序列模型預測中的應用[J]．統(tǒng)計與決策，2017（13）：82-84．

[10] 南亞翔，李紅利，修春波，等．基于卡爾曼濾波的空氣質量指數(shù)預測方法[J]．環(huán)境科學導刊，2016（3）：80-84．

[11] 許哲明．基于卡爾曼濾波的地表移動變形預測[J]．沈陽工業(yè)大學學報，2017（5）：557-561．

[12] 張恒德，咸云浩，謝永華，等．基于時間序列分析和卡爾曼濾波的霾預報技術[J]．計算機應用，2017（11）：3311-3316．

[13] 張智勇，張丹丹，賈建，等．基于改進卡爾曼濾波的軌道交通站臺短時客流預測[J]．武漢理工大學學報（交通科學與工程版），2017（6）：974-977．

[14] 彭湃，程漢湘，陳杏燦，等．基于自適應卡爾曼濾波的鋰電池SOC估計[J]．電源技術，2017（11）：1541-1544．