人民币汇率变动对我国进出口贸易的影响论文(修改).doc

1、湖北汽车工业学院大学生数学建模邀请赛 承诺书我们仔细阅读了第二届华中地区大学生数学建模邀请赛的竞赛细则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们的参赛报名号为：05004参赛队员（签名）：队员1:T933-2杨进队员2:T843

2、-3李志刚队员3:T923-1叶胜峰0第二届华中地区大学生数学建模邀请赛编号专用页选择的题号：A参赛的编号：05004（以下内容参赛队伍不需要填写）竞赛评阅编号:ii探讨人民币汇率变动对我国进由口贸易的影响【摘要】近几十年来，世界经济朝全球化发展，资本流动自由化、国际金融、贸易的发展在 不同方向上增加或降低汇率波动对国际贸易产生了深远的影响。我国经济一直保持着良好的发展势头，特别是出口贸易的增长让一些发达国家羡慕。人们对中国的经济前景比 较乐观，人民币的国际地位越来越高，经济实力的提升客观上推动了人民币升值。为此，我们针对人民币汇率的不断提高对我国进出口商品贸易的影响，各个产业结 构的变化以及

3、由此对我国国内居民消费水平的影响做了定性和定量的分析。首先，我们用灰色系统从长期的角度分析了我国未来 5年的人民币汇率发展趋势， 及短期几天内人民币汇率发展趋势， 这样使得我们可以对人民币汇率短期的变化情况又 一个更加直观的了解。结果表明未来我国人民币升值是不可阻挡的趋势。接着，我们通过对各个产业的进出口贸易情况做分析处理，得出历年来我国各个产业的贸易特化系数.通过贸易特化系数的变化来更加直观的看出人民币汇率的改变对我 国各个产业的进出口贸易的影响.分析之后发现人民币汇率的改变对我国大部分产业都 有较大的影响.之后我们以各个产业的贸易特化系数为自变量，人民币汇率为因变量建立每个产业与人民币汇率

4、之间的多元回归分析 ，定量分析了人民币汇率波动对各类商品 影响程度。最后通过整体上分析我国进口贸易额、出口贸易额增长率与人民币汇率增长率、居 民消费指数增长率的关系，以进出口贸易额增长率为自变量，人民币汇率增长率和居民消费指数增长率为因变量建立典型相关分析， 结果表明了人民币汇率增长率对进口贸易 额增长率起到决定性的作用，对出口贸易额也有明显的影响。通过以上人民币汇率波动对我国进出口贸易影响的分析及当前我国人民币汇率面 临的形势，为我国经济快速平稳的发展人民币汇率应作如何调整提出一定的建议和意 见。【关键词】人民币汇率灰色系统贸易特化系数实证分析多元回归模型 居民消费价格指数典型相关模型一、问

5、题的提出2011年8月24日人民币对美元汇率中间价突破 6.4 ,再创2005年人民币汇率改革 以来新高。关于人民币升值的利弊人们褒贬不一。有的人为人民币升值痛心疾首，有的无动于 衷。也有人认为人民币升值现在弊大于利，未来利大于弊，一时间，争论不断。因此， 研究人民币升值对我国的影响是个热门问题。作为大学生我们应当学会关注我国经济的发展，树立强烈的主人翁意识。针对人民 币汇率升高的现象，我们从定性和定量的角度分析去考虑人民币汇率的波动对我国进出 口贸易的影响。我们提出如下问题：1、对我国未来几年的人民币对美元的汇率走势作出预测。2、分析一下人民币汇率的波动对我国进出口商品的影响。3、分析一下人

6、民币汇率、我国居民消费 CPI指数变化率与我国贸易进出口额变化率的 关系。4、结合上面的分析，请说明一下我国政府应该如何调整人民币汇率政策。二、问题的分析人民币升值对我国国民经济及生活各个方面带来多方面的影响。基于这样的前提 下，我们通过所学的知识定量的分析了人民币升值对我国进出口贸易的影响。对于第一问，我们从国家统计局网站收集到我国人民币历年的年均汇率，及从国家 外汇管理局收集了我国近一个月的人民币汇率中间价。对以上的数据分别用灰色系统模型分析，对我国人民币汇率在未来 5年、未来10天的波动状况。针对第二问，我们从国家统计局网站收集到我国近几年以来的贸易进出口分类总 额，对初级产品、工业制品

7、、其他三大类分别计算了其下的各个分支产业的贸易特化系 数，然后做出各个产业的贸易特化系数波动曲线，通过实证分析了人民币汇率波动对各 个不同产业的影响。然后根据利用各求解出的历年来产业的贸易特化系数与人民币汇率 建立线性回归模型，分析人民币汇率的变化对各个产业的影响关系。对于第三问，根据从国家统计局网站得到居民消费水平、货物进出口贸易额、年均 汇率等相关数据，对居民消费价格指数的增长百分比、年均汇率的增长百分比与进口、 出口贸易额增长率建立了典型相关分析，来研究其间的关系。对于第四问，我们结合自己所了解的知识，及上面的分析，对我国汇率应该如何调 整提出了一些建议。三、模型假设与符号说明3.1模型

8、假设：(1)假设数据来源真实可靠。(2)假设选取的指标有一定的代表性。(3)假设进出口贸易影响中只考虑供我国居民的消费指数、人民币汇率的影响(4)假设人民币汇率和商品的贸易特化系数间存在线性关系。(5)假设典型分析的两组数据都是连续变量，均服从多元正态分布。3.2符号说明:符号符号说明i代表产业编号k代表时间C灰色模型后验差X代表产业的贸易特化系数X(n) (k)利率的n次累加序列(100美元折算)x代表居民消费指数增长百分比yi代表进口贸易增长百分比y2代表出口贸易增长百分比z代表年均人民币汇率增长百分比注：其他符号见文中说明四、模型建立与求解4.1 问题一4.1.1 基于灰色系统预测未来5

9、年人民币对美元汇率的数学模型1:分析建模给定原始时间1994 2009年100美元兑换成人民币的序列X(0)(k),对X(0)(k)生成1-AGO序列X(k),另外可得Yn见下表：k123141516X(0)861.87835.10831.42760.40694.51683.10X861.871696.972528.3 911528.4812222.9912906. 09Yn1696.972528.3 911528.4812222.9912906. 09其中，X0=AGOX1),即 kX(k)八 x(0) (i)i 1Yn= X (2) , X(3),，X(15) , X(1)(16) T对于

10、上述的X(0) (k)的GM(1,1)参数按照下式进行最小二乘估计得一-z(1)(2)1-z(3)-z(1)(4) -z(5)B= 11111-z(22)-z(23)-z(24) ")(25)1111-0.5X x(2)1| _0.5X(1)(2)+x(1)(3)JI _0.5X(1)(3)+x(1)(4) _0儿(4)+x(1)(5).0.5k(12)+x(13)】-0.5x(1)(13) +x(1) (14)】 l-0.5x(1)(14) +x(1) (15)】-0.5x(1)(15) +x(1) (16)】11111I .1111将B,Yn代入辨识算式，有:得GM(1,1 )模

11、型为:1)灰微分方程：X(k)-a z(k)=b2)(1)白化方程：d- adtz (k)=b3)白化方程的时间相应式：x(1) (t+1)= x(0)一be力+2 a a4) GM(1,1)模型的检验：残差序列的均方差:( 其中名=:一)0 20原序列的均方差：S2 = X0(k)-X 2'、X(0)(k)(其中 X =)n -1n后验差检验：C=SS2由结果可知：C=0.3889<0.5,故可知预测结果合格。利用MATLA酎程见附录7.2.1 ;运行结果见附录7.2.1结果。其中预测值如下:880人民币汇率年变化及五年预测值86084082080078076074072070

12、0680 19941996199820002002200420062008201020122014预测未来5年人民币对美元汇率表格如下年份20102011201220132014预测值P 738.9522731.4055723.9359716.5426709.2249对结果的评价：由结果可以看出100美元兑成人民币的数量呈递减趋势，即人民币 与美元的汇率越来越小，人民币升值越来越严重。虽然该预测与现在的实际情况不符， 但其趋势是相同的。其可能原因为：以年为时间单位误差太大、预测的时间跨度太大或 原始数据过去的时间太久对现在的影响很小。为此，我们建立如下模型来预测短期情况。4.1.2 基于灰色系

13、统预测未来10天人民币对美元汇率的数学模型：分析建模给定原始时间2011.7.1 -2011.7.29内100美元兑换成人民币3的序列X(0)(k),对X(0)(k)生成1-AGOff列X(k),另外可得Yn见下表：k123192021X646.85646.61646.5644.26644.38644.42X646.851293.461939.9612280.6412924.6 213569.04Yn1293.461939.9612280.6412924.6 213569.04理论知识与4.1.1中的一样，在这里不再累赘。由结果可知：C=0.3089<0.35,故可知预测结果优。利用MA

14、TLA酎程见附录7.2.1 ;运行结果见附录7.2.1结果。其中预测值如下:日 起7/307/318/18/28/38/48/58/68/78/8预 测 值644.40644.2 4644.08643.9 1643.7 5643.5 9643.43643.26643.1 0642.94（备注：预测值为100美元可以兑换人民币数量（元）对结果的评价：由预测结果可以看出人民币兑美元的汇率虽然存在大小波动，但其整体 趋势是不断减小的。预测值也与现实相符，说明人民币升值进一步加剧。4.2 问题二4.2.1 首先根据10个产业的贸易进出口额状况，作出各个产业的贸易特化系数（ TS。, 然后作出其相应的曲

15、线变化图，对其进行一定的说明。问题分析：将我国基础产业分为：初级产品、工业制品、其他。这里只对初级产品、工 业制品分析（收集各产品的进出口额）。出口贸易商品结构有不同层次，在这里我们采用了贸易特化系数的变化趋势来阐述某种贸易商品或贸易产业的结构优化。 贸易特化系 数反应的是一国某贸易产业或贸易产品的净出口额与该产品的进出口总额之比。 其计算 公式为：TSC=(X-M)/(X+M)-1<TSGC 1当TSC指标值介于-1到-0.3时，则代表该国产业的进口值远大于出口值，故该产业属 于“进口强势”型产业，即该产业在国际市场上占有绝对的比较优势； 若TSC指标值介于-0.3与0.3之间，表示其

16、进口与出口值差异不大，代表该产业在国际 上水平较高，属于“水平竞争”型产业； 当TSC指标值介于0.3 1时，表示其出口值远大于其进口值，故该产业在该国处于绝 对的劣势，属于“出口强势”型产业。因此，从各种贸易行业贸易特化系数的变化我们可以大概地看出我国贸易结构调整的方 向。表1:对外贸协商品结构组成产拓名称代号产品类型食品及主要供食用的活动物SITCO初级产品饮料及烟类SITC1初级产品燃料以外的非食用粗原料SITC2初级产品:矿物燃料，润泄油及有关原料SITC3初级产品动植物油脂及蜡SITC4初级产拈化学品及有关产品SITC5资本与技术密集型产品轻纺产品、橡胶制品、矿冶产品及其制品SITC

17、6资源与劳动密集型产品机械和运输设饴SITC7资本与技术密集型产品杂项制品SITCS资源与劳动密集型产品:未分类的其他商品SITC9未分类将所有进出口产品分成两类：初级产品：食品及主要供食用的活动物、饮料及烟类、非食用原料、矿物燃料、动植物、 油脂及蜡工业制成品：化学品及有关产品、轻纺产品、橡胶制品、矿冶产品及其制品、机 械及 运输设备、杂项制品MATLA褊程见附录7.2.1 ,结果如下:1）初级产品及其组成项的图像:初级产品贸易特化系数0.60.40.2-0.2-0.4-0.6-0.810152025从上面可以看出我国初级产品贸易特化曲线持续下降，说明了该产业主要是通过进口，属于强进口型，其

18、对国外市场依赖过强，因此人民币汇率升值是有利于该产业发展 的。食品及主要供食用贸易特化系数0.60.50.40.30.211从上面可以看出该产业的贸易特化曲线基本处于正的，也就是说，主要是向国外出 口，属于强出口口型，人民币升值会不利于此类商品在国外的发展。-0.4饮料及烟类贸易特化系数0.80.60.40.20-0.210152025从上面曲线可以看出特化曲线后一阶段处于不断下降的趋势，该产业正从水平竞争型向进口型转化的趋势，人民币升值对其发展是有利。非食品原料贸易特化系数0-0.1-0.2-0.3-0.4-0.5-0.6-0.7-0.8-0.9矿物，润滑油相关原料贸易特化系数1,0.8.-

19、0.6-r0.4-0.2-r0-_-0.2Lr-0.4 - X-0.6.-0.8 :IIC0510152025从此图可以看出非食品原料贸易特化系数最后将会向 -1趋近，说明了非食品原料贸 易几乎全部依靠进口来提供，若人民币升值，会有利于我国进口非食品原料。上图反映了矿物润滑油相关原料的特化系数逐渐逼近 -1,说明了我国的矿物润滑油 相关原料主要依靠进口，人民币升值将有利于我国对此产品的进口与需求。从上面可以看出我国动植物油脂及蜡的贸易特化系数逐步趋向-1,说明了其进口趋势显著加强，人民币升值有利于其发展。由结果得出如下分析：初级产品、饮料及烟类、非食用原料、矿物燃料、动植物油脂及蜡的贸易特化系

20、数均随着人民币的升值而减小， 食品及主要供食用物的贸易特化系数虽 然有上下波动但整体比较稳定，由此可见人民币升值对我国的进出口贸易有着密切影 响。使我国的初级产品类产业的进口值远大于出口值， 即“进口强势”型产业越来越多 而食品的波动比较平稳，可能受国家相关政策的影响，使我国的食品价格保持平稳。有 利于人民的生活和社会的稳定。2）工业制成品及其组成项的图像:工业制成品贸易特化系数0.30.20.10-0.1-0.2-0.3-0.4-0.50510152025从上面图像可以看出我国工业制成品贸易特化曲线基本处于-0.3到0.3之间，说明了此类产业属于水平竞争型，人民币汇率的升高对其影响不大。但是

21、我们也可以看到在汇率 升高的同时其贸易特化系数也是呈上升趋势，如果人民币持续升值，该产业又转化成为 出口强势型的产业。化学品贸易特化系数-0.2-0.25-0.3-0.35-0.4-0.45-0.5-0.55从上图可以看出，化学品相关产业的贸易特化系数基本上处于 -0.3 一下，即便是又 波动影响也不大，可以认为该产业是进口强型产业，在人民币汇率升高有利于该产业的 发展。从下图可以看出：我国的轻纺橡胶、矿冶贸易特化曲线，在前期处于进口强势到了 后期逐渐转变为出口强势产业，人民币的汇率提高不利于其再国际市场的竞争力。轻纺.橡胶.矿冶贸易特化系数0.5(0.4 -1az" jF10.3

22、Lr0.2 -n/0.1 -.J I/I A,1J 0 -0.1 -1 I 1 J- 0.2 -1 V- 0.3 -一II- 0.4 -c L 'rrrr- 0.5 05101520250.4机械与运输设备贸易特化系数0.20-0.2-0.4-0.6-0.8-1101520从此图可以看出我国机械运输设备贸易特化曲线一直处于上升的趋势, 会向出口型发展，人民币的升值将会不利于其发展。25说明了该产业将杂项制品贸易特化系数0.90.80.70.60.50.40.30.2 01015202513由结果得出如下分析：工业制成品、化学品及有关产品、轻纺产品、橡胶制品、矿 冶产品及其制品、机械及运

23、输设备的贸易特化系数均随着人民币的升值而增大，杂项制 品的贸易特化系数虽然有上下波动但整体比较稳定，由此可见人民币升值对我国的进出口贸易有着密切影响。是我国的工业类产业的进口值远小于出口值，即“出口弱势”型 产业越来越多。对比初级产品 和工业制成品可知：我国的工业水品低下，且只能用缺少 技术含量的初级产品和廉价的劳动力来弥补进出口差额，国家应加大科技的发展，尽快 使工业水平赶上来，加大工业制成品的出口 ，减少因人民币升值而带来的影响。4.2.2多元回归分析建模过程：1）模型的建立设自变量x1，X2,X3.Xp的观测值Xii,Xi2.Xip及因变量y对应的观测值y i满足关系y 1 ='

24、+£ PjXj +易 i=1,2,3- n； j=1,2,3,- p其中，当，,.wn是相互独立且都服从正态分布的随机变量 由最小二乘法，得P =(XTX)(XTY)2）回归方程的显著性检验：利用 F分布函数对变量前的系数检验3）回归系数的显著性检验：利用t分布函数检验。最后得到结果如下：产业名称编R方程显著性 检验（F）参数显著 性检验（t）回归方程备注食品及主 要供食用 的活动物10.02880.01080.0288y - 773.15277 xxx + 309.7116汇率对其 影响显著饮料及烟 酒20.1440（不合格）汇率与此 行业没有 显著关系非食用原 料30.00590

25、.0010.0059y 二 一419392239 X xa + 40C.43293汇率对其 影响显著矿物燃料4<0.0001<0.0001<0.0001y - -309.74464 xx+ + 657.01952汇率对其 影响显著动植物油 脂50.000960.00420.00096y =-458.67428 xxs + 3821132C汇率对其 影响显著化学品60.1256（不合格）汇率对其 影口何不显 著轻纺橡胶7不合格汇率对其 影口何不显 著机械及运 输8<0.0001<0.0001<0.0001y = 406.07219 xxfl + 791,929

26、68汇率对其 影响显著杂项90.2222（不合格）汇率对其 影口何不显 者（备注：y代表年均人民币汇率，100美元兑换的人民币数量（元），Xi代表贸易特化系数）由上表看出，人民币汇率的变化对我国的大部分产业都有显著的影响。也就是说， 人民币持续升值对我国各行各业都会产生决定性改变，因此，为了稳定中国国民经济稳 步发展，在人民币升值过程中国家需要进行干预，严格控制人民币的升值幅度，这样才 能保证中国经济快速稳步发展，中国人民能更美好的生活。4.3 问题三 典型相关分析思想：考察一组随机变量和另一组随机变量之间的相关性，可以用他们之间的相关性。考虑到 两组变量的线性组合：u1 = a11xx1+a

27、21xx2+-.+anlxxnv1 = b11xy1+b21xy2 + -.+amlxyin并研究他们之间的相关系数pvj ,在所有的线性组合中找一对相关系数最大的线 性组合，用这两个组合的单相关系数来表示两组变量的相关性，这两个线性组合叫做一 对典型变量。下一步，再在两组比昂两的与 u, v不相关的线性组合中，找一对相关系数最大的线性 组合，设为：u2=a12xx1+iZ2xxI+_+anlxxnv2 = b12xy1+b22xy2 + -,+am2xyin且P（%，vJ就是第二个典型相关系数，这样下去可以得到若干对典型变量，从而提取 出两组变量间的全部信息。就本题而言，用到典型相关分析。利

28、用 SAS编程得到了：人民币汇率增长率、居民消费指数增长率与进出口贸易额增长率间的两个典型相关系数 及特征值如下表：厅P典型相关系数标准误差特征值方差比率10.9521250.0199259.70010.935720.6324310.1279270.66660.0643从上表可以看出：第一对典型变量的相关系数较高，说明了这一对相关变量间的关系是 密切的。对相关系数检验如下表:厅PF计算值自由度F检验的显著性概 率130.624<0.001213.3310.0016这两对典型相关变量，进行F统计检验，均通过了检验，说明了这两对典型变量间的相 关性显著，可以用来解释产出变量。对冗余度分析得到

29、下表:典型变量的解释能力厅P产出组方差 被影响组解 释比例对产出组的 解释能力产出方差被 典型变量解 释比例对影响组的 解释能力影响组方差 被典型变量 的解释比例.10.906510.95190.86290.87600.7941 120.40000.4810.01930.12400.0496从此表可以看出：这两对典型相关变量的解释能力均较强，特别是第一对相关变量对产 出组和影响组均有很高的解释能力；第一对典型变量对产出组方差有95%Z上的解释能力，第二对典型相关变量较弱些，达到了 48位右。最后得出的典型变量模型：厅P典型相关模型1u1 = 0.0S46 xx + 0.8714 X 2wT =

30、 0.7245 xyt+0.8619 X %u2 = 04305 X x + 0.8869 X 2w2 = 0,8048 X 为+0.8816 X y2（x代表居民消费指数增长百分比、z代表人民币汇率增长百分比；y1代表出口贸易额 增长率、y2代表进口贸易额的增长率；）结果分析：人民币汇率增长率对进口贸易有显著的作用。从上面第一对典型变量分析结果可以看出：人民币汇率增长百分比在典型变量 u1占得比重达到了 0.8714,出口贸易额增长率 在w1中占的比重达到0.8619.由于u1,w1由前面分析知道其相关系数达到了 0.95以上, 可以认为人民币汇率的增长率进口贸易额增长率变化起到决定性的作用

31、。从第二对典型相关变量分析知道：人民币汇率的增长率、居民消费指数增长率对进 出口贸易额增长率都一定的影响，不过期间的影响程度不是那么显著。综上所述：人民币汇率的增长率对我国进口贸易额有着决定性的因素，同时它和居 民消费指数增长率对我国进出口额增长率都有密切的关联。4.4 问题四通过上面的分析在人民币不断升值的条件下， 我国政府应该在一定保证人民币汇率 在小范围的增长。对于那些出口型的产业特别是工业制品，我们应该提高产品的技术优 势，而不是通过廉价的劳动力去争取国际市场的优势。对于那些进口型的产业要从提高 国内资源的利用率和经济效率的角度出发，以减少对国际市场的依赖程度，特别是那些 矿物润油类。

32、最后要不断加强我国产业结构的调整经济制度的调整，利用人民币汇率的 小范围的增长来优化我国产业，同时来促进我国企业竞争力。五、模型评价与推广第一问用了灰色预测，对于人民币汇率的影响的因素太多，而之间的确切关系又无 法确定。因此用灰色系统可以对那些只知道一部分，而对整个其它不甚了解的状态下， 是一种比较好的预测方法。第二问用了典型变量分析的方法，对于多组变量间的关系是错综复杂的，通过其间 的组合关系可以更清楚解释各个组变量间的关联程度。第三问应用了多元回归，可以用来研究多个变量对一个变量的影响程度，对于生活 中的一些存在线性关系的变量间的关系讨论，是一种典型的方法。六、参考文献1汪晓银 周保平，数

33、学建模与数学实验，北京：科学出版社，2010。2中华人民共和国国家统计局网站 3国家外汇管理局网站 七、附录7.1 附录程序:7.1.1 问题一相关程序:基于灰色系统预测未来21年汇率的程序：function X=RmbExRata()clear all;clc;人民币汇率(年平均价)','Sheet1','B16:B31');n=length(X0);for i=1:nX0(i)=X0(i);endk=5;X=GM11(X0,k);figure;plot(1994:2009,X0,'g*-');title('人民币汇率年变化及五

34、年预测值');hold onplot(1994:2014,X);end%function X,c,error1,error2=GM11(X0,k)% 建立函数X,c,error1,error2=example9_3_2_3(X0,k)%其中X0为输入序列，k为预测长度，% X为预测输出序列，c为后验差检验数，error1为残差，error2为相对误差 format;n=length(X0);X1=;X1(1)=X0(1);for i=2:nX1(i)=X1(i-1)+X0(i);%计算累加生成序列endfor i=1:n-1B(i,1)=-0.5*(X1(i)+X1(i+1); %计算

35、 B, YnB(i,2)=1;Y(i)=X0(i+1);endB；Y;20alpha=(B'*B)A(-1)*B'*Y'%a=alpha(1,1); b=alpha(2,1); d=b/a;%c=X1-d; X2(1)=X0(1); X(1)=X0(1); for i=1:n-1X2(i+1)=c*exp(-a*i)+d;X(i+1)=X2(i+1)-X2(i);%end for i=(n+1):(n+k)X2(i)=c*exp(-a*(i-1)+d;%X(i)=X2(i)-X2(i-1);endfor i=1:nerror(i)=X(i)-X0(i);error1(i

36、)=abs(error(i); % error2(i)=error1(i)/X0(i); % end c=std(error1)/std(X0) % end做最小二乘估计计算时间响应函数参数计算预测序列计算预测序列计算残差计算相对误差计算后验差检验数基于灰色系统预测未来31天汇率的程序：function X=RmbExRata1()clear all;clc;最近一个月人民币汇率中间价,Sheet1','B24:B4');n=length(X0);for i=1:n/2temp=X0(i);X0(i)=X0(n-i+1);X0(n-i+1)=temp;endY=1 4

37、5 6 7 8 11 12 13 14 15 18 19 20 21 22 25 26 27 28 29;k=10;X=GM11(X0,k);figure;plot(Y,X0,'g*-');title('人民币汇率月变化及十天预测值);hold onY=Y 30 31 32 33 34 35 36 37 38 39;plot(Y,X);end%function X,c,erro1error2=GM11(X0,k)% 建立函数X,c,erro1error2=example9_3_2_3(X0,k)%其中X0为输入序列，k为预测长度，%X为预测输出序列，c为后验差检验数，

38、format; n=length(X0);X1=;X1(1)=X0(1); for i=2:nX1(i)=X1(i-1)+X0(i);%end for i=1:n-1 B(i,1)=-0.5*(X1(i)+X1(i+1); %B(i,2)=1;Y(i)=X0(i+1); end B；Y; alpha=(B'*B)A(-1)*B'*Y'%a=alpha(1,1); b=alpha(2,1); d=b/a;%c=X1-d; X2(1)=X0(1); X(1)=X0(1); for i=1:n-1X2(i+1)=c*exp(-a*i)+d;X(i+1)=X2(i+1)-X2(

39、i);%end for i=(n+1):(n+k)X2(i)=c*exp(-a*(i-1)+d;%X(i)=X2(i)-X2(i-1); end for i=1:nerror(i)=X(i)-X0(i);error1(i)=abs(error(i); % error2(i)=error1(i)/X0(i); % end c=std(error1)/std(X0) % enderrorl 为残差，error2为相对误差计算累加生成序列计算B, Yn做最小二乘估计计算时间响应函数参数计算预测序列计算预测序列计算残差计算相对误差计算后验差检验数7.1.2 问题二相关程序: function TSCc

40、lear all;clc;%J级产品的处理历年出口货物分类总额 历年进口货物分类总额 n=length(CJCKZE);%M级产品贸易特化系数的计算for i=1:nCJCK(i)=CJCKZE(i)-CJJKZE(i);CJZE(i)=CJCKZE(i)+CJJKZE(i); CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title('初级产品贸易特化系数)；%t品及主要供食用的产品的处理历年出口货物分类总额 历年进口货物分类总额 n=length(SPCKZE);%t品及主要供食用的贸易特化系数的计算for i=1:nSPCK

41、(i)=SPCKZE(i)-SPJKZE(i);SPZE(i尸SPCKZE(i)+SPJKZE(i); SPTSC(i)=SPCK(i)/SPZE(i);endfigure;plot(1:n,SPTSC);title('食品及主要供食用贸易特化系数);%须料及烟类的产品的处理历年出口货物分类总额 历年进口货物分类总额 n=length(YYCKZE);%须料及烟类的贸易特化系数的计算for i=1:nYYCK(i)=YYCKZE(i)-YYJKZE(i);YYZE(i)=YYCKZE(i)+YYJKZE(i); YYTSC(i)=YYCK(i)/YYZE(i);endfigure;&#

42、39;,'Sheet1','C9:C30');','Sheet1','C9:C30');','Sheet1','D9:D30');','Sheet1','D9:D30');','Sheet1','E9:E30');','Sheet1','E9:E30');plot(1:n,YYTSC);title('饮料及烟类贸易特化系数)；%E食用原料的产品的处理历年出口

43、货物分类总额历年进口货物分类总额n=length(FSPCKZE);%E食品原料的贸易特化系数的计算for i=1:nFSPCK(i)=FSPCKZE(i)-FSPJKZE(i);FSPZE(i)=FSPCKZE(i)+FSPJKZE(i);FSPTSC(i)=FSPCK(i)/FSPZE(i);endfigure;plot(1:n,FSPTSC);title(非食品原料贸易特化系数);%T物，润滑油相关原料的产品的处理历年出口货物分类总额历年进口货物分类总额n=length(KRCKZE);%T物，润滑油相关原料的贸易特化系数的计算for i=1:nKRCK(i尸KRCKZE(i)-KRJK

44、ZE(i);KRZE(i)=KRCKZE(i)+KRJKZE(i);KRTSC(i)=KRCK(i)/KRZE(i);endfigure;plot(1:n,KRTSC);title('矿物，润滑油相关原料贸易特化系数);%痴植物油脂及蜡相关原料的产品的处理历年出口货物分类总额历年进口货物分类总额n=length(DZYCKZE);%痴植物油脂及蜡的贸易特化系数的计算for i=1:nDZYCK(i)=DZYCKZE(i)-DZYJKZE(i);DZYZE(i)=DZYCKZE(i)+DZYJKZE(i);DZYTSC(i)=DZYCK(i)/DZYZE(i);endfigure;plo

45、t(1:n,DZYTSC);title('动植物油脂及蜡贸易特化系数);','Sheet1','F9:F30');','Sheet1','F9:F30');','Sheet1','G9:G30');','Sheet1','G9:G30');','Sheet1','H9:H30');','Sheet1','H9:H30');end22工业制成品的贸易特

46、化系数（TS。的程序:function TSC1clear all;clc;历年出口货物分类总额历年进口货物分类总额n=length(CJCKZE);for i=1:nCJCK(i)=CJCKZE(i)-CJJKZE(i);CJZE(i)=CJCKZE(i)+CJJKZE(i); CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title('工业制成品贸易特化系数)；历年出口货物分类总额历年进口货物分类总额n=length(CJCKZE);for i=1:nCJCK(i)=CJCKZE(i)-CJJKZE(i);CJZE(i)=CJ

47、CKZE(i)+CJJKZE(i); CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title('化学品贸易特化系数);历年出口货物分类总额历年进口货物分类总额n=length(CJCKZE);for i=1:nCJCK(i)=CJCKZE(i)-CJJKZE(i);CJZE(i)=CJCKZE(i)+CJJKZE(i); CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title(' 轻纺.橡月里.矿冶贸易特化系数);历年出口货物分类总额历年进口货物分类总额n=len

48、gth(CJCKZE);','Sheet1','I9:I30');','Sheet1','I9:I30');','Sheet1','J9:J30');','Sheet1','J9:J30');','Sheet1','K9:K30');','Sheet1','K9:K30');','Sheet1','L9:L30');&

49、#39;,'Sheet1','L9:L30');for i=1:n23-0.90920.2940478.320.32920.3707 -0.07460.6092 -0.7183-0.28120.7156532.330.44160.4513 -0.17870.3846 -0.6548-0.41690.7441551.460.45070.5016 -0.29510.1359 -0.5813-0.4391-0.40630.7194576.20.58400.5724 -0.2810-0.1722-0.4201-0.3546-0.49320.7131861.870.523

50、00.8729 -0.28620.0042 -0.5703-0.3209-0.40310.7613835.10.23760.5533 -0.39800.0196 -0.7032-0.31090.7369831.420.28670.4595 -0.4513-0.0739-0.6372-0.3420-0.21590.7385-0.0190-0.45180.1708 -0.50180.1588 -0.4655-0.0886-0.2702-0.09490.0568 -0.2527-0.0483CJCK(i)=CJCKZE(i)-CJJKZE(i);CJZE(i)=CJCKZE(i)+CJJKZE(i)

51、;CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title('机械与运输设备贸易特化系数);历年出口货物分类总额,Sheet1','M9:M30');历年进口货物分类总额,Sheet1','M9:M30');n=length(CJCKZE);for i=1:nCJCK(i)=CJCKZE(i)-CJJKZE(i);CJZE(i)=CJCKZE(i)+CJJKZE(i);CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title('杂项制品贸易特化系数);end 建立多元