2011经济增长问题的数学建模论文6

PAGE2GDP与一二三产业的散点图道格拉斯生产函数资本与劳动力讨省内生产总值增长与资本及劳动之间的关系多元线性回归显著性检验道格拉斯函数数据拟合2012高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：聊城大学参赛队员(打印并签名)：1.孔鹏2.郑国梁3.吴其诚指导教师或指导教师组负责人(打印并签名)：日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：假定在一段不太长的时间内技术水平不变.定义与符号说明山东省省内生产总值复相关系数统计量随机误差置信水平工业产值建筑业产值农林渔业产值省内生产总值增长劳动投入量资本投入数量劳动对产出的贡献程度资本对产出的贡献程度投资金额指数就业人数指数省内生产总值指数道格拉斯函数常数模型的建立与求解5.1问题一模型的建立与求解：回归分析方法是统计分析的重要组成部分，用回归分析方法来研究自变量与变量的关系函数是一种常用的有效方法.我们通过回归模型的建立，定量预测了未来经济的发展.5.1.1GDP与工业值、建筑业及农林渔业产值数量模型：通过在互联网上搜集到1981年～2010年，山东省GDP与工业值、建筑业及农林渔业产值的数据（见附表1），可以定性的看出GDP与工业、建筑业及农林渔业产值为整体上升的趋势.为了大致分析GDP与工业值、建筑业及农林渔业产值关系，我们首先利用附录数据做出了GDP与工业产值的关系散点图（如图1）.图1工业产值与GDP散点图 从图可以发现，随着工业产值的增加，GDP的值有比较明显的线性增长趋势.图中的直线是用线性模型拟合的.同理我们也分别作出了建筑业产值与GDP的关系散点图（图2）、农林渔业产值与GDP的关系散点图（图3）.图2建筑业产值与GDP散点图图3农林渔业产值与GDP散点图 通过图2，图3可以看出建筑业产值与农林渔业产值同样有很强的线性关系，同样也分别用直线模型对其拟合.建筑业产值与GDP线性模型农林渔业产值与GDP线性模型因此，综上所述四者之间有很强的线性关系，可建立多元线性回归模型在模型中除了工业，建筑业，农林渔业外，影响国内生产总值的其他因素的作用都包含在随机误差内，这里假设相互独立，且服从均值为零的正态分布，.对模型直接利用matlab统计工具箱求解，得到回归系数估计值及其置信区间（置信水平=0.05），检验统计量，，的结果见表1.参数参数估计置信区间-80.7854[-265.3513,103.7804]1.1521[1.0356,1.2685]5.5154[4.3792,6.6515]1.1139[0.6934,1.5345]表1模型的计算结果5.1.2结果分析:表1显示，指因变量（国内生产总值）的99.97%可由模型确定，值远远超过检验的临界值，p远小于，因而模型从整体上来看是可用的.表1的回归系数给出了模型中，，，的估计值，即，，.检查它们的置信区间发现与的置信区间都包含零点，这表明回归模型常数与回归变量对模型的影响不太显著.这也符合这一事实，农林渔业产值对GDP的影响较小，工业与建筑业对GDP的影响较大.但是一般情况下常数的值都保留在模型中，不剔除.回归变量系数区间右端点相对距零点较远，所以我们也保留在模型中.因此最终确定的模型为：5.2问题二模型的建立与求解：在问题二中为了讨论省内生产总值增长与资本及劳动之间的关系，通过查阅相关资料，我们了解到省内生产总值通常取决于相关的生产资料和劳动力等相关重要因素.于是我们通过建立柯布—道格拉斯生产函数，来探讨省内生产总值增长与资本及劳动之间的关系，进而利用已有的数据验证其结果.5.2.1模型的建立:在经济学的分析中，为了简化分析，通常假定生产中只有劳动和资本这两种生产要素.若以表示劳动投入量，以表示资本投入数量，则生产函数可以写为：生产函数表示生产中的投入量和产出量之间的依存关系，这种关系普遍存在于各种生产过程中.一家工厂必然具有一个生产函数，一家饭店也是如此，甚至一所学校或者医院同样会存在着各自的生产函数，产品可能是实实在在的有形产品，也可能是无形产品比如服务.估计和研究生产函数，对于经济理论和实践经验都具有一定意义.柯布—道格拉斯（Cobb-Dauglas）生产函数是由数学家柯布和经济学家道格拉斯于20世纪30年代初一起提出来的.柯布—道格拉斯生产函数被认为是一种很有用的生产函数，因为该函数以极简单的形式描述了经济学家所关心的一些性质，它在经济理论的分析和实证研究中都有具有一定意义.柯布—道格拉斯生产函数的函数表达式如下：其中，代表产出量,代表资本投入量,代表劳动投入量,、、为未知参数.表示技术或管理等参数对经济增长的影响系数,和分别表示劳动和资本对产出的贡献程度,且.对该生产函数取对数得:由于柯布―道格拉斯生产函数假设技术、管理水平不变,即是一个常数,在此可以忽略的影响.所以，可简化为:求出道格拉斯函数以后，我们通过道格拉斯函数可以预测出来每一年的GDP总产值，然后利用题目所提供的数据进行检验，可以发现道格拉斯很好的表示出来了省内生产总值的增长和投资与劳动之间的关系.5.2.2模型的求解:为了求解上述模型，通过分析题目所给数据和从网上查找相关数据，我们列出了我省从1981年到2003的GDP总值，投资金额总和和我省就业人数的表2如下：表2年份GDP值投资金额就业人数1981346.5779.602885.41982395.3885.002948.61983459.8396.463029.71984581.56140.153140.41985680.46194.333276.81986742.05223.083425.01987892.29297.773561.519881117.66369.823680.119891293.94305.543803.419901511.19335.663927.419911810.54439.824219.319922196.53601.504302.619932770.37892.484379.319943844.51108.004382.119954953.351320.975207.419965883.81558.015227.919976537.071792.225256.019987021.352056.975287.619997493.842222.175314.720008337.472542.655441.820019195.042807.795475.3200210275.53509.295527.0200312078.155328.445620.62010GDP:亿元；投资金额：亿元；就业人数：万人；在实际生产中，人们关心的往往是生产的增长量，而不是绝对量，因此定义投资金额指数，就业人数指数和省内生产总值指数分别为利用上述定义的三个指数公式，通过使用matlab软件计算出表3中1981年到2003年的省内生产总值指数，投资金额指数，和就业人数指数的一组数据，取1990年为基年,则t=0.t-90.21270.69690.2612-80.27240.71600.2848-70.31660.73030.3198-60.40570.74410.3870-50.56300.77030.4830-40.69080.81250.5489-30.83960.79200.6438-21.03020.83910.8033-10.97630.85450.9083011111.23841.01151.16621.78861.02171.439132.89371.03181.883743.94621.04182.570154.54341.05123.195365.17071.06493.747375.59151.07834.170386.35681.09094.435596.58651.10264.7269107.32941.11335.2355118.40891.12785.7734129.97161.13896.36251312.97541.14957.2215表3从表中可知，在正常的经济发展过程中（除个别年份外），上述三个指标都是随时间增长的，但是很难直接从表中发现具体的经济规律.为了定量分析，定义两个新的变量根据表中数据，在直角坐标系上做出的散点图，发现基本上成正比例关系（散点位于一条直线的附近），如图5图5散点图我们可以用MATLAB软件中的curvefit()函数来作数据拟合，即寻求函数Q(K，L)中的未知参数,,，使这个函数尽量逼近表5-2-2所给出的统计数据.则可以得到：,,.于是公式变为：这就是产值随资金K、劳动力L的变化规律.5.2.3模型的检验：为了对所建立的模型进行检验，我们利用得出的道格拉斯函数对每年的GDP指数做出了预测，结果如下表4所示，然后利用已有的GDP指数进行比较，最后得出所建立的道格拉斯函数是有意义的，可以正确表示出省内生产总值与资金与劳动之间的关系.t预测值实际值-90.27610.2612-80.33620.2848-70.37930.3198-60.46150.3870-50.59910.4830-40.70610.5489-30.82620.6438-20.97420.8033-10.93820.908300.9889111.16971.166021.55801.439132.26551.883742.88622.570153.22503.195363.57453.747373.80794.170384.21594.435594.34334.7269104.72765.2355115.27275.7734126.02896.3625137.40637.2215表4为了形象的表示出预测值与实际值之间的关系，我们做出了下图6，通过图表可以发现，道格拉斯函数已经可以很精确的表示来省内生产总值的变化趋势：图6省内生产总值的预测与实际比较图增加生产、发展经济所依靠的主要因素有增加投资、增加劳动力以及技术革新等，在研究省民经济产值与这些因素的数量关系时，由于技术水平不像资金、劳动力那样容易定量化，作为初步的模型，可认为技术水平不变，只讨论产值和资金、劳动力之间的关系.在科学发展不快时，如资本主义经济发展的前期，这种模型是有意义的.从而可以说明省内生产总值增长与资本及劳动之间满足柯布—道格拉斯（Cobb-Dauglas）生产函数的关系.模型的评价与推广6.1模型的优点：在问题一中，多元回归模型，因变量省内生产总值的99.89%可由模型确定，说明模型从整体上来看是可用的.在问题二中，运用了柯布—道格拉斯生产函数，使该模型的建立有理论依据作支撑，且有助于对模型的结果进行分析.在分析省内生产总值与投资和劳动力关系是，忽略其他因素，从而简化了模型，便于大概的预测.6.2模型的缺点：问题一中，由于省内生产总值受国际经济、政府政策、自然灾害等因素的影响，所以某一时期GDP波动幅度较大，因此影响了模型整体预测的准确性.问题二中，忽略其他因素对省内生产总值的影响，和实际问题存在的误差.一定历史时期的生产函数是反映当时的社会生产力水平的.6.3模型的推广与改进：推广：模型一是一类基于统计分析的随机模型，因此适用于大量数据的随机现象.如经济增长，灾害预测等.模型二中，在信息经济时代，所投入的生产要素的核心成分从资本、劳动力逐渐转变为以信息技术为代表的高新技术.当信息资源应用于生产中时，对生产人员、资本、流程等形成革命性的影响作用，极大地提高了生产要素生产率，促进了经济发展.综合上述原因，需要对柯布——道格拉斯生产函数做出了一定的修正，使之适用于信息时代的生产力发展水平.改进：模型一中参数与的置信区间包含零点，说明模型中还存在缺点，变量之间很可能存在交互作用.因此应在模型中加入交互项，改进原有的模型，建立新的回归模型.模型二较原来的模型增加了信息技术设备的资本投入和信息技术的劳动力投入后，得到使得模型成为更贴近时代的生产模型，改进后的柯布—道格拉斯生产函数是在现代信息工业经济时代构造出的反映了现代信息工业经济时代生产力特征的函数模型.改进后的柯布—道格拉斯生产函数模型更具有时代特色，适用性更广.参考文献[1]姜启源，谢金星，数学模型，高等教育出版社，2007.[2]韩中庚，数学建模方法及其应用，高等教育出版社，2006.[3]周品，赵新芬，MATLAB数学建模与仿真，2009.[4]王兵团，数学建模基础，2005.[5]齐微，柯布—道格拉斯生产函数模型，中国科技论文在线.附录附表1：2-11952-2010年地区生产总值单位:亿元年份地区

生产总值第一产业第二产业第三产业人均地区

生产总值

(元)工业建筑业#交通运输仓储

邮电通信业#批发零售

贸易餐饮业1981346.57132.21155.41138.0917.3258.9513.1816.564721982395.38154.07166.05147.1018.9575.2614.1022.395311983459.83185.57178.75159.1519.6095.5117.1232.076111984581.56222.13239.27214.2025.07120.1622.2537.977651985680.46235.96293.07259.4233.65151.4328.2246.148871986742.05252.73313.21274.8038.41176.1134.0349.859561987892.29287.31384.57341.3143.26220.4147.4560.58113119881117.66331.94497.10435.5161.59288.6255.5086.71139519891293.94359.14579.65513.9765.68355.1566.55108.51159519901511.19425.29635.98568.2567.73449.9281.87132.83181519911810.54521.85745.90663.9082.00542.7998.97161.00212219922196.53534.62999.11889.59109.52662.80120.18200.93255619932770.37596.631355.711201.67154.04818.03141.09239.98321219943844.50775.031891.431692.10199.331178.04213.58346.88444119954953.351010.132355.782098.06257.731587.44296.08471.95570119965883.801200.172784.092475.99308.101899.54362.05574.97674619976537.071195.003147.372796.02351.352194.70420.49664.92746119987021.351215.813408.063008.45399.612397.49437.37735.63796819997493.841221.003644.323197.16447.162628.52484.65787.90848320008337.471268.574164.453665.74498.712904.45545.13856.94932620019195.041359.494556.014004.09551.923279.53657.57972.3310195200210275.501390.005184.984518.87666.113700.52655.641142.5411340200312078.151480.676485.055706.71778.344112.43710.181283.7013268200415021.841778.458478.697576.12902.574764.70969.131431.5816413200518366.871963.5110478.629418.581060.045924.7419934200621900.192138.9012574.0311378.821195.217187.2623603200725776.912509.1414647.5313283.721363.818620.2427604200830933.283002.6517571.9815894.951677.0310358.6432936200933896.653226.6418901.8316896.142005.6911768.1835894201039169.923588.2821238.4918861.452377.0414343.1441106附表2：年份投资资金来源山东省预算省内贷款利用外资自筹和内资金其他资金总量(亿元)1981269.8122.036.453.91982279.3136.160.574.51983339.7145.566.684.31984421.0158.570.7108.71985407.8160.391.5133.61986455.6188.5107.3189.21987496.6192.0112.0221.11988432.0217.8145.3298.71989366.1223.0161.1290.31990393.0245.5194.6294.41991380.4254.7208.9350.41992347.5264.0218.7540.01993483.7272.0234.3562.41994529.6297.6269.0631.01995621.1298.7275.9709.21996625.9317.7276.6712.41997696.7352.6283.9726.519981197.4382.9297.0739.619991852.1455.9306.8759.720002109.5727.3396.3827.420012546.4729.8430.7830.020023161.0859.1485.0841.920032687.81044.4499.4854.820043254.91188.0585.7936.320053554.31319.0648.81138.720063672.01590.5684.31160.220074257.12044.2782.712769.720084524.82643.7791.913204.92010构成(%)198128.112.73.855.4198222.714.34.958.1198323.812.34.759.2198423.014.13.959.0198516.020.13.660.3198614.621.14.459.9198713.123.04.859.119889.321.05.963.819898.317.36.667.819908.719.66.365.419916.823.55.764.019924.327.45.862.519933.723.57.365.519943.022.49.964.719953.020.511.265.319962.719.611.866.019972.818.910.667.719984.219.39.167.419996.219.26.767.820006.420.35.168.220016.719.14.669.620027.019.74.668.720034.620.54.470.520044.418.54.472.720054.417.34.274.120063.916.53.676.020073.915.33.477.420084.314.52.978.3附表3：3-2历年人口出生率、死亡率、自然增长率年份出生率

(‰)死亡率

(‰)自然增长率

(‰)出生人口数

(万人)死亡人口数

(万人)自然增长人数

(万人)1981-16.48(6.41)(10.07)(121)(47)(74)1982-17.05(6.10)(10.95)(127)(45)(82)198315.10(12.76)6.73(5.87)8.37(6.89)114(96)51(44)63(52)198413.80(12.99)5.80(6.03)8.00(6.96)104(99)44(46)60(53)198515.12(11.75)6.64(5.90)8.48(5.85)116(90)51(45)65(45)198619.90(14.71)7.28(5.86)12.62(8.85)156(114)57(46)99(68)198723.35(17.43)7.07(5.64)16.28(11.79)184(137)56(44)128(93)198817.54(17.95)6.04(5.95)11.50(12.00)140(143)48(47)92(96)198916.88(18.87)5.70(5.51)11.18(13.36)137(153)46(45)91(108)199018.21(26.10)6.96(6.02)11.25(20.08)152(217)58(50)94(167)199115.40(16.39)6.54(5.73)8.86(10.66)131(139)56(49)75(90)199211.43(10.95)6.88(6.02)4.55(4.93)98(94)59(52)39(42)199310.49(9.47)6.76(5.84)3.73(3.63)90(81)58(50)32(31)19949.69(9.31)6.67(5.99)3.02(3.32)84(80)58(52)26(28)19959.82(9.66)6.47(5.83)3.35(3.83)85(84)56(51)29(33)199610.60(10.33)6.76(6.04)3.84(4.29)92(90)59(53)33(37)199711.28(10.84)6.65(5.90)4.63(4.94)99(95)58(52)41(43)199811.58(11.52)6.12(5.95)5.46(5.57)102(102)54(53)48(49)199911.08(10.23)6.27(5.72)4.81(4.51)98(91)55(51)43(40)200010.75(11.38)6.29(6.70)4.46(4.68)97(102)56(60)40(42)200111.12(9.93)6.24(5.46)4.88(4.47)100(89)56(49)44(40)200211.17(10.20)6.62(5.86)4.55(4.34)101(92)60(53)41(39)200311.42(9.31)6.64(6.07)4.78(3.24)104(85)61(55)43(30)200412.50(10.59)6.49(5.60)6.01(4.99)114(97)59(51)55(46)200512.14(10.17)6.31(5.85)5.83(4.32)112(94)58(54)54(40)200611.60(9.59)6.10(5.62)5.50(3.97)108(89)57(52)51(37)200711.11(10.05)6.11(6.47)5.00(3.58)104(94)57(60)47(33)200811.25(10.13)6.16(6.81)5.09(3.32)106((95)58(64)48(31)200911.70(10.96)6.08(6.11)5.62(4.86)110(103)57(58)53(46)201011.61(15.82)6.24(8.58)5.34(7.24)111(150)60(81)51(69)模型二计算程序：Q=[1117.661193.941511.191810.542196.532770.373844.504953.355883.806537.077021.357493.848337.47...9195.0410275.5012078.1515021.8418366.8721900.1925776.9130933.2833896.0539169.92];IQ=Q/11795.0K=[421.0407.8455.6496.6432.0366.1393.0380.4347.5483.7529.6621.1625.9696.71197.41852.1...2109.52546.43161.02687.83254.93554.33672.0];IK=K/1334.5L=[2885.42948.63029.73140.43276.83425.03561.53680.13803.43927.44219.34302.64379.34382.15207.4...5227.95256.05287.65314.75441.85475.35527.05620.6];IL=L/4318.3Et=zeros(23);Et=Et(1,1:23);fort=1:1:23;Et(t)=log(IL(t)/IK(t));endEt;Wt=Et;fort=1:1:23;Wt(t)=log(IQ(t)/IK(t));endWt;x=Et;y=Wt;plot(x,y,'*');xlabel('E');ylabel('W');a=[0.26120.28480.31980.38700.48300.54890.64380.8033...0.90831.00001.16601.43911.88372.57013.19533.7473...4.17034.43554.72695.23555.77346.36257.2215];y=[0.21270.27240.31660.40570.56300.69080.83931.0302...0.97631.00001.23841.78862.89373.94624.54345.1707...5.59156.35686.58657.32948.40899.971612.9754;0.69690.71600.73030.74410.77030.79200.81520.8391...0.85451.00001.01151.02171.03181.04181.05121.0649...1.07831.09091.10261.11331.12781.13891.1495];curvefun=inline('x(1)*(y(1,:).^x(2)).*(y(2,:).^x(3))','x','y')x0=[0.1,0.1,0.2];x=lsqcurvefit(curve