数学建模B题优秀论文.pdf
2010高教社杯全国大学生数学建模竞赛 承 诺 书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。我们参赛选择的题号是(从A/B/C/D 中选择一项填写):B 我们的参赛报名号为(如果赛区设置报名号的话):所属学校(请填写完整的全名):参赛队员(打印并签名):1.王静茹 2.杨曼 3.朱元霞 指导教师或指导教师组负责人 (打印并签名):日期:年 月 日 赛区评阅编号(由赛区组委会评阅前进行编号):2010 高教社杯全国大学生数学建模竞赛 编 号 专 用 页 赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评 阅 人 评 分 备 注 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):2010 年上海世博会经济影响力的定量评估 摘要 本文选取 2010 年上海世博会对上海经济的影响作为研究对象,首先,我们选择了五届影响力较大的世博会与上海世博会进行了定量的纵向评估。利用互联网的相关数据,运用层次分析法确定了各级评价指标的相对权重,然后利用模糊综合评判法给这六届世博会的经济影响力进行了定量评估,利用MATLAB 计算出了 1933 年芝加哥世博会以来六届综合性世博会的经济影响力的综合评分依次为75.12、80.01、80、11、77.35、79.35、80.75,由表我们可以肯定上海世博会的经济影响力是继 1851 年伦敦世博会以来较强的。其次我们采用投入产出模型模型的核心思想,以年份与GDP的对数值的二次相关关系和上海市社会固定资产总投入与GDP的对数值的线性关系,利用上海统计年鉴发布的数据,分别建立无世博影响的表达式iiixxxeQ21210904.01117.00032.06278.81,与有世博影响的表达式iiixxxeQ21212955.00176.00019.01211.82,两式的预测误差均在 1.1%以内。与2008 年真实值比较,用表达式1Q预测 2008 年的GDP的值可以得出世博会对 2008 年上海市经济贡献率达到 20.9%。并且在得知申办世博会后第i年上海市固定投入总额的前提下由%100212QQQ可求出世博会对上海地区经济的持续性积极影响。如假设2011 年市固定资产总投资为 5600 亿元,则世博会对上海经济有16%的积极影响。最后,经过对 2010 年上海世博会的经济影响力的两方面的评估,我们得知上海世博会在历届世博会的经济影响力的综合评分中是最高的。由此得出,上海世博会对上海经济的影响力是非常大的,此次世博会除了对上海的直接收益影响明显外,世博会对上海地区经济的持续性积极影响。关键词:层次分析 模糊综合评判 投入产出模型 回归模型 一、问题重述 2010 年上海世博会是首次在中国举办的世界博览会。从 1851 年伦敦的“万国工业博览会”开始,世博会正日益成为各国人民交流历史文化、展示科技成果、体现合作精神、展望未来发展等的重要舞台。请你们选择感兴趣的某个侧面,建立数学模型,利用互联网数据,定量评估 2010 年上海世博会的影响力。二、问题分析 本文选择了经济影响力这一侧面。通过上海世博会与前六届世博会对经济的影响进行对比,应用层次分析法和模糊数学评判法来评价上海世博会对经济的影响。“影响力”仅仅是一个宽泛的概念,因此在实际生活中需要将影响力转化为有参考价值的数据时,往往需要利用数学建模的方法对其进行定性与定量的分析。本文在评估 2010年上海世博会的影响力时,正是应用了这种方法。于是建立模型二,是通过举办世博会所带来的GDP与不举办世博会GDP进行对比,建立一套投入产出模型,得出影响力的一个参考系数。考察世博会自筹办以来对上海经济增加的拉动作用,以及对上海经济发展的促进作用。世博效应的作用使得世博会结束后可能出现经济回落幅度较小,而后续效应时间周期很长,至少有 10 到 50 年,之后还将影响很长时间,因此我们将世博会对上海的影响分为直接影响和间接影响分别进行评判。三、模型假设 假设确定的各级评价指标的各个因素是合理的且并不相互影响;假设专家对每个因素的评价是合乎实际的;假设在构造成对比较矩阵时对各因素的权重赋值是合理的;四、符号说明 1,2ii表示两个主要因素的分类指标的权重向量;1,2iR i表示第i因素的模糊评判矩阵;1,2iB i表示一级评判矩阵的元素(1,2)iZ i表示第i因素的综合评分;Z表示世博总体的综合评分;,A B分别表示两主要影响因素的成对比较矩阵;,AB分别表示两成对比较矩阵的最大特征根;p表示两主要因素的权重相对于目标层的权重向量;nx1表示从 1997 年起第n年 nx2为从 1997 年起第n年的上海市固定投入总量 1Q从1997 年至 2002 年上海市GDP总量;2Q从 2003 年至 2008 年上海市GDP总量 五、模型建立和求解 4.1 纵向比较历届世博会对经济的影响力 世博会本身是用来展示人类文明成就的,但对于主办国或主办城市来说,争办世博会是有其明确的自身目标利益,不同的主办国与主办城市有不同的目标利益,但有一点是共同的,通过举办世博会赋予主办国和主办城市新的价值、新的地位。世博会的举办将在经济、文化、科技、政治地位等方面为其带来巨大的影响。下面将利用层次分析法以及模糊评判法纵向比较历届世博会,从而得出上海世博会对经济影响力。4.2 模型的建立 4.2.1 运用层次分析法确定评价指标权重.层次分析法是把专家的经验认识与理性认识结合起来,以科学的解决确定加权系数的问题。首先我们逐一判断世博会每一层次上指标的相对重要程度,并将两两比较判断的结果按给定的比率标度定量化,从而构成成对比较矩阵,通过计算矩阵的最大特征值极其相应的特征向量,最终得出该层次指标权重系数。我们将决策问题分解为3个层次,最上层为目标层,即为世博会经济影响力的综合定量评估,最下层为方案层,1,26iP i表示六届世博会经济影响力的综合评分的最终排序,中间层为准则层,有直接影响和间接影响两个一级准则,下设旅游经济、洒店经济投资融资、就业率、提长升国际地位、带动区域经济、产业结构调整七个二级准则,具体关系如图1 图 1 上海世博会经济影响力综合评估的层次结构图 (1)构建成对比较矩阵。成对比较矩阵元素的数值反映了各元素的相对重要程度,采用1-10及其倒数来进行标度(见表1)。表 1:评价尺度赋值表 重要性 同等 重要 过渡值 稍微 重要 过渡值 明显 重要 过渡值 重要 得多 过渡值 绝对 重要 过渡值 赋值 1 2 3 4 5 6 7 8 9 10(2)采用MATLAB编程的方法来求解所得比较矩阵的最大特征根及其对应当的归一特征 向量。具体过程见程序 1 和程序 2(3)进行一致性检验 一致性指标CI可由1nCIn得出;平均随机性指标RI可由表 2 得出。表2:多阶判断矩阵的RI值 阶数 1 2 3 4 5 6 7 8 9 RI 0.00 0.00 0.58 0.90 1.12 1.24 1.32 1.41 1.45 当随机一致性比率0.10CICRRI时,认为层次分析法的结果有满意的一致性,即权重的分配是合理的。否则,要调整判断矩阵的元素取值,重新分配权重的值。4.2.2 运用模糊分析法对六届世博会进行综合评价.模糊综合评判的一般提法:设12,nUu uu为研究对象的n因素(或指标),称之为因素集(或指标集)。12,mVv vv为诸因素(或指标)的种评判所构成的评判集(或称评语集,评价集,决策集等),它们的元素个数和名称均可根据实际问题的需要和决策人主观确定.实际中,很多问题的因素评判集都是模糊的,因此,综合评判应该是V上的一个模糊子集 12,()mBb bbF V,其中kb为评判kv对模糊子集B的隶属度:()(1,2,)Bkkvb km,即反映了第k种评判kv在综合评价中所起的作用.综合评判B依赖于各因素的权重,即它应该是U上的模糊子集12,()nAa aaF U,且11niia,其中ia表示第i种因素的权重.于是,当权重A给定以后,则相应地就可以给定一个综合评判B.依据层次分析法得出的各级指标权重,建立历届世博会模糊定量评价尺度模型。一级指标集12,BB B,相应的权重集为12,bb b,12,B B分别代表直接影响和间接影响。二级指标集127,CC CC,相应的权重集为127,cc cc,127,c cc分别表示旅游经济,酒店经济产业结构调整。世博会经济影响力综合定量评估 直接影间接影旅游经济酒店经济投资融资交通经济提升 国际 地位 带动区域 经济增产业结构定义评语集为12345,UU U U U U 非常好,较好,好,一般,不太好,对其赋值为:95,85,75,65,55U。利用专家或群众对世博会融入举办地经济社会的情况进行评价,得到从二级指标到评语集的模糊矩阵。根据评语集划分的五个等级,ijr=评价第i个因素为第j个等级的人数/评委的总人数,且其表示第i个二级因素第j个等级的隶属度。于是可以得到二个主要因素的模糊评判矩阵12,R R;根据层次分析法确定权重向量:12,确定一级评价指标第四步:求一级评判矩阵12BKB 111*BwR,222*BwR,各主要因素的权重B 相对于目标层的权重向量为:p 由此便可得评判矩阵 设分数值95,85,75,65,55U;则可得到两主要影响因素的综合评分为:/iiijZBUB,1,2i 其中总体综合评分为:4.3 模型的求解 在构造上海世博会综合评分的成对比较矩阵和确定各评价指标的相对重要程度时,我们综合分析了一个世纪来影响力较大的六届世博会(如表 3)的相关数据,同时权衡了世博会本身各评价指标的相对重要程度,然后确定成对比较矩阵。表 3:历届世博会相关数据 年份 举办城市 经济影响力 主题 1893 美国 芝加哥 被称为“改变美国的一届世博会”宣告美国跨入世界一流强国的行列 并从此一跃成为 GDP 世界排名第一的国家。进步的世纪 1967 加拿大 蒙特利尔 加速建造了新的地铁系统 扩充了城郊高速公路网 人类与世界 1970 日本大阪 上世纪六七十年代日本 1970 年 GDP 达到最高值21.2%在全日本产生了 15600 亿日元收入和 5700 亿日元附加值 人类的进步与和谐 1992 西班牙 塞维亚 投资约 100 亿美元用于园内(14%)园外(86%)建设 建成了全长 470 公里的马德里塞尔维亚高速火车 改善和扩建了室内机场、道路、桥梁等基础设施 发现的时代 2000 德国 汉诺威 直接指出为 102 亿马克 在德国全国 GPD 的效应接近 0.26%人类 给当地政府带来了 20 亿马克的财政税收 2010 中国上海 市场开发良好,运营期间收入 16.91 亿元 基础设施投资约 300 亿元 入园人数预测达7000千万,创1851伦敦世博会以来最高 城市,让生活更美好 图 3:六届世博会参观人数与举办天数的柱形图 例如对直接影响因素中“旅游经济”这一评价指标,我们通过这六届世博会的总参观人数(如图 3,其中上海的数据是官方预测得到)来确定旅游经济这一评价指标的权重 大小,由下表我们可以根据参观人数的情况依次给美国芝加哥世博会,加拿大蒙特利尔 世博会中国上海世博会赋值:3、5、6、4、3、7。下面我们以上海世博会为例,对直接影响的四个评价指标和间接影响的三个评价指标分别构造成对比较矩阵为:根据和法求出其最大特征根和特征向量分别为:02.4A,48.0,08.0,15.0,29.01,02.3B,35.0,45.0,20.02;对以上两个成对比较矩阵进行一致检验:对于A有:0067.01nnCIA,而4,0.90nRI 从而我们可以得到1.0007.0RICICR,因此A的不一致性在容许的范围内,可用其特征向量1作为其权向量 对于B有:01.01nnCIB,而3,0.58nRI 从而我们可以得到1.0017.0RICICR,因此B的不一致性在容许的范围内,可用其特征向量2作为其权向量.取两个主要因素的权向量 通过参考上海世博会的各项数据和专家和群众对各项指标的评价意见,我们得到:由层次分析法我们得到:20.19,0.45,0.36;从而我们可以得到(程序见附录 1):根据以上方法,我们同理可以得到上述六届世博会的综合评分如表 4:表 4:六届世博会的综合评分汇总表 举办时间 1933 1967 1970 1992 2000 2010 举办地 芝加哥 蒙特利尔 大阪 塞尔维亚 汉诺威 上海 综合评为 75.12 80.01 80.11 77.35 79.85 80.81 由上表可知,通过对历届世博会经济效益的各项评价指标的综合评估,上海世博会的经济效益的影响力是历届最强的,这也从一个方面说明上海世博会必将对上海的经济发展乃至全国的经济增长产生积极的推动作用,世界各地的游客,商人,政客涌入上海,必将促使上海成为国际化的政治和金融中心。模型二 本模型中,先以上海不申办世博会为假设,建立出一套投入产出模型后,在此前提下预计出 2008 年上海GDP总额,与 08 年上海实际GDP总额进行对比,从而可得到世博会的影响力。同时,以该模型为基础,可以分别算出世博会第n年后,是否举办世博会对该地区所产生的持续性影响。4.1 模型建立与求解 根据 2009 年上海统计年鉴发布的数据 年 份 固定投资合 计 上海市 生产总值 年 份 固定投 资合 计 上海市 生产总值(亿元)(亿元)1978 27.91 272.81 1993 653.91 1519.23 1979 35.58 286.43 1994 1123.29 1990.86 1980 45.43 311.89 1995 1601.79 2499.43 1981 54.6 324.76 1996 1952.05 2957.55 1982 71.34 337.07 1997 1977.59 3438.79 1983 75.94 351.81 1998 1964.83 3801.09 1984 92.3 390.85 1999 1856.72 4188.73 1985 118.56 466.75 2000 1869.67 4771.17 1986 146.93 490.83 2001 1994.73 5210.12 1987 186.3 545.46 2002 2187.06 5741.03 1988 245.27 648.3 2003 2452.11 6694.23 1989 214.76 696.54 2004 3084.66 8072.83 1990 227.08 781.66 2005 3542.55 9164.10 1991 258.3 893.77 2006 3925.09 10366.37 1992 357.38 1114.32 2007 4458.61 12188.85 1993 653.91 1519.23 2008 4829.45 13698.15 由以上数据经过初步处理后,可得出以下结论 进入 90 年代后上海GDP的呈指数型增长(计 1978 年为第 1 年,2008 年为第 30 年)GDP总量的对数值可以认为与与时间增长有着二次相关关系(见程序 4)(以 1978 年为第一年)GDP总量的对数值与投入资本呈一次线性关系。(见程序 5)(计 1978 年为第一年)为了使得出结果更有代表性我们截取(世博申办成功前)1997 年至 2002 年的数据 年份 1997 1998 1999 2000 2001 2002 固定投入 1977.59 1964.83 1856.72 1869.67 1994.73 2187.06 GDP 总量 3438.79 3801.09 4188.73 4771.17 5210.12 5741.03 根据上述结论,建立GDP总量与时间、投入关系的 的投入产出函数模型 其中1Q为从 1997 年起第nx1年的上海市GDP总量 其中nx2为从 1997 年起第nx1年的上海市固定投入总量 将等式两边取对数,变形为 即将问题转化为一个二元二次回归问题(求解过程见程序 6)。在 MATLAB 求解过程中,取前五组数据解得投入与产出模型的解:并将第六组固定资产投入带入后进行检验,可得 2002 年的GDP预测总量为:5803.9 亿元,与实际值比较,误差为 1.1%,符合误差要求。因此,在不申办世博的假设下,我们对 2008 年上海GDP总量进行预测 1082582945.40904.081117.080032.06278.812eQ(单位:亿元)再与申办世博后上海市 2008 年实际 GDP 总量 13698.15 亿元进行对比 即可以认为在 2008 年世博会对上海市GDP存在 20.9%的正面影响作用。再以上述模型为基础,建立举办世博会的前提下的投入产出模型 年份 2003 2004 2005 2006 2007 2008 固定投入 2452.11 3084.66 3542.55 3925.09 4458.61 4829.45 GDP 总量 6694.23 8072.83 9164.10 10366.37 12188.85 13698.15 同样以前五组数据为基础,在 MATLAB 拟合后得如下解(见程序 7)。并将第六组固定资产投入带入后进行检验,可得 2008 年的GDP预测总量为:13753 亿元,与实际值比较,误差为 0.4%,符合误差要求。因此在已知第n年上海市固定投入总额的前提下 由%100212QQQ 可以得出世博会后,世博期间投入对上海地区的持续性影响。如:假设 2011 上海市固定投资为 5600 亿元,计算得:%18 则 2011 年时,世博会对上海经济有18%的积极影响。五、模型的评价 模型一的优点:该模型采用了模糊判别矩阵,能很好的将各种二级因素通过运算得到了各因素的评分和综合评分,这些都表明了所述因素都是符合客观实际的,这一模型是可以应用于实际的,而且具有相当广泛的通用性。模型一的缺点:在模型中关于各因素的权重主要是依据相关数据主观上给定的,可能带有片面性,实际中可以根据侧重点适当调整,不影响整体效果。如果能详细的分析各因素之间的关系,通过两两比较矩阵确定权向量,结果可能会更客观,但工作量可能会相对更大。模型二的优点:1、模型一重点考虑世博会对经济的影响,从世博申办前上海经济发展状况与申办之后上海经济发展状况的角度,采用层层递进以及图像拟合的方法,最终找出了上海市GDP 的对数值与年份和社会固定投资的二元二次关系。2、我们利用建立的数学模型对2009 年上海市经济进行预测,发现误差率为3%,可以说明模型是可以用来预测2010 年上海市 GDP 的。3、上海市自从2003年申办到世博会举办权,经济开始突增,所以建立两个表达式,用表达式一模拟 2003 年到 2010 年 GDP,来和实际 GDP 进行比较,凸显世博会所带来的前期经济效益。思路清晰明了,简单易行。模型二的缺点:1、考虑因素全面而且深刻,采用主成分分析的方法,找出对经济贡献率较大的三大因素。另外我们在讨论经济的同时,也考虑到社会效益,使用定性与定量相结合的处理方法。2、对模型二进行灵敏度分析以及稳定性分析,另外采用对比的方法来刻画影响力的程度,方法科学严谨。参考文献 1姜启源,数学模型(第三版),高等教育出版社,2003 2Mark M.Meerschaert,Mathematical Modeling(Third Edition),机械工业出版社,2009 3於世为,高等教育经济效益模糊综合评价改进模型,科技进步与对策,2005 年 02期 4王建伟,MATLAB7.X 程序设计,中国水利水电出版社,2007 5曲雪莲,2010 年世博会对上海经济拉动效应研究,6肖道刚,基于成本收益理论的上海世博经济影响研究,7上海市统计局,6 月份世博园消费更趋活跃 ges/00001652.gif?/2005shtj/tjfx/jdxx/images/00001653.gif 8罗燕倩,星级酒店客房出租率达 80%0910/09/6G7A2FAC00014AED.html 1 姜启源,数学模型(第三版)M,北京:高等教育出版社,2003。2 凤凰网,日本人难忘的大阪世博会,2010 年 9 月 10 日 3 中华网,?大阪世博会,强力推动日本经济,2010 年 9 月 10 日 4 人民网,看大阪世博会如何助推日本经济,2010 年 9 月 10 日 5 上海统计局,上海统计年鉴 2009,附录 程序 1:b=1,2,4,1/2;1/2,1,2,1/3;1/4,1/2,1,1/5;2,3,5,1 b=1.0000 2.0000 4.0000 0.5000 0.5000 1.0000 2.0000 0.3333 0.2500 0.5000 1.0000 0.2000 2.0000 3.0000 5.0000 1.0000 x,lumda=eig(b)x=0.4932 0.1700-0.4437i 0.1700+0.4437i 0.5345 0.2635 0.0940+0.0255i 0.0940-0.0255i -0.8018 0.1385 0.0506+0.1117i 0.0506-0.1117i 0.2673 0.8174 -0.8659 -0.8659 -0.0000 lumda=4.0211 0 0 0 0 -0.0106+0.2913i 0 0 0 0 -0.0106-0.2913i 0 0 0 0 0.0000 r=abs(sum(lumda)r=4.0211 0.2915 0.2915 0.0000 n=find(r=max(r)n=1 max_lumda_b=x(:,n)max_lumda_b=0.4932 0.2635 0.1385 0.8174 w=x(:,n)/sum(x(:,n)w=0.2880 0.1539 0.0809 0.4773 max_lumda_b=lumda(n,n)max_lumda_b=4.0211 程序 2:a=1 1/2 1/2;2 1 3/2;2 2/3 1 a=1.0000 0.5000 0.5000 2.0000 1.0000 1.5000 2.0000 0.6667 1.0000 x,lumda=eig(a)x=0.3280 0.1640+0.2841i 0.1640-0.2841i 0.7510 -0.7510 -0.7510 0.5731 0.2865-0.4963i 0.2865+0.4963i lumda=3.0183 0 0 0 -0.0091+0.2348i 0 0 0 -0.0091-0.2348i r=abs(sum(lumda)r=3.0183 0.2350 0.2350 n=find(r=max(r)n=1 max_lumda_a=lumda(n,n)max_lumda_a=3.0183 w=x(:,n)/sum(x(:,n)w=0.1985 0.4546 0.3469 程序 3:d1=0.29,0.15,0.08,0.48 d1=0.2900 0.1500 0.0800 0.4800 d2=0.20,0.45,0.35 d2=0.2000 0.4500 0.3500 a4=0.5,0.3,0.1,0.05,0.05;0.4,0.2,0.2,0.1,0.1;0.2,0.3,0.4,0.05,0.05;0.3,0.3,0.2,0.1,0.1 a4=0.5000 0.3000 0.1000 0.0500 0.0500 0.4000 0.2000 0.2000 0.1000 0.1000 0.2000 0.3000 0.4000 0.0500 0.0500 0.3000 0.3000 0.2000 0.1000 0.1000 b4=0.4,0.3,0.1,0.1,0.1;0.3,0.3,0.1,0.2,0.1;0.2,0.3,0.1,0.3,0.1 b4=0.4000 0.3000 0.1000 0.1000 0.1000 0.3000 0.3000 0.1000 0.2000 0.1000 0.2000 0.3000 0.1000 0.3000 0.1000 p=0.4,0.6 p=0.4000 0.6000 m1=d1*a4 m1=0.3650 0.2850 0.1870 0.0815 0.0815 m2=d2*b4 m2=0.2850 0.3000 0.1000 0.2150 0.1000 R=p*m1;m2 R=0.3170 0.2940 0.1348 0.1616 0.0926 z1=m1*95;85;75;65;55 z1=82.7050 z2=m2*95;85;75;65;55 z2=79.5500 z3=R*95;85;75;65;55 z3=80.8120 程序 4:t=3 4 5 6 7 8 s=8.8090 8.9963 9.1230 9.2463 9.4083 9.5250 T=ones(6,1)t(t.2);b,bint,r,rint,stats=regress(s,T);b,stats t=3 4 5 6 7 8 s=8.8090 8.9963 9.1230 9.2463 9.4083 9.5250 b=8.3051 0.1827 -0.0038 stats=0.9972 529.1794 0.0002 0.0003 plot(t,s,r)hold on fplot(8.3051+0.1827*t-0.0038*t2,3 8)hold on title(上海 GDP 对数值与年份关系);xlabel(年份-2000),ylabel(上海 GDP 对数值)legend(实际数据,拟合曲线)程序 5:t=2452.11 3084.66 3542.55 3925.09 4458.61 4829.45/1000 s=8.8090 8.9963 9.1230 9.2463 9.4083 9.5250 plot(t,s,r)t=2.4521 3.0847 3.5426 3.9251 4.4586 4.8294 s=8.8090 8.9963 9.1230 9.2463 9.4083 9.5250 plot(t,s,r)hold on fplot(8.0662+0.3010*t,2 5)hold on title(上海 GDP 对数值与投资总额关系);xlabel(投资总额/1000),ylabel(上海 GDP 对数值)legend(实际数据,拟合曲线)程序 6:Q1 参数确定:x1=1997-2000 1998-2000 1999-2000 2000-2000 2001-2000 x2=1977.59/1000 1964.83/1000 1856.72/1000 1869.67/1000 1994.73/1000 y=8.1429 8.2430 8.3402 8.4703 8.5584;X=ones(5,1)x1 x2(x1.2);b,bint,r,rint,stats=regress(y,X);b,stats x1=-3 -2 -1 0 1 x2=1.9776 1.9648 1.8567 1.8697 1.9947 b=8.6278 0.1117 -0.0904 0.0032 stats=0.9974 130.2157 0.0643 0.0003 y=8.6278+0.0032*4+0.1117*2-0.0904*2.18706 y=8.6663 exp(y)ans=5.8039e+003 (y-5741.03)/5741.03 ans=-0.9985 (exp(y)-5741.03)/5741.03 ans=0.0110 程序 7:Q2 参数确定:x1=2003-2000 2004-2000 2005-2000 2006-2000 2007-2000 x2=2452.11/1000 3084.66/1000 3542.55/1000 3925.09/1000 4458.61/1000 y=log(6694.23 8072.83 9164.10 10366.37 12188.85);X=ones(5,1)x1 x2(x1.2);b,bint,r,rint,stats=regress(y,X);b,stats x1=3 4 5 6 7 x2=2.4521 3.0847 3.5426 3.9251 4.4586 b=8.1211 -0.0176 0.2955 0.0019 stats=1.0e+003*0.0010 1.7622 0.0000 0.0000