1、湖北汽車工業(yè)學(xué)院大學(xué)生數(shù)學(xué)建模邀請(qǐng)賽 承諾書我們仔細(xì)閱讀了第二屆華中地區(qū)大學(xué)生數(shù)學(xué)建模邀請(qǐng)賽的競(jìng)賽細(xì)則。我們完全明白，在競(jìng)賽開(kāi)始后參賽隊(duì)員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊(duì)外的任何人（包括指導(dǎo)教師）研究、討論與賽題有關(guān)的問(wèn)題。我們知道，抄襲別人的成果是違反競(jìng)賽規(guī)則的，如果引用別人的成果或其他公開(kāi)的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻(xiàn)的表述方式在正文引用處和參考文獻(xiàn)中明確列出。我們鄭重承諾，嚴(yán)格遵守競(jìng)賽規(guī)則，以保證競(jìng)賽的公正、公平性。如有違反競(jìng)賽規(guī)則的行為，我們將受到嚴(yán)肅處理。我們的參賽報(bào)名號(hào)為：05004參賽隊(duì)員（簽名）：隊(duì)員1:T933-2楊進(jìn)隊(duì)員2:T843

2、-3李志剛隊(duì)員3:T923-1葉勝峰0第二屆華中地區(qū)大學(xué)生數(shù)學(xué)建模邀請(qǐng)賽編號(hào)專用頁(yè)選擇的題號(hào)：A參賽的編號(hào)：05004（以下內(nèi)容參賽隊(duì)伍不需要填寫）競(jìng)賽評(píng)閱編號(hào):ii探討人民幣匯率變動(dòng)對(duì)我國(guó)進(jìn)由口貿(mào)易的影響【摘要】近幾十年來(lái)，世界經(jīng)濟(jì)朝全球化發(fā)展，資本流動(dòng)自由化、國(guó)際金融、貿(mào)易的發(fā)展在 不同方向上增加或降低匯率波動(dòng)對(duì)國(guó)際貿(mào)易產(chǎn)生了深遠(yuǎn)的影響。我國(guó)經(jīng)濟(jì)一直保持著良好的發(fā)展勢(shì)頭，特別是出口貿(mào)易的增長(zhǎng)讓一些發(fā)達(dá)國(guó)家羨慕。人們對(duì)中國(guó)的經(jīng)濟(jì)前景比 較樂(lè)觀，人民幣的國(guó)際地位越來(lái)越高，經(jīng)濟(jì)實(shí)力的提升客觀上推動(dòng)了人民幣升值。為此，我們針對(duì)人民幣匯率的不斷提高對(duì)我國(guó)進(jìn)出口商品貿(mào)易的影響，各個(gè)產(chǎn)業(yè)結(jié) 構(gòu)的變化以及

3、由此對(duì)我國(guó)國(guó)內(nèi)居民消費(fèi)水平的影響做了定性和定量的分析。首先，我們用灰色系統(tǒng)從長(zhǎng)期的角度分析了我國(guó)未來(lái) 5年的人民幣匯率發(fā)展趨勢(shì)， 及短期幾天內(nèi)人民幣匯率發(fā)展趨勢(shì)， 這樣使得我們可以對(duì)人民幣匯率短期的變化情況又 一個(gè)更加直觀的了解。結(jié)果表明未來(lái)我國(guó)人民幣升值是不可阻擋的趨勢(shì)。接著，我們通過(guò)對(duì)各個(gè)產(chǎn)業(yè)的進(jìn)出口貿(mào)易情況做分析處理，得出歷年來(lái)我國(guó)各個(gè)產(chǎn)業(yè)的貿(mào)易特化系數(shù).通過(guò)貿(mào)易特化系數(shù)的變化來(lái)更加直觀的看出人民幣匯率的改變對(duì)我 國(guó)各個(gè)產(chǎn)業(yè)的進(jìn)出口貿(mào)易的影響.分析之后發(fā)現(xiàn)人民幣匯率的改變對(duì)我國(guó)大部分產(chǎn)業(yè)都 有較大的影響.之后我們以各個(gè)產(chǎn)業(yè)的貿(mào)易特化系數(shù)為自變量，人民幣匯率為因變量建立每個(gè)產(chǎn)業(yè)與人民幣匯率

4、之間的多元回歸分析 ，定量分析了人民幣匯率波動(dòng)對(duì)各類商品 影響程度。最后通過(guò)整體上分析我國(guó)進(jìn)口貿(mào)易額、出口貿(mào)易額增長(zhǎng)率與人民幣匯率增長(zhǎng)率、居 民消費(fèi)指數(shù)增長(zhǎng)率的關(guān)系，以進(jìn)出口貿(mào)易額增長(zhǎng)率為自變量，人民幣匯率增長(zhǎng)率和居民消費(fèi)指數(shù)增長(zhǎng)率為因變量建立典型相關(guān)分析， 結(jié)果表明了人民幣匯率增長(zhǎng)率對(duì)進(jìn)口貿(mào)易 額增長(zhǎng)率起到?jīng)Q定性的作用，對(duì)出口貿(mào)易額也有明顯的影響。通過(guò)以上人民幣匯率波動(dòng)對(duì)我國(guó)進(jìn)出口貿(mào)易影響的分析及當(dāng)前我國(guó)人民幣匯率面 臨的形勢(shì)，為我國(guó)經(jīng)濟(jì)快速平穩(wěn)的發(fā)展人民幣匯率應(yīng)作如何調(diào)整提出一定的建議和意 見(jiàn)。【關(guān)鍵詞】人民幣匯率灰色系統(tǒng)貿(mào)易特化系數(shù)實(shí)證分析多元回歸模型 居民消費(fèi)價(jià)格指數(shù)典型相關(guān)模型一、問(wèn)

5、題的提出2011年8月24日人民幣對(duì)美元匯率中間價(jià)突破 6.4 ,再創(chuàng)2005年人民幣匯率改革 以來(lái)新高。關(guān)于人民幣升值的利弊人們褒貶不一。有的人為人民幣升值痛心疾首，有的無(wú)動(dòng)于 衷。也有人認(rèn)為人民幣升值現(xiàn)在弊大于利，未來(lái)利大于弊，一時(shí)間，爭(zhēng)論不斷。因此， 研究人民幣升值對(duì)我國(guó)的影響是個(gè)熱門問(wèn)題。作為大學(xué)生我們應(yīng)當(dāng)學(xué)會(huì)關(guān)注我國(guó)經(jīng)濟(jì)的發(fā)展，樹立強(qiáng)烈的主人翁意識(shí)。針對(duì)人民 幣匯率升高的現(xiàn)象，我們從定性和定量的角度分析去考慮人民幣匯率的波動(dòng)對(duì)我國(guó)進(jìn)出 口貿(mào)易的影響。我們提出如下問(wèn)題：1、對(duì)我國(guó)未來(lái)幾年的人民幣對(duì)美元的匯率走勢(shì)作出預(yù)測(cè)。2、分析一下人民幣匯率的波動(dòng)對(duì)我國(guó)進(jìn)出口商品的影響。3、分析一下人

6、民幣匯率、我國(guó)居民消費(fèi) CPI指數(shù)變化率與我國(guó)貿(mào)易進(jìn)出口額變化率的 關(guān)系。4、結(jié)合上面的分析，請(qǐng)說(shuō)明一下我國(guó)政府應(yīng)該如何調(diào)整人民幣匯率政策。二、問(wèn)題的分析人民幣升值對(duì)我國(guó)國(guó)民經(jīng)濟(jì)及生活各個(gè)方面帶來(lái)多方面的影響。基于這樣的前提 下，我們通過(guò)所學(xué)的知識(shí)定量的分析了人民幣升值對(duì)我國(guó)進(jìn)出口貿(mào)易的影響。對(duì)于第一問(wèn)，我們從國(guó)家統(tǒng)計(jì)局網(wǎng)站收集到我國(guó)人民幣歷年的年均匯率，及從國(guó)家 外匯管理局收集了我國(guó)近一個(gè)月的人民幣匯率中間價(jià)。對(duì)以上的數(shù)據(jù)分別用灰色系統(tǒng)模型分析，對(duì)我國(guó)人民幣匯率在未來(lái) 5年、未來(lái)10天的波動(dòng)狀況。針對(duì)第二問(wèn)，我們從國(guó)家統(tǒng)計(jì)局網(wǎng)站收集到我國(guó)近幾年以來(lái)的貿(mào)易進(jìn)出口分類總 額，對(duì)初級(jí)產(chǎn)品、工業(yè)制品

7、、其他三大類分別計(jì)算了其下的各個(gè)分支產(chǎn)業(yè)的貿(mào)易特化系 數(shù)，然后做出各個(gè)產(chǎn)業(yè)的貿(mào)易特化系數(shù)波動(dòng)曲線，通過(guò)實(shí)證分析了人民幣匯率波動(dòng)對(duì)各 個(gè)不同產(chǎn)業(yè)的影響。然后根據(jù)利用各求解出的歷年來(lái)產(chǎn)業(yè)的貿(mào)易特化系數(shù)與人民幣匯率 建立線性回歸模型，分析人民幣匯率的變化對(duì)各個(gè)產(chǎn)業(yè)的影響關(guān)系。對(duì)于第三問(wèn)，根據(jù)從國(guó)家統(tǒng)計(jì)局網(wǎng)站得到居民消費(fèi)水平、貨物進(jìn)出口貿(mào)易額、年均 匯率等相關(guān)數(shù)據(jù)，對(duì)居民消費(fèi)價(jià)格指數(shù)的增長(zhǎng)百分比、年均匯率的增長(zhǎng)百分比與進(jìn)口、 出口貿(mào)易額增長(zhǎng)率建立了典型相關(guān)分析，來(lái)研究其間的關(guān)系。對(duì)于第四問(wèn)，我們結(jié)合自己所了解的知識(shí)，及上面的分析，對(duì)我國(guó)匯率應(yīng)該如何調(diào) 整提出了一些建議。三、模型假設(shè)與符號(hào)說(shuō)明3.1模型

8、假設(shè)：(1)假設(shè)數(shù)據(jù)來(lái)源真實(shí)可靠。(2)假設(shè)選取的指標(biāo)有一定的代表性。(3)假設(shè)進(jìn)出口貿(mào)易影響中只考慮供我國(guó)居民的消費(fèi)指數(shù)、人民幣匯率的影響(4)假設(shè)人民幣匯率和商品的貿(mào)易特化系數(shù)間存在線性關(guān)系。(5)假設(shè)典型分析的兩組數(shù)據(jù)都是連續(xù)變量，均服從多元正態(tài)分布。3.2符號(hào)說(shuō)明:符號(hào)符號(hào)說(shuō)明i代表產(chǎn)業(yè)編號(hào)k代表時(shí)間C灰色模型后驗(yàn)差X代表產(chǎn)業(yè)的貿(mào)易特化系數(shù)X(n) (k)利率的n次累加序列(100美元折算)x代表居民消費(fèi)指數(shù)增長(zhǎng)百分比yi代表進(jìn)口貿(mào)易增長(zhǎng)百分比y2代表出口貿(mào)易增長(zhǎng)百分比z代表年均人民幣匯率增長(zhǎng)百分比注：其他符號(hào)見(jiàn)文中說(shuō)明四、模型建立與求解4.1 問(wèn)題一4.1.1 基于灰色系統(tǒng)預(yù)測(cè)未來(lái)5

9、年人民幣對(duì)美元匯率的數(shù)學(xué)模型1:分析建模給定原始時(shí)間1994 2009年100美元兌換成人民幣的序列X(0)(k),對(duì)X(0)(k)生成1-AGO序列X(k),另外可得Yn見(jiàn)下表：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對(duì)于

10、上述的X(0) (k)的GM(1,1)參數(shù)按照下式進(jìn)行最小二乘估計(jì)得一-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代入辨識(shí)算式，有:得GM(1,1 )模

11、型為:1)灰微分方程：X(k)-a z(k)=b2)(1)白化方程：d- adtz (k)=b3)白化方程的時(shí)間相應(yīng)式：x(1) (t+1)= x(0)一be力+2 a a4) GM(1,1)模型的檢驗(yàn)：殘差序列的均方差:( 其中名=:一)0 20原序列的均方差：S2 = X0(k)-X 2'、X(0)(k)(其中 X =)n -1n后驗(yàn)差檢驗(yàn)：C=SS2由結(jié)果可知：C=0.3889<0.5,故可知預(yù)測(cè)結(jié)果合格。利用MATLA酎程見(jiàn)附錄7.2.1 ;運(yùn)行結(jié)果見(jiàn)附錄7.2.1結(jié)果。其中預(yù)測(cè)值如下:880人民幣匯率年變化及五年預(yù)測(cè)值86084082080078076074072070

12、0680 19941996199820002002200420062008201020122014預(yù)測(cè)未來(lái)5年人民幣對(duì)美元匯率表格如下年份20102011201220132014預(yù)測(cè)值P 738.9522731.4055723.9359716.5426709.2249對(duì)結(jié)果的評(píng)價(jià)：由結(jié)果可以看出100美元兌成人民幣的數(shù)量呈遞減趨勢(shì)，即人民幣 與美元的匯率越來(lái)越小，人民幣升值越來(lái)越嚴(yán)重。雖然該預(yù)測(cè)與現(xiàn)在的實(shí)際情況不符， 但其趨勢(shì)是相同的。其可能原因?yàn)椋阂阅隇闀r(shí)間單位誤差太大、預(yù)測(cè)的時(shí)間跨度太大或 原始數(shù)據(jù)過(guò)去的時(shí)間太久對(duì)現(xiàn)在的影響很小。為此，我們建立如下模型來(lái)預(yù)測(cè)短期情況。4.1.2 基于灰色系

13、統(tǒng)預(yù)測(cè)未來(lái)10天人民幣對(duì)美元匯率的數(shù)學(xué)模型：分析建模給定原始時(shí)間2011.7.1 -2011.7.29內(nèi)100美元兌換成人民幣3的序列X(0)(k),對(duì)X(0)(k)生成1-AGOff列X(k),另外可得Yn見(jiàn)下表：k123192021X646.85646.61646.5644.26644.38644.42X646.851293.461939.9612280.6412924.6 213569.04Yn1293.461939.9612280.6412924.6 213569.04理論知識(shí)與4.1.1中的一樣，在這里不再累贅。由結(jié)果可知：C=0.3089<0.35,故可知預(yù)測(cè)結(jié)果優(yōu)。利用MA

14、TLA酎程見(jiàn)附錄7.2.1 ;運(yùn)行結(jié)果見(jiàn)附錄7.2.1結(jié)果。其中預(yù)測(cè)值如下:日 起7/307/318/18/28/38/48/58/68/78/8預(yù) 測(cè) 值644.40644.2 4644.08643.9 1643.7 5643.5 9643.43643.26643.1 0642.94（備注：預(yù)測(cè)值為100美元可以兌換人民幣數(shù)量（元）對(duì)結(jié)果的評(píng)價(jià)：由預(yù)測(cè)結(jié)果可以看出人民幣兌美元的匯率雖然存在大小波動(dòng)，但其整體 趨勢(shì)是不斷減小的。預(yù)測(cè)值也與現(xiàn)實(shí)相符，說(shuō)明人民幣升值進(jìn)一步加劇。4.2 問(wèn)題二4.2.1 首先根據(jù)10個(gè)產(chǎn)業(yè)的貿(mào)易進(jìn)出口額狀況，作出各個(gè)產(chǎn)業(yè)的貿(mào)易特化系數(shù)（ TS。, 然后作出其相應(yīng)的曲

15、線變化圖，對(duì)其進(jìn)行一定的說(shuō)明。問(wèn)題分析：將我國(guó)基礎(chǔ)產(chǎn)業(yè)分為：初級(jí)產(chǎn)品、工業(yè)制品、其他。這里只對(duì)初級(jí)產(chǎn)品、工 業(yè)制品分析（收集各產(chǎn)品的進(jìn)出口額）。出口貿(mào)易商品結(jié)構(gòu)有不同層次，在這里我們采用了貿(mào)易特化系數(shù)的變化趨勢(shì)來(lái)闡述某種貿(mào)易商品或貿(mào)易產(chǎn)業(yè)的結(jié)構(gòu)優(yōu)化。 貿(mào)易特化系 數(shù)反應(yīng)的是一國(guó)某貿(mào)易產(chǎn)業(yè)或貿(mào)易產(chǎn)品的凈出口額與該產(chǎn)品的進(jìn)出口總額之比。 其計(jì)算 公式為：TSC=(X-M)/(X+M)-1<TSGC 1當(dāng)TSC指標(biāo)值介于-1到-0.3時(shí)，則代表該國(guó)產(chǎn)業(yè)的進(jìn)口值遠(yuǎn)大于出口值，故該產(chǎn)業(yè)屬 于“進(jìn)口強(qiáng)勢(shì)”型產(chǎn)業(yè)，即該產(chǎn)業(yè)在國(guó)際市場(chǎng)上占有絕對(duì)的比較優(yōu)勢(shì)； 若TSC指標(biāo)值介于-0.3與0.3之間，表示其

16、進(jìn)口與出口值差異不大，代表該產(chǎn)業(yè)在國(guó)際 上水平較高，屬于“水平競(jìng)爭(zhēng)”型產(chǎn)業(yè)； 當(dāng)TSC指標(biāo)值介于0.3 1時(shí)，表示其出口值遠(yuǎn)大于其進(jìn)口值，故該產(chǎn)業(yè)在該國(guó)處于絕 對(duì)的劣勢(shì)，屬于“出口強(qiáng)勢(shì)”型產(chǎn)業(yè)。因此，從各種貿(mào)易行業(yè)貿(mào)易特化系數(shù)的變化我們可以大概地看出我國(guó)貿(mào)易結(jié)構(gòu)調(diào)整的方 向。表1:對(duì)外貿(mào)協(xié)商品結(jié)構(gòu)組成產(chǎn)拓名稱代號(hào)產(chǎn)品類型食品及主要供食用的活動(dòng)物SITCO初級(jí)產(chǎn)品飲料及煙類SITC1初級(jí)產(chǎn)品燃料以外的非食用粗原料SITC2初級(jí)產(chǎn)品:礦物燃料，潤(rùn)泄油及有關(guān)原料SITC3初級(jí)產(chǎn)品動(dòng)植物油脂及蠟SITC4初級(jí)產(chǎn)拈化學(xué)品及有關(guān)產(chǎn)品SITC5資本與技術(shù)密集型產(chǎn)品輕紡產(chǎn)品、橡膠制品、礦冶產(chǎn)品及其制品SITC

17、6資源與勞動(dòng)密集型產(chǎn)品機(jī)械和運(yùn)輸設(shè)飴SITC7資本與技術(shù)密集型產(chǎn)品雜項(xiàng)制品SITCS資源與勞動(dòng)密集型產(chǎn)品:未分類的其他商品SITC9未分類將所有進(jìn)出口產(chǎn)品分成兩類：初級(jí)產(chǎn)品：食品及主要供食用的活動(dòng)物、飲料及煙類、非食用原料、礦物燃料、動(dòng)植物、 油脂及蠟工業(yè)制成品：化學(xué)品及有關(guān)產(chǎn)品、輕紡產(chǎn)品、橡膠制品、礦冶產(chǎn)品及其制品、機(jī) 械及 運(yùn)輸設(shè)備、雜項(xiàng)制品MATLA褊程見(jiàn)附錄7.2.1 ,結(jié)果如下:1）初級(jí)產(chǎn)品及其組成項(xiàng)的圖像:初級(jí)產(chǎn)品貿(mào)易特化系數(shù)0.60.40.2-0.2-0.4-0.6-0.810152025從上面可以看出我國(guó)初級(jí)產(chǎn)品貿(mào)易特化曲線持續(xù)下降，說(shuō)明了該產(chǎn)業(yè)主要是通過(guò)進(jìn)口，屬于強(qiáng)進(jìn)口型，其

18、對(duì)國(guó)外市場(chǎng)依賴過(guò)強(qiáng)，因此人民幣匯率升值是有利于該產(chǎn)業(yè)發(fā)展 的。食品及主要供食用貿(mào)易特化系數(shù)0.60.50.40.30.211從上面可以看出該產(chǎn)業(yè)的貿(mào)易特化曲線基本處于正的，也就是說(shuō)，主要是向國(guó)外出 口，屬于強(qiáng)出口口型，人民幣升值會(huì)不利于此類商品在國(guó)外的發(fā)展。-0.4飲料及煙類貿(mào)易特化系數(shù)0.80.60.40.20-0.210152025從上面曲線可以看出特化曲線后一階段處于不斷下降的趨勢(shì)，該產(chǎn)業(yè)正從水平競(jìng)爭(zhēng)型向進(jìn)口型轉(zhuǎn)化的趨勢(shì)，人民幣升值對(duì)其發(fā)展是有利。非食品原料貿(mào)易特化系數(shù)0-0.1-0.2-0.3-0.4-0.5-0.6-0.7-0.8-0.9礦物，潤(rùn)滑油相關(guān)原料貿(mào)易特化系數(shù)1,0.8.-

19、0.6-r0.4-0.2-r0-_-0.2Lr-0.4 - X-0.6.-0.8 :IIC0510152025從此圖可以看出非食品原料貿(mào)易特化系數(shù)最后將會(huì)向 -1趨近，說(shuō)明了非食品原料貿(mào) 易幾乎全部依靠進(jìn)口來(lái)提供，若人民幣升值，會(huì)有利于我國(guó)進(jìn)口非食品原料。上圖反映了礦物潤(rùn)滑油相關(guān)原料的特化系數(shù)逐漸逼近 -1,說(shuō)明了我國(guó)的礦物潤(rùn)滑油 相關(guān)原料主要依靠進(jìn)口，人民幣升值將有利于我國(guó)對(duì)此產(chǎn)品的進(jìn)口與需求。從上面可以看出我國(guó)動(dòng)植物油脂及蠟的貿(mào)易特化系數(shù)逐步趨向-1,說(shuō)明了其進(jìn)口趨勢(shì)顯著加強(qiáng)，人民幣升值有利于其發(fā)展。由結(jié)果得出如下分析：初級(jí)產(chǎn)品、飲料及煙類、非食用原料、礦物燃料、動(dòng)植物油脂及蠟的貿(mào)易特化系

20、數(shù)均隨著人民幣的升值而減小， 食品及主要供食用物的貿(mào)易特化系數(shù)雖 然有上下波動(dòng)但整體比較穩(wěn)定，由此可見(jiàn)人民幣升值對(duì)我國(guó)的進(jìn)出口貿(mào)易有著密切影 響。使我國(guó)的初級(jí)產(chǎn)品類產(chǎn)業(yè)的進(jìn)口值遠(yuǎn)大于出口值， 即“進(jìn)口強(qiáng)勢(shì)”型產(chǎn)業(yè)越來(lái)越多 而食品的波動(dòng)比較平穩(wěn)，可能受國(guó)家相關(guān)政策的影響，使我國(guó)的食品價(jià)格保持平穩(wěn)。有 利于人民的生活和社會(huì)的穩(wěn)定。2）工業(yè)制成品及其組成項(xiàng)的圖像:工業(yè)制成品貿(mào)易特化系數(shù)0.30.20.10-0.1-0.2-0.3-0.4-0.50510152025從上面圖像可以看出我國(guó)工業(yè)制成品貿(mào)易特化曲線基本處于-0.3到0.3之間，說(shuō)明了此類產(chǎn)業(yè)屬于水平競(jìng)爭(zhēng)型，人民幣匯率的升高對(duì)其影響不大。但是

21、我們也可以看到在匯率 升高的同時(shí)其貿(mào)易特化系數(shù)也是呈上升趨勢(shì)，如果人民幣持續(xù)升值，該產(chǎn)業(yè)又轉(zhuǎn)化成為 出口強(qiáng)勢(shì)型的產(chǎn)業(yè)。化學(xué)品貿(mào)易特化系數(shù)-0.2-0.25-0.3-0.35-0.4-0.45-0.5-0.55從上圖可以看出，化學(xué)品相關(guān)產(chǎn)業(yè)的貿(mào)易特化系數(shù)基本上處于 -0.3 一下，即便是又 波動(dòng)影響也不大，可以認(rèn)為該產(chǎn)業(yè)是進(jìn)口強(qiáng)型產(chǎn)業(yè)，在人民幣匯率升高有利于該產(chǎn)業(yè)的 發(fā)展。從下圖可以看出：我國(guó)的輕紡橡膠、礦冶貿(mào)易特化曲線，在前期處于進(jìn)口強(qiáng)勢(shì)到了 后期逐漸轉(zhuǎn)變?yōu)槌隹趶?qiáng)勢(shì)產(chǎn)業(yè)，人民幣的匯率提高不利于其再國(guó)際市場(chǎng)的競(jìng)爭(zhēng)力。輕紡.橡膠.礦冶貿(mào)易特化系數(shù)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機(jī)械與運(yùn)輸設(shè)備貿(mào)易特化系數(shù)0.20-0.2-0.4-0.6-0.8-1101520從此圖可以看出我國(guó)機(jī)械運(yùn)輸設(shè)備貿(mào)易特化曲線一直處于上升的趨勢(shì), 會(huì)向出口型發(fā)展，人民幣的升值將會(huì)不利于其發(fā)展。25說(shuō)明了該產(chǎn)業(yè)將雜項(xiàng)制品貿(mào)易特化系數(shù)0.90.80.70.60.50.40.30.2 01015202513由結(jié)果得出如下分析：工業(yè)制成品、化學(xué)品及有關(guān)產(chǎn)品、輕紡產(chǎn)品、橡膠制品、礦 冶產(chǎn)品及其制品、機(jī)械及運(yùn)

23、輸設(shè)備的貿(mào)易特化系數(shù)均隨著人民幣的升值而增大，雜項(xiàng)制 品的貿(mào)易特化系數(shù)雖然有上下波動(dòng)但整體比較穩(wěn)定，由此可見(jiàn)人民幣升值對(duì)我國(guó)的進(jìn)出口貿(mào)易有著密切影響。是我國(guó)的工業(yè)類產(chǎn)業(yè)的進(jìn)口值遠(yuǎn)小于出口值，即“出口弱勢(shì)”型 產(chǎn)業(yè)越來(lái)越多。對(duì)比初級(jí)產(chǎn)品 和工業(yè)制成品可知：我國(guó)的工業(yè)水品低下，且只能用缺少 技術(shù)含量的初級(jí)產(chǎn)品和廉價(jià)的勞動(dòng)力來(lái)彌補(bǔ)進(jìn)出口差額，國(guó)家應(yīng)加大科技的發(fā)展，盡快 使工業(yè)水平趕上來(lái)，加大工業(yè)制成品的出口 ，減少因人民幣升值而帶來(lái)的影響。4.2.2多元回歸分析建模過(guò)程：1）模型的建立設(shè)自變量x1，X2,X3.Xp的觀測(cè)值Xii,Xi2.Xip及因變量y對(duì)應(yīng)的觀測(cè)值y i滿足關(guān)系y 1 ='

24、+£ PjXj +易 i=1,2,3- n； j=1,2,3,- p其中，當(dāng)，,.wn是相互獨(dú)立且都服從正態(tài)分布的隨機(jī)變量 由最小二乘法，得P =(XTX)(XTY)2）回歸方程的顯著性檢驗(yàn)：利用 F分布函數(shù)對(duì)變量前的系數(shù)檢驗(yàn)3）回歸系數(shù)的顯著性檢驗(yàn)：利用t分布函數(shù)檢驗(yàn)。最后得到結(jié)果如下：產(chǎn)業(yè)名稱編R方程顯著性 檢驗(yàn)（F）參數(shù)顯著 性檢驗(yàn)（t）回歸方程備注食品及主 要供食用 的活動(dòng)物10.02880.01080.0288y - 773.15277 xxx + 309.7116匯率對(duì)其 影響顯著飲料及煙 酒20.1440（不合格）匯率與此 行業(yè)沒(méi)有 顯著關(guān)系非食用原 料30.00590

25、.0010.0059y 二 一419392239 X xa + 40C.43293匯率對(duì)其 影響顯著礦物燃料4<0.0001<0.0001<0.0001y - -309.74464 xx+ + 657.01952匯率對(duì)其 影響顯著動(dòng)植物油 脂50.000960.00420.00096y =-458.67428 xxs + 3821132C匯率對(duì)其 影響顯著化學(xué)品60.1256（不合格）匯率對(duì)其 影口何不顯 著輕紡橡膠7不合格匯率對(duì)其 影口何不顯 著機(jī)械及運(yùn) 輸8<0.0001<0.0001<0.0001y = 406.07219 xxfl + 791,929

26、68匯率對(duì)其 影響顯著雜項(xiàng)90.2222（不合格）匯率對(duì)其 影口何不顯 者（備注：y代表年均人民幣匯率，100美元兌換的人民幣數(shù)量（元），Xi代表貿(mào)易特化系數(shù)）由上表看出，人民幣匯率的變化對(duì)我國(guó)的大部分產(chǎn)業(yè)都有顯著的影響。也就是說(shuō)， 人民幣持續(xù)升值對(duì)我國(guó)各行各業(yè)都會(huì)產(chǎn)生決定性改變，因此，為了穩(wěn)定中國(guó)國(guó)民經(jīng)濟(jì)穩(wěn) 步發(fā)展，在人民幣升值過(guò)程中國(guó)家需要進(jìn)行干預(yù)，嚴(yán)格控制人民幣的升值幅度，這樣才 能保證中國(guó)經(jīng)濟(jì)快速穩(wěn)步發(fā)展，中國(guó)人民能更美好的生活。4.3 問(wèn)題三 典型相關(guān)分析思想：考察一組隨機(jī)變量和另一組隨機(jī)變量之間的相關(guān)性，可以用他們之間的相關(guān)性。考慮到 兩組變量的線性組合：u1 = a11xx1+a

27、21xx2+-.+anlxxnv1 = b11xy1+b21xy2 + -.+amlxyin并研究他們之間的相關(guān)系數(shù)pvj ,在所有的線性組合中找一對(duì)相關(guān)系數(shù)最大的線 性組合，用這兩個(gè)組合的單相關(guān)系數(shù)來(lái)表示兩組變量的相關(guān)性，這兩個(gè)線性組合叫做一 對(duì)典型變量。下一步，再在兩組比昂兩的與 u, v不相關(guān)的線性組合中，找一對(duì)相關(guān)系數(shù)最大的線性 組合，設(shè)為：u2=a12xx1+iZ2xxI+_+anlxxnv2 = b12xy1+b22xy2 + -,+am2xyin且P（%，vJ就是第二個(gè)典型相關(guān)系數(shù)，這樣下去可以得到若干對(duì)典型變量，從而提取 出兩組變量間的全部信息。就本題而言，用到典型相關(guān)分析。利

28、用 SAS編程得到了：人民幣匯率增長(zhǎng)率、居民消費(fèi)指數(shù)增長(zhǎng)率與進(jìn)出口貿(mào)易額增長(zhǎng)率間的兩個(gè)典型相關(guān)系數(shù) 及特征值如下表：廳P典型相關(guān)系數(shù)標(biāo)準(zhǔn)誤差特征值方差比率10.9521250.0199259.70010.935720.6324310.1279270.66660.0643從上表可以看出：第一對(duì)典型變量的相關(guān)系數(shù)較高，說(shuō)明了這一對(duì)相關(guān)變量間的關(guān)系是 密切的。對(duì)相關(guān)系數(shù)檢驗(yàn)如下表:廳PF計(jì)算值自由度F檢驗(yàn)的顯著性概 率130.624<0.001213.3310.0016這兩對(duì)典型相關(guān)變量，進(jìn)行F統(tǒng)計(jì)檢驗(yàn)，均通過(guò)了檢驗(yàn)，說(shuō)明了這兩對(duì)典型變量間的相 關(guān)性顯著，可以用來(lái)解釋產(chǎn)出變量。對(duì)冗余度分析得到

29、下表:典型變量的解釋能力廳P產(chǎn)出組方差 被影響組解 釋比例對(duì)產(chǎn)出組的 解釋能力產(chǎn)出方差被 典型變量解 釋比例對(duì)影響組的 解釋能力影響組方差 被典型變量 的解釋比例.10.906510.95190.86290.87600.7941 120.40000.4810.01930.12400.0496從此表可以看出：這兩對(duì)典型相關(guān)變量的解釋能力均較強(qiáng)，特別是第一對(duì)相關(guān)變量對(duì)產(chǎn) 出組和影響組均有很高的解釋能力；第一對(duì)典型變量對(duì)產(chǎn)出組方差有95%Z上的解釋能力，第二對(duì)典型相關(guān)變量較弱些，達(dá)到了 48位右。最后得出的典型變量模型：廳P典型相關(guān)模型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代表居民消費(fèi)指數(shù)增長(zhǎng)百分比、z代表人民幣匯率增長(zhǎng)百分比；y1代表出口貿(mào)易額 增長(zhǎng)率、y2代表進(jìn)口貿(mào)易額的增長(zhǎng)率；）結(jié)果分析：人民幣匯率增長(zhǎng)率對(duì)進(jìn)口貿(mào)易有顯著的作用。從上面第一對(duì)典型變量分析結(jié)果可以看出：人民幣匯率增長(zhǎng)百分比在典型變量 u1占得比重達(dá)到了 0.8714,出口貿(mào)易額增長(zhǎng)率 在w1中占的比重達(dá)到0.8619.由于u1,w1由前面分析知道其相關(guān)系數(shù)達(dá)到了 0.95以上, 可以認(rèn)為人民幣匯率的增長(zhǎng)率進(jìn)口貿(mào)易額增長(zhǎng)率變化起到?jīng)Q定性的作用

31、。從第二對(duì)典型相關(guān)變量分析知道：人民幣匯率的增長(zhǎng)率、居民消費(fèi)指數(shù)增長(zhǎng)率對(duì)進(jìn) 出口貿(mào)易額增長(zhǎng)率都一定的影響，不過(guò)期間的影響程度不是那么顯著。綜上所述：人民幣匯率的增長(zhǎng)率對(duì)我國(guó)進(jìn)口貿(mào)易額有著決定性的因素，同時(shí)它和居 民消費(fèi)指數(shù)增長(zhǎng)率對(duì)我國(guó)進(jìn)出口額增長(zhǎng)率都有密切的關(guān)聯(lián)。4.4 問(wèn)題四通過(guò)上面的分析在人民幣不斷升值的條件下， 我國(guó)政府應(yīng)該在一定保證人民幣匯率 在小范圍的增長(zhǎng)。對(duì)于那些出口型的產(chǎn)業(yè)特別是工業(yè)制品，我們應(yīng)該提高產(chǎn)品的技術(shù)優(yōu) 勢(shì)，而不是通過(guò)廉價(jià)的勞動(dòng)力去爭(zhēng)取國(guó)際市場(chǎng)的優(yōu)勢(shì)。對(duì)于那些進(jìn)口型的產(chǎn)業(yè)要從提高 國(guó)內(nèi)資源的利用率和經(jīng)濟(jì)效率的角度出發(fā)，以減少對(duì)國(guó)際市場(chǎng)的依賴程度，特別是那些 礦物潤(rùn)油類。

32、最后要不斷加強(qiáng)我國(guó)產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整經(jīng)濟(jì)制度的調(diào)整，利用人民幣匯率的 小范圍的增長(zhǎng)來(lái)優(yōu)化我國(guó)產(chǎn)業(yè)，同時(shí)來(lái)促進(jìn)我國(guó)企業(yè)競(jìng)爭(zhēng)力。五、模型評(píng)價(jià)與推廣第一問(wèn)用了灰色預(yù)測(cè)，對(duì)于人民幣匯率的影響的因素太多，而之間的確切關(guān)系又無(wú) 法確定。因此用灰色系統(tǒng)可以對(duì)那些只知道一部分，而對(duì)整個(gè)其它不甚了解的狀態(tài)下， 是一種比較好的預(yù)測(cè)方法。第二問(wèn)用了典型變量分析的方法，對(duì)于多組變量間的關(guān)系是錯(cuò)綜復(fù)雜的，通過(guò)其間 的組合關(guān)系可以更清楚解釋各個(gè)組變量間的關(guān)聯(lián)程度。第三問(wèn)應(yīng)用了多元回歸，可以用來(lái)研究多個(gè)變量對(duì)一個(gè)變量的影響程度，對(duì)于生活 中的一些存在線性關(guān)系的變量間的關(guān)系討論，是一種典型的方法。六、參考文獻(xiàn)1汪曉銀 周保平，數(shù)

33、學(xué)建模與數(shù)學(xué)實(shí)驗(yàn)，北京：科學(xué)出版社，2010。2中華人民共和國(guó)國(guó)家統(tǒng)計(jì)局網(wǎng)站 3國(guó)家外匯管理局網(wǎng)站 七、附錄7.1 附錄程序:7.1.1 問(wèn)題一相關(guān)程序:基于灰色系統(tǒng)預(yù)測(cè)未來(lái)21年匯率的程序：function X=RmbExRata()clear all;clc;人民幣匯率(年平均價(jià))','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、年預(yù)測(cè)值');hold onplot(1994:2014,X);end%function X,c,error1,error2=GM11(X0,k)% 建立函數(shù)X,c,error1,error2=example9_3_2_3(X0,k)%其中X0為輸入序列，k為預(yù)測(cè)長(zhǎng)度，% X為預(yù)測(cè)輸出序列，c為后驗(yàn)差檢驗(yàn)數(shù)，error1為殘差，error2為相對(duì)誤差 format;n=length(X0);X1=;X1(1)=X0(1);for i=2:nX1(i)=X1(i-1)+X0(i);%計(jì)算累加生成序列endfor i=1:n-1B(i,1)=-0.5*(X1(i)+X1(i+1); %計(jì)算

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做最小二乘估計(jì)計(jì)算時(shí)間響應(yīng)函數(shù)參數(shù)計(jì)算預(yù)測(cè)序列計(jì)算預(yù)測(cè)序列計(jì)算殘差計(jì)算相對(duì)誤差計(jì)算后驗(yàn)差檢驗(yàn)數(shù)基于灰色系統(tǒng)預(yù)測(cè)未來(lái)31天匯率的程序：function X=RmbExRata1()clear all;clc;最近一個(gè)月人民幣匯率中間價(jià),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('人民幣匯率月變化及十天預(yù)測(cè)值);hold onY=Y 30 31 32 33 34 35 36 37 38 39;plot(Y,X);end%function X,c,erro1error2=GM11(X0,k)% 建立函數(shù)X,c,erro1error2=example9_3_2_3(X0,k)%其中X0為輸入序列，k為預(yù)測(cè)長(zhǎng)度，%X為預(yù)測(cè)輸出序列，c為后驗(yàn)差檢驗(yàn)數(shù)，

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為相對(duì)誤差計(jì)算累加生成序列計(jì)算B, Yn做最小二乘估計(jì)計(jì)算時(shí)間響應(yīng)函數(shù)參數(shù)計(jì)算預(yù)測(cè)序列計(jì)算預(yù)測(cè)序列計(jì)算殘差計(jì)算相對(duì)誤差計(jì)算后驗(yàn)差檢驗(yàn)數(shù)7.1.2 問(wèn)題二相關(guān)程序: function TSCc

40、lear all;clc;%J級(jí)產(chǎn)品的處理歷年出口貨物分類總額 歷年進(jìn)口貨物分類總額 n=length(CJCKZE);%M級(jí)產(chǎn)品貿(mào)易特化系數(shù)的計(jì)算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('初級(jí)產(chǎn)品貿(mào)易特化系數(shù))；%t品及主要供食用的產(chǎn)品的處理歷年出口貨物分類總額 歷年進(jìn)口貨物分類總額 n=length(SPCKZE);%t品及主要供食用的貿(mào)易特化系數(shù)的計(jì)算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('食品及主要供食用貿(mào)易特化系數(shù));%須料及煙類的產(chǎn)品的處理歷年出口貨物分類總額 歷年進(jìn)口貨物分類總額 n=length(YYCKZE);%須料及煙類的貿(mào)易特化系數(shù)的計(jì)算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('飲料及煙類貿(mào)易特化系數(shù))；%E食用原料的產(chǎn)品的處理歷年出口

43、貨物分類總額歷年進(jìn)口貨物分類總額n=length(FSPCKZE);%E食品原料的貿(mào)易特化系數(shù)的計(jì)算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(非食品原料貿(mào)易特化系數(shù));%T物，潤(rùn)滑油相關(guān)原料的產(chǎn)品的處理歷年出口貨物分類總額歷年進(jìn)口貨物分類總額n=length(KRCKZE);%T物，潤(rùn)滑油相關(guān)原料的貿(mào)易特化系數(shù)的計(jì)算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('礦物，潤(rùn)滑油相關(guān)原料貿(mào)易特化系數(shù));%癡植物油脂及蠟相關(guān)原料的產(chǎn)品的處理歷年出口貨物分類總額歷年進(jìn)口貨物分類總額n=length(DZYCKZE);%癡植物油脂及蠟的貿(mào)易特化系數(shù)的計(jì)算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('動(dòng)植物油脂及蠟貿(mào)易特化系數(shù));','Sheet1','F9:F30');','Sheet1','F9:F30');','Sheet1','G9:G30');','Sheet1','G9:G30');','Sheet1','H9:H30');','Sheet1','H9:H30');end22工業(yè)制成品的貿(mào)易特

46、化系數(shù)（TS。的程序:function TSC1clear all;clc;歷年出口貨物分類總額歷年進(jìn)口貨物分類總額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('工業(yè)制成品貿(mào)易特化系數(shù))；歷年出口貨物分類總額歷年進(jìn)口貨物分類總額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('化學(xué)品貿(mào)易特化系數(shù));歷年出口貨物分類總額歷年進(jìn)口貨物分類總額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(' 輕紡.橡月里.礦冶貿(mào)易特化系數(shù));歷年出口貨物分類總額歷年進(jìn)口貨物分類總額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('機(jī)械與運(yùn)輸設(shè)備貿(mào)易特化系數(shù));歷年出口貨物分類總額,Sheet1','M9:M30');歷年進(jìn)口貨物分類總額,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('雜項(xiàng)制品貿(mào)易特化系數(shù));end 建立多元