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

《人民币汇率变动对我国进出口贸易的影响论文(修改)(27页DOC).docx》由会员分享，可在线阅读，更多相关《人民币汇率变动对我国进出口贸易的影响论文(修改)(27页DOC).docx（33页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

2、3-2 杨进 队员2： T843-3 李志刚 队员3： T923-1 叶胜峰 第二届华中地区大学生数学建模邀请赛编 号 专 用 页选择的题号： A 参赛的编号： 05004 （以下内容参赛队伍不需要填写） 竞赛评阅编号： 最新精品资料整理推荐，更新于二二年十二月二十八日2020年12月28日星期一21:15:39探讨人民币汇率变动对我国进出口贸易的影响【摘 要】近几十年来，世界经济朝全球化发展，资本流动自由化、国际金融、贸易的发展在不同方向上增加或降低汇率波动对国际贸易产生了深远的影响。我国经济一直保持着良好的发展势头，特别是出口贸易的增长让一些发达国家羡慕。人们对中国的经济前景比较乐观，人民

3、币的国际地位越来越高，经济实力的提升客观上推动了人民币升值。为此，我们针对人民币汇率的不断提高对我国进出口商品贸易的影响，各个产业结构的变化以及由此对我国国内居民消费水平的影响做了定性和定量的分析。首先，我们用灰色系统从长期的角度分析了我国未来5年的人民币汇率发展趋势，及短期几天内人民币汇率发展趋势，这样使得我们可以对人民币汇率短期的变化情况又一个更加直观的了解。结果表明未来我国人民币升值是不可阻挡的趋势。接着，我们通过对各个产业的进出口贸易情况做分析处理,得出历年来我国各个产业的贸易特化系数.通过贸易特化系数的变化来更加直观的看出人民币汇率的改变对我国各个产业的进出口贸易的影响.分析之后发现

4、人民币汇率的改变对我国大部分产业都有较大的影响.之后我们以各个产业的贸易特化系数为自变量,人民币汇率为因变量建立每个产业与人民币汇率之间的多元回归分析,定量分析了人民币汇率波动对各类商品影响程度。最后通过整体上分析我国进口贸易额、出口贸易额增长率与人民币汇率增长率、居民消费指数增长率的关系，以进出口贸易额增长率为自变量,人民币汇率增长率和居民消费指数增长率为因变量建立典型相关分析，结果表明了人民币汇率增长率对进口贸易额增长率起到决定性的作用，对出口贸易额也有明显的影响。通过以上人民币汇率波动对我国进出口贸易影响的分析及当前我国人民币汇率面临的形势，为我国经济快速平稳的发展人民币汇率应作如何调整

5、提出一定的建议和意见。【关键词】 人民币汇率 灰色系统 贸易特化系数 实证分析 多元回归模型居民消费价格指数 典型相关模型一、问题的提出2011年8月24日人民币对美元汇率中间价突破6.4，再创2005年人民币汇率改革以来新高。关于人民币升值的利弊人们褒贬不一。有的人为人民币升值痛心疾首，有的无动于衷。也有人认为人民币升值现在弊大于利，未来利大于弊，一时间，争论不断。因此，研究人民币升值对我国的影响是个热门问题。作为大学生我们应当学会关注我国经济的发展，树立强烈的主人翁意识。针对人民币汇率升高的现象，我们从定性和定量的角度分析去考虑人民币汇率的波动对我国进出口贸易的影响。我们提出如下问题：1、

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

7、10天的波动状况。针对第二问，我们从国家统计局网站收集到我国近几年以来的贸易进出口分类总额，对初级产品、工业制品、其他三大类分别计算了其下的各个分支产业的贸易特化系数，然后做出各个产业的贸易特化系数波动曲线，通过实证分析了人民币汇率波动对各个不同产业的影响。然后根据利用各求解出的历年来产业的贸易特化系数与人民币汇率建立线性回归模型，分析人民币汇率的变化对各个产业的影响关系。对于第三问，根据从国家统计局网站得到居民消费水平、货物进出口贸易额、年均汇率等相关数据，对居民消费价格指数的增长百分比、年均汇率的增长百分比与进口、出口贸易额增长率建立了典型相关分析，来研究其间的关系。对于第四问，我们结合自

8、己所了解的知识，及上面的分析，对我国汇率应该如何调整提出了一些建议。三、模型假设与符号说明3.1模型假设：（1） 假设数据来源真实可靠。（2） 假设选取的指标有一定的代表性。（3） 假设进出口贸易影响中只考虑供我国居民的消费指数、人民币汇率的影响。（4） 假设人民币汇率和商品的贸易特化系数间存在线性关系。（5） 假设典型分析的两组数据都是连续变量，均服从多元正态分布。3.2符号说明： 符号符号说明k代表产业编号代表时间灰色模型后验差Xi代表产业的贸易特化系数X(n)(k)利率的n次累加序列（100美元折算）x代表居民消费指数增长百分比y1代表进口贸易增长百分比y2z代表出口贸易增长百分比代表年

9、均人民币汇率增长百分比注：其他符号见文中说明 四、模型建立与求解4.1 问题一4.1.1基于灰色系统预测未来5年人民币对美元汇率的数学模型1：分析建模 给定原始时间19942009年100美元兑换成人民币的序列X(k),对X(k)生成1-AGO序列X(k)，另外可得Y见下表：k123141516X861.87835.10831.42760.40694.51683.10X861.871696.972528.3911528.4812222.9912906.09Y1696.972528.3911528.4812222.9912906.09其中，X=AGO（X），即 X(k)= Y= X(2)，X(3

10、)，X(15)，X(16)对于上述的X(k)的GM(1,1)参数按照下式进行最小二乘估计得B=将B,Y代入辨识算式,有: =(BB) B Y=得GM（1,1）模型为:1) 灰微分方程： X(k)-a =b2) 白化方程： - a =b3) 白化方程的时间相应式：(t+1)= e+4) GM(1,1)模型的检验:残差序列的均方差： (其中)原序列的均方差： S (其中)后验差检验：C=由结果可知：C=0.38890.5,故可知预测结果合格。利用MATLAB编程见附录7.2.1；运行结果见附录7.2.1结果。其中预测值如下：预测未来5年人民币对美元汇率表格如下:年份2010201120122013

11、2014预测值738.9522731.4055723.9359716.5426709.2249对结果的评价:由结果可以看出100美元兑成人民币的数量呈递减趋势，即人民币与美元的汇率越来越小，人民币升值越来越严重。虽然该预测与现在的实际情况不符，但其趋势是相同的。其可能原因为：以年为时间单位误差太大、预测的时间跨度太大或原始数据过去的时间太久对现在的影响很小。为此，我们建立如下模型来预测短期情况。4.1.2基于灰色系统预测未来10天人民币对美元汇率的数学模型：分析建模 给定原始时间2011.7.12011.7.29内 100美元兑换成人民币3的序列X(k),对X(k)生成1-AGO序列X(k)，

12、另外可得Y见下表：k123192021X646.85646.61646.5644.26644.38644.42X646.851293.461939.9612280.6412924.6213569.04Y1293.461939.9612280.6412924.6213569.04理论知识与4.1.1中的一样，在这里不再累赘。由结果可知：C=0.30890.35,故可知预测结果优。利用MATLAB编程见附录7.2.1；运行结果见附录7.2.1结果。其中预测值如下：日起7/30 7/318/18/28/38/48/58/68/78/8预测值644.40644.24644.08643.91643.75

13、643.59643.43643.26643.10642.94（备注：预测值为100美元可以兑换人民币数量（元）对结果的评价:由预测结果可以看出人民币兑美元的汇率虽然存在大小波动，但其整体趋势是不断减小的。预测值也与现实相符，说明人民币升值进一步加剧。4.2 问题二4.2.1首先根据10个产业的贸易进出口额状况，作出各个产业的贸易特化系数（TSC），然后作出其相应的曲线变化图，对其进行一定的说明。问题分析：将我国基础产业分为：初级产品、工业制品、其他。这里只对初级产品、工业制品分析（收集各产品的进出口额）。出口贸易商品结构有不同层次，在这里我们采用了贸易特化系数的变化趋势来阐述某种贸易商品或贸易

14、产业的结构优化。贸易特化系数反应的是一国某贸易产业或贸易产品的净出口额与该产品的进出口总额之比。其计算公式为： TSC=(X-M)/(X+M) -1TSC1当TSC指标值介于-1到-0.3时，则代表该国产业的进口值远大于出口值，故该产业属于“进口强势”型产业，即该产业在国际市场上占有绝对的比较优势；若TSC指标值介于-0.3与0.3之间，表示其进口与出口值差异不大，代表该产业在国际上水平较高，属于“水平竞争”型产业；当TSC指标值介于0.31时，表示其出口值远大于其进口值，故该产业在该国处于绝对的劣势，属于“出口强势”型产业。因此，从各种贸易行业贸易特化系数的变化我们可以大概地看出我国贸易结构

15、调整的方向。将所有进出口产品分成两类：初级产品：食品及主要供食用的活动物、饮料及烟类、非食用原料、矿物燃料、动植物、油脂及蜡工业制成品：化学品及有关产品、轻纺产品、橡胶制品、矿冶产品及其制品、机 械及运输设备、杂项制品MATLAB编程见附录7.2.1，结果如下：1）初级产品及其组成项的图像：从上面可以看出我国初级产品贸易特化曲线持续下降，说明了该产业主要是通过进口，属于强进口型，其对国外市场依赖过强，因此人民币汇率升值是有利于该产业发展的。从上面可以看出该产业的贸易特化曲线基本处于正的，也就是说，主要是向国外出口，属于强出口口型，人民币升值会不利于此类商品在国外的发展。从上面曲线可以看出特化曲

16、线后一阶段处于不断下降的趋势，该产业正从水平竞争型向进口型转化的趋势，人民币升值对其发展是有利。从此图可以看出非食品原料贸易特化系数最后将会向-1趋近，说明了非食品原料贸易几乎全部依靠进口来提供，若人民币升值，会有利于我国进口非食品原料。上图反映了矿物润滑油相关原料的特化系数逐渐逼近-1，说明了我国的矿物润滑油相关原料主要依靠进口，人民币升值将有利于我国对此产品的进口与需求。 从上面可以看出我国动植物油脂及蜡的贸易特化系数逐步趋向-1，说明了其进口趋势显著加强，人民币升值有利于其发展。由结果得出如下分析：初级产品、饮料及烟类、非食用原料、矿物燃料、动植物油脂及蜡的贸易特化系数均随着人民币的升值

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

18、学品相关产业的贸易特化系数基本上处于-0.3一下，即便是又波动影响也不大，可以认为该产业是进口强型产业，在人民币汇率升高有利于该产业的发展。从下图可以看出：我国的轻纺橡胶、矿冶贸易特化曲线，在前期处于进口强势到了后期逐渐转变为出口强势产业，人民币的汇率提高不利于其再国际市场的竞争力。从此图可以看出我国机械运输设备贸易特化曲线一直处于上升的趋势，说明了该产业将会向出口型发展，人民币的升值将会不利于其发展。由结果得出如下分析：工业制成品、化学品及有关产品、轻纺产品、橡胶制品、矿冶产品及其制品、机械及运输设备的贸易特化系数均随着人民币的升值而增大，杂项制品的贸易特化系数虽然有上下波动但整体比较稳定，

19、由此可见人民币升值对我国的进出口贸易有着密切影响。是我国的工业类产业的进口值远小于出口值，即“出口弱势”型产业越来越多。对比初级产品和工业制成品可知：我国的工业水品低下，且只能用缺少技术含量的初级产品和廉价的劳动力来弥补进出口差额，国家应加大科技的发展，尽快使工业水平赶上来，加大工业制成品的出口,减少因人民币升值而带来的影响。4.2.2 多元回归分析建模过程：1） 模型的建立:设自变量的观测值及因变量y对应的观测值y满足关系式： y= i=1,2,3n；j=1,2,3,p其中，是相互独立且都服从正态分布的随机变量。由最小二乘法，得 2）回归方程的显著性检验：利用F分布函数对变量前的系数检验。3

20、）回归系数的显著性检验： 利用t分布函数检验。最后得到结果如下：产业名称编号方程显著性检验（F）参数显著性检验（t）回归方程备注食品及主要供食用的活动物10.02880.01080.0288y=773.15277x1+389.7118汇率对其影响显著饮料及烟酒20.1440(不合格)-汇率与此行业没有显著关系非食用原料30.00590.0010.0059y=-4193.92239x3+480.43293汇率对其影响显著矿物燃料40.00010.00010.0001y=-309.74464x4+657.01952汇率对其影响显著动植物油脂50.000960.00420.00096y=-458.6

21、7428x5+382.11320汇率对其影响显著化学品60.1256(不合格)-汇率对其影响不显著轻纺橡胶7不合格-汇率对其影响不显著机械及运输80.00010.00010.0001y=486.07219x8+791.92968汇率对其影响显著杂项90.2222(不合格)-汇率对其影响不显著（备注：y代表年均人民币汇率，100美元兑换的人民币数量(元)，xi代表贸易特化系数）由上表看出，人民币汇率的变化对我国的大部分产业都有显著的影响。也就是说，人民币持续升值对我国各行各业都会产生决定性改变，因此，为了稳定中国国民经济稳步发展，在人民币升值过程中国家需要进行干预，严格控制人民币的升值幅度，这样

22、才能保证中国经济快速稳步发展，中国人民能更美好的生活。4.3 问题三典型相关分析思想：考察一组随机变量和另一组随机变量之间的相关性，可以用他们之间的相关性。考虑到两组变量的线性组合：u1=a11x1+a21x2+an1xnv1=b11y1+b21y2+am1ym并研究他们之间的相关系数（u1，v1），在所有的线性组合中找一对相关系数最大的线性组合，用这两个组合的单相关系数来表示两组变量的相关性，这两个线性组合叫做一对典型变量。下一步，再在两组比昂两的与u，v不相关的线性组合中，找一对相关系数最大的线性组合，设为：u2=a12x1+a22x2+an2xnv2=b12y1+b22y2+am2ym且

23、（u2，v2）就是第二个典型相关系数，这样下去可以得到若干对典型变量，从而提取出两组变量间的全部信息。就本题而言，用到典型相关分析。利用SAS编程得到了：人民币汇率增长率、居民消费指数增长率与进出口贸易额增长率间的两个典型相关系数及特征值如下表： 序号典型相关系数标准误差特征值方差比率10.9521250.0199259.70010.935720.6324310.1279270.66660.0643从上表可以看出：第一对典型变量的相关系数较高，说明了这一对相关变量间的关系是密切的。对相关系数检验如下表：序号F计算值自由度F检验的显著性概率130.6240.001213.3310.0016这两对

24、典型相关变量，进行F统计检验，均通过了检验，说明了这两对典型变量间的相关性显著，可以用来解释产出变量。对冗余度分析得到下表：典型变量的解释能力序号产出组方差被影响组解释比例对产出组的解释能力产出方差被典型变量解释比例对影响组的解释能力影响组方差被典型变量的解释比例10.90650.95190.86290.87600.794120.40000.4810.01930.12400.0496从此表可以看出：这两对典型相关变量的解释能力均较强，特别是第一对相关变量对产出组和影响组均有很高的解释能力；第一对典型变量对产出组方差有95%以上的解释能力，第二对典型相关变量较弱些，达到了48%左右。最后得出的典

25、型变量模型：序号典型相关模型1u1=0.0546x+0.8714zw1=0.7245y1+0.8619y22u2=0.4305x+0.8869zw2=0.8048y1+0.8816y2（x代表居民消费指数增长百分比、z代表人民币汇率增长百分比；y1代表出口贸易额增长率、y2代表进口贸易额的增长率；）结果分析：人民币汇率增长率对进口贸易有显著的作用。从上面第一对典型变量分析结果可以看出：人民币汇率增长百分比在典型变量u1占得比重达到了0.8714,出口贸易额增长率在w1中占的比重达到0.8619.由于u1,w1由前面分析知道其相关系数达到了0.95以上，可以认为人民币汇率的增长率进口贸易额增长率

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

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

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

29、014,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:n X1(i)=X1(i-1)+X0(i); %计算累加生成序列endfor i=1:n-1B(i,1)=-0.5*(X1(i)+X1(i+1); %计算B，Yn B(i,2)=1; Y(i)=X0(i+1)

30、;endB;Y;alpha=(B*B)(-1)*B*Y; %做最小二乘估计a=alpha(1,1);b=alpha(2,1);d=b/a; %计算时间响应函数参数c=X1(1)-d;X2(1)=X0(1);X(1)=X0(1);for i=1:n-1 X2(i+1)=c*exp(-a*i)+d; X(i+1)=X2(i+1)-X2(i); %计算预测序列endfor i=(n+1):(n+k) X2(i)=c*exp(-a*(i-1)+d; %计算预测序列 X(i)=X2(i)-X2(i-1);endfor i=1:n error(i)=X(i)-X0(i); error1(i)=abs(er

31、ror(i); %计算残差 error2(i)=error1(i)/X0(i); %计算相对误差endc=std(error1)/std(X0) %计算后验差检验数end 基于灰色系统预测未来31天汇率的程序：function X=RmbExRata1()clear all;clc;X0=xlsread(E:最近一个月人民币汇率中间价,Sheet1,B24:B4);n=length(X0);for i=1:n/2 temp=X0(i); X0(i)=X0(n-i+1); X0(n-i+1)=temp;endY=1 4 5 6 7 8 11 12 13 14 15 18 19 20 21 22

32、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,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(

33、1)=X0(1);for i=2:n X1(i)=X1(i-1)+X0(i); %计算累加生成序列endfor i=1:n-1B(i,1)=-0.5*(X1(i)+X1(i+1); %计算B，Yn B(i,2)=1; Y(i)=X0(i+1);endB;Y;alpha=(B*B)(-1)*B*Y; %做最小二乘估计a=alpha(1,1);b=alpha(2,1);d=b/a; %计算时间响应函数参数c=X1(1)-d;X2(1)=X0(1);X(1)=X0(1);for i=1:n-1 X2(i+1)=c*exp(-a*i)+d; X(i+1)=X2(i+1)-X2(i); %计算预测序列e

34、ndfor i=(n+1):(n+k) X2(i)=c*exp(-a*(i-1)+d; %计算预测序列 X(i)=X2(i)-X2(i-1);endfor i=1:n error(i)=X(i)-X0(i); error1(i)=abs(error(i); %计算残差 error2(i)=error1(i)/X0(i); %计算相对误差endc=std(error1)/std(X0) %计算后验差检验数end7.1.2 问题二相关程序: 初级产品的贸易特化系数（TSC）的程序：function TSCclear all;clc;%初级产品的处理CJCKZE=xlsread(E:历年出口货物分类

35、总额,Sheet1,C9:C30);CJJKZE=xlsread(E:历年进口货物分类总额,Sheet1,C9:C30);n=length(CJCKZE);%初级产品贸易特化系数的计算for i=1:n CJCK(i)=CJCKZE(i)-CJJKZE(i); CJZE(i)=CJCKZE(i)+CJJKZE(i); CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title(初级产品贸易特化系数);%食品及主要供食用的产品的处理SPCKZE=xlsread(E:历年出口货物分类总额,Sheet1,D9:D30);SPJKZE=xlsrea

36、d(E:历年进口货物分类总额,Sheet1,D9:D30);n=length(SPCKZE);%食品及主要供食用的贸易特化系数的计算for i=1:n SPCK(i)=SPCKZE(i)-SPJKZE(i); SPZE(i)=SPCKZE(i)+SPJKZE(i); SPTSC(i)=SPCK(i)/SPZE(i);endfigure;plot(1:n,SPTSC);title(食品及主要供食用贸易特化系数);%饮料及烟类的产品的处理YYCKZE=xlsread(E:历年出口货物分类总额,Sheet1,E9:E30);YYJKZE=xlsread(E:历年进口货物分类总额,Sheet1,E9:

37、E30);n=length(YYCKZE);%饮料及烟类的贸易特化系数的计算for i=1:n YYCK(i)=YYCKZE(i)-YYJKZE(i); YYZE(i)=YYCKZE(i)+YYJKZE(i); YYTSC(i)=YYCK(i)/YYZE(i);endfigure;plot(1:n,YYTSC);title(饮料及烟类贸易特化系数);%非食用原料的产品的处理FSPCKZE=xlsread(E:历年出口货物分类总额,Sheet1,F9:F30);FSPJKZE=xlsread(E:历年进口货物分类总额,Sheet1,F9:F30);n=length(FSPCKZE);%非食品原料

38、的贸易特化系数的计算for i=1:n FSPCK(i)=FSPCKZE(i)-FSPJKZE(i); FSPZE(i)=FSPCKZE(i)+FSPJKZE(i); FSPTSC(i)=FSPCK(i)/FSPZE(i);endfigure;plot(1:n,FSPTSC);title(非食品原料贸易特化系数);%矿物,润滑油相关原料的产品的处理KRCKZE=xlsread(E:历年出口货物分类总额,Sheet1,G9:G30);KRJKZE=xlsread(E:历年进口货物分类总额,Sheet1,G9:G30);n=length(KRCKZE);%矿物,润滑油相关原料的贸易特化系数的计算f

39、or i=1:n KRCK(i)=KRCKZE(i)-KRJKZE(i); KRZE(i)=KRCKZE(i)+KRJKZE(i); KRTSC(i)=KRCK(i)/KRZE(i);endfigure;plot(1:n,KRTSC);title(矿物,润滑油相关原料贸易特化系数);%动植物油脂及蜡相关原料的产品的处理DZYCKZE=xlsread(E:历年出口货物分类总额,Sheet1,H9:H30);DZYJKZE=xlsread(E:历年进口货物分类总额,Sheet1,H9:H30);n=length(DZYCKZE);%动植物油脂及蜡的贸易特化系数的计算for i=1:n DZYCK(

40、i)=DZYCKZE(i)-DZYJKZE(i); DZYZE(i)=DZYCKZE(i)+DZYJKZE(i); DZYTSC(i)=DZYCK(i)/DZYZE(i);endfigure;plot(1:n,DZYTSC);title(动植物油脂及蜡贸易特化系数);end 工业制成品的贸易特化系数（TSC）的程序：function TSC1clear all;clc;CJCKZE=xlsread(E:历年出口货物分类总额,Sheet1,I9:I30);CJJKZE=xlsread(E:历年进口货物分类总额,Sheet1,I9:I30);n=length(CJCKZE);for i=1:n C

41、JCK(i)=CJCKZE(i)-CJJKZE(i); CJZE(i)=CJCKZE(i)+CJJKZE(i); CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title(工业制成品贸易特化系数);CJCKZE=xlsread(E:历年出口货物分类总额,Sheet1,J9:J30);CJJKZE=xlsread(E:历年进口货物分类总额,Sheet1,J9:J30);n=length(CJCKZE);for i=1:n CJCK(i)=CJCKZE(i)-CJJKZE(i); CJZE(i)=CJCKZE(i)+CJJKZE(i); CJTSC(i)=CJCK(i)/CJZE(i);endfigure;plot(1:n,CJTSC);title(化学品贸易特化系数);