《经济增长数学管理与财务知识分析建模论文.docx》由会员分享,可在线阅读,更多相关《经济增长数学管理与财务知识分析建模论文.docx(22页珍藏版)》请在三一办公上搜索。
1、2011数学建模第二次模拟承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题.我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出.我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性.如有违反竞赛规则的行为,我们将受到严肃处理.我们参赛选择的题号是(从A/B/C/D中选择一项填写): 我们的参赛报名号为(如果赛区设置报名号的话):
2、所属学校(请填写完整的全名):河南理工大学万方科技学院 参赛队员 (打印并签名) :1. 关海超 2. 刘源 3. 冯艳伟 指导教师或指导教师组负责人 (打印并签名): 日期: 2011 年 8 月 21 日赛区评阅编号(由赛区组委会评阅前进行编号):2011数学建模第二次模拟编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):经济增长问题摘 要国内生产总值(GDP)常被公认为衡量国家经济状况的最佳指标.它不但可反映一个国家的经济发展情况,更可以反
3、映一国的国力与财富.因此分析各产业对于GDP的影响,并研究GDP的增长规律是具有现实意义的.在问题一中,我们分别做出了GDP与工业、建筑业及农林渔业产值关系的散点图,分析得出GDP的值与各产业之间存在明显的线性关系. 回归分析是统计分析的重要组成部分,用回归分析方法来研究自变量与因变量的关系函数是一种常用的有效方法.因此我们建立起了多元线性回归模型,用MATLAB计算得到的模型为.在对该模型进行显著性检验中,我们对各参数进行了显著性分析,得到模型的复相关系数=0.999,统计量=30900. 统计量的值远超过检验的临界值,因此可以验证模型是可用的.最后,我们利用所建立的模型对20102014年
4、的GDP值做出了预测,分析了各产业对GDP的影响.通过处理预测的数据,我们得出平均每年GDP的增长率为10%左右,其中建筑业与工业对GDP的影响较大,而农林渔业对GDP的影响较小,这也符合中国的产业结构与经济发展情况.在问题二中,为了讨论国内生产总值增长与资本及劳动之间的关系,我们通过分析数据、查阅相关资料,了解到了国内生产总值的大小通常取决于相关的生产资料和劳动力等相关重要因素.于是,我们通过建立柯布道格拉斯生产函数,定义了三个指数分别为:投资金额指数,就业人数指数和国内生产总值指数.利用定义的三个指数公式,计算出1981年到2008年的国内生产总值指数,投资金额指数和就业人数指数的一组数据
5、,并探讨国内生产总值增长与资本及劳动之间的关系.但是上述三个指标都是随时间增长的,很难直接从表中发现具体的经济规律.为了定量分析,我们定义两个新的变量分别,通过做散点图发现这两个变量基本上成正比例关系.我们用MATLAB软件中的curvefit()函数来作数据拟合,求得函数中的未知参数,,通过检验进而得处道格拉斯生产函数为,这就是产值随资金K、劳动力L的变化规律.为了验证第二问中的结果,我们用求得的道格拉斯生产函数来预测每年的国内生产总值,然后与题目提供的数据进行比较来进行检验.通过检验可以发现预测值的误差很小,因此道格拉斯生产函数可以表示出国内生产总值增长与资本及劳动之间的关系.关键字:多元
6、线性回归 显著性检验 道格拉斯函数 数据拟合1. 问题重述国内生产总值(Gross Domestic Product,简称GDP)是指在一定时期内(一个季度或一年),一个国家或地区的经济中所生产出的全部最终产品和服务的价值,常被公认为衡量国家经济状况的最佳指标.它不但可反映一个国家的经济表现,更可以反映一国的国力与财富.(1) 建立国内生产总值与工业值、建筑业及农林渔业产值之间的数量模型,利用数据对未来经济做出预测;(2) 讨论国内生产总值增长与资本及劳动之间的关系,利用数据验证其结果.2. 问题分析2.1 建立GDP与工业值、建筑业及农林渔业产值之间的数量模型,并对未来经济做出预测.在问题一
7、中,我们通过分析材料得出这是研究对象的内在特性和各个因素间关系的问题,即研究GDP与工业值、建筑业及农林渔业产值关系.一般用机理分析的方法建立数学模型.由于经济问题是一种随机的问题,所以通常的办法是搜集大量的数据,基于对数据的统计分析去建立模型.因为影响GDP的因素有三个,即工业值、建筑业及农林渔业产值,且各个产业与GDP都为线性关系.所以我们建立起一个多元线性回归模型,并检验模型显著性,通过对模型的反复修改与检验,建立更合理的模型.2.2 讨论国内生产总值增长与资本及劳动之间的关系,并验证其结果.在问题二中,为了讨论国内生产总值增长与资本及劳动之间的关系,通过查阅相关资料,我们了解到国内生产
8、总值通常取决于相关的生产资料和劳动力等相关重要因素. 要建立道格拉斯生产函数,我们只需要讨论产值和资金,劳动之间的关系,从而达到我们的目的.这样处理不仅能简化问题,而且是合理的在生产产值上的预测,柯布-道格拉斯(Cobb-Douglas)生产函数预测的结果近似就是准确生产值.于是我们通过建立柯布道格拉斯生产函数,来探讨国内生产总值增长与资本及劳动之间的关系,进而利用已有的数据验证其结果.3. 模型假设1. 假设所统计的数据都在误差允许的范围之内;2. 忽略由于非正常条件下的引起的数据的巨大波动;3. 假设在短期内国内生产总值只取决于投资和劳动力因素;4. 假定在一段不太长的时间内技术水平不变.
9、4. 定义与符号说明国内生产总值复相关系数统计量随机误差置信水平工业产值建筑业产值农林渔业产值国内生产总值增长劳动投入量资本投入数量劳动对产出的贡献程度资本对产出的贡献程度投资金额指数就业人数指数国内生产总值指数道格拉斯函数常数5. 模型的建立与求解5.1 问题一模型的建立与求解:回归分析方法是统计分析的重要组成部分,用回归分析方法来研究自变量与变量的关系函数是一种常用的有效方法.我们通过回归模型的建立,定量预测了未来经济的发展.5.1.1 GDP与工业值、建筑业及农林渔业产值数量模型: 通过在互联网上搜集到1978年2009年,中国GDP与工业值、建筑业及农林渔业产值的数据(见附表1),可以
10、定性的看出GDP与工业、建筑业及农林渔业产值为整体上升的趋势.为了大致分析GDP与工业值、建筑业及农林渔业产值关系,我们首先利用附录数据做出了GDP与工业产值的关系散点图(如图1).图1 工业产值与GDP散点图 从图可以发现,随着工业产值的增加,GDP的值有比较明显的线性增长趋势.图中的直线是用线性模型 拟合的.同理我们也分别作出了建筑业产值与GDP的关系散点图(图2)、农林渔业产值与GDP的关系散点图(图3).图2 建筑业产值与GDP散点图图3 农林渔业产值与GDP散点图通过图2,图3可以看出建筑业产值与农林渔业产值同样有很强的线性关系,同样也分别用直线模型对其拟合.建筑业产值与GDP线性模
11、型农林渔业产值与GDP线性模型因此,综上所述四者之间有很强的线性关系,可建立多元线性回归模型在模型中除了工业,建筑业,农林渔业外,影响国内生产总值的其他因素的作用都包含在随机误差内,这里假设相互独立,且服从均值为零的正态分布,.对模型直接利用matlab统计工具箱求解,得到回归系数估计值及其置信区间(置信水平= 0.05),检验统计量,的结果见表1.参数参数估计置信区间732.3-589.6 2054.21.91.7 2.04.02.9 5.10.03-0.03 0.4 表1 模型的计算结果5.1.2 结果分析: 表1显示,指因变量(国内生产总值)的99.89%可由模型确定,值远远超过检验的临
12、界值,p远小于,因而模型从整体上来看是可用的.表1的回归系数给出了模型中,的估计值,即,.检查它们的置信区间发现与的置信区间都包含零点,这表明回归模型常数与回归变量对模型的影响不太显著.这也符合这一事实,农林渔业产值对GDP的影响较小,工业与建筑业对GDP的影响较大.但是一般情况下常数的值都保留在模型中,不剔除.回归变量系数区间右端点相对距零点较远,所以我们也保留在模型中.因此最终确定的模型为:5.1.3 未来经济预测将回归系数的估计值代入模型,即可预测未来的GDP情况.代入得到的模型即 只要能预测未来的工业值、建筑业及农林渔业产值即能预测未来GDP的产值.根据统计数据(见附表1),我们用ma
13、tlab计算得到各产值的平均年增长率的中位数,工业值增长率13.5%,建筑业产值增长率16.8%,农林渔业产值10.7%.我们对GDP的值做短期的预测,预测未来五年的GDP情况见下图4.图4 2010-2014中国GDP情况预测从模型中看,工业值、建筑业及农林渔业产值的预测直接影响GDP预测的准确性.其中工业与建筑业对预测的影响较大,农林渔业影响较小.从预测结果中看中国GDP总值成上升趋势,从数值中计算得平均值为12.3%.当然GDP的增长率的计算还应除去通货膨胀、消费指数等因素的影响,所以实际中应小一些.2010年的GDP为39.789万亿,这也与预测结果相符合,说明模型的合理性.5.2 问
14、题二模型的建立与求解:在问题二中为了讨论国内生产总值增长与资本及劳动之间的关系,通过查阅相关资料,我们了解到国内生产总值通常取决于相关的生产资料和劳动力等相关重要因素.于是我们通过建立柯布道格拉斯生产函数,来探讨国内生产总值增长与资本及劳动之间的关系,进而利用已有的数据验证其结果.5.2.1 模型的建立:在经济学的分析中,为了简化分析,通常假定生产中只有劳动和资本这两种生产要素.若以表示劳动投入量,以表示资本投入数量,则生产函数可以写为:生产函数表示生产中的投入量和产出量之间的依存关系,这种关系普遍存在于各种生产过程中.一家工厂必然具有一个生产函数,一家饭店也是如此,甚至一所学校或者医院同样会
15、存在着各自的生产函数,产品可能是实实在在的有形产品,也可能是无形产品比如服务.估计和研究生产函数,对于经济理论和实践经验都具有一定意义.柯布道格拉斯(Cobb-Dauglas)生产函数是由数学家柯布和经济学家道格拉斯于20世纪30年代初一起提出来的.柯布道格拉斯生产函数被认为是一种很有用的生产函数,因为该函数以极简单的形式描述了经济学家所关心的一些性质,它在经济理论的分析和实证研究中都有具有一定意义.柯布道格拉斯生产函数的函数表达式如下:其中,代表产出量,代表资本投入量,代表劳动投入量,、为未知参数.表示技术或管理等参数对经济增长的影响系数, 和分别表示劳动和资本对产出的贡献程度,且.对该生产
16、函数取对数得:由于柯布道格拉斯生产函数假设技术、管理水平不变,即是一个常数, 在此可以忽略的影响.所以,可简化为:求出道格拉斯函数以后,我们通过道格拉斯函数可以预测出来每一年的GDP总产值,然后利用题目所提供的数据进行检验,可以发现道格拉斯很好的表示出来了国内生产总值的增长和投资与劳动之间的关系. 5.2.2 模型的求解:为了求解上述模型,通过分析题目所给数据和从网上查找相关数据,我们列出了我国从1981年到2003的GDP总值,投资金额总和和我国就业人数的表2如下:年份GDP总值投资金额就业人数19814889.5961.04.512619825330.51230.44.6358198359
17、85.61430.14.728619847243.81832.94.817919859040.72543.24.9873198610274.43120.65.1282198712050.63791.75.2783198815036.84653.85.4334198917000.94410.45.5329199018718.34517.56.4749199121826.25594.56.5491199226937.38080.16.6152199335260.013072.36.6808199448108.517827.16.7455199559810.520524.96.80651996701
18、42.523358.66.8951199778060.825259.76.9821199883024.328716.97.0637199988479.229754.67.1394200098000.533110.47.20852001108068.237987.07.30252002119095.745046.97.37412003135174.058616.37.4432表2GDP:亿元;投资金额:亿元;就业人数:亿人;在实际生产中,人们关心的往往是生产的增长量,而不是绝对量,因此定义投资金额指数,就业人数指数和国内生产总值指数分别为利用上述定义的三个指数公式,通过使用matlab软件计算出
19、表3中1981年到2008年的国内生产总值指数,投资金额指数,和就业人数指数的一组数据,取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.8837
20、43.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从表中可知,在正常的经济发展过程中(除个别年份外),上述三个指标都是随时间增长的,但是很难直接从表中发现具体的经济规律.为了定量分析,定义两个新的变量 根据表中数据,在直角坐标系上做出的散点图,发现基本上成
21、正比例关系(散点位于一条直线的附近),如图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
22、-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为了形象的表示出预测值与实际值之间的关
23、系,我们做出了下图6,通过图表可以发现,道格拉斯函数已经可以很精确的表示来国内生产总值的变化趋势:图6 国内生产总值的预测与实际比较图增加生产、发展经济所依靠的主要因素有增加投资、增加劳动力以及技术革新等,在研究国民经济产值与这些因素的数量关系时,由于技术水平不像资金、劳动力那样容易定量化,作为初步的模型,可认为技术水平不变,只讨论产值和资金、劳动力之间的关系.在科学发展不快时,如资本主义经济发展的前期,这种模型是有意义的.从而可以说明国内生产总值增长与资本及劳动之间满足柯布道格拉斯(Cobb-Dauglas)生产函数的关系.6. 模型的评价与推广6.1 模型的优点:在问题一中,多元回归模型,
24、因变量国内生产总值的99.89%可由模型确定,说明模型从整体上来看是可用的.在预测2010-2014年的GDP的值时.我们计算得中国平均年GDP的增长量为10%左右,这也完全符合中国的经济发展情况.在问题二中,运用了柯布道格拉斯生产函数,使该模型的建立有理论依据作支撑,且有助于对模型的结果进行分析.在分析国内生产总值与投资和劳动力关系是,忽略其他因素,从而简化了模型,便于大概的预测.6.2 模型的缺点:问题一中,由于国内生产总值受国际经济、政府政策、自然灾害等因素的影响,所以某一时期GDP波动幅度较大,因此影响了模型整体预测的准确性.问题二中,忽略其他因素对国内生产总值的影响,和实际问题存在的
25、误差.一定历史时期的生产函数是反映当时的社会生产力水平的.6.3 模型的推广与改进:推广:模型一是一类基于统计分析的随机模型,因此适用于大量数据的随机现象.如经济增长,灾害预测等.模型二中,在信息经济时代,所投入的生产要素的核心成分从资本、劳动力逐渐转变为以信息技术为代表的高新技术.当信息资源应用于生产中时,对生产人员、资本、流程等形成革命性的影响作用,极大地提高了生产要素生产率,促进了经济发展.综合上述原因,需要对柯布道格拉斯生产函数做出了一定的修正,使之适用于信息时代的生产力发展水平.改进:模型一中参数与的置信区间包含零点,说明模型中还存在缺点,变量之间很可能存在交互作用.因此应在模型中加
26、入交互项,改进原有的模型,建立新的回归模型.模型二较原来的模型增加了信息技术设备的资本投入和信息技术的劳动力投入后,得到使得模型成为更贴近时代的生产模型,改进后的柯布道格拉斯生产函数是在现代信息工业经济时代构造出的反映了现代信息工业经济时代生产力特征的函数模型.改进后的柯布道格拉斯生产函数模型更具有时代特色,适用性更广.7. 参考文献1 姜启源,谢金星,数学模型,高等教育出版社,2007.2 韩中庚,数学建模方法及其应用,高等教育出版社,2006.3 周品,赵新芬,MATLAB数学建模与仿真,2009.4 王兵团,数学建模基础,2005.5 齐微,柯布道格拉斯生产函数模型,中国科技论文在线.8
27、. 附录附表1:年份GDP值工业建筑业农林渔业19783645.21607.0138.21027.519794062.61769.7143.81270.219804545.61996.5195.51371.619814889.52048.4207.11559.519825330.52162.3220.71777.419835985.62375.6270.61978.419847243.82789.0316.72316.119859040.73448.7417.92564.4198610274.43967.0525.72788.7198712050.64585.8665.83233.019881
28、5036.85777.2810.03865.4198917000.96484.0794.04265.9199018718.36858.0859.45062.0199121826.28087.11015.15342.2199226937.310284.51415.05866.6199335260.014188.02266.56963.8199448108.519480.72964.79572.7199559810.524950.63728.812135.8199670142.529447.64387.414015.4199778060.832921.44621.614441.9199883024
29、.334018.44985.814817.6199988479.235861.55172.114770.0200098000.540033.65522.314944.72001108068.243580.65931.715781.32002119095.747431.36465.516537.02003135174.054945.57490.817381.72004159586.765210.08694.321412.72005185808.677230.810367.322420.02006217522.791310.912408.624040.02007267763.7110534.915
30、296.528627.02008316228.8130260.218743.233702.02009343464.7135239.922398.835226.0附表2:年 份投资资金来源国家预算国内贷款利用外资自 筹 和内 资 金其他资金总量 (亿元)1981269.8122.036.4532.91982279.3176.160.5714.51983339.7175.566.6848.31984421.0258.570.71082.71985407.8510.391.51533.61986455.6658.5137.31869.21987496.6872.0182.02241.11988432
31、.0977.8275.32968.71989366.1763.0291.12990.31990393.0885.5284.62954.41991380.41314.7318.93580.41992347.52214.0468.75050.01993483.73072.0954.38562.41994529.63997.61769.011531.01995621.14198.72295.913409.21996(629.7)(4576.5)(2747.4)(15465.4)625.94573.72746.615412.41997696.74782.62683.917096.519981197.4
32、5542.92617.019359.619991852.15725.92006.820169.720002109.56727.31696.322577.420012546.47239.81730.726470.020023161.08859.12085.030941.920032687.812044.42599.441284.820043254.913788.03285.754236.320054154.316319.03978.870138.720064672.019590.54334.390360.220075857.123044.25132.7116769.720087954.82644
33、3.75311.9143204.9构成(%)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.3
34、19962.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:人口出生率、死亡率和自然增长率单位:年 份出生率死亡率自然增长率年 份出生率死亡率自然增长率197818.256.25
35、12.00199517.126.5710.55198018.216.3411.87199616.986.5610.42198120.916.3614.55199716.576.5110.06198222.286.6015.68199815.646.509.14198320.196.9013.29199914.646.468.18198419.906.8213.08200014.036.457.58198521.046.7814.26200113.386.436.95198622.436.8615.57200212.866.416.45198723.336.7216.61200312.416.4
36、06.01198822.376.6415.73200412.296.425.87198921.586.5415.04200512.406.515.89199021.066.6714.39200612.096.815.28199119.686.7012.98200712.106.935.17199218.246.6411.60200812.147.065.08199318.096.6411.45200912.137.085.05199417.706.4911.21模型二计算程序:Q=4889.5 5330.5 5985.6 7243.8 9040.7 10274.4 12050.6 15036.
37、8 17000.9 18718.3 21826.2 26937.3 35260.0 . 48108.5 59810.5 70142.5 78060.8 83024.3 88479.2 98000.5 108068.2 119095.7 135174.0;IQ=Q/18718.3K=961 1230.4 1430.1 1832.9 2543.2 3120.6 3791.7 4653.8 4410.4 4517.5 5594.5 8080.1 13072.3 17827.1 20524.9 . 23358.6 25259.7 28716.9 29754.6 33110.4 37987.0 4504
38、6.9 58616.3;IK=K/4517.5L=4.5126 4.6358 4.7286 4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749 6.5491 6.6152 6.6808 6.7455. 6.8065 6.8951 6.9821 7.0637 7.1394 7.2085 7.3025 7.3741 7.4432;IL=L/6.4749Et=zeros(23);Et=Et(1,1:23);for t=1:1:23; Et(t)=log(IL(t)/IK(t);endEt;Wt=Et;for t=1:1:23; Wt(t)=log(IQ(
39、t)/IK(t);endWt;x=Et;y=Wt;plot(x,y,*);xlabel(E);ylabel(W);a=0.2612 0.2848 0.3198 0.3870 0.4830 0.5489 0.6438 0.8033. 0.9083 1.0000 1.1660 1.4391 1.8837 2.5701 3.1953 3.7473. 4.1703 4.4355 4.7269 5.2355 5.7734 6.3625 7.2215;y=0.2127 0.2724 0.3166 0.4057 0.5630 0.6908 0.8393 1.0302. 0.9763 1.0000 1.2384 1.7886 2.8937 3.9462 4.5434 5.1707. 5.5915 6.3568 6.5865