人民币汇率变动对我国进出口贸易的影响论文

1、巴湖北汽车工业学拜院大白学生数学建模邀拌请赛扒承 诺 书斑我们仔细阅读了佰第二届华中地区安大学生数学建模傲邀请赛稗的竞赛细则。懊我们完全明白，板在竞赛开始后参岸赛队员不能以任绊何方式（包括电爸话、电子邮件、跋网上咨询等）与稗队外的任何人（办包括指导教师）拌研究、讨论与赛肮题有关的问题。蔼我们知道，抄袭邦别人的成果是违白反竞赛规则的,跋 如果引用别人鞍的成果或其他公傲开的资料（包括吧网上查到的资料凹），必须按照规巴定的参考文献的拌表述方式在正文瓣引用处和参考文拌献中明确列出。吧我们郑重承诺，肮严格遵守竞赛规挨则，以保证竞赛暗的公正、公平性奥。如有违反竞赛懊规则的行为，我坝们将受到严肃处靶理。俺我

2、们的参赛报名案号为：懊 0500唉4 昂参赛队员 (伴签名) ：霸队员1：把 跋T933-2 邦杨进 疤队员2：昂 翱T843-哀3 搬李志刚般 挨队员3：傲 矮T923-1 背叶胜峰 搬第二届华中地区挨大学生数学建模耙邀请赛安编 号 专 用敖 页肮选择的题号：案 鞍 A 伴参赛的编号：奥 碍05004凹 八 埃 坝 爸 坝 捌 敖 敖 办 瓣（以下内容参赛跋队伍不需要填写搬）版 昂 竞安赛评阅埃编号：碍 颁 暗探讨人民币汇率耙变动对我国蔼进出口贸易奥的影响【摘 要】八近几十年来，世吧界经济朝全球化翱发展，资本流动扒自由化、国际金昂融、贸易的发展伴在不同方向上增颁加或降低汇率波阿动对国际贸易产

3、摆生了深远的影响颁。我国经济一直伴保持着良好的发俺展势头，特别是澳出口贸易的增长按让一些发达国家罢羡慕。人们对中埃国的经济前景比摆较乐观，人民币傲的国际地位越来按越高，经济实力班的提升客观上推颁动了人民币升值艾。叭为此，我们针对氨人民币汇率的不柏断提高对我国进胺出口商品贸易的案影响摆，各个产业结构芭的变化背以及由此对我国唉国内居民消费水碍平的霸影响鞍做了定性柏和定量耙的分析。耙首先，我们用灰癌色系统从长期的办角度分析了我国芭未来5年的人民鞍币汇率发展趋势案，及短期几天内肮人民币汇率发展瓣趋势，柏这样使得我们可矮以对人民币汇率拔短期的变化情况板又一个更加直观啊的了解。伴结果表明未来我暗国人民币升

4、值是翱不可阻挡的趋势罢。叭接着，八我们通过对各个版产业的进出口贸颁易情况做分析处凹理,得出历年来岸我国各个产业的背贸易特化系数.叭通过贸易特化系绊数的变化来更加坝直观的看出人民熬币汇率的改变对拔我国各个产业的唉进出口贸易的影八响.分析之后发鞍现人民币汇率的吧改变对我国大部笆分产业都有较大碍的影响.之后我耙们瓣以百各个产业的稗贸易特化系数罢为自变量,人民办币汇率为因变量巴建立每个产业与拜人民币汇率之间哎的多元回归分析癌,安定量分析了人民把币汇率波动对各败类商品影响程度按。癌最后通过整体上摆分析我国进坝口贸易额、癌出口贸易板额哎增长率与人民币埃汇率增长率奥、居民消费指数岸增长率爸的关系，稗以进出口

5、贸易额艾增长率为自变量案,人民币汇率增拜长率和居民消费奥指数增长率为因斑变量建立氨典型哎相关凹分析扒，盎结果疤表明了人民币汇胺率增长率对进口跋贸易额增长率起半到决定性的作用挨，对出口贸易额罢也有明显的影响耙。氨通过以上坝人民俺币汇率波动对我瓣国进出口贸易影氨响的分析及当前跋我国人民币汇率岸面临的形势哀，岸为我国经济快速矮平稳的发展人民唉币汇率应作如何柏调整提出一定的奥建议和意见稗。鞍【邦关键词俺】哀 爱人民币汇率 懊灰色系统 昂贸易特化系数 俺 实证分析 跋多元回归模型暗居民消费价格指鞍数 典型相关哀模型一、问题的提出懊2011年8月疤24日办人民币对美元汇般率中间价突破6爱.4，再创20把0

6、5年人民币汇靶率改革以来新高哀。澳关于人民币升值安的利弊人们褒贬耙不一。懊有的人般为人民币升值痛半心疾首，有的无柏动于衷。也有人叭认为人民币升值盎现在弊大于利，吧未来利大于弊，唉一时间，争论不邦断。因此，研究坝人民币升值对我哎国的影响是个热疤门问题。霸作为大学生我们扒应当学会关注我鞍国经济的发展，版树立强烈的主人版翁意识。针对人耙民币汇率升高的俺现象，我们从定班性和定量的角度版分析去考虑人民霸币汇率的波动对芭我国进出口贸易爸的影响。八我们提出如下问扳题：鞍对我国未来几年白的人民币对美元邦的汇率走势作出熬预测傲。敖分析一下人民币盎汇率的波动对我跋国进出口商品的笆影响霸。败分析一下人民币俺汇率、我

7、国居民暗消费CPI指数岸变化率与我国贸佰易进出口额变化佰率的关系罢。坝结合上面的分析稗，请说明一下我艾国政府应该如何板调整人民币汇率耙政策。扒 碍二、问题的分析般人民币升值对我般国国民经济及生昂活各个方面带来鞍多方面的影响。癌基于这样的前提阿下，我们通过所扳学的知识定量的败分析了人民币升巴值对我国进出口拌贸易的影响。绊对于第一问，我八们从国家统计局矮网站收集到我国拌人民币历年的年班均汇率，及从国昂家外汇管理局收鞍集了我国近一个俺月的人民币汇率伴中间价。对以上巴的数据分别用灰懊色系统模型分析败，对我国人民币板汇率在未来5年巴、未来10天的案波动状况。疤针对第二问，我安们从案国家统计局网站鞍收集到

8、我国近几背年以来的贸易进熬出口分类总额，白对初级产品、工绊业制品、其他三拌大类分别计算了拌其下的各个分支盎产业的贸易特化碍系数，然后做出俺各个产业的贸易拔特化系数波动曲白线，通过实证分拌析了人民币汇率坝波动对各个不同斑产业的影响。然叭后根据利用各求安解出的历年来产稗业的贸易特化系版数与人民币汇率盎建立线性回归模昂型，分析人民币按汇率的变化对各蔼个产业的影响关翱系。熬对于第三问，根凹据从吧国家统计局网站绊得到居民消费水叭平、货物进出口背贸易额、年均汇办率等相关数据，邦对居民消费价格耙指数的增长百分颁比、年均汇率的氨增长百分比与进奥口、出口贸易额背增长率建立了典爸型相关分析，来般研究其间的关系懊。

9、碍对于第四问，我绊们结合自己所了罢解的知识，及上俺面的分析，对我罢国汇率应该如何斑调整提出了一些绊建议。巴三、俺模型假设背与芭符号说明啊3.1模型胺假设：埃（1）哎 假设数据来源矮真实可靠。熬（2）凹 假设选取的指盎标有一定的代表癌性。半（3）般 假设进出口贸白易影响中只考虑氨供我国居民的消斑费指数、人民币搬汇率的影响。凹（4）澳 瓣假设袄人民币汇率和商凹品的贸易特化系氨数间存在线性关暗系。版（5）背 耙假设典型分析的袄两组数据都是连阿续变量，均服从啊多元正态分布。白3.2芭符号说明： 捌符号凹符号说明肮k芭代表产业编号版代表时间袄灰色模型后验差暗X袄i斑代表产业的贸易搬特化系数皑X安(n)鞍

10、(k)靶利率的n次累加扳序列（100美邦元折算）傲x笆代表居民消费指碍数增长百分比挨y矮1熬代表进口贸易增敖长百分比啊y巴2袄z霸代表出口贸易增肮长百分比芭代表年均人民币般汇率增长百分比背注：其他符号见氨文中说明拔 四霸、模型建立与求隘解4.1 问题一疤4.1.1绊基于灰色系统爱预测未来坝5翱年人民币对美元办汇率拌的数学模型颁1把：蔼分析建模 给袄定原始时间19霸9癌4啊瓣200癌9案年靶100美元兑换蔼成人民币的半序列X柏(k),对X百(k)生成1-袄AGO序列X邦(k)，另外可跋得Y隘见下表：爸k白1凹2按3按板14跋15昂16捌X翱861.87凹835.10癌831.42暗瓣760.40

11、霸694.51疤683.10捌X班861.87拜1696.97扮2528.39熬坝11528.4胺8埃12222.9唉9绊12906.0蔼9败Y澳捌1696.97哀2528.39佰凹11528.4奥8吧12222.9癌9般12906.0版9白其中，X版=AGO（X澳），即昂 暗 挨 阿 X斑(k)=靶 拜 碍 靶 Y霸= X搬(2)，X般(3)，澳熬，X胺(案15阿)，X翱(肮16傲)八对于上述的X扮(k)的GM(袄1,1)参数按皑照下式进行最小佰二乘估计得B=按将B,败Y把代入辨识算式,坝有:巴 翱 半=胺=(B案B)隘 B岸 Y巴=矮得GM（1,1半）模型为:跋灰微分方程： 佰 X耙(k)

12、-a案 袄=b疤白化方程： 巴- a扳 啊=b吧白化方程的时间半相应式蔼：案(t+1)=吧 佰般e暗+叭GM(1,1)矮模型的检验袄:佰残差序列的均方扳差：哎 霸 (百其中把)按原序列的均方差芭：叭 S哀 (其背中蔼)拜后验差检验芭：百C=版由结果可知：C搬=0.3889敖0.5,故可癌知预测结果合格芭。败利用MATLA搬B编程见附录傲7.2.1扒；运行结果见附哎录7.2.1结哀果。其中预测值百如下：背预测未来懊5盎年人民币对美元爱汇率扮表格如下:佰年份绊2010阿2011哀2012白2013矮2014稗预测值皑738.952半2拌731.405澳5瓣723.935芭9扮716.542吧6皑7

13、09.224半9八对结果的评价:斑由结果可以看出按100美元兑成拔人民币的数量呈碍递减趋势，即人扳民币与美元的汇罢率越来越小，人笆民币升值越来越哎严重。虽然该预埃测与现在的实际跋情况不符，但其班趋势是相同的。挨其可能原因为：捌以年为时间单位般误差太大版、般预测的时间跨度伴太大或原始数据巴过去的时间太久坝对现在的影响很爸小。为此，我们吧建立如下模型懊来预测短期情况懊。败4.1.搬2安基于灰色系统邦预测未来啊10霸天人民币对美元爱汇率霸的数学模型：吧分析建模 给敖定原始时间柏2011.7.斑1搬白2011.7.啊29内 靶100美元兑换版成人民币啊绊3澳八的搬序列X哀(k),对X半(k)生成1-胺

14、AGO序列X搬(k)，另外可捌得Y唉见下表：袄k哎1胺2安3蔼扮19般20板21背X案646.85靶646.61爸646.5拌搬644.26艾644.38爸644.42耙X哎646.85斑1293.46版1939.96案叭12280.6跋4皑12924.6胺2笆13569.0阿4般Y背疤1293.46拔1939.96般暗12280.6扒4芭12924.6百2耙13569.0跋4隘理论知识与傲4.1.1胺中的一样，在这吧里不再累赘。阿由结果可知：C罢=0.3板0拌890.伴3捌5,故可知预测阿结果邦优霸。懊利用MATLA癌B编程见附录败7.2.1稗；运行结果见附盎录7.2.1结邦果。其中预测值摆

15、如下：把日起爱7唉/30 柏7氨/31暗8爱/扳1拜8坝/鞍2芭8案/斑3班8唉/皑4疤8班/5蔼8背/6敖8癌/7版8背/8澳预测值摆644.40鞍644.2拜4皑644.0霸8爸643.91哀643.75懊643.5背9柏643.4摆3哀643.26扮643.10罢642.94把（备注：预测值柏为100美元可胺以兑换人民币数颁量（元）按对结果的评价:肮由预测结果可以颁看出岸人民币兑美元的氨汇率虽然存在大般小波动，但其整扒体趋势是不断减霸小的。预测值也半与现实相符，说伴明人民币升值进按一步加剧。4.2 问题二盎4.2.1癌首先根据10个八产业的贸易进出耙口额状况，作出盎各个产业的贸易唉特化系

16、数（TS绊C），然后作出耙其相应的曲线变坝化图，对其进行哎一定的说明稗。岸问题分析：鞍将我国基础产业扒分为：初级产品柏、工业制品、其班他。这里只对初跋级白产品、工业制品唉分析背（收集各产品的熬进出口额）埃。出口贸易商品摆结构有不同层次瓣，在这里我们采斑用了贸易特化系拔数的变化趋势来安阐述某种贸易商罢品或贸易产业的瓣结构优化。贸易靶特化系数反应的挨是一国某贸易埃产业或贸易产品敖的净出口额与该埃产品的进出口总版额之比。其计算颁公式为：唉 哎 拌 T吧SC=(爱X-M)/(X罢+M) 吧 -1熬矮TSC蔼艾1白当TSC指标值半介于-1拌到-0.3时，柏则代表该国产业岸的进口值远大于氨出口值，故该产阿

17、业属于霸“瓣进口强势爸”爱型产业，即该产班业在国际市场上蔼占有绝对的比较把优势；芭若TSC指标值瓣介于-0.3与版0.3之间，表哀示其进口与出口般值差异不大，代翱表该产业在国际艾上水平较高，属八于吧“般水平竞争疤”盎型产业；笆当TSC指标值暗介于0.3板败1时，表示其出绊口值远大于其进爸口值，故该产业阿在该国处于绝对胺的劣势，属于袄“巴出口强势艾”疤型产业。啊因此，从各种贸邦易行业贸易特化芭系数的变化我们哎可以大概地看出板我国贸易结构调伴整的方向。安将所有进出口产澳品分成两类：爸初级产品拌：熬食品及主要供食背用的昂活动物、饮料及癌烟类、非食用原爱料、矿物燃料、爸动鞍植物、油脂及蜡班工业制成品芭

18、：吧化学品及有关产翱品、轻纺产品、氨橡胶制品、矿冶按产品及其制品、伴机 挨械及运输设备、稗杂项制品绊MATLAB编安程见附录爸7.2.1邦，结果如下：昂1）胺初级产品捌及其组成项的图瓣像：胺从上面可以看出哀我国初级产品贸澳易特化曲线持续绊下降，说明了该板产业主要是通过暗进口百，属于强懊进口哀型，其对国外市颁场依赖过强，因瓣此人民币汇率升皑值是罢有利于该产业发跋展的。败从上面可以看出俺该产业的贸易特背化曲线基本处于拜正的，也就是说阿，主要是隘向背国外疤出口扒，属于强柏出口胺口型，人民币升凹值会扳不邦利于案此类商品在国外拜的发展霸。办从上面曲线可以八看出特化曲线后扒一阶段处于不断办下降的趋势，该瓣

19、产业正从水平竞佰争型向搬进百口型转化的趋势傲，人民币升值对啊其发展是搬有利柏。半从此图可以看出坝非食品原料贸易办特化系数最后将爸会向-1趋近，芭说明了非食品原埃料贸易几乎全部白依靠进口来提供摆，若人民币升值阿，会有利于我国八进口非食品原料吧。蔼上图反映埃了矿物润滑油相百关原料的特化系罢数逐渐逼近-1哎，说明了我国的罢矿物润滑油相关瓣原料主要依靠进八口，人民币升值爱将有利于版我国对此产品的埃进口与需求。 挨从上面可以看出疤我国动植物油脂半及蜡的贸易特化瓣系数逐步趋向-版1，说明了其进唉口趋势显著加强坝，人民币升值有暗利于其发展。摆由安结果办得出如下分析：八初级产品鞍、矮饮料及烟类、非胺食用原料、

20、矿物皑燃料、拜动植物拌油脂及蜡矮的贸易特化系数笆均随着人民币的皑升值而减小，搬食品及主要供食罢用熬物胺的贸易特化系数袄虽然有上下波动伴但整体比较稳定搬，搬由此可见扮人民币升值对我暗国的进出口贸易班有着密切影响。岸使八我国的笆初级产品绊类板产业的进口值远唉大于出口值，即斑“疤进口强势笆”佰型产业越来越多盎。而食品的波动澳比较平稳，可能邦受国家相关政策啊的影响，使我国邦的食品价格保持败平稳。有利于人败民的生活和社会颁的稳定板。盎2）捌工业制成品鞍及其组成项的图瓣像：啊从上面图像可以坝看出笆我国工业制成品斑贸易特化曲线基靶本处于-0.3摆到0.3之间，跋说明了此类产业般属于水平竞争型瓣，人民币汇率的

21、挨升高对其影响不邦大。但是我们也斑可以看到在汇率佰升高的同时其贸罢易特化系数也是翱呈上升趋势，如敖果人民币持续升皑值，该产业又转拜化成为出口强势鞍型的产业。把从上图可以看出版，化学品相关产爱业的贸易特化系挨数基本上处于-颁0.3一下，即佰便是又波动影响凹也不大，可以认袄为该产业是白进口强稗型案产业，在人民币版汇率升高有利于奥该产业的发展。斑从下图可以看出蔼：我国的轻纺橡袄胶、矿冶贸易特颁化曲线，在前期百处于进口强势到澳了后期逐渐转变巴为出口强势产业把，人民币的汇率皑提高不利于其再哎国际市场的竞争懊力。胺从此图可以看出氨我国机械运输设奥备贸易特化曲线敖一直处于上升的搬趋势，说明了该把产业将会向胺

22、出口型发展，人凹民币的升值将会挨不利于其发展。皑由扮结果埃得出如下分析：版工业制成品邦、化学品及有关熬产品、轻纺产品吧、橡胶制品、矿伴冶产品及其制品罢、机械及运输设扒备碍的贸易特化系数把均随着人民币的鞍升值而增大，杂把项制品的贸易特翱化系数虽然有上奥下波动但整体比跋较稳定，由此可绊见人民币升值对把我国的进出口贸叭易有着密切影响瓣。摆是癌我国的盎工业类摆产业的进口值远叭小拜于出口值，即埃“碍出口弱势瓣”唉型产业越来越多俺。对比背初级产品耙和笆工业制成品袄可知：我国的工昂业水品低下，且碍只能用缺少技术班含量的初级产品败和廉价的劳动力版来弥补进出口差瓣额，国家应加大挨科技的发展，尽绊快使工业水平赶傲

23、上来，加大背工业制成品扮的澳出口,减少因人哎民币升值而带来鞍的影响。熬4.2.2隘 多元回归分析建模过程：模型的建立:阿设自变量挨的观测值敖及因变量y对应矮的观测值y鞍满足关系式：巴 般 矮 y翱=背 i=肮1,2,3颁唉n；j=1,2暗,3,翱爸p懊其中，拌是相互独立且都啊服从正态分布的扮随机变量。唉由最小二乘法，熬得俺 暗 拔 肮2）般回归方程的显著吧性检验：利用F袄分布函数对变量笆前的系数检验。奥3）矮回归系数的显著败性检验： 利用坝t分布函数检验安。俺最后得到结果如扳下：癌产业名称扳编号扮方程显著性检验扳（F）哀参数显著鞍性检验（t）按回归方程蔼备注傲食品及主要供食败用的活动物芭1拌0

24、.0288伴0.0108靶0.0288y=773.15277x1+389.7118吧汇率对其影响显吧著柏饮料及烟酒隘2昂0.1440暗(不合格)傲-癌-八-唉-哀汇率与此行业没八有显著关系百非食用原料肮3阿0.0059熬0.001熬0.0059y=4193.92239x3+480.43293伴汇率对其影响显百著版矿物燃料版4胺0.0001胺0.0001捌0.0001y=309.74464x4+657.01952白汇率对其影响显伴著鞍动植物油脂把5阿0.00096稗0.0042鞍0.00096y=458.67428x5+382.11320拜汇率对其影响显隘著扒化学品扒6瓣0.1256阿(不合格)

25、拜-吧-疤-坝-拜汇率对其影响不柏显著胺轻纺橡胶敖7啊不合格瓣-扒-熬-皑-佰汇率对其影响不挨显著背机械及运输埃8百0.0001安0.0001安0.0001y=486.07219x8+791.92968岸汇率对其影响显俺著坝杂项哎9叭0.2222袄(不合格)盎-胺-安-肮-敖-皑汇率对其影响不扳显著蔼（备注：y代表伴年均人民币汇率敖，100美元兑坝换的人民币数量皑(元)凹，x搬i半代表贸易特化系敖数翱）按由上表看出，人罢民币汇率的变化坝对我国的大部分败产业都有显著的奥影响。也就是说板，人民币持续升岸值对我国各行各叭业都会产生决定班性改变，因此，皑为了稳定中国国巴民经济稳步发展矮，在人民币升值背

26、过程中国家需要佰进行干预，严格鞍控制人民币的升版值幅度，这样才矮能保证中国经济肮快速稳步发展，艾中国人民能更美按好的生活。4.3 问题三阿典型相关分析思柏想暗：绊考察一组随机变阿量和另一组随机把变量之间的相关办性，可以用他们坝之间的相关性。吧考虑到两组变量盎的线性组合：u1=a11x1+a21x2+an1xnv1=b11y1+b21y2+am1ym矮并研究他们之间啊的相关系数（u1，v1）澳，在所有的线性凹组合中找一对相隘关系数最大的线笆性组合，用这两熬个组合的单相关碍系数来表示两组袄变量的相关性，肮这两个线性组合拌叫做一对典型变半量。半下一步，再在两罢组比昂两的与u板，v不相关的线拔性组合中

27、，找一矮对相关系数最大岸的线性组合，设案为：u2=a12x1+a22x2+an2xnv2=b12y1+b22y2+am2ym矮且（u2，v2）八就是第二个典型颁相关系数，这样爸下去可以得到若白干对典型变量，背从而提取出两组盎变量间的全部信瓣息。俺就本题而言，用俺到典型相关分析把。利用SAS编邦程得到了：皑人民币汇率增长叭率、居民消费指鞍数增长率与进出把口贸易额增长率碍间的两个典型相绊关系数及特征值矮如下表： 凹序号颁典型相关系数半标准误差胺特征值昂方差比率隘1啊0.95212笆5拌0.01992蔼5鞍9.7001暗0.9357吧2佰0.63243坝1鞍0.12792澳7佰0.6666拜0.06

28、43安从上表可以看出案：第一对典型变拔量的相关系数较扮高，说明了这一扒对相关变量间的蔼关系是密切的。佰对傲相矮关系数检验按如下表：暗序号柏F计算值斑自由度矮F检验的显著性按概率拔1哎30.62案4捌0.001半2啊13.33柏1暗0.0016半这两对典型相关耙变量，进行F统皑计检验，均通过爸了检验，说明了案这两对典型变量稗间的相关性显著半，可以用来解释氨产出变量。八对冗余度分析得坝到下表氨：按典型变量的解释把能力癌序号埃产出组方差被影啊响组解释比例班对产出组的解释隘能力柏产出方差被典型柏变量解释比例白对影响组的解释八能力敖影响组方差被典蔼型变量的解释比安例懊1斑0.9065背0.9519拜0.

29、8629疤0.8760扮0.7941碍2巴0.4000芭0.481巴0.0193翱0.1240疤0.0496拜从此表可以看出挨：这两对典型相澳关变量的解释能八力均较强，特别拜是第一对相关变哀量对产出组和影瓣响组均有很高的霸解释能力；第一矮对典型变量对产碍出组方差有95版%以上的解释能安力，第二对典型癌相关变量较弱些巴，达到了48%疤左右。奥最后得出的芭典型变量模型柏：把序号搬典型相关模型百1u1=0.0546x+0.8714zw1=0.7245y1+0.8619y2岸2u2=0.4305x+0.8869zw2=0.8048y1+0.8816y2笆（版x代表居民消费澳指数增长百分比瓣、z代表人民

30、币澳汇率增长百分比爸；y1代表出口拜贸易额增长率、翱y2代表进口贸蔼易额的增长率；摆）结果分析：败人民币汇率增长俺率对进口贸易有板显著的作用。从哎上面第一对典型安变量分析结果可败以看出：人民币啊汇率增长百分比颁在典型变量u1岸占得比重达到了哀0.8714,岸出口贸易额增长绊率在w1中占的熬比重达到0.8啊619.由于u隘1,w1由前面啊分析知道其相关盎系数达到了0.按95以上，可以叭认为人民币汇率疤的增长率进口贸板易额增长率变化邦起到决定性的作罢用。胺从第二对典型相奥关变量分析知道胺：人民币汇率的霸增长率、居民消版费指数增长率对稗进出口贸易额增颁长率都一定的影叭响，不过期间的搬影响程度不是那氨

31、么显著。败综上所述：人民柏币汇率的增长率鞍对我国进口贸易挨额有着决定性的安因素，同时它和氨居民消费指数增办长率对我国进出邦口额增长率都有把密切的关联。4.4 问题四百通过上面的分析斑在人民币不断升霸值的条件下，我昂国政府应该在一癌定保证人民币汇罢率在小范围的增爸长。对于那些出半口型的产业特别百是工业制品，我笆们应该提高产品矮的技术优势，而疤不是通过廉价的翱劳动力去争取国啊际市场的优势。佰对于那些进口型霸的产业要从提高俺国内资源的利用懊率和经济效率的罢角度出发，以减扳少对国际市场的俺依赖程度，特别爱是那些矿物润油澳类。最后要不断搬加强我国产业结叭构的调整经济制案度的调整，利用颁人民币汇率的小拜范

32、围的增长来优哎化我国产业，同芭时来促进我国企笆业竞争力。癌五凹、模型评价与推阿广办第一问用了芭灰色预测，对于隘人民币汇率的影癌响的因素太多，氨而之间的确切关哎系又无法确定。胺因此用灰色系统巴可以对那些只知坝道一部分，而对挨整个其它不甚了艾解的状态下，是胺一种比较好的预百测方法。佰第二问用了典型背变量分析的方法耙，对于多组变量碍间的关系是错综班复杂的，通过其埃间的组合关系可捌以更清楚解释各巴个组变量间的关百联程度疤。捌第三问应用了多氨元回归，可以用芭来研究多个变量罢对一个变量的影爱响程度，对于生爱活中的一些存在鞍线性关系的变量鞍间的关系讨论，皑是一种典型的方癌法。六、参考文献扳1敖 柏汪晓银 周

33、保平叭，数学建模与数办学实验，北京：百科学出版社，2疤010。靶2疤 柏中华人民共和国氨国家统计局网站败http:/把www.sta昂.唉cn/tjsj暗/ndsj/挨3 国家外八汇管理局网站h百ttp:/w扒ww.safe鞍.氨/model_白safe/in肮dex.htm哀l 七、附录澳7.1 附录程办序:扮7.1.1蔼 问题一相关程敖序:傲基于背灰色系统预测邦未来21年汇率霸的八程序：稗functio哎n X=R拜mbExRat般a()邦clear a罢ll;clc;扒X0=xlsr版ead(E:凹人民币汇率(矮年平均价),爸Sheet1巴,B16:扮B31);拔n=lengt霸h(X0

34、);凹for i=1哀:n胺 X0(啊i)=X0(i拜);endk=5;傲X=GM1背1(X0,k)唉;figure;芭plot(19伴94:2009败,X0,g*芭-);办title(稗人民币汇率年变爱化及五年预测值艾);安hold on稗 般plot(19百94:2014案,X);end%伴functio案n X,c,俺error1,阿error2八=GM11(X埃0,k)疤% 建立函数百X,c,err办or1,err唉or2=ex白ample9_袄3_2_3(X挨0,k)版% 其中X0为隘输入序列，k为班预测长度，跋% X为预测输凹出序列，c为后般验差检验数，e盎rror1为残柏差，er

35、ror绊2为相对误差format;败n=lengt埃h(X0);X1=;败X1(1)=X吧0(1);扮for i=2翱:n哎 X1(摆i)=X1(i瓣-1)+X0(鞍i); 罢 拔 %计算累加隘生成序列end败for i=1搬:n-1稗B(i,1)=阿-0.5*(X办1(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

36、)搬-d;矮X2(1)=X熬0(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); 啊 艾 %计算预测序拜列end按for i=(哎n+1):(n俺+k)靶 X2(奥i)=c*ex哎p(-a*(i捌-1)+d;肮 耙%计算预测序列岸 X(i罢)=X2(i)班-X2(i-1班);end办for i=1鞍:n半 err安or(i)=X耙(i)-X0(稗i);耙 err艾or1(i)=背abs(err败or(i);斑 把%计算残差颁 err八or2(i)=傲error1(半i)/X0(i

37、扒); 碍%计算相对误差end巴c=std(e摆rror1)/笆std(X0)暗 敖 %计算后癌验差检验数end肮基于矮灰色系统预测霸未来31天汇率哎的蔼程序：碍functio八n X=R哀mbExRat靶a1()跋clear a唉ll;clc;敖X0=xlsr艾ead(E:碍最近一个月人哀民币汇率中间价熬,Shee敖t1,B2搬4:B4);唉n=lengt捌h(X0);背for i=1癌:n/2蔼 tem背p=X0(i)碍;半 X0(唉i)=X0(n阿-i+1);阿 X0(肮n-i+1)=白temp;end碍Y=1 4 扒5 6 7 8捌 11 12 案13 14 1爱5 18 19吧 20

38、 21 安22 25 2俺6 27 28罢 29;k=10;捌X=GM1把1(X0,k)拌;figure;矮plot(Y,扮X0,g*-皑);半title(稗人民币汇率月变碍化及十天预测值挨);hold on蔼Y=Y 30芭 31 32 霸33 34 3挨5 36 37扮 38 39氨;霸plot(Y,扮X);end%般functio般n X,c,矮error1,笆error2胺=GM11(X阿0,k)颁% 建立函数拌X,c,err板or1,err艾or2=ex拔ample9_伴3_2_3(X挨0,k)昂% 其中X0为版输入序列，k为爸预测长度，邦% X为预测输皑出序列，c为后唉验差检验数，e

39、笆rror1为残扒差，error爱2为相对误差format;伴n=lengt绊h(X0);X1=;佰X1(1)=X爱0(1);蔼for i=2癌:n澳 X1(霸i)=X1(i胺-1)+X0(翱i); 拜 拔 %计算累加耙生成序列end耙for i=1鞍:n-1疤B(i,1)=邦-0.5*(X挨1(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; 捌 翱 啊 案 伴 %计佰算时间响

40、应函数皑参数俺c=X1(1)氨-d;白X2(1)=X白0(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); 拔 霸 %计算预测序案列end皑for i=(澳n+1):(n斑+k)版 X2(爸i)=c*ex班p(-a*(i啊-1)+d;版 傲%计算预测序列吧 X(i疤)=X2(i)白-X2(i-1蔼);end按for i=1搬:n埃 err昂or(i)=X案(i)-X0(板i);把 err懊or1(i)=班abs(err白or(i);暗 皑%计算残差搬 err爱or2(i)=懊er

41、ror1(罢i)/X0(i背); 俺%计算相对误差end蔼c=std(e爱rror1)/白std(X0)跋 伴 %计算后袄验差检验数end白7.1.2搬 问题二相关程搬序:版懊 拌初级产品埃的贸易特化系数挨（TSC）的程靶序：班functio巴n TSC斑clear a哎ll;clc;稗%初级产品的岸处理坝CJCKZE=斑xlsread霸(E:历年霸出口货物分类总癌额,She八et1,C背9:C30)板;伴CJJKZE=坝xlsread胺(E:历年肮进口货物分类总柏额,She俺et1,C矮9:C30)啊;败n=lengt哀h(CJCKZ蔼E);盎%初级产品贸吧易特化系数的计捌算氨for i=1

42、扒:n爸 CJC佰K(i)=CJ啊CKZE(i)扮-CJJKZE罢(i);埃 CJZ扮E(i)=CJ霸CKZE(i)把+CJJKZE隘(i);吧 CJT拔SC(i)=C哎JCK(i)/伴CJZE(i)八;endfigure;爱plot(1:埃n,CJTSC捌);班title(捌初级产品贸易特爸化系数);艾%食品及主要靶供食用的产品的昂处理笆SPCKZE=岸xlsread斑(E:历年挨出口货物分类总邦额,She白et1,D坝9:D30)俺;搬SPJKZE=办xlsread哎(E:历年坝进口货物分类总拜额,She澳et1,D昂9:D30)吧;班n=lengt跋h(SPCKZ吧E);百%食品及主要败

43、供食用的贸易特岸化系数的计算疤for i=1板:n罢 SPC盎K(i)=SP澳CKZE(i)岸-SPJKZE巴(i);板 SPZ案E(i)=SP盎CKZE(i)袄+SPJKZE阿(i);把 SPT坝SC(i)=S叭PCK(i)/挨SPZE(i)叭;endfigure;啊plot(1:拜n,SPTSC癌);败title(矮食品及主要供食摆用贸易特化系数傲);绊%饮料及烟类扳的产品的处理吧YYCKZE=翱xlsread昂(E:历年吧出口货物分类总安额,She办et1,E颁9:E30)安;扒YYJKZE=板xlsread瓣(E:历年凹进口货物分类总拜额,She傲et1,E艾9:E30)拜;捌n=le

44、ngt癌h(YYCKZ伴E);拌%饮料及烟类巴的贸易特化系数稗的计算版for i=1昂:n肮 YYC败K(i)=YY吧CKZE(i)柏-YYJKZE绊(i);捌 YYZ懊E(i)=YY搬CKZE(i)傲+YYJKZE矮(i);扳 YYT隘SC(i)=Y哎YCK(i)/班YYZE(i)蔼;endfigure;邦plot(1:捌n,YYTSC绊);安title(敖饮料及烟类贸易岸特化系数);拜%非食用原料颁的产品的处理傲FSPCKZE挨=xlsrea岸d(E:历埃年出口货物分类搬总额,Sh伴eet1,艾F9:F30扳);氨FSPJKZE坝=xlsrea班d(E:历班年进口货物分类按总额,Sh哀ee

45、t1,八F9:F30肮);靶n=lengt摆h(FSPCK坝ZE);班%非食品原料板的贸易特化系数昂的计算懊for i=1哎:n拔 FSP肮CK(i)=F叭SPCKZE(岸i)-FSPJ佰KZE(i);半 FSP蔼ZE(i)=F霸SPCKZE(办i)+FSPJ坝KZE(i);熬 FSP昂TSC(i)=袄FSPCK(i靶)/FSPZE白(i);endfigure;唉plot(1:翱n,FSPTS八C);邦title(暗非食品原料贸易班特化系数);爱%矿物,润滑阿油相关原料的产瓣品的处理阿KRCKZE=扮xlsread笆(E:历年耙出口货物分类总坝额,She坝et1,G稗9:G30)扳;爸KRJK

46、ZE=扳xlsread佰(E:历年稗进口货物分类总败额,She唉et1,G捌9:G30)哀;矮n=lengt懊h(KRCKZ跋E);班%矿物,润滑埃油相关原料的贸蔼易特化系数的计板算凹for i=1昂:n斑 KRC背K(i)=KR斑CKZE(i)瓣-KRJKZE哎(i);捌 KRZ盎E(i)=KR暗CKZE(i)搬+KRJKZE芭(i);靶 KRT阿SC(i)=K昂RCK(i)/背KRZE(i)班;endfigure;吧plot(1:罢n,KRTSC挨);拜title(耙矿物,润滑油相按关原料贸易特化伴系数);矮%动植物油脂按及蜡相关原料的霸产品的处理昂DZYCKZE绊=xlsrea斑d(E:

47、历挨年出口货物分类碍总额,Sh斑eet1,八H9:H30爸);暗DZYJKZE阿=xlsrea胺d(E:历阿年进口货物分类癌总额,Sh跋eet1,巴H9:H30胺);暗n=lengt懊h(DZYCK叭ZE);岸%动植物油脂爸及蜡的贸易特化阿系数的计算八for i=1奥:n唉 DZY扮CK(i)=D吧ZYCKZE(澳i)-DZYJ搬KZE(i);板 DZY疤ZE(i)=D胺ZYCKZE(案i)+DZYJ摆KZE(i);胺 DZY皑TSC(i)=阿DZYCK(i伴)/DZYZE蔼(i);endfigure;爸plot(1:般n,DZYTS爸C);半title(啊动植物油脂及蜡澳贸易特化系数挨);e

48、nd矮背 扮工业制成品的贸俺易特化系数（T坝SC）的程序：翱functio盎n TSC1澳clear a安ll;clc;按CJCKZE=唉xlsread岸(E:历年靶出口货物分类总凹额,She暗et1,I蔼9:I30)傲;俺CJJKZE=皑xlsread爱(E:历年碍进口货物分类总巴额,She肮et1,I岸9:I30)疤;跋n=lengt艾h(CJCKZ袄E);背for i=1爱:n班 CJC版K(i)=CJ八CKZE(i)绊-CJJKZE安(i);绊 CJZ跋E(i)=CJ唉CKZE(i)背+CJJKZE胺(i);班 CJT柏SC(i)=C傲JCK(i)/昂CJZE(i)拌;endfigur

49、e;白plot(1:败n,CJTSC拜);奥title(稗工业制成品贸易傲特化系数);拜CJCKZE=斑xlsread哎(E:历年扒出口货物分类总翱额,She白et1,J隘9:J30)扒;摆CJJKZE=艾xlsread隘(E:历年澳进口货物分类总挨额,She班et1,J隘9:J30)佰;邦n=lengt斑h(CJCKZ芭E);隘for i=1唉:n佰 CJC俺K(i)=CJ傲CKZE(i)肮-CJJKZE办(i);芭 CJZ坝E(i)=CJ捌CKZE(i)案+CJJKZE芭(i);斑 CJT佰SC(i)=C拔JCK(i)/袄CJZE(i)爱;endfigure;昂plot(1:叭n,CJTS

50、C般);肮title(氨化学品贸易特化俺系数);挨CJCKZE=霸xlsread笆(E:历年绊出口货物分类总胺额,She皑et1,K稗9:K30)颁;胺CJJKZE=板xlsread笆(E:历年按进口货物分类总柏额,She奥et1,K扮9:K30)胺;跋n=lengt败h(CJCKZ艾E);胺for i=1叭:n版 CJC坝K(i)=CJ耙CKZE(i)摆-CJJKZE皑(i);捌 CJZ绊E(i)=CJ版CKZE(i)敖+CJJKZE懊(i);笆 CJT癌SC(i)=C版JCK(i)/板CJZE(i)败;endfigure;稗plot(1:办n,CJTSC肮);凹title(懊轻纺.橡胶.矿

51、霸冶贸易特化系数半);凹CJCKZE=巴xlsread伴(E:历年安出口货物分类总昂额,She柏et1,L败9:L30)巴;坝CJJKZE=捌xlsread挨(E:历年败进口货物分类总阿额,She柏et1,L坝9:L30)巴;胺n=lengt哀h(CJCKZ哀E);昂for i=1碍:n傲 CJC熬K(i)=CJ扳CKZE(i)罢-CJJKZE奥(i);霸 CJZ邦E(i)=CJ鞍CKZE(i)绊+CJJKZE凹(i);板 CJT爱SC(i)=C败JCK(i)/颁CJZE(i)颁;endfigure;安plot(1:疤n,CJTSC半);拌title(耙机械与运输设备肮贸易特化系数懊);白CJ

52、CKZE=肮xlsread般(E:历年捌出口货物分类总凹额,She颁et1,M碍9:M30)蔼;办CJJKZE=板xlsread哀(E:历年矮进口货物分类总埃额,She岸et1,M敖9:M30)颁;靶n=lengt拔h(CJCKZ班E);跋for i=1艾:n版 CJC啊K(i)=CJ爱CKZE(i)颁-CJJKZE邦(i);颁 CJZ板E(i)=CJ笆CKZE(i)柏+CJJKZE靶(i);爸 CJT疤SC(i)=C吧JCK(i)/白CJZE(i)颁;endfigure;板plot(1:唉n,CJTSC吧);胺title(扮杂项制品贸易特蔼化系数);end摆建立多元回归关拔系的程序扒(SAS

53、)班：稗data ra埃ta;傲input y跋 x1-x9唉;cards;暗157.7捌百0.0098 佰澳0.3684 鞍八-0.3500把 盎傲0.9094 吧邦-0.5987爱 把按-0.4440靶 凹白-0.0190案 矮肮-0.7172捌 百板0.6791 邦293.66搬蔼0.4201 岸版-0.3248暗 摆搬-0.0990熬 奥鞍0.9529 安艾0.0506 邦澳-0.5339瓣 耙搬-0.4518般 胺跋-0.9092敖 隘肮0.2940 蔼478.32岸笆0.3292 笆俺0.3707 背奥-0.0746败 哎癌0.6092 半办-0.7183啊 安俺-0.2812碍

54、爸艾0.1708 版按-0.5018傲 袄碍0.7156 芭532.33鞍背0.4416 伴扳0.4513 稗斑-0.1787败 拔奥0.3846 扮奥-0.6548翱 敖隘-0.4169邦 靶澳0.1588 啊耙-0.4655扳 蔼奥0.7441 爱551.46爸疤0.4507 靶搬0.5016 俺坝-0.2951百 罢吧0.1359 挨矮-0.5813背 岸扮-0.4391芭 案扮-0.0886胺 澳吧-0.4063罢 百懊0.7194 岸576.2柏懊0.5840 奥把0.5724 笆伴-0.2810哎 案柏-0.1722艾 半凹-0.4201拌 袄霸-0.3546安 蔼敖-0.2702

55、柏 败懊-0.4932按 颁盎0.7131 叭861.87稗肮0.5230 靶霸0.8729 捌稗-0.2862碍 笆拜0.0042 疤澳-0.5703埃 靶爱-0.3209扒 埃版-0.0949板 熬拜-0.4031岸 蔼蔼0.7613 颁835.1罢八0.2376 拜颁0.5533 板笆-0.3980稗 颁敖0.0196 哀拔-0.7032耙 爱八-0.3109傲 拌巴0.0568 艾耙-0.2527八 笆八0.7369 班831.42澳隘0.2867 伴胺0.4595 叭隘-0.4513暗 瓣扮-0.0739按 盎颁-0.6372俺 稗凹-0.3420班 罢哎-0.0483办 霸澳-0.

56、2159啊 碍吧0.7385 岸828.98懊啊0.4403 矮案0.5325 俺扳-0.4821叭 扒芭-0.1919暗 癌艾-0.4449捌 阿罢-0.3072板 办啊0.0332 芭佰-0.0940蔼 柏伴0.7836 啊827.91搬佰0.4702 稗昂0.6898 隘敖-0.5056扳 盎凹-0.1340瓣 瓣暗-0.6585佰 霸疤-0.3227邦 靶拔0.0221 爱敖-0.0619叭 瓣邦0.7850 扒827.83巴巴0.4858 氨吧0.5751 芭跋-0.5293扒 奥皑-0.3134白 跋吧-0.8239八 败颁-0.3970岸 拌罢-0.0156艾 半皑-0.0828稗