【正文】 ,[0:n1])%求已知數(shù)據(jù)的預(yù)測(cè)值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運(yùn)算，還原數(shù)據(jù)epsilon=x039。a39。 %累加運(yùn)算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。t39。 %求微分方程的符號(hào)解x=subs(x,{39。) %計(jì)算級(jí)比的范圍x1=cumsum(x0)。 %代入估計(jì)參數(shù)值和初始值yuce1=subs(x,39。)。lamda=x0(1:n1)./x0(2:n) %計(jì)算級(jí)比range=minmax(lamda39。},{u(1),u(2),x0(1)})。x(0)=x039。 %注意這里為列向量n=length(x0)。x039。,39。 %計(jì)算級(jí)比偏差值，u(1)=a程序結(jié)果：yuce = +003 * Columns 1 through 5 Columns 6 through 10 epsilon = Columns 1 through 5 0 Columns 6 through 10 delta = Columns 1 through 5 0 Columns 6 through 10 rho = Columns 1 through 5 Columns 6 through 9 程序7天然氣消費(fèi)量clc,clearx0=[]。,39。Dx+a*x=b39。) %計(jì)算相對(duì)誤差rho=1(*u(1))/(1+*u(1))*lamda39。b39。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。endAmax=max(max(B))。[x,y]=size(A)for j=1:yD(:,j)=A(:,j)./A(:,1)。A=[8936 6900 7319 8833 7338 8471 8585 6548 9371 12 8571 6791 10049 11 8118 69778714 7098 10695 8770 6763 11765 8963 7097 12879 9663 7678 14369 10849 8525 16743 11061 715913272 11020 8438 5198 20031 13652 16779 14496 10525 24616 16895 18428 15811 12283 28526 18742 19952 17948 14471 33837 20604 20737 19467 16993 11889 40014 21776 21003 18706 18566 13629 44253 23709 23100 21655 21754 17131 52840 25774 ]。[J].山東：山東大學(xué)學(xué)報(bào)(理學(xué)版)，2007[3]在各種預(yù)測(cè)領(lǐng)域都有著廣泛的運(yùn)用，是處理小樣本預(yù)測(cè)問(wèn)題的有效工具；4　 多目標(biāo)規(guī)劃模型也就是優(yōu)化模型能更好地兼顧統(tǒng)籌處理多種目標(biāo)的關(guān)系，在整個(gè)碳排放約束下，考慮節(jié)能目標(biāo)和江蘇省經(jīng)濟(jì)發(fā)展目標(biāo)求得更切合實(shí)際要求的解；5　 關(guān)聯(lián)度分析發(fā)是以發(fā)展態(tài)勢(shì)為立足點(diǎn)，對(duì)樣本量的多少?zèng)]有過(guò)分的要求，也不需要典型的分布規(guī)律，為我們的操作減輕了負(fù)擔(dān)；216。3　 保護(hù)環(huán)境，走可持續(xù)發(fā)展道路江蘇省政府可通過(guò)制定相關(guān)的條款來(lái)制約企業(yè)發(fā)展所帶來(lái)的環(huán)境污染，如誰(shuí)污染，誰(shuí)治理。b=21634第三產(chǎn)業(yè)煤炭消耗量根據(jù)表21中第三產(chǎn)業(yè)煤炭消耗量的數(shù)據(jù)，一次累加得到取代入④得到運(yùn)用最小二乘法⑥得=[],a=。 若為參數(shù)列, 且 ⑤則GM (1 , 1) 模型的最小二乘估計(jì)參數(shù)列滿足： ⑥求灰色微分倒數(shù) ⑦求解, 從而獲得時(shí)間響應(yīng)函數(shù)： ⑧獲取時(shí)間響應(yīng)序列為： ⑨ 模型的求解根據(jù)題目要求將江蘇省主要能源消費(fèi)量匯總到同一張表中，并取2000年到2010年的數(shù)據(jù)，見(jiàn)表51 表51 江蘇省主要能源（煤炭、石油、天然氣等）的消費(fèi)量年份煤炭消費(fèi)量石油消費(fèi)量天然氣消費(fèi)量12345678910煤炭消費(fèi)量 根據(jù)表51中煤炭消費(fèi)量的數(shù)據(jù)，依次累加得到取代入④得到 運(yùn)用最小二乘法⑥得 =[], a=。根據(jù)灰色加權(quán)關(guān)聯(lián)度的大小，對(duì)各評(píng)價(jià)對(duì)象進(jìn)行排序，可建立評(píng)價(jià)對(duì)象的關(guān)聯(lián)序，關(guān)聯(lián)度越大，其影響程度越大?；疑P(guān)聯(lián)度分析具體步驟如下：（1） 由表一知，確定的因素有m個(gè)，所選數(shù)據(jù)有n個(gè)，參考數(shù)列為，比較數(shù)列為 。做出折線圖41：由上折線圖可知從1995年至2010年人均地區(qū)生產(chǎn)總值逐年遞增，是影響該地區(qū)生產(chǎn)總值的主要因素。短序列灰色增量模型還可以對(duì)全國(guó)人口的檢驗(yàn)與預(yù)測(cè)。請(qǐng)對(duì)預(yù)測(cè)模型和預(yù)測(cè)結(jié)果的合理性進(jìn)行檢驗(yàn)和說(shuō)明。題目要求我們找出影響江蘇省煤炭消費(fèi)總量以及第一、二、三產(chǎn)業(yè)煤炭消費(fèi)量的指標(biāo)及其影響情況，江蘇省煤炭消費(fèi)總量以及第一、二、三產(chǎn)業(yè)煤炭消費(fèi)量，并在此基礎(chǔ)上考慮節(jié)能目標(biāo)和江蘇省經(jīng)濟(jì)發(fā)展目標(biāo)，預(yù)測(cè)和優(yōu)化未來(lái)十年江蘇省主要能源消費(fèi)的結(jié)構(gòu)。 關(guān)鍵字：灰色關(guān)聯(lián)分析法、灰色系統(tǒng)理論GM(1,1)模型、低碳經(jīng)濟(jì)、 煤炭消費(fèi)量預(yù)測(cè)、關(guān)聯(lián)度分析、優(yōu)化、Excel、matlab在選擇數(shù)據(jù)的過(guò)程中，我們有意識(shí)只保留了影響較大的指標(biāo)（即1995年以后的數(shù)據(jù)）。因此，解決好未來(lái)江蘇巨大的煤炭供需缺口，分析預(yù)測(cè)江蘇未來(lái)的，可以為江蘇戰(zhàn)略性能源開(kāi)發(fā)供應(yīng)提供依據(jù)，從而保證全省經(jīng)濟(jì)社會(huì)的綠色發(fā)展。 針對(duì)問(wèn)題四，我們根據(jù)前幾道題的結(jié)果，從產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)、環(huán)境保護(hù)等方面提出了自己的意見(jiàn)和建議。江蘇省為響應(yīng)國(guó)際組織的號(hào)召，制訂了碳排放約束指標(biāo)，見(jiàn)表11. 表11 江蘇省碳排放約束指標(biāo)碳排放約束指標(biāo)名稱(chēng)單位地區(qū)生產(chǎn)總值能耗降低（%）單位地區(qū)生產(chǎn)總值二氧化碳排放減少（%）非化石能源占一次能源消費(fèi)比重（％）指標(biāo)值2015年18197左右2020年50注：2015年指標(biāo)值是以2010年為基期，2020年指標(biāo)值是以2005年為基期。根據(jù)題目的相關(guān)信息，我們覺(jué)得用最近十年的相關(guān)信息去預(yù)測(cè)未來(lái)十年更具有代表性，于是我們舍棄了2000年以前的相關(guān)數(shù)據(jù)。 對(duì)問(wèn)題4：我們根據(jù)前幾道題的結(jié)果，從產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)、環(huán)境保護(hù)等方面提出了自己的意見(jiàn)和建議。；、石油、天然氣這些燃料燃燒排放；4． 符號(hào)說(shuō)明序號(hào)符號(hào)（可帶單位）符號(hào)說(shuō)明1 分別表示、第二產(chǎn)業(yè)產(chǎn)量、總?cè)丝跀?shù)、能源消費(fèi)總量、煤炭生產(chǎn)量、煤炭消費(fèi)量,2 表示、第二產(chǎn)業(yè)產(chǎn)量、總?cè)丝跀?shù)、能源消費(fèi)總量、煤炭生產(chǎn)量、煤炭消費(fèi)量的系數(shù)3 表示、第二產(chǎn)業(yè)產(chǎn)量、總?cè)丝跀?shù)、能源消費(fèi)總量、煤炭生產(chǎn)量、煤炭消費(fèi)量的灰色關(guān)聯(lián)系數(shù)4為第i個(gè)評(píng)價(jià)對(duì)象對(duì)理想對(duì)象的灰色加權(quán)關(guān)聯(lián)度5參數(shù)列6,B矩陣5． 模型的建立與求解 對(duì)問(wèn)題1的分析與求解 我們首先通過(guò)觀察數(shù)據(jù)之間的不同關(guān)系，從中選擇出我們認(rèn)為有用的數(shù)據(jù)，在選擇數(shù)據(jù)的過(guò)程中，我們有意識(shí)的刪除了一些不重要的數(shù)據(jù)和影響較小的數(shù)據(jù)，然后得到影響較大的指標(biāo)，將這些指標(biāo)匯總到一起，利用灰色關(guān)聯(lián)分析法得到它們的關(guān)聯(lián)系數(shù)，并計(jì)算灰色加權(quán)關(guān)聯(lián)度，根據(jù)灰色加權(quán)關(guān)聯(lián)度的大小，對(duì)各評(píng)價(jià)對(duì)象進(jìn)行排序，建立評(píng)價(jià)對(duì)象的關(guān)聯(lián)序，關(guān)聯(lián)度越大，其評(píng)價(jià)結(jié)果越好。為了使資源的平衡利用，應(yīng)該使用天然氣代替部分煤炭資源，天然氣既經(jīng)濟(jì)又環(huán)保，是清潔的綠色能源，值得推廣。（4） 計(jì)算灰色加權(quán)關(guān)聯(lián)度。 預(yù)測(cè)模型建立與求解（ 模型） 由已知數(shù)據(jù)，對(duì)于20002010 年某項(xiàng)指標(biāo)記為 ② 其中序列數(shù)據(jù)累加生成。 根據(jù)題目要求將江蘇省煤炭消費(fèi)總量以及第一、二、三產(chǎn)業(yè)煤炭消費(fèi)量匯總到同一張表中，并取2000年到2010年的數(shù)據(jù)，見(jiàn)表54表54 江蘇省煤炭消費(fèi)總量以及第一、二、三產(chǎn)業(yè)煤炭消費(fèi)量年份能源第一產(chǎn)業(yè)煤炭消耗量第二產(chǎn)業(yè)煤炭消耗量第三產(chǎn)業(yè)煤炭消耗量123460567638923709631025774能源根據(jù)表21中能源的數(shù)據(jù)，依次累加得到取代入④得到 運(yùn)用最小二乘法⑥得=[],a=。逐步減小煤炭能源的比重，大幅增加其他能源（如石油、天然氣、風(fēng)能，水能、非化石能源等）的比重，從理論上講，這樣既可以保證經(jīng)濟(jì)增長(zhǎng)，又可以實(shí)現(xiàn)節(jié)能減排的目標(biāo)，達(dá)到約束條件。 優(yōu)點(diǎn)：1　 計(jì)算方便，借助excel，matlab軟件使得數(shù)據(jù)處理簡(jiǎn)便、快捷、準(zhǔn)確性好。短序列灰色增量模型還可以對(duì)全國(guó)人口的檢驗(yàn)與預(yù)測(cè)。姜啟源，謝金星，葉俊，數(shù)學(xué)模型（第四版），北京高等教育出版社，2011，1，第4版[6]for i=1:y1B(:,i)=abs(A(:,1)A(:,i+1))。A=[ 6900 7319 7338 8471 6548 9371 12 53 6791 10049 11 8118 6977 7098 10695 6763 11765 7097 12879 7678 14369 8525 16743 11061 7159 11020 8438 5198 20031 13652 14496 10525 24616 16895 15811 12283 28526 18742 17948 14471 33837 20604 19467 16993 11889 40014 21776 18706 18566 13629 44253 23709 21655 21754 17131 52840 25774 ]。) %計(jì)算級(jí)比的范圍x1=cumsum(x0)。 %求微分方程的符號(hào)解x=subs(x,{39。t39。 %累加運(yùn)算B=[*(x1(1:n