【毕业学位论文】（Word原稿）粮食安全问题研究——以安徽省为例-统计教育学

粮食安全问题研究 以安徽省为例 摘 要 安徽省在全国粮食生产中占有特殊的地位， 对于稳定我国粮食生产起着重要作用。本文以安徽粮食安全发展状况为背景，通过建立适当的经济模型和数学模型，解决安徽粮食安全相关问题。 首先，利用主成分方法建立粮食安全综合指标。然后利用灰色关联分析方法对众多影响粮食安全因素进行排序，得出以下结论：影响程度从大到小依次为 人口数量 、 粮食作物播种面积 、 农业劳动人口 、 乡城收入比 、 单位面积粮食产量 、 受灾面积 、 城市化率 、人均国民生产总值增长率 。 其次，通过建立向量自回归模型分析各因素对粮 食安全的影响，并且进行了协整分析、脉冲响应分析。结果表明： 粮食作物播种面积、单位面积粮食产量、乡城收入比、人均国民生产总值增长率 、城市化率对粮食安全起到正向作用，受灾面积、农业劳动人口、人口数量对粮食安全起到反向作用。 接着，结合粮食安全指数及其各影响因素之间的关系建立起的协整方程，从机制设计理论的角度出发，为解决我国粮食安全问题提出关键影响环节，并就关键影响环节展开分析，提出合理的解决粮食安全的机制解决办法。 最后，文章在以上研究的基础上，采用神经网络预测方法对安徽省 2008预测，预测结果显示： 2008 年和 2011年安徽省粮食安全程度较高；2009年和 2010年粮食安全程度较低，尤其是 2010 年不容乐观，文章据此 对安徽省粮食安全的发展状况 提出相关政策建议。 关键词： 粮食安全；主成分； 冲响应；机制设计；神经网络 一、基本现状与问题的提出 民以食为天，国以粮为本。 粮食是关系国计民生的重要战略物资，粮食安全与社会的和谐、政治的稳定、经济的持续发展息息相关。 保障粮食安全，对实现全面建设小康社会的目标、构建社会主义和谐社会和推进社会主义新农村建设具有十分重要的意义。当前 我国粮食安全形势总体是好的，粮食综合生产能力稳步提高，食物供给日益丰富，供需基本平衡。但我国人口众多，对粮食的需求量大，粮食安全的基础比较脆弱。从今后发展趋势看，随着工业化、城镇化的发展以及人口增加和人民生活水平提高，粮食消费需求将呈刚性增长，而耕地减少、水资源短缺、气候变化等对粮食生产的约束日益突出。我国粮食的供需将长期处于紧平衡状态，保障粮食安全面临严峻挑战。 安徽省在全国粮食生产中占有特殊的地位， 对于稳定我国粮食生产起着重要作用。一是因为安徽是一个传统的农业大省和产量大省，是全国 13 个粮食主产省和全 国 6 个粮食调出省份之一，年均调出粮食 500 万吨 。二是因为安徽是全国农村改革的先行省份，1978年凤阳小岗村率先实行大包干， 2001年全国农业税费改革在安徽省全省试点， 2005年在全国率先取消了农业税。三是安徽的农业发展状况在全过具有代表性。 因此，如何构建适应社会主义市场经济发展要求和符合安徽省省情的粮食安全保障体系，对于维持社会稳定与国家自立具有重要意义。安徽省基本概括如下： 全省耕地面积 408万公顷，土地肥沃，适宜各种农作物生长。长江和淮河自西向东横贯全境，其中：淮河流域 长江流域 新安江 流域 安徽农业的大发展源于 30年前发生在江淮大地上的一场变革。 1978年夏天，安徽大旱。夏收后，凤阳县小岗村 18户农民签下了分田到户的契约，拉开了大包干的序幕。第二年，粮食产量便达到了 66吨，相当于 1966年 1970年的总和。大包干的成功为安徽的粮食生产开辟出一片新天地。农村改革 30年来，安徽粮食生产的总趋 势是上升的（见图1），粮食产量由 幅达 年平均递增 安徽实现了由农业省向农业大省、产粮大省的转变，粮食年产量常年居全国第 6位。随着粮食生产水平的提升，安徽对全国粮食安全的贡献越来越大，不仅解决了本省人民的吃饭问题，而且每年调出粮食约 500万吨。安徽省的粮食生产还是全国粮食增产的助推器，2006年安徽粮食增产 到了全国总增量的 40%。 2008年，安徽省夏粮总产 1170万吨，创历史最好水平，增量占全国的 20%。 但是 ，在取得这些成就 的同时，安徽省的粮食安全也存在一些潜在的威胁。 自改革开放以来 ， 随着安徽经济快速发展 、 人口增长、工业化和城市化进程的加快 ， 耕地资源日益受到工业和城市土地利用的经济竞争 ， 部分耕地非农利用的趋势不可逆转 ； 同时 ,由于水土流失、工业废物排放、农药化肥的大量使用 ， 不合理的灌溉方式及自然灾害等导致耕地质量下降 ， 制约着安徽省的粮食生产与粮食安全。 图 1 安徽省 1978二、理论背景和文献综述 在中国，粮食安全是指主要粮食品种的产和销、供与求动态的平衡，及应对突发事件的实力和 能力。 粮食安全的最终目标 是 ，确保所有的人在任何时候既买得到又买得起他们所需的基本食品。 ” 这个概念包括三个具体目标： 确保生产足够数量的粮食； 最大限度地稳定粮食供应； 确保所有需要粮食的人都能获得粮食。 对于粮食安全，应该从以下几个方面来理解 第一，粮食安全包括宏观安全和微观安全。宏观安全就是从总体、总量上看粮食供求是均衡的微观安全是从局部、个体上看粮食供求是均衡的，任何人都能买得到并且能够买得起为维持生存和健康所必需 的足够食品。只有同时实现了宏观安全和微观安全才算实现了粮食安全。 第二，粮食安全包括短期安全和长期安全。即保证不仅在短期内任何人都能买得到又能买得起为维持生存和健康所必需的足够食品，更重要的是在长期内，在任何时候，任何人都能买得到又能买得起为维持生存和健康所必需的足够粮食。粮食安全的目标应该是实现长期粮食安全。 第三，粮食安全包括量的安全和质的安全。也就是说，粮食安全的目标不仅仅是维持生存还要保证健康不仅仅是从数量上有足够的粮食，还要确保粮食食用卫生、安全。 本文主要是从粮食的长期安全以及粮食量的安全角度出 发，研究安徽省的宏观粮食安全问题。 最早对中国的粮食安全问题进行研究的是美国的布朗， 1994年美国世界观察研究所所长莱斯特布朗（ 华盛顿邮报上发表“中国能使世界挨饿吗？ 它的崛起正在消费掉全球粮食供给”。他在 1995 年又进一步认为，到 2030 年中国由于耕地面积的减少等因素，粮食将减少 20%，如果不考虑膳食结构的改善，中国将进口相当于现在世界粮食贸易总额的粮食，达 这以后，世界上许多组织或个人都对此进行了有针对性的研究评价。美国农业部费 里德克科鲁克认为，布朗博士在长期预测中没有把市场经济的自我矫正机制因素考虑进去，中国政府、生产者和消费者不会对经济环境的变化无动于衷。日本农林水产省白石和良认为，布朗的预测依据不十分明确，过分强调特定条件，因而产生了不符合实际的结论。世界银行专门进行了关于中国长期粮食安全的课题研究，认为中国在 2020年粮食进口量约为 6000万吨，主要粮食出口国完全可以在不大幅度提价的情况下供应这些粮食。关于世界粮食安全形势的分析研究，联合国粮农组织（ 国际谷物理事会（ 对 1997世界粮食形势作 了测算，认为供给形势有所改善，安全系数有所提高。 1999 年在英国举办了关于未来世界粮食形势的研讨会，会后出版了一个名为养活 80 多亿人，对科学的挑战论文集，其中大量篇目充斥着悲观的论调，只有极少数作者持相反态度。 改革开放以来，尤其是布朗博士对中国未来粮食供给表示深深的忧虑以来，我国很多学者都对中国粮食安全状况进行了不同程度的评估。主要分为以下几种： 关于粮食安全测度基准与方法的研究 朱泽（ 1998）通过比较中国与世界重要粮食生产与消费国在粮食总产量 波动系数、粮食自给率、粮食库存水平、人均粮食占有量与低收入人口粮食供应水平等指标方面的差异性，以及粮食安全系数的计算与比较，得出了中国粮食安全整体水平仅次于加拿大、法国、美国与澳大利亚等四国，高于世界平均水平的结论。徐奉贤 (1999)在四指标中加上了贫困者的粮食保障水平，通过五项指标的简单平均计算了粮食安全系数。马九杰等(2001)则认为粮食安全综合指数是食物和膳食能量供求平衡指数、粮食生产波动指数、粮食储备一需求比率、粮食国际贸易依存度系数、粮食及食物市场稳定性等各指标得分值的加权平均值。刘晓梅 (2004)构建了一个四指标加权平均安全系数，认为我国的粮食安全系数是 世界范围来看是比较高的。吕耀等在中国食物保障可持续性及其评价一文中，突破了许多人仅仅或主要从粮食产量与人口这两个因素来考虑粮食供应保障的情况，而着力于从农业生产及其资源生产效率、经济效益、资源利用率以及资源与环境质量等方面来构建中国食物保障可持续性评价指标体系，并依据中国自然资源区划的九大农业区内个县市的有关资料进行了具体评价，并得出了中国区际食物保障可持续水平有较大差异、样本间的食物保障可持续指数差异很大与各区域影响食物保障可持续 性的主要限制因子不同等重要结论。 对中国粮食供求的研究 在粮食生产与供给方面， 1995年莱斯特布朗指出，由于耕地将被工业和城市用地挤占，中国谷物产量将会从 1990年的 020年的 吨，而人口又大量增加，所以世界无法养活中国；让邓肯 （ 1997） 1997年预测中国谷物产量 2010年将达到 吨，其预测结果与莱斯特布朗完全不一样；朱玫（ 1997）对各个阶段影响中国粮食波动因素的粮食产量增长贡献作了实证分析，结论是：单产贡献率是主 要的。 在粮食需求方面，常根发（ 2000）从微观农户消费行为探讨了农村口粮消费的实现模式，农户口粮消费与粮食生产具有阶段性和不一致性，二者受农户消费效用和收入的双重约束；朱彬（ 2001）认为，尽管高收入会带来肉类食品的高消费，但由于消费结构向低粮耗肉类倾斜，故而饲料粮需求的增长速度将被减缓。 对各省份粮食安全状况的研究 目前关于省份地区的粮食安全问题，已有相当多学者对福建、浙江、江苏、陕西、山东、湖北、吉林、河北、湖南等省从影响因素、预警系统的设计等角度进行了研究，如宁自军（ 2003）就粮食购销市场化条件下浙江省粮食安全预警系统的设计及其指标体系的构建进行了初步研究。门可佩等（ 2008）以江苏省粮食生产情况为例，分析粮食安全相关因素，江苏省粮食安全问题存在隐患的结论。李梦觉（ 2009）通过对“粮食安全”内涵和评价指标的研究，分析了湖南省安全状况及其面临的问题。 而关于粮食生产大省 安徽省的粮食安全问题研究寥寥无几。金宝石（ 2005） 在持下 ,建立了 1980、 2002 年两个时期安徽省各县域耕地资源、粮食生产以及影响粮食生产的地形地貌、土壤、光热水等自然条件的空间数据库 和多年社会属性数据库。本文主要从影响粮食安全的社会经济因素角度进行研究，研究各影响因素对安徽省粮食安全的关系，并对未来安徽省粮食安全状况进行预测，分析安徽省粮食安全可能存在的粮食危机，并结合实际提出相关政策建议。 三、数据的收集与指标的选取 根据上文粮食安全的定义，本文将粮食安全指标化为粮食安全指数 Y ，它由人均粮食播种面积影响指数、人均粮食占有量影响指数、粮食总产波动率影响指数及粮食自给率影响指数（设定安徽粮食自给率的安全底线定为 92）用主成分方法求得，旨在衡量一定时期内一国或地区面对世界粮食危机时的自我承受能力。影响一国粮食安全的社会经济因素很多，基于安徽当前社会经济发展的实际，本文分别从粮食的供给和需求的角度选择以下因素作为影响粮食安全的关联因子。 1X ：粮食作物播种面积，反映一定时期内谷物、豆类、薯类总的播种面积，是粮食供给变化的最直接影响因素； 2X ：单位面积粮食产量，是技术进步、投入变化、制度创新、灾害气候等各种因素效应的综合反映；3X：受灾面积，是影响农作物生产的直接因素，短期内对粮食供应产生直接的影响； 4X ：农业劳动人口，反映了一定时期内投入到农业生产中的劳动力数量，是粮食供给的直接影响因素；5X：乡城收入比，反映了农业劳动者从业的机会成本，是农业劳动者固守土地的激励因素，可以间接影响到粮食的供给；6X：人口数量，是粮食需求的直接影响因素；7X：人均国民生产总值增长率 ，指一定时期分摊到每个国民身上的国民生产总值的变化率，用以反映一个国家的经济发展程度，是粮食安全的需求影响因素；8X：城市化率，反映一定时期内城镇人口占总人口的变化情况，该指标通常体现了一国或地区的经济发展水平，是粮食需求的间接影响因素（见 表 1）。 表 1 安徽省各指标数据 影响因素指标 粮食安全综合指标 年份 1X 2X 35 人均粮食占有量 粮食总产波动率 粮食自给率 1995 780 4956 000 996 6029 968 054 997 984 109 998 5991 904 4986 152 999 860 5001 205 6 000 442 048 278 8 235 001 068 325 266 002 155 5080 369 431 003 091 410 2 190 004 119 461 501 005 351 148 516 474 006 429 160 593 639 007 479 208 676 716 、实证分析 安徽当前的粮食安全系统，包含着各种各样的影响因素，多种因素共同作用的结果决定了该系统的发展态势。而哪些是主要因素，哪些是次要因素，则是系统分析中相关人士普遍关心的问题。传统的数理统计中的回归分析、方差分析等都经常被用来进行系统的影响分析，然而这些方法的正常运用都必须具备一定的统计条件，如要求具备相对研究对象个数而言的大量样本数据，要求数据样本服从或近似服从某个典型的概率分布，还要求有相关理论作研究支撑等。尽管已有众多的文献对粮食安全系统进行了各种角度、各种层次的数理研究，但局限数理统计分析方法 的缺陷，一般仅有少数几个研究对象，且对统计数据的性质并没有得到真正严格的要求。从此角度来看，灰色关联分析方法对样本量的多少和样本有无规律都同样适用，也不会出现量化结果与定性分析结果不符的情况，这些在一定程度上弥补了采用传统数理统计方法作系统分析所导致的缺憾。因此，运用灰色关联理论对安徽近年的粮食安全形势进行相关性分析和研究，至少是研究方法上的一种积极尝试。 主成分分析的主要目的是希望用较少的变量去解释原来资料中的大部分变异，将许多相关性很高的变量转化成彼此相互独立或不相关的变量。通 常是选出比原始变量个数少，能解释大部分资料中的变异的几个新变量，即所谓主成分，并用以解释资料的综合性指标 ，所以说 主成分分析实际上是一种降维方法。 在实证问题研究中，为了全面、系统地分析问题，我们必须考虑众多影响因素。这些涉及的因素一般称为指标，在多元统计分析中也称为变量。 本文中人均粮食播种面积影响指数、人均粮食占有量影响指数、粮食总产波动率影响指数及粮食自给率影响指数等都是变量， 因为每个变量都在不同程度上反映了所研究问题的某些信息，并且指标之间彼此有一定的相关性，因而所得的统计数据反映的信息在一定程度上有重 叠。在用统计方法研究多变量问题时，变量太多会增加计算量和增加分析问题的复杂性，人们希望在进行定量分析的过程中，涉及的变量较少，得到的信息量较多。主成分分析 正好可以满足要求。 主成分分析法是一种数学变换的方法 , 它把给定的一组相关变量通过线性变换转成另一组不相关的变量，这些新的变量按照方差依次递减的顺序排列。在数学变换中保持变量的总方差不变，使第一变量具有最大的方差，称为第一主成分，第二变量的方差次大，并且和第一变量不相关，称为第二主成分。依次类推， i 个变量就有 i 个主成分。 其中p 维正交化向量（ ），i 称 的第 i 个主成分。设 X 的协方差矩阵为 ，则 必为半正定对称矩阵，求特征值i（按从大到小排序）及其特征向量，可以证明，i所对应的正交化特征向量，即为第 i 个主成分 /， 通常要求提取的主成分的数量 k 满足 i 。 主成分的 分析 步骤为 ： 利用初始化方法对 数据 进行 标准化 处理； 求相关系数矩阵 ； 经过 一系列正交变换，使非对角线上的数置 0，加到主对角上 ； 得特征根i（即相应那个主成分引起变异的方差 )， 并按照从大到小的顺序把特征根排列 ； 求各个特征根对应的特征向 量 用下式计算每个特征根的贡献率iv /； 根据 提取的主成分的数量 k 满足 i 确定提取的主成分个数，并写出主成分表达式。 本文中提取 3个主成分，分别为 9192121111 9292221212 9392321313 其中 ,(321921 ),(321321 本文要研究粮食安全与其影响因素之间 的关系 ，由于衡量粮食安全的因素众多，所以任何一个粮食安全指标不足以表明粮食安全情况，因此，需要考虑综合指标。这里，我们利用前面主成分分析得到的结果，对 4个粮食安全指标进行加权综合，得到一个综合指标 H ，来很好的表示粮食安全。 综合指标表达式为332211 ，其中21 , , 表 2 1978年份 综合指数 年份 综合指数 年份 综合指数 1978 988 998 979 989 999 980 990 000 981 991 001 982 992 002 983 993 003 984 994 004 985 995 005 986 996 006 987 997 007 灰色系统理论是由著名学者邓聚龙教授首创的一种系统科学理论 (其中的灰色关联分析是根据各因素变化曲线几何形状的相似程度，来判断因素之间关联程度的方法。此方法通过对动态过程发展态势的量化分析，完成对系统内时间序列有关统计数据几何关系的比较，求出参考数列与各比较数列之间的灰色关联度。与参考数列关联度越大的比较数列，其发展方向和速率与参考数列越接近，与参考数列的关系越紧密。灰色关联分析方法要求样本容量可以少到 4个，对数据无规律同样适用，不会出现量化结果与定性分析结果不符的情况。 其基本思想是将评价指标原始观测数进行无量纲化处理，计算关联系数、关联度以及根据关联度的大小对待评指标进行排序。灰色关联度的应用涉及社会科学和自然科学的各个领域，尤其在社会经济领域，如国民经济各部门投资收益、区域经济优势分析、产业结构调整等方面，都取得较好的应用效果。 关联度有绝对关联度和相对关联度之分，绝对关联度采用初始点零化法进行初值化处理，当分析的因素差异较大时，由于变量间的量纲不一致，往往影响分析，难以得出合理的结果。而相对关联度用相对量进行分析，计算结果仅与序列相对于初始点的变化速率有关，与各观测数 据大小无关，这在一定程度上弥补了绝对关联度的缺陷。 灰色关联分析的步骤为： 第一步， 确定分析数列 设参考数列（又称母序列）为 ,2,1)( 比较数列（又称子序列） ,2,1)(i i ,2, 。 第二步， 变量的无量纲化 由于系统中各因素列中的数据可能因量纲不同，不便于比较或在比较时难以得到正确的结论。因此在进行灰色关联度分析时，一般都要进行数据的无量纲化处理。 ,2,1,0;,2,1,)1( )()( 第三步， 计算关联系数 )(0 )(关联系数 )()(m a xm a x)()()()(m a xm a x)()(m i nm i n)( 记 )()()( ，则 )(m a xm a x)()(m a xm a x)(m ),0( ，称为分辨系数。 越小，分辨力越大，一般 的取值区间为 1,0 ，具体取值可视情况而定。当 时，分辨力最好，通常取 。 第 四 步， 计算关联度 参考数列与比较数列的关联度为 nk ii )(1 ， .,2,1 第 五 步， 关联度排序 关联度按大小排序，如果 21 ，则参考数列 y 与比较数列 2x 更相似。 在算出 )(列的关联系数后，计算各类关联系数的平均值，平均值( )(结果，如表 3所示： 表 3 1978 2007年粮食安全与各影响因素之间的关联度 人口数量 人） 食作物播种面积 公顷） 业劳动人口 人） 城收入比 位面积粮食产量 斤 /公顷） 灾面积 公顷） 市化率 %） 均国民生产总值增长率 据上述灰色关联模型，选取 1978 2007 年安徽粮食安全综合指数作为系统参考序列，另分别选取安徽 1978 2007 年粮食作物播种面积、单位面积粮食产量、受灾面积农业劳动人口、乡城收入比、人口数量、人均国民生产总值增长率、城市化率作为系统的比较序列，分别计算得各关联度及排序。 从整体来看， 1978 2007 年，各因素对粮食安全的影响程度相差无几。但本文为了进一步研究，还是对影响因素关联度的大小进行了对比。在现有的影响安徽粮食安全的因子中，人口数量因子的影响最大，其次分别是粮食作物播种面积、城乡收入比、农业劳动人口、单位面积粮食产量、受灾面积、城市化率、人均国民生产总值增长率。安徽是一个人口大省，在“民以食为天”的前提下，尽管安省的人口增速有所放慢，但庞大的人口基数显然仍然是影响粮食安全的的第一因素。单位面积粮食产量、受灾面积分别位列第五位、第六位，说明二者对安徽粮食安全并非主要影 响因素。事实上安徽的粮食作物受灾情况虽然每年呈现不同的形式，但总受灾面积一般比较稳定。而安徽的单位面积粮食产量也一直处于全国前列，在没有重大生物技术突破的情况下，可挖掘的空间并不大。也就是说，在此期间，二者几乎可以当作每年都产生的固定因素。粮食作物播种面积、农业劳动人口在安徽粮食安全影响因子中分别排序为第二位、第三位，说明二者是影响安徽粮食安全的主要因素。粮食作物播种面积、农业劳动人口对安徽粮食安全的影响主要作用于粮食的供给方面，结合上文对受灾面积、单位面积粮食产量因子的分析，显然当短期内无法显著提高单位面 积粮食产量情况下，扩大粮食作物播种面积，加大农业劳动的投入才能显著地提高粮食的供给，而现实情况是粮食作物播种面积因为工业化和城市化进程而不断减少，农业劳动力由于机会成本所致而蜂拥流入城市，因此二者成为安徽粮食安全的重要影响因子是由安徽现实的经济发展水平和二元经济结构决定的。另外，城市化率、人均国民生产总值增长率从粮食的需求方面也影响到安徽的粮食安全，二者的关联度均达到 78%以上，但相比其它影响因子，其影响较小。 食安全 因素 分析 为了进一步研究 影响 粮食安全 的因素 ，对数据利用向量自回归模型进行了协整分 析、脉冲响应分析。 时间序列的平稳性检验 对多个时间序列进行协整分析的第一步就是确定每个序列是否平稳，而判断序列平稳性，通常借助单位根检验。常用的单位根检验主要有 于 R(1)过程的单位根检验，因此，在这里我们采用 用 中所涉及的所有变量进行平稳性检验。 表 4 序列 检验类型 (c,t,k) 验值 1%显著 水平 5%显著 水平 10%显著 水平 判断 结论 c,0,1) 平稳 0,t,2) 平稳 c,0,1) 平稳 0,t,3) 平稳 c,0,1) 平稳 0,t,4) 平稳 c,0,1) 平稳 0,t,1) 平稳 Y (c,0,2) 平稳 0,t,1) 稳 c,0,3) 稳 0,t,1) 稳 0,0,5) 稳 c,t,5) 稳 0,0,1) 稳 0,t,2) 稳 c,0,1) 稳 Y (0,t,4) 稳 经检验，时间序列 ,过一阶差分之后都平稳 ,因而是一阶单整的 )1(I 序列。 协整分析 协整检验主要是判断变量间是否存在长期均衡关系，目前主要有两种方法：一是 1987 年提出的基于协整回归残差的两步法检验；二是整检验 ,由 988、 1991)和 990)提出的基于 量自回归 )方法的协整系统检验。需指出的是， 变量协整关系检验通常采用 整检验法。为了克服多变量和小样本条件文采用 法进行协整检验。 表 5 协整检验 原假设 特征值 似然函数 5%显著性水平 1%显著性水平 无 多一个 多二个 多三个 多四个 表 5，从零假设 r=0 开始，似然函数的值为 过 5%显著性水平的临界值 明应拒绝零假设的 r=0，接受 r=1 备择假设，说明变量之间存在长期均衡关系。依据检验结果，可建立 8个解释变量与 Y 之间的协整方程 : 9 5 2 9 5 . 5 8 8 5 从上式可知，各影响因素指标与之间在长期内呈显著相关关系。其中，粮食作物播种面积、单位面积粮食产量、乡城收入比、 人均国民生产总值增长率 、城市化率对粮食安全起到正向作用，受灾面积、农业劳动人口、人口数量对粮食安全起到反向作用。 建立 型 由于上述变量间存在长期均衡关系 ,彼此间相互影响，相互制约，可建立向量自回归模型，并可在向量自回归模型的基础上用脉冲响应函数和方差分解来研究它的动态性质。 一般模型数学表达式为： 111 其中，m 维内生变量向量 ,d 维外生变量向量 , 1A , ,B , , 待估计的参数矩阵，内生变量和外生变量分别有 P 和 r 阶滞后期。t是随机扰动项，与同时刻的元素可以彼此相关，但不能与自身滞后 值和模型右边的变量相关。在实际应用中，通常希望滞后期 P 和 r 足够大，以完整地反映所构造 动态特征。但滞后期越长，模型中待估计的参数就越多，自由度就越小。因此，应在滞后期与自由度之间寻求一种均衡状态，一般根据 息量取值最小的准则确定模型的阶数。经过多次试验，由 型整体检验结果表明，滞后期 P 为 1 时效果最好。每个方程仅有一部分滞后项是显著的，但在建立 型时一般不进行筛选，各个滞后变量仍被保留。 型可以反映粮食安全与各影响因素之间的相互影响程度。 表 6 1）模型结果 Y 2 4 6 8 Y(1(2(3(4(5(6(7(8( C 型 协整分析结果表明了各变量之 间的长期关系及其变动趋势，但上述关系是否构成因果关系还需进一步检验。这里主要通过构建 型和 果检验。这样不仅反映了传统的因果检验所考虑的变量间短期因果关系，而且能够解释两个或两个以上变量间的长期因果关系。由于误差修正模型的滞后期是无约束 型一阶差分变量的滞后期，因此，对应的 后期应为 1，采用协整方程有截距但没有确定趋势的形式，构造 型，各参数估计见表 7。 表 7 型估计结果 D(Y) D(D(D(D(D(D(D(D(Y( (1) (1) (1) (1) (1) (1) (1) (1) C 脉冲响应函数 下图是基于 )和渐近解析法模拟的脉冲响应函数曲线 ,横轴代表响应函数的追踪期数 ,纵轴代表因变量对解释变量的响应程度。下图中实线为响应函数的计算值 ,在模型中 ,我们将响应函数的追踪期数设定为十年。 图 2 图 3 图 4 图 5 图 6 图 7 图 8 图 9 图 10 从图 2中，可以观察到粮食安全对其自身的一个标准差新息立刻有较强的反应，在第 2期变为负值，第 3 期变为正值，第 4期变为负值，之后较为稳定；从图 3可知，粮食安全对单位面积粮食产量标准新息的脉冲反应在前三期稳定，之后上下波动，一直持续到结束。从图 4来看，粮食安全对受灾面积脉冲反应在第一期至第三期为负值，第五期之后较稳定。从图 5 中，农业劳动人口在短期内对粮食安全产生负面影响，但从长期来看，影响是上下波动的，这与经 济理论是一致的。从图 6、图 7、图 8、图 9、图 10可知，在长期内，乡城收入比、人口数量、人均国民生产总值增长率、城市化率对粮食安的影响较为稳定，不会发生较大的变化。 五、粮食安全机制设计 由第四部分中的协整方程可以得到，单位面积粮食产量、城市化率对粮食安全指数有正向影响，而农业劳动人口则对粮食安全指数有反向的影响。单位面积粮食产量增加，粮食安全指数就会增加，这说明技术因素对粮食安全问题至关重要，在技术密集