11年中国粮食生产周期的提出及预测实践

1、  11年中国粮食生产周期的提出及预测实践   1988年3月至6月，我曾在吉林大学欧阳植教授指导下， 应用 时间序列 分析 方法 ，进行过太阳黑子相对数预测预报 研究 工作。在这项研究过程中，我大量阅读到了有关气候、农业、疾病等与太阳黑子变化有一定关系的专著和论文，对太阳黑子相对数11年周期有了进一步的了解；但同时也发现：缺少更有说服力的证据和定量分析结果证明这一观点。我们曾企图进一步研究太阳黑子相对数同气候或疾病之间变化的关系，因缺少相关资料只好做罢。 1988年8月，我来到甘肃省统计局从事农产量抽样调查工作，使我有机会全面了解农业生产的全

2、貌和获得比较完整的粮食产量 历史 资料。1989年3月至6月，在进行甘肃省粮食产量周期研究与预测的过程中，我应用了时间序列分析的boxjenkins模型；并发现尽管甘肃省粮食产量历史数据的平稳性较差，但所做的拟合预测具有一定的可信度。这一结果直接鼓励我将研究对象，转到了 中国 粮食生产周期的研究与预测；并于1989年7月完成了中国粮食生产周期的研究与预测一文，提出了11年中国粮食生产周期的观点。在当时对1989年粮食产量普遍预计预测减产的情况下，根据11年周期和预测结果，大胆提出了1989、 1990年将连续创历史最高水平的预测预报。同年8月，受到了当时农业部长何康的肯定；同时，论文也分别在统

3、计与决策1989·5， 发展 1989·6上发表，并且前者加一编者按。1990年6月，又将当年的预测结果寄呈何康部长。19911992年进行了跟踪预测。 1993年7月，又一次在普遍预计预测全年减产的形势下，给有关领导人寄呈了1993年全国粮食产量将创历史最高水平的预测报告。 1994年6月又撰写了1994年及未来我国粮食产量的预测与分析一文，并将此文上报国家统计局、农调总队、农业部等有关部门及领导人。1994年7月，此文和中国粮食生产周期的研究与预测作为“全国第三届市场 经济 与统计改革研讨会”的征文参加研讨，后一篇论文还被评为优秀论文二等奖。从过去已知的预测结果和实践看

4、：此预测预报模型不仅准确地预测了未来中国粮食生产的增减趋势，而且19891993年预测值与实际数的误差在-0.370.95之间，都非常小，具有很好的可信度。 现在将这几年的研究成果，分以下四个部分：(1)中国粮食生产周期的研究与预测的主要 内容 及预测结  果；(2)预测结果与实际数的对比及相关指标的预测；(3)未来中国粮食产量的预测与分析；(4)结束语，向诸位做一介绍。 一、中国粮食生产周期的研究与预测的主要内容及预测结果  (一)方法简介 一个经济时间序列xt(t1,2,n)，通常认为由三种因素组成，即长期趋势、季节因素(周期因素)、随机因

5、素。本文选用自相关分析图来判别序列的平稳性与周期，并且通过自相关和偏相关分析图，确定arima模型的自回归阶与动平阶。 自相关系数 n-kt=1(xt-x)xt+k-x)  rknt=1(xt-x)2其中x为xt(t1,2,n)的平均值， rk为滞后k期的自相关系数。(1)平稳性识别如果rj(j1,k)随着j增大而迅速靠近零，或散乱地分布在零点周围，则认为序列平稳；否则非平稳。对于非平稳序列，通过差分，消除其趋势。(2)周期识别对于一平稳序列，观察其自相关分析图，如果每隔时间t，自相关系数显著偏高，可以认为该序列具有周期t；否则，无周期(无季节性)。偏自相关系数在

6、已知自相关系数的条件下，解如下一系列方程组：1 r1 rk-1k1 r1r1 1 rk-2k2r2   rk-1    1 kk rk得到偏自相关系数11，22，kk。然后根据自相关系数和偏自相关系数的截尾与拖尾确定自回归阶p与动平均阶q。参数估计。 (二) 周期的确立“经济周期也叫商业循环或经济波动，这种周期在每一次重复出现中，周期的长度和振幅都不同，完全相似的经济周期是不存在的”，“只有经济从一个高峰到另一个高峰，或者从一个底谷到另一个底谷，波动的时间在15个月以上，才可以算作一个经济周期”。这种对经济周期的定义，笔者认为着重于对经济的事后分析而不是

7、事前预测，并且忽略了两个极其重要的事实：即经济数据的离散性和由于经济增长趋势的 影响 ，可能使后期的底谷值大于前期的高峰值，这样就使对底谷的预见性大幅度降低，造成了对经济周期研究困难性进一步加大。可以考虑通过差分方法剔除趋势项，并把根据自相关分析图得到的周期视为经济周期。它如同一个数学函数如 y=sinx周期(最小正周期2)，并不规定周期有某种特定的始末点(如高峰点或低谷点)，仅仅要求它使函数满足sin(2+x)sinx(sin(2k+x)sinx，k=±1,)，或具有统计意义上的类似关系。当然一个周期是否有效或有意义，就看采用这个周期所做的预测模型，是否能够拟合和反映现实。经过差分

8、、 计算 、分析、筛选和拟合预测结果及对比(见表1)，最后确定中国粮食生产周期t11年。表1 19851988年粮食产量拟合预测及误差 单位：万吨、   预测值  误差年份   实测值    t=5年    t=11年 t=5年 t=11年1985 38060.638458.8038885.80 1.36 2.161986 39285.138614.3538689.17-1.71 -1.521987 40297.742271.6340352.51 4.89 0.131988 39401.042903.56

9、39814.49 8.89 1.05资料来源:19491984年中国农业光辉成就，全国农业 统计资料(续编)19841987，1988年是公布的实测数。(三) 粮食产量预测根据19501988年粮食产量数据，对19892000年产量进行多步预测，分别选用arima(2,2,3)、arima(2,2,2)模型，预测方程如下：(1-1b-pbp)·(1-b)dxt=(1-1b-qbq)t其中：p=2，d=2，q=3，21989年粮食产量将创历史最高实测水平预测表明，1989年实测粮食产量最低限度将达到 40561.8841136.33万吨，比1988年增加1160.88  17

10、35.33万吨，增长2.954.40；比1984年增加222.78797.20万吨，增长0.551.97，将创历史最高实测水平，也有很大希望突破历史最高年报数40730.5万吨。19902000年的趋势预测从预测结果看(见表2)，1990年将继续保持增长趋势，重写 历史 最高水平，并且在19891995年可能形成类似于19781984年的高速 发展 时期；同样，高速增长之后出现徘徊的可能性依然存在。从现实和长期发展的角度出发，为使19952000年不重蹈19851988年之辙，我们必须在争取高速发展态势的情况下尽可能避免较大波动和徘徊，下面着重根据这一态势的发展趋势进行讨论(略)。表2 中国

11、粮食产量的多步预测值 单位:万吨   预 测 值 d=2 t=11年年份 实际数*  p=2 q=3 p=2 q=2 p=2 q=21989 40561.88 40684.48* 41136.33 40754.91990 42456.92 43067.69 43917.83 44624.31991 40502.72 41882.25 42648.74 43529.31992 40693.84 41985.94 42862.56 44258.01993 42153.14 44471.46 45604.66 45643.51994 44522.33 47360.75

12、48662.98 1995 45202.52 48557.13 50017.03 1996 41967.14 45784.26 47457.03 1997 42211.14 46480.00 48378.87 1998 42218.57 46940.52 49063.70 1999 40292.38 45450.66 47814.33 2000 40399.42 46098.91 49172.96 * 1989年的预测值为修正值，19902000年为多步预测值。实际数栏为此次所加，实际数为公布的实测数。 (四) 太阳黑子 理论各国 科学 家一直致力于太阳活动(其主要标志是太阳黑子相对数)同地球

13、大气、气候、磁砀等一系列物理现象变化关系的 研究 ，证实了太阳对这些变化 影响 的存在。同时直接对太阳黑子相对数同粮食产量进行相关研究发现，太阳活动变化(通过气候、降雨量和 自然 条件变化的作用)经常最终导致粮食单产以520的幅度增减。但是由于太阳异常活动在不同地域范围常常导致了互异的结果，太阳黑子相对数同粮食产量又不存在整体上一致的正相关或负相关，也没有发现它们周期之间有某种确定的关系。所以，尽管早在1884年英国 经济 学家威廉·斯坦利·杰文斯(william stanley jenvons)就提出太阳黑子同经济特别是农业存在着一定的联系，但并没有引起经济学界应有的关注

14、和认可。然而，现在令人感到惊奇、兴奋和有趣的是，中国粮食生产同太阳黑子相对数变化存在着相同的11年周期，而 应用 这个周期对我国最近四年粮食产量的预测又十分有效。虽然我们不能由此得出太阳活动 规律 性变化决定了我国粮食生产周期，但至少可以说太阳活动的变化同其它因素的共同作用，使中国粮食生产选择并呈现出11年的周期。同时，非常有意思的是，最近两次太阳黑子大爆炸分别始于1978年和1988年，间隔时间恰好是11年的周期，根据预测， 19881992年将是太阳活动的22号峰年期，使19881995年太阳黑子相对数呈现出相似于19781984年的规律，从而可能导致形成类似的自然气候条件和发展农业的外部

15、机遇。二、预测结果与实际数的对比及相关指标的预测现在，从以下几个方面分别考察11年中国粮食生产周期和预测模型的有效性及预测值的准确性：从表2多步预测值可以看出：1989、1990年将连续创造历史最高水平，并且1991、1992年将处于一种调整期，而1993年将再创历史最高纪录；19891993年粮食产量的实际增减趋势和运动轨迹完全证明了这一趋势预测的正确性。19891993年一步预测值与实际粮食产量的比较。从表3可见，此模型不仅准确地预测了未来中国粮食生产的趋势，而且预测误差在-0.370.95之间，都非常小，具有很好的可信度。同时，以后逐年所做的多步预测值和1989年所做的预测值的趋势相同。

16、19751993年，相隔11年两年间的粮食播种面积、单产和粮食产量增减趋势表现出惊人的相似性和一致，特别是单产、总产量之间的增减趋势完全一致，进一步证实了11年中国粮食生产周期的存在，见表3、表4。同时，从表3还可以看出：1978年以后，尽管中国 农村 实行了家庭联产承包责任制，七十年代末八十年代初，粮食生产呈现出高速非常规增长态势。但19891993年中国粮食产量的实际增减趋势变化，可以说是19781982年粮食产量增长变化的翻版和重现1978、1989产量都比上年增长；1979、1990都比上年有大幅度增长(9.0、 9.5)，分别跨上3.3亿吨、4.4亿吨两个新台阶；1980、1981年

17、处于一种调整期；1991、1992年同样处于一种调整期；1982、1993都分别创造了(当时的和现在的)历史最高水平，运动规律完全相同。表3中国粮食总产量、一步预测值及误差 单位:万吨年份   实 际 数 年份 实 际 数预 测 值误 差1978 30476.5 1989 40754.941136.33 0.951979 33211.5 1990 44624.344534.68-0.201980 32055.5 1991 43529.343556.48 0.061981 32502.0 1992 44258.044095.87-0.371982 35450.0 1993 4

18、5643.546030.54 0.85表4 全国粮食播种面积、单产、总产量增长率增减表年 份 75 76 77 78 79 80 81 82 83 84 85增长率 + -+ - +* -* - -+ -* +* -* -=年份 86 87 88 89 90 91 92 93 94 95 96增长率 + + - + +* - -+-+ -+ -+ -注：+代表正增长、*增长率>4，-负增长、=负增长率<-4；第一、二、三三个符号分别代表面积、单产、总产量的增长率。9496年为预测趋势。 粮食播种面积、单产、分季作物产量的预测1992年，在得到完整的中国粮食播种面积、单产、分季作物粮

19、食产量历史资料之后，我又进一步对它们分别进行了 计算 、拟合、 分析 和预测。得到如下结论：尽管它们都同样表现出11年周期性变化，但应用11年周期，它们的预测误差依次呈递增趋势。因为，面积具有平稳性；单产波动幅度大，平稳性差；而分季作物产量之间具有互补性；同时，面积因增减幅度小，导致预测趋势有时与实际相反。 三、未来粮食产量预测与分析现在，根据19501993年中国粮食产量历史数据，继续运用11年周期和boxjenkins模型，对1994年产量进行预测；同时，把1994年的预测值看做实际数，对1995-2000年进行多步预测。1994年粮食产量的预测预测结果表明：1994年粮食总产量将达到45305.846433.0万吨，较大可能达到45616.6万吨左右，将1993年(45643.5万吨的水平)持平。从现有的情况分析，由于今年粮食种植结构的调整，粮食播种面积减少；再加上干旱、洪涝等灾情又交替出现，导致夏粮减产(-4.2)，对秋粮生产也不利，从而使全年呈现出减产趋势。但从种种预测结果看，既使全年减产，最低限度也将达到 44435.6万吨以上，更大幅度减产的可能性非常小。同时，(一旦减产)这将成为自1975年以来，第一次打破相隔11年粮食单产、总产量