《人民币汇率变动对我国进出口贸易的影响论文(修改)41500.docx》由会员分享，可在线阅读，更多相关《人民币汇率变动对我国进出口贸易的影响论文(修改)41500.docx（44页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、湖北汽车工业学学院大学生数数学建模邀请请赛承 诺诺 书我们仔细阅读了了第二届华中中地区大学生生数学建模邀邀请赛的竞赛赛细则。我们完全明白，在在竞赛开始后后参赛队员不不能以任何方方式（包括电电话、电子邮邮件、网上咨咨询等）与队队外的任何人人（包括指导导教师）研究究、讨论与赛赛题有关的问问题。我们知道，抄袭袭别人的成果果是违反竞赛赛规则的, 如果引用别别人的成果或或其他公开的的资料（包括括网上查到的的资料），必必须按照规定定的参考文献献的表述方式式在正文引用用处和参考文文献中明确列列出。我们郑重承诺，严严格遵守竞赛赛规则，以保保证竞赛的公公正、公平性性。如有违反反竞赛规则的的行为，我们们将受到严肃
2、肃处理。我们的参赛报名名号为： 050004 参赛队员 (签名) ：队员1： T933-2 杨进 队员2： T843-3 李志刚 队员3： T923-1 叶胜峰峰 第二届华中地区区大学生数学学建模邀请赛赛编 号 专专 用 页选择的题号： A 参赛的编号： 055004 （以下内容参赛赛队伍不需要要填写） 竞赛评阅阅编号： 探讨人民币汇率率变动对我国国进出口贸易易的影响【摘 要】近几十年来，世世界经济朝全全球化发展，资资本流动自由由化、国际金金融、贸易的的发展在不同同方向上增加加或降低汇率率波动对国际际贸易产生了了深远的影响响。我国经济济一直保持着着良好的发展展势头，特别别是出口贸易易的增长让一
3、一些发达国家家羡慕。人们们对中国的经经济前景比较较乐观，人民民币的国际地地位越来越高高，经济实力力的提升客观观上推动了人人民币升值。为此，我们针对对人民币汇率率的不断提高高对我国进出出口商品贸易易的影响，各各个产业结构构的变化以及及由此对我国国国内居民消消费水平的影影响做了定性性和定量的分析析。首先，我们用灰灰色系统从长长期的角度分分析了我国未未来5年的人人民币汇率发发展趋势，及及短期几天内内人民币汇率率发展趋势，这样使得我们可以对人民币汇率短期的变化情况又一个更加直观的了解。结果表明未来我国人民币升值是不可阻挡的趋势。接着，我们通过过对各个产业业的进出口贸贸易情况做分分析处理,得得出历年来我
4、我国各个产业业的贸易特化化系数.通过过贸易特化系系数的变化来来更加直观的的看出人民币币汇率的改变变对我国各个个产业的进出出口贸易的影影响.分析之之后发现人民民币汇率的改改变对我国大大部分产业都都有较大的影影响.之后我我们以各个产业的的贸易特化系系数为自变量量,人民币汇汇率为因变量量建立每个产产业与人民币币汇率之间的的多元回归分分析,定量分分析了人民币币汇率波动对对各类商品影影响程度。最后通过整体上上分析我国进进口贸易额、出出口贸易额增长率与人人民币汇率增增长率、居民民消费指数增增长率的关系系，以进出口口贸易额增长长率为自变量量,人民币汇汇率增长率和和居民消费指指数增长率为为因变量建立立典型相关
5、分析，结果表明了人人民币汇率增增长率对进口口贸易额增长长率起到决定定性的作用，对对出口贸易额额也有明显的的影响。通过以上人民币币汇率波动对对我国进出口口贸易影响的的分析及当前前我国人民币币汇率面临的的形势，为我国经济济快速平稳的的发展人民币币汇率应作如如何调整提出出一定的建议议和意见。【关键词】 人民币汇率率 灰色系系统 贸易易特化系数 实证分析析 多元回回归模型居民消费价格指指数 典型相相关模型一、问题的提出出2011年8月月24日人民民币对美元汇汇率中间价突突破6.4，再再创20055年人民币汇汇率改革以来来新高。关于人民币升值值的利弊人们们褒贬不一。有的人为人民币升值痛心疾首，有的无动于
6、衷。也有人认为人民币升值现在弊大于利，未来利大于弊，一时间，争论不断。因此，研究人民币升值对我国的影响是个热门问题。作为大学生我们们应当学会关关注我国经济济的发展，树树立强烈的主主人翁意识。针针对人民币汇汇率升高的现现象，我们从从定性和定量量的角度分析析去考虑人民民币汇率的波波动对我国进进出口贸易的的影响。我们们提出如下问问题：1、 对我国未来几年年的人民币对对美元的汇率率走势作出预预测。2、 分析一下人民币币汇率的波动动对我国进出出口商品的影影响。3、 分析一下人民币币汇率、我国国居民消费CCPI指数变变化率与我国国贸易进出口口额变化率的的关系。4、 结合上面的分析析，请说明一一下我国政府府
7、应该如何调调整人民币汇汇率政策。 二、问题的的分析人民币升值对我我国国民经济济及生活各个个方面带来多多方面的影响响。基于这样样的前提下，我我们通过所学学的知识定量量的分析了人人民币升值对对我国进出口口贸易的影响响。对于第一问，我我们从国家统统计局网站收收集到我国人人民币历年的的年均汇率，及及从国家外汇汇管理局收集集了我国近一一个月的人民民币汇率中间间价。对以上上的数据分别别用灰色系统统模型分析，对对我国人民币币汇率在未来来5年、未来来10天的波波动状况。针对第二问，我我们从国家统统计局网站收收集到我国近近几年以来的的贸易进出口口分类总额，对对初级产品、工工业制品、其其他三大类分分别计算了其其下
8、的各个分分支产业的贸贸易特化系数数，然后做出出各个产业的的贸易特化系系数波动曲线线，通过实证证分析了人民民币汇率波动动对各个不同同产业的影响响。然后根据据利用各求解解出的历年来来产业的贸易易特化系数与与人民币汇率率建立线性回回归模型，分分析人民币汇汇率的变化对对各个产业的的影响关系。对于第三问，根根据从国家统统计局网站得得到居民消费费水平、货物物进出口贸易易额、年均汇汇率等相关数数据，对居民民消费价格指指数的增长百百分比、年均均汇率的增长长百分比与进进口、出口贸贸易额增长率率建立了典型型相关分析，来来研究其间的的关系。对于第四问，我我们结合自己己所了解的知知识，及上面面的分析，对对我国汇率应应
9、该如何调整整提出了一些些建议。三、模型假设与与符号说明3.1模型假设设：（1） 假设数数据来源真实实可靠。（2） 假设选选取的指标有有一定的代表表性。（3） 假设进进出口贸易影影响中只考虑虑供我国居民民的消费指数数、人民币汇汇率的影响。（4） 假设人人民币汇率和和商品的贸易易特化系数间间存在线性关关系。（5） 假设典典型分析的两两组数据都是是连续变量，均均服从多元正正态分布。3.2符号说明明： 符号符号说明k代表产业编号代表时间灰色模型后验差差Xi代表产业的贸易易特化系数X(n)(k)利率的n次累加加序列（1000美元折算算）x代表居民消费指指数增长百分分比y1代表进口贸易增增长百分比y2z代
10、表出口贸易增增长百分比代表年均人民币币汇率增长百百分比注：其他符号见见文中说明 四、模型型建立与求解解4.1 问题一一4.1.1基于于灰色系统预预测未来5年人民币对对美元汇率的的数学模型1：分析建模 给给定原始时间间19942009年100美元元兑换成人民民币的序列XX(k),对对X(k)生生成1-AGGO序列X(k)，另外外可得Y见下下表：k123141516X861.87835.10831.42760.40694.51683.10X861.871696.9772528.39911528.44812222.99912906.009Y1696.9772528.39911528.44812222
11、.99912906.009其中，X=AGGO（X），即即 X(k)= Y= X(2)，XX(3)，X(15)，X(16)对于上述的X(k)的GMM(1,1)参数按照下下式进行最小小二乘估计得得B=将B,Y代入辨辨识算式,有有: =(BB) B Y=得GM（1,11）模型为:1) 灰微分方程： X(k)-a =b2) 白化方程： - a =b3) 白化方程的时间间相应式：(t+1)= e+4) GM(1,1)模型的检验验:残差序列的均方方差： (其中)原序列的均方差差： S (其中)后验差检验：CC=由结果可知：CC=0.388890.5,故可知知预测结果合合格。利用MATLAAB编程见附附录7
12、.2.1；运行结结果见附录77.2.1结结果。其中预预测值如下：预测未来5年人人民币对美元元汇率表格如下下:年份20102011201220132014预测值738.95222731.40555723.93559716.54226709.22449对结果的评价:由结果可以以看出1000美元兑成人人民币的数量量呈递减趋势势，即人民币币与美元的汇汇率越来越小小，人民币升升值越来越严严重。虽然该该预测与现在在的实际情况况不符，但其其趋势是相同同的。其可能能原因为：以以年为时间单单位误差太大大、预测的时间间跨度太大或或原始数据过过去的时间太太久对现在的的影响很小。为为此，我们建建立如下模型型来预测短期
13、期情况。4.1.2基于于灰色系统预预测未来100天人民币对对美元汇率的数学模模型：分析建模 给给定原始时间间2011.7.12011.7.29内内 100美美元兑换成人人民币3的序列X(k),对X(kk)生成1-AGO序列列X(k)，另另外可得Y见见下表：k123192021X646.85646.61646.5644.26644.38644.42X646.851293.4661939.96612280.66412924.66213569.004Y1293.4661939.96612280.66412924.66213569.004理论知识与4.1.1中的的一样，在这这里不再累赘赘。由结果可知：
14、CC=0.300890.35,故可知知预测结果优优。利用MATLAAB编程见附附录7.2.1；运行结结果见附录77.2.1结结果。其中预预测值如下：日起7/30 7/318/18/28/38/48/58/68/78/8预测值644.40644.24644.08643.91643.75643.59643.43643.26643.10642.94（备注：预测值值为100美美元可以兑换换人民币数量量（元）对结果的评价:由预测结果果可以看出人人民币兑美元元的汇率虽然然存在大小波波动，但其整整体趋势是不不断减小的。预预测值也与现现实相符，说说明人民币升升值进一步加加剧。4.2 问题二二4.2.1首先先根
15、据10个个产业的贸易易进出口额状状况，作出各各个产业的贸贸易特化系数数（TSC），然然后作出其相相应的曲线变变化图，对其其进行一定的的说明。问题分析：将我我国基础产业业分为：初级级产品、工业业制品、其他他。这里只对对初级产品、工工业制品分析析（收集各产产品的进出口口额）。出口口贸易商品结结构有不同层层次，在这里里我们采用了了贸易特化系系数的变化趋趋势来阐述某某种贸易商品品或贸易产业业的结构优化化。贸易特化化系数反应的的是一国某贸贸易产业或贸贸易产品的净净出口额与该该产品的进出出口总额之比比。其计算公公式为： TSSC=(X-M)/(XX+M) -1TSC1当TSC指标值值介于-1到到-0.3时
16、时，则代表该该国产业的进进口值远大于于出口值，故故该产业属于于“进口强势”型产业，即即该产业在国国际市场上占占有绝对的比比较优势；若TSC指标值值介于-0.3与0.33之间，表示示其进口与出出口值差异不不大，代表该该产业在国际际上水平较高高，属于“水平竞争”型产业；当TSC指标值值介于0.331时，表示示其出口值远远大于其进口口值，故该产产业在该国处处于绝对的劣劣势，属于“出口强势”型产业。因此，从各种贸贸易行业贸易易特化系数的的变化我们可可以大概地看看出我国贸易易结构调整的的方向。将所有进出口产产品分成两类类：初级产品：食品品及主要供食食用的活动物物、饮料及烟烟类、非食用用原料、矿物物燃料、
17、动植物、油脂脂及蜡工业制成品：化化学品及有关关产品、轻纺纺产品、橡胶胶制品、矿冶冶产品及其制制品、机 械械及运输设备备、杂项制品品MATLAB编编程见附录77.2.1，结结果如下：1）初级产品及及其组成项的的图像：从上面可以看出出我国初级产产品贸易特化化曲线持续下下降，说明了了该产业主要要是通过进口口，属于强进口口型，其对国国外市场依赖赖过强，因此此人民币汇率率升值是有利利于该产业发发展的。从上面可以看出出该产业的贸贸易特化曲线线基本处于正正的，也就是是说，主要是是向国外出口，属于强强出口口型，人人民币升值会会不利于此类商品品在国外的发发展。从上面曲线可以以看出特化曲曲线后一阶段段处于不断下下
18、降的趋势，该该产业正从水水平竞争型向向进口型转化的的趋势，人民民币升值对其其发展是有利利。从此图可以看出出非食品原料料贸易特化系系数最后将会会向-1趋近近，说明了非非食品原料贸贸易几乎全部部依靠进口来来提供，若人人民币升值，会会有利于我国国进口非食品品原料。上图反映了矿物物润滑油相关关原料的特化化系数逐渐逼逼近-1，说说明了我国的的矿物润滑油油相关原料主主要依靠进口口，人民币升升值将有利于于我国对此产产品的进口与与需求。 从上面可以看出出我国动植物物油脂及蜡的的贸易特化系系数逐步趋向向-1，说明明了其进口趋趋势显著加强强，人民币升升值有利于其其发展。由结果得出如下下分析：初级级产品、饮料及烟类
19、类、非食用原原料、矿物燃燃料、动植物物油脂及蜡的贸贸易特化系数数均随着人民民币的升值而而减小，食品品及主要供食食用物的贸易特化化系数虽然有有上下波动但但整体比较稳稳定，由此可可见人民币升升值对我国的的进出口贸易易有着密切影影响。使我国的初级产产品类产业的进口口值远大于出出口值，即“进口强势”型产业越来来越多。而食食品的波动比比较平稳，可可能受国家相相关政策的影影响，使我国国的食品价格格保持平稳。有有利于人民的的生活和社会会的稳定。2）工业制成品品及其组成项项的图像：从上面图像可以以看出我国工工业制成品贸贸易特化曲线线基本处于-0.3到00.3之间，说说明了此类产产业属于水平平竞争型，人人民币汇
20、率的的升高对其影影响不大。但但是我们也可可以看到在汇汇率升高的同同时其贸易特特化系数也是是呈上升趋势势，如果人民民币持续升值值，该产业又又转化成为出出口强势型的的产业。从上图可以看出出，化学品相相关产业的贸贸易特化系数数基本上处于于-0.3一一下，即便是是又波动影响响也不大，可可以认为该产产业是进口强强型产业，在人人民币汇率升升高有利于该该产业的发展展。从下图可以看出出：我国的轻轻纺橡胶、矿矿冶贸易特化化曲线，在前前期处于进口口强势到了后后期逐渐转变变为出口强势势产业，人民民币的汇率提提高不利于其其再国际市场场的竞争力。从此图可以看出出我国机械运运输设备贸易易特化曲线一一直处于上升升的趋势，说
21、说明了该产业业将会向出口口型发展，人人民币的升值值将会不利于于其发展。由结果得出如下下分析：工业业制成品、化化学品及有关关产品、轻纺纺产品、橡胶胶制品、矿冶冶产品及其制制品、机械及及运输设备的的贸易特化系系数均随着人人民币的升值值而增大，杂杂项制品的贸贸易特化系数数虽然有上下下波动但整体体比较稳定，由由此可见人民民币升值对我我国的进出口口贸易有着密密切影响。是是我国的工业类类产业的进口口值远小于出口值，即即“出口弱势”型产业越来来越多。对比比初级产品和工业制成品品可知：我国国的工业水品品低下，且只只能用缺少技技术含量的初初级产品和廉廉价的劳动力力来弥补进出出口差额，国国家应加大科科技的发展，尽
22、尽快使工业水水平赶上来，加加大工业制成成品的出口,减少少因人民币升升值而带来的的影响。4.2.2 多多元回归分析析建模过程：1） 模型的建立:设自变量的观测测值及因变量量y对应的观观测值y满足足关系式： yy= i=1,22,3n；j=11,2,3,p其中，是相互独独立且都服从从正态分布的的随机变量。由最小二乘法，得得 2）回归方程的的显著性检验验：利用F分分布函数对变变量前的系数数检验。3）回归系数的的显著性检验验： 利用tt分布函数检检验。最后得到结果如如下：产业名称编号方程显著性检验验（F）参数显著性检验（t）回归方程备注食品及主要供食食用的活动物物10.02880.01080.0288
23、y=773.15277x1+389.7118汇率对其影响显显著饮料及烟酒20.1440(不合格)-汇率与此行业没没有显著关系系非食用原料30.00590.0010.0059y=-4193.92239x3+480.43293汇率对其影响显显著矿物燃料40.000110.000110.00011y=-309.74464x4+657.01952汇率对其影响显显著动植物油脂50.0009660.00420.000966y=-458.67428x5+382.11320汇率对其影响显显著化学品60.1256(不合格)-汇率对其影响不不显著轻纺橡胶7不合格-汇率对其影响不不显著机械及运输80.000110.
24、000110.00011y=486.07219x8+791.92968汇率对其影响显显著杂项90.2222(不合格)-汇率对其影响不不显著（备注：y代表表年均人民币币汇率，1000美元兑换换的人民币数数量(元)，xxi代表贸易特特化系数）由上表看出，人人民币汇率的的变化对我国国的大部分产产业都有显著著的影响。也也就是说，人人民币持续升升值对我国各各行各业都会会产生决定性性改变，因此此，为了稳定定中国国民经经济稳步发展展，在人民币币升值过程中中国家需要进进行干预，严严格控制人民民币的升值幅幅度，这样才才能保证中国国经济快速稳稳步发展，中中国人民能更更美好的生活活。4.3 问题三三典型相关分析思思
25、想：考察一组随机变变量和另一组组随机变量之之间的相关性性，可以用他他们之间的相相关性。考虑虑到两组变量量的线性组合合：u1=a11x1+a21x2+an1xnv1=b11y1+b21y2+am1ym并研究他们之间间的相关系数数（u1，v1），在所有的的线性组合中中找一对相关关系数最大的的线性组合，用用这两个组合合的单相关系系数来表示两两组变量的相相关性，这两两个线性组合合叫做一对典典型变量。下一步，再在两两组比昂两的的与u，v不不相关的线性性组合中，找找一对相关系系数最大的线线性组合，设设为：u2=a12x1+a22x2+an2xnv2=b12y1+b22y2+am2ym且（u2，v2）就是第
26、二个典典型相关系数数，这样下去去可以得到若若干对典型变变量，从而提提取出两组变变量间的全部部信息。就本题而言，用用到典型相关关分析。利用用SAS编程程得到了：人民币汇率增长长率、居民消消费指数增长长率与进出口口贸易额增长长率间的两个个典型相关系系数及特征值值如下表： 序号典型相关系数标准误差特征值方差比率10.95212250.01992259.70010.935720.63243310.12792270.66660.0643从上表可以看出出：第一对典典型变量的相相关系数较高高，说明了这这一对相关变变量间的关系系是密切的。对相关系数检验验如下表：序号F计算值自由度F检验的显著性性概率130.6
27、240.001213.3310.0016这两对典型相关关变量，进行行F统计检验验，均通过了了检验，说明明了这两对典典型变量间的的相关性显著著，可以用来来解释产出变变量。对冗余度分析得得到下表：典型变量的解释释能力序号产出组方差被影影响组解释比比例对产出组的解释释能力产出方差被典型型变量解释比比例对影响组的解释释能力影响组方差被典典型变量的解解释比例10.90650.95190.86290.87600.794120.40000.4810.01930.12400.0496从此表可以看出出：这两对典典型相关变量量的解释能力力均较强，特特别是第一对对相关变量对对产出组和影影响组均有很很高的解释能能力；
28、第一对对典型变量对对产出组方差差有95%以以上的解释能能力，第二对对典型相关变变量较弱些，达达到了48%左右。最后得出的典型型变量模型：序号典型相关模型1u1=0.0546x+0.8714zw1=0.7245y1+0.8619y22u2=0.4305x+0.8869zw2=0.8048y1+0.8816y2（x代表居民消消费指数增长长百分比、zz代表人民币币汇率增长百百分比；y11代表出口贸贸易额增长率率、y2代表表进口贸易额额的增长率；）结果分析：人民币汇率增长长率对进口贸贸易有显著的的作用。从上上面第一对典典型变量分析析结果可以看看出：人民币币汇率增长百百分比在典型型变量u1占占得比重达到
29、到了0.87714,出口口贸易额增长长率在w1中中占的比重达达到0.86619.由于于u1,w11由前面分析析知道其相关关系数达到了了0.95以以上，可以认认为人民币汇汇率的增长率率进口贸易额额增长率变化化起到决定性性的作用。从第二对典型相相关变量分析析知道：人民民币汇率的增增长率、居民民消费指数增增长率对进出出口贸易额增增长率都一定定的影响，不不过期间的影影响程度不是是那么显著。综上所述：人民民币汇率的增增长率对我国国进口贸易额额有着决定性性的因素，同同时它和居民民消费指数增增长率对我国国进出口额增增长率都有密密切的关联。4.4 问题四四通过上面的分析析在人民币不不断升值的条条件下，我国国政
30、府应该在在一定保证人人民币汇率在在小范围的增增长。对于那那些出口型的的产业特别是是工业制品，我我们应该提高高产品的技术术优势，而不不是通过廉价价的劳动力去去争取国际市市场的优势。对对于那些进口口型的产业要要从提高国内内资源的利用用率和经济效效率的角度出出发，以减少少对国际市场场的依赖程度度，特别是那那些矿物润油油类。最后要要不断加强我我国产业结构构的调整经济济制度的调整整，利用人民民币汇率的小小范围的增长长来优化我国国产业，同时时来促进我国国企业竞争力力。五、模型评价与与推广第一问用了灰色色预测，对于于人民币汇率率的影响的因因素太多，而而之间的确切切关系又无法法确定。因此此用灰色系统统可以对那
31、些些只知道一部部分，而对整整个其它不甚甚了解的状态态下，是一种种比较好的预预测方法。第二问用了典型型变量分析的的方法，对于于多组变量间间的关系是错错综复杂的，通通过其间的组组合关系可以以更清楚解释释各个组变量量间的关联程程度。第三问应用了多多元回归，可可以用来研究究多个变量对对一个变量的的影响程度，对对于生活中的的一些存在线线性关系的变变量间的关系系讨论，是一一种典型的方方法。六、参考文献1 汪晓银银 周保平，数数学建模与数数学实验，北北京：科学出出版社，20010。2 中华人人民共和国国国家统计局网网站httpp:/3 国家外外汇管理局网网站httpp:/ 七、附录录7.1 附录程程序:7.
32、1.1 问问题一相关程程序: 基于灰色系统预预测未来211年汇率的程程序：functioon X=RmbEExRataa()clear aall;clc;X0=xlsrread(E:人民民币汇率(年年平均价),Sheeet1,B16:B31);n=lengtth(X0);for i=11:n X0(i)=X00(i);endk=5;X=GM111(X0,k);figure;plot(19994:20009,X00,g*-);title(人民币汇率率年变化及五五年预测值);hold onn plot(19994:20014,X);end%functioon X,c,errror1,eerror2
33、2=GM111(X0,k)% 建立函数X,c,eerror11,erroor2=eexamplle9_3_2_3(XX0,k)% 其中X0为为输入序列，kk为预测长度度，% X为预测输输出序列，cc为后验差检检验数，errror1为为残差，errror2为为相对误差format;n=lengtth(X0);X1=;X1(1)=XX0(1);for i=22:n X1(i)=X11(i-1)+X0(ii); %计算累加加生成序列endfor i=11:n-1B(i,1)=-0.5*(X1(ii)+X1(i+1); %计算B，YYn B(ii,2)=11; Y(ii)=X0(i+1);endB;Y
34、;alpha=(B*B)(-1)*B*YY; %做最小二乘乘估计a=alphaa(1,1);b=alphaa(2,1);d=b/a; %计计算时间响应应函数参数c=X1(1)-d;X2(1)=XX0(1);X(1)=X00(1);for i=11:n-1 X2(i+1)=c*expp(-a*ii)+d; X(ii+1)=XX2(i+11)-X2(i); %计算预测序序列endfor i=(n+1):(n+k) X2(i)=c*exp(-a*(i-1)+dd; %计计算预测序列列 X(ii)=X2(i)-X22(i-1);endfor i=11:n errror(i)=X(i)-X0(ii);
35、errror1(ii)=abss(erroor(i); %计计算残差 errror2(ii)=errror1(ii)/X0(i); %计计算相对误差差endc=std(eerror11)/stdd(X0) %计算后验差差检验数end 基于灰色系统预预测未来311天汇率的程程序：functioon X=RmbEExRataa1()clear aall;clc;X0=xlsrread(E:最近近一个月人民民币汇率中间间价,SSheet11,B224:B4);n=lengtth(X0);for i=11:n/2 temmp=X0(i); X0(i)=X00(n-i+1); X0(n-i+11)=te
36、mmp;endY=1 4 5 6 77 8 111 12 113 14 15 188 19 220 21 22 255 26 227 28 29;k=10;X=GM111(X0,k);figure;plot(Y,X0,gg*-);title(人民币汇率率月变化及十十天预测值);hold onnY=Y 300 31 332 33 34 355 36 337 38 39;plot(Y,X);end%functioon X,c,errror1,eerror22=GM111(X0,k)% 建立函数X,c,eerror11,erroor2=eexamplle9_3_2_3(XX0,k)% 其中X0为为输
37、入序列，kk为预测长度度，% X为预测输输出序列，cc为后验差检检验数，errror1为为残差，errror2为为相对误差format;n=lengtth(X0);X1=;X1(1)=XX0(1);for i=22:n X1(i)=X11(i-1)+X0(ii); %计算累加加生成序列endfor i=11:n-1B(i,1)=-0.5*(X1(ii)+X1(i+1); %计算B，YYn B(ii,2)=11; Y(ii)=X0(i+1);endB;Y;alpha=(B*B)(-1)*B*YY; %做最小二乘乘估计a=alphaa(1,1);b=alphaa(2,1);d=b/a; %计计算时
38、间响应应函数参数c=X1(1)-d;X2(1)=XX0(1);X(1)=X00(1);for i=11:n-1 X2(i+1)=c*expp(-a*ii)+d; X(ii+1)=XX2(i+11)-X2(i); %计算预测序序列endfor i=(n+1):(n+k) X2(i)=c*exp(-a*(i-1)+dd; %计计算预测序列列 X(ii)=X2(i)-X22(i-1);endfor i=11:n errror(i)=X(i)-X0(ii); errror1(ii)=abss(erroor(i); %计计算残差 errror2(ii)=errror1(ii)/X0(i); %计计算相对
39、误差差endc=std(eerror11)/stdd(X0) %计算后验差差检验数end7.1.2 问问题二相关程程序: 初级产品的的贸易特化系系数（TSCC）的程序：functioon TSCCclear aall;clc;%初级产品的的处理CJCKZE=xlsreead(EE:历年出出口货物分类类总额,Sheett1,CC9:C300);CJJKZE=xlsreead(EE:历年进进口货物分类类总额,Sheett1,CC9:C300);n=lengtth(CJCCKZE);%初级产品贸贸易特化系数数的计算for i=11:n CJCCK(i)=CJCKZZE(i)-CJJKZZE(i);
40、CJZZE(i)=CJCKZZE(i)+CJJKZZE(i); CJTTSC(i)=CJCKK(i)/CCJZE(ii);endfigure;plot(1:n,CJTTSC);title(初级产品贸贸易特化系数数);%食品及主要要供食用的产产品的处理SPCKZE=xlsreead(EE:历年出出口货物分类类总额,Sheett1,DD9:D300);SPJKZE=xlsreead(EE:历年进进口货物分类类总额,Sheett1,DD9:D300);n=lengtth(SPCCKZE);%食品及主要要供食用的贸贸易特化系数数的计算for i=11:n SPCCK(i)=SPCKZZE(i)-SPJ
41、KZZE(i); SPZZE(i)=SPCKZZE(i)+SPJKZZE(i); SPTTSC(i)=SPCKK(i)/SSPZE(ii);endfigure;plot(1:n,SPTTSC);title(食品及主要要供食用贸易易特化系数);%饮料及烟类类的产品的处处理YYCKZE=xlsreead(EE:历年出出口货物分类类总额,Sheett1,EE9:E300);YYJKZE=xlsreead(EE:历年进进口货物分类类总额,Sheett1,EE9:E300);n=lengtth(YYCCKZE);%饮料及烟类类的贸易特化化系数的计算算for i=11:n YYCCK(i)=YYCKZZE(i)-YYJKZZE(i); YYZZE(i)=YYCKZZE(i)+YYJKZZE(i); YYTTSC(i)=YYCKK(i)/YYYZE(ii);endfigure;plot(1: