基于最小一乘准则的中国粮食产量与影响因素的相关性分析_精品.doc

基于最小一乘准则的中国粮食产量与影响因素的相关性分析 Chinese with English abstract 0 引 言 最小一乘准则是 1757 年数学家 Boscovitch 和 1793 年 Laplace 在研究拟合直线时考虑和研究的问 题，比最小二乘法还早了 40 多年。最小一乘法由 于存在着绝对值而不便于计算，这使得它的发展长 期处于停滞状态，但对它的研究一直没有中断。 1987年9月瑞士召开了有关这个准则研究的国际会 议；陈希孺[1]系统提出了最小一乘线性回归问题； Gu Lemin[2]提出了线性函数类最小一乘问题的一 般解；谢开贵等[3]对最小一乘的回归模型进行了研 究；冯守平[4]提出了全最小一乘法；顾乐民[5-6]给出 了曲线拟合最小一乘法的一般解并提出了更具实收稿日期：2013-05-01 修订日期：2013-05-26 作者简介：顾乐民（1952－），男，主要从事最佳逼近问题研究与应用。 上海 同济大学材料科学与工程学院，201804。 Email：gulemin@//0. 用的广义最小一乘法；Gu Lemin[7]对广义最小一乘 法做了补充，并在生物学、生命科学中作了应用。 由于最小一乘解的问题近年来有大的突破，使得非 线性最小一乘问题的实现成为了可能。 影响粮食产量的因素很多，对于改革开放 30 a 来影响粮食产量的主要因素，姜长云等[8]从 4 个不 同时间阶段角度；刘东等[9]从土地资源承载力角度； 李茂松等[10]自然灾害影响角度；刘忠等[11]从土地利 用和人口密度角度对此进行了分析。 Cobb Douglas生产函数，简称为C-D生产函数， 是一个典型的多因素数学模型[12]。宋廷山等[13-15] 通过 C-D 生产函数模型，将粮食产量与主要影响因 素建立关联并进行相关性分析。在粮食产量预测方 面，梁仕莹等[16]用 3 种不同的方法；季翔等[17]用内 涵解析方法对中国今后的粮食产量进行了预测。 大量数据处理表明，在各种影响因素中，有些 是呈幂指数关系，而有些呈指数关系，因此构成一 ?综合研究? 农业工程学报2013 年2 种混合型的生产函数关系。混合型的生产函数模型 是对纯幂指数的 C-D 生产函数模型的推广，构成一 种广义的生产函数模型，由于获得更小的拟合误 差，从而具有更好的描述性，并构成了一个在理论 上和应用中都值得探讨的课题。 本文主要从数据处理的角度，将影响粮食产量 的 5 个主要因素（化肥施用量、粮食播种面积、成 灾面积、农业机械总动力、第一产业就业人数）与 中国粮食产量建立关联；对生产函数数学模型进行 拟合；在获得最小一乘最佳逼近结果的基础上，对 数据结论进行解释和分析。 1 模型与方法 1.1 函数定义 1）广义的生产函数模型（定义 1） 大量数据处理结果表明，各种投入与产出之间 的函数关系，用指数与幂指数的混合型数学模型， 即广义的生产函数模型能取得更好的结果。 设投入量总数有n 个，用 xx1，x2,…，xn表示； 参数 aa0，a1，a2，…，an为 1n + 个不全为零的实 数；广义的生产函数 fx,a是 x 和 a 的函数，为指数 和幂指数混合型，用数 j 作为二者间量化的区分，即 11 1 2 2 0 1 1 0 1 , e e. e , 1 1 n n k k k a aa x a x n n j n a x a k k k j f a x x a x j n ? +∏ ∏ x a ≤ ≤ （1） 当 j1 时 1 1 e k k j a x k ∏ 项不存在，式（1）为 C-D 生 产函数数学模型 1 2 0 1 2 0 1 , . n k n a aa a n k k f a x x x a x ? ∏x a（2） 当 jn+1 时 k n a k k j x ∏ 项不存在，式（1）为指数型 生产函数模型 1 1 2 2 0 0 1 , e e. e en n k k n a x a xa x a x k f a a ∏x a （3） 文中采用的是 n5 的指数型生产函数模型。 2）指数型生产函数模型（定义 2） 3 3 5 51 1 2 2 4 4 3 3 5 51 1 2 2 4 4 0 5 0 0 1 , e e e e e e e k k a x a xa x a x a x a x a x a xa x a x a x k f a a a++ + + ? ∏ x a （4） 在 x 上取 m 个点，构成 i1,2,…，m，m≥n+1 个离散数据点 1, 2, 3, 4, 5, 6, , , , , , i i i i i i ix x x x x xx （5）