经济增长问题的数学建模论文6实用精品资料(00001)

1、1. GDP与一二三产业的散点图2. 道格拉斯生产函数-资本与劳动力讨省内生产总值增长与资本及劳动之间的关系多元线性回归显著性检验 道格拉斯函数 数据拟合2012高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛

2、的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）： 聊城大学 参赛队员 (打印并签名) ：1. 孔鹏 2. 郑国梁 3. 吴其诚 指导教师或指导教师组负责人 (打印并签名)： 日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：题目经济增长问题摘 要山东省生产总值（GDP）常被公认为衡量省内经济状况的最佳指标.它不但可反映一个省的经济发展情况，更可以反映一省的省力与财富.因此分析各产业对于GDP的影响，并研究GDP的增长规律是具有现实意义的.在问题一

3、中，我们分别做出了GDP与工业、建筑业及农林渔业产值关系的散点图，分析得出GDP的值与各产业之间存在明显的线性关系. 回归分析是统计分析的重要组成部分，用回归分析方法来研究自变量与因变量的关系函数是一种常用的有效方法.因此我们建立起了多元线性回归模型，用MATLAB计算得到的模型为.在对该模型进行显著性检验中，我们对各参数进行了显著性分析，得到模型的复相关系数=0.998，统计量=29666 统计量的值远超过检验的临界值，因此可以验证模型是可用的.分析了各产业对GDP的影响.通过处理预测的数据，我们得出平均每年GDP的增长率为9%左右，其中建筑业与工业对GDP的影响较大，而农林渔业对GDP的影

4、响较小，这也符合山东省的产业结构与经济发展情况.在问题二中，为了讨论省内生产总值增长与资本及劳动之间的关系，我们通过分析数据、查阅相关资料，了解到了省内生产总值的大小通常取决于相关的生产资料和劳动力等相关重要因素.于是，我们通过建立柯布道格拉斯生产函数，定义了三个指数分别为：投资金额指数，就业人数指数和省内生产总值指数.利用定义的三个指数公式，计算出1981年到2010年的省内生产总值指数，投资金额指数和就业人数指数的一组数据，并探讨省内生产总值增长与资本及劳动之间的关系.但是上述三个指标都是随时间增长的，很难直接从表中发现具体的经济规律.为了定量分析，我们定义两个新的变量分别,通过做散点图发

5、现这两个变量基本上成正比例关系.我们用MATLAB软件中的curvefit()函数来作数据拟合，求得函数中的未知参数,，通过检验进而得处道格拉斯生产函数为，这就是产值随资金K、劳动力L的变化规律.为了验证第二问中的结果，我们用求得的道格拉斯生产函数来预测每年的山东省生产总值，然后与题目提供的数据进行比较来进行检验.通过检验可以发现预测值的误差很小，因此道格拉斯生产函数可以表示出省内生产总值增长与资本及劳动之间的关系.关键字：多元线性回归显著性检验 道格拉斯函数 数据拟合1. 问题描述山东省生产总值(Gross Domestic Product,简称GDP)是指在一定时期内(一个季度或一年)，一

6、个省或地区的经济中所生产出的全部最终产品和服务的价值，常被公认为衡量一个省经济状况的最佳指标.它不但可反映一个省的经济表现，更可以反映一省的省力与财富.（1） 建立省内生产总值与工业值、建筑业及农林渔业产值之间的数量模型。（2） 讨论省内生产总值增长与资本及劳动之间的关系，利用数据验证其结果.2. 问题分析2.1 建立GDP与工业值、建筑业及农林渔业产值之间的数量模型，并对未来经济做出预测.在问题一中，我们通过分析材料得出这是研究对象的内在特性和各个因素间关系的问题，即研究GDP与工业值、建筑业及农林渔业产值关系.一般用机理分析的方法建立数学模型.由于经济问题是一种随机的问题，所以通常的办法是

7、搜集大量的数据，基于对数据的统计分析去建立模型.因为影响GDP的因素有三个，即工业值、建筑业及农林渔业产值，且各个产业与GDP都为线性关系.所以我们建立起一个多元线性回归模型，并检验模型显著性，通过对模型的反复修改与检验，建立更合理的模型.2.2 讨论省内生产总值增长与资本及劳动之间的关系，并验证其结果.在问题二中，为了讨论省内生产总值增长与资本及劳动之间的关系，通过查阅相关资料，我们了解到省内生产总值通常取决于相关的生产资料和劳动力等相关重要因素. 要建立道格拉斯生产函数，我们只需要讨论产值和资金，劳动之间的关系，从而达到我们的目的.这样处理不仅能简化问题，而且是合理的在生产产值上的预测，柯

8、布-道格拉斯（Cobb-Douglas）生产函数预测的结果近似就是准确生产值.于是我们通过建立柯布道格拉斯生产函数，来探讨省内生产总值增长与资本及劳动之间的关系，进而利用已有的数据验证其结果.3. 模型假设1. 假设所统计的数据都在误差允许的范围之内；2. 忽略由于非正常条件下的引起的数据的巨大波动；3. 假设在短期内省内生产总值只取决于投资和劳动力因素;4. 假定在一段不太长的时间内技术水平不变.4. 定义与符号说明山东省省内生产总值复相关系数统计量随机误差置信水平工业产值建筑业产值农林渔业产值省内生产总值增长劳动投入量资本投入数量劳动对产出的贡献程度资本对产出的贡献程度投资金额指数就业人数

9、指数省内生产总值指数道格拉斯函数常数5. 模型的建立与求解5.1 问题一模型的建立与求解：回归分析方法是统计分析的重要组成部分，用回归分析方法来研究自变量与变量的关系函数是一种常用的有效方法.我们通过回归模型的建立，定量预测了未来经济的发展.5.1.1GDP与工业值、建筑业及农林渔业产值数量模型： 通过在互联网上搜集到1981年2010年，山东省GDP与工业值、建筑业及农林渔业产值的数据（见附表1），可以定性的看出GDP与工业、建筑业及农林渔业产值为整体上升的趋势.为了大致分析GDP与工业值、建筑业及农林渔业产值关系，我们首先利用附录数据做出了GDP与工业产值的关系散点图（如图1）.图1 工业

10、产值与GDP散点图 从图可以发现，随着工业产值的增加，GDP的值有比较明显的线性增长趋势.图中的直线是用线性模型拟合的.同理我们也分别作出了建筑业产值与GDP的关系散点图（图2）、农林渔业产值与GDP的关系散点图（图3）.图2 建筑业产值与GDP散点图图3 农林渔业产值与GDP散点图通过图2，图3可以看出建筑业产值与农林渔业产值同样有很强的线性关系，同样也分别用直线模型对其拟合.建筑业产值与GDP线性模型农林渔业产值与GDP线性模型因此，综上所述四者之间有很强的线性关系，可建立多元线性回归模型在模型中除了工业，建筑业，农林渔业外，影响国内生产总值的其他因素的作用都包含在随机误差内，这里假设相互

11、独立，且服从均值为零的正态分布，.对模型直接利用matlab统计工具箱求解，得到回归系数估计值及其置信区间（置信水平= 0.05），检验统计量，的结果见表1.参数参数估计置信区间-80.7854-265.3513,103.78041.15211.0356,1.26855.51544.3792,6.65151.11390.6934,1.5345 表1 模型的计算结果5.1.2 结果分析: 表1显示，指因变量（国内生产总值）的99.97%可由模型确定，值远远超过检验的临界值，p远小于，因而模型从整体上来看是可用的.表1的回归系数给出了模型中，的估计值，即，.检查它们的置信区间发现与的置信区间都包含

12、零点，这表明回归模型常数与回归变量对模型的影响不太显著.这也符合这一事实，农林渔业产值对GDP的影响较小，工业与建筑业对GDP的影响较大.但是一般情况下常数的值都保留在模型中，不剔除.回归变量系数区间右端点相对距零点较远，所以我们也保留在模型中.因此最终确定的模型为：5.2 问题二模型的建立与求解：在问题二中为了讨论省内生产总值增长与资本及劳动之间的关系，通过查阅相关资料，我们了解到省内生产总值通常取决于相关的生产资料和劳动力等相关重要因素.于是我们通过建立柯布道格拉斯生产函数，来探讨省内生产总值增长与资本及劳动之间的关系，进而利用已有的数据验证其结果.5.2.1 模型的建立:在经济学的分析中

13、，为了简化分析，通常假定生产中只有劳动和资本这两种生产要素.若以表示劳动投入量，以表示资本投入数量，则生产函数可以写为：生产函数表示生产中的投入量和产出量之间的依存关系，这种关系普遍存在于各种生产过程中.一家工厂必然具有一个生产函数，一家饭店也是如此，甚至一所学校或者医院同样会存在着各自的生产函数，产品可能是实实在在的有形产品，也可能是无形产品比如服务.估计和研究生产函数，对于经济理论和实践经验都具有一定意义.柯布道格拉斯（Cobb-Dauglas）生产函数是由数学家柯布和经济学家道格拉斯于20世纪30年代初一起提出来的.柯布道格拉斯生产函数被认为是一种很有用的生产函数，因为该函数以极简单的形

14、式描述了经济学家所关心的一些性质，它在经济理论的分析和实证研究中都有具有一定意义.柯布道格拉斯生产函数的函数表达式如下：其中，代表产出量,代表资本投入量,代表劳动投入量,、为未知参数.表示技术或管理等参数对经济增长的影响系数, 和分别表示劳动和资本对产出的贡献程度,且.对该生产函数取对数得:由于柯布道格拉斯生产函数假设技术、管理水平不变,即是一个常数, 在此可以忽略的影响.所以，可简化为:求出道格拉斯函数以后，我们通过道格拉斯函数可以预测出来每一年的GDP总产值，然后利用题目所提供的数据进行检验，可以发现道格拉斯很好的表示出来了省内生产总值的增长和投资与劳动之间的关系. 5.2.2 模型的求解

15、:为了求解上述模型，通过分析题目所给数据和从网上查找相关数据，我们列出了我省从1981年到2003的GDP总值，投资金额总和和我省就业人数的表2如下： 表2年份GDP值投资金额就业人数1981346.5779.60 2885.4 1982395.3885.00 2948.6 1983459.8396.46 3029.7 1984581.56140.15 3140.4 1985680.46194.33 3276.8 1986742.05223.08 3425.0 1987892.29297.77 3561.5 19881117.66369.82 3680.1 19891293.94305.54

16、3803.4 19901511.19335.66 3927.4 19911810.54439.82 4219.3 19922196.53601.50 4302.6 19932770.37892.48 4379.3 19943844.51108.00 4382.1 19954953.351320.97 5207.4 19965883.81558.01 5227.9 19976537.071792.22 5256.0 19987021.352056.97 5287.6 19997493.842222.17 5314.7 20008337.472542.65 5441.8 20019195.0428

17、07.79 5475.3 200210275.53509.29 5527.0 200312078.155328.44 5620.6 2010GDP:亿元；投资金额：亿元；就业人数：万人；在实际生产中，人们关心的往往是生产的增长量，而不是绝对量，因此定义投资金额指数，就业人数指数和省内生产总值指数分别为利用上述定义的三个指数公式，通过使用matlab软件计算出表3中1981年到2003年的省内生产总值指数，投资金额指数，和就业人数指数的一组数据，取1990年为基年,则t=0.t-90.21270.69690.2612-80.27240.71600.2848-70.31660.73030.3198

18、-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.32

19、941.11335.2355118.40891.12785.7734129.97161.13896.36251312.97541.14957.2215表3从表中可知，在正常的经济发展过程中（除个别年份外），上述三个指标都是随时间增长的，但是很难直接从表中发现具体的经济规律.为了定量分析，定义两个新的变量根据表中数据，在直角坐标系上做出的散点图，发现基本上成正比例关系（散点位于一条直线的附近），如图5图5 散点图我们可以用MATLAB软件中的curvefit()函数来作数据拟合，即寻求函数Q(K，L)中的未知参数,，使这个函数尽量逼近表5-2-2所给出的统计数据.则可以得到：,.于是公式变为：这

20、就是产值随资金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

21、.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 省内生产总值的预测与实际比较图增加生产、发展经济所依靠的主要因素有增加投资、增加劳动力以及技术革新等，在研究省民经济产值与这些因素的数量关

22、系时，由于技术水平不像资金、劳动力那样容易定量化，作为初步的模型，可认为技术水平不变，只讨论产值和资金、劳动力之间的关系.在科学发展不快时，如资本主义经济发展的前期，这种模型是有意义的.从而可以说明省内生产总值增长与资本及劳动之间满足柯布道格拉斯（Cobb-Dauglas）生产函数的关系.6. 模型的评价与推广6.1 模型的优点：在问题一中，多元回归模型，因变量省内生产总值的99.89%可由模型确定，说明模型从整体上来看是可用的.在问题二中，运用了柯布道格拉斯生产函数，使该模型的建立有理论依据作支撑，且有助于对模型的结果进行分析.在分析省内生产总值与投资和劳动力关系是，忽略其他因素，从而简化了

23、模型，便于大概的预测.6.2 模型的缺点：问题一中，由于省内生产总值受国际经济、政府政策、自然灾害等因素的影响，所以某一时期GDP波动幅度较大，因此影响了模型整体预测的准确性.问题二中，忽略其他因素对省内生产总值的影响，和实际问题存在的误差.一定历史时期的生产函数是反映当时的社会生产力水平的.6.3 模型的推广与改进：推广：模型一是一类基于统计分析的随机模型，因此适用于大量数据的随机现象.如经济增长，灾害预测等.模型二中，在信息经济时代，所投入的生产要素的核心成分从资本、劳动力逐渐转变为以信息技术为代表的高新技术.当信息资源应用于生产中时，对生产人员、资本、流程等形成革命性的影响作用，极大地提

24、高了生产要素生产率，促进了经济发展.综合上述原因，需要对柯布道格拉斯生产函数做出了一定的修正，使之适用于信息时代的生产力发展水平.改进：模型一中参数与的置信区间包含零点，说明模型中还存在缺点，变量之间很可能存在交互作用.因此应在模型中加入交互项，改进原有的模型，建立新的回归模型.模型二较原来的模型增加了信息技术设备的资本投入和信息技术的劳动力投入后，得到使得模型成为更贴近时代的生产模型，改进后的柯布道格拉斯生产函数是在现代信息工业经济时代构造出的反映了现代信息工业经济时代生产力特征的函数模型.改进后的柯布道格拉斯生产函数模型更具有时代特色，适用性更广.7. 参考文献1 姜启源，谢金星，数学模型

25、，高等教育出版社，2007.2 韩中庚，数学建模方法及其应用，高等教育出版社，2006.3 周品，赵新芬，MATLAB数学建模与仿真，2009.4 王兵团，数学建模基础，2005.5 齐微，柯布道格拉斯生产函数模型，中国科技论文在线.8. 附录附表1：2-1 1952-2010年地区生产总值单位:亿元年 份地 区生产总值第一产业第二产业第三产业人均地区生产总值(元)工 业建筑业#交通运输仓储 邮电通信业#批发零售 贸易餐饮业1981346.57 132.21 155.41 138.09 17.32 58.95 13.18 16.56 472 1982395.38 154.07 166.05 1

26、47.10 18.95 75.26 14.10 22.39 531 1983459.83 185.57 178.75 159.15 19.60 95.51 17.12 32.07 611 1984581.56 222.13 239.27 214.20 25.07 120.16 22.25 37.97 765 1985680.46 235.96 293.07 259.42 33.65 151.43 28.22 46.14 887 1986742.05 252.73 313.21 274.80 38.41 176.11 34.03 49.85 956 1987892.29 287.31 384.5

27、7 341.31 43.26 220.41 47.45 60.58 1131 19881117.66 331.94 497.10 435.51 61.59 288.62 55.50 86.71 1395 19891293.94 359.14 579.65 513.97 65.68 355.15 66.55 108.51 1595 19901511.19 425.29 635.98 568.25 67.73 449.92 81.87 132.83 1815 19911810.54 521.85 745.90 663.90 82.00 542.79 98.97 161.00 2122 199221

28、96.53 534.62 999.11 889.59 109.52 662.80 120.18 200.93 2556 19932770.37 596.63 1355.71 1201.67 154.04 818.03 141.09 239.98 3212 19943844.50 775.03 1891.43 1692.10 199.33 1178.04 213.58 346.88 4441 19954953.35 1010.13 2355.78 2098.06 257.73 1587.44 296.08 471.95 5701 19965883.80 1200.17 2784.09 2475.

29、99 308.10 1899.54 362.05 574.97 6746 19976537.07 1195.00 3147.37 2796.02 351.35 2194.70 420.49 664.92 7461 19987021.35 1215.81 3408.06 3008.45 399.61 2397.49 437.37 735.63 7968 19997493.84 1221.00 3644.32 3197.16 447.16 2628.52 484.65 787.90 8483 20008337.47 1268.57 4164.45 3665.74 498.71 2904.45 54

30、5.13 856.94 9326 20019195.04 1359.49 4556.01 4004.09 551.92 3279.53 657.57 972.33 10195 200210275.50 1390.00 5184.98 4518.87 666.11 3700.52 655.64 1142.54 11340 200312078.15 1480.67 6485.05 5706.71 778.34 4112.43 710.18 1283.70 13268 200415021.84 1778.45 8478.69 7576.12 902.57 4764.70 969.13 1431.58

31、 16413 200518366.87 1963.5110478.62 9418.58 1060.04 5924.7419934 200621900.19 2138.90 12574.03 11378.82 1195.21 7187.26 23603 200725776.91 2509.14 14647.53 13283.72 1363.81 8620.24 27604 200830933.28 3002.65 17571.98 15894.95 1677.03 10358.64 32936 200933896.65 3226.64 18901.83 16896.14 2005.69 1176

32、8.18 35894 201039169.92 3588.28 21238.49 18861.45 2377.04 14343.14 41106 附表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.32

33、98.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.420012

34、546.4729.8430.7830.020023161.0859.1485.0841.920032687.81044.4499.4854.820043254.91188.0585.7 936.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.

35、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.

36、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

37、(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)1

38、40(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

39、)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(

40、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)2005

41、12.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.

42、34(7.24)111(150)60(81)51(69)模型二计算程序：Q=1117.661193.941511.191810.542196.532770.373844.504953.355883.80 6537.07 7021.35 7493.84 8337.47 .9195.04 10275.50 12078.15 15021.84 18366.87 21900.19 25776.91 30933.28 33896.05 39169.92;IQ=Q/11795.0K=421.0 407.8 455.6 496.6 432.0 366.1 393.0 380.4 347.5 483.7 52

43、9.6 621.1 625.9 696.7 1197.4 1852.1 .2109.5 2546.4 3161.0 2687.8 3254.9 3554.3 3672.0;IK=K/1334.5L=2885.4 2948.6 3029.7 3140.4 3276.8 3425.0 3561.5 3680.1 3803.4 3927.4 4219.3 4302.6 4379.3 4382.1 5207.4.5227.9 5256.0 5287.6 5314.7 5441.8 5475.3 5527.0 5620.6 ;IL=L/4318.3Et=zeros(23);Et=Et(1,1:23);f

44、or t=1:1:23; Et(t)=log(IL(t)/IK(t);endEt;Wt=Et;for t=1:1:23; Wt(t)=log(IQ(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 7.3294 8.4089 9.9716 12.9754; 0.6969 0.7160 0