【正文】 參考數(shù)列為，比較數(shù)列為 。越小，分辨率越小。根據(jù)灰色加權(quán)關(guān)聯(lián)度的大小，對各評價對象進行排序，可建立評價對象的關(guān)聯(lián)序，關(guān)聯(lián)度越大，其影響程度越大。特別是依據(jù)目前已有的數(shù)據(jù)對未來的發(fā)展趨勢做出預(yù)測分析。 若為參數(shù)列, 且 ⑤則GM (1 , 1) 模型的最小二乘估計參數(shù)列滿足： ⑥求灰色微分倒數(shù) ⑦求解, 從而獲得時間響應(yīng)函數(shù)： ⑧獲取時間響應(yīng)序列為： ⑨ 模型的求解根據(jù)題目要求將江蘇省主要能源消費量匯總到同一張表中，并取2000年到2010年的數(shù)據(jù)，見表51 表51 江蘇省主要能源（煤炭、石油、天然氣等）的消費量年份煤炭消費量石油消費量天然氣消費量12345678910煤炭消費量 根據(jù)表51中煤炭消費量的數(shù)據(jù)，依次累加得到取代入④得到 運用最小二乘法⑥得 =[], a=。 對問題3的分析與求解 與第二小問一樣，我們采用灰色預(yù)測GM（1,1）對未來十年江蘇省煤炭消費總量以及第一、二、三產(chǎn)業(yè)煤炭消費量進行預(yù)測，并用模型對原始數(shù)據(jù)進行估計和驗證，求出誤差，找出模型的靈敏度。b=21634第三產(chǎn)業(yè)煤炭消耗量根據(jù)表21中第三產(chǎn)業(yè)煤炭消耗量的數(shù)據(jù)，一次累加得到取代入④得到運用最小二乘法⑥得=[],a=。所以調(diào)整優(yōu)化能源結(jié)構(gòu)是江蘇省在未來幾年中需要重視的一項工作。3　 保護環(huán)境，走可持續(xù)發(fā)展道路江蘇省政府可通過制定相關(guān)的條款來制約企業(yè)發(fā)展所帶來的環(huán)境污染，如誰污染，誰治理。 模型的分析216。在各種預(yù)測領(lǐng)域都有著廣泛的運用，是處理小樣本預(yù)測問題的有效工具；4　 多目標(biāo)規(guī)劃模型也就是優(yōu)化模型能更好地兼顧統(tǒng)籌處理多種目標(biāo)的關(guān)系，在整個碳排放約束下，考慮節(jié)能目標(biāo)和江蘇省經(jīng)濟發(fā)展目標(biāo)求得更切合實際要求的解；5　 關(guān)聯(lián)度分析發(fā)是以發(fā)展態(tài)勢為立足點，對樣本量的多少沒有過分的要求，也不需要典型的分布規(guī)律，為我們的操作減輕了負擔(dān)；216。它不受傳統(tǒng)模型對原始數(shù)據(jù)種種要求的約束，具有實用性強、預(yù)測性能好的優(yōu)點是由于疾病預(yù)測，以探求一種較為實用的傳染病定量預(yù)測方法。[J].山東：山東大學(xué)學(xué)報(理學(xué)版)，2007[3][5]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 ]。A=[ 6900 7319 150 7338 8471 6548 9371 12 6791 10049 11 8118 6977 7098 10695 6763 11765 7097 12879 7678 14369 8525 16743 11061 715960 11020 8438 5198 20031 13652 14496 10525 24616 16895 15811 12283 28526 18742 63 17948 14471 33837 20604 19467 16993 11889 40014 21776 63 18706 18566 13629 44253 23709 21655 21754 17131 52840 25774 ][x,y]=size(A)for j=1:yD(:,j)=A(:,j)./A(:,1)endA=D。[x,y]=size(A)for j=1:yD(:,j)=A(:,j)./A(:,1)。C=(Amin+*Amax)./(B+*Amax)%關(guān)聯(lián)系數(shù)的求解for i=1:y1a(i)=(1/x)*sum(C(:,i))%灰色關(guān)聯(lián)度數(shù)End程序結(jié)果：C = Columns 1 through 5 Columns 6 through 8 a = Columns 1 through 5 Columns 6 through 8 程序4第三產(chǎn)業(yè)煤炭消耗量clear all。endAmax=max(max(B))。lamda=x0(1:n1)./x0(2:n) %計算級比range=minmax(lamda39。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。)。b39。 %代入估計參數(shù)值和初始值yuce1=subs(x,39。) %計算相對誤差rho=1(*u(1))/(1+*u(1))*lamda39。) %計算級比的范圍x1=cumsum(x0)。Dx+a*x=b39。 %求微分方程的符號解x=subs(x,{39。,39。t39。 %計算級比偏差值，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天然氣消費量clc,clearx0=[]。 %累加運算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。,39。a39。x039。,[0:n1])%求已知數(shù)據(jù)的預(yù)測值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運算，還原數(shù)據(jù)epsilon=x039。 %注意這里為列向量n=length(x0)。Y=x0(2:n)。x(0)=x039。,39。},{u(1),u(2),x0(1)})。yuce %計算殘差delta=abs(epsilon./x039。lamda=x0(1:n1)./x0(2:n) %計算級比range=minmax(lamda39。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。)。b39。 %代入估計參數(shù)值和初始值yuce1=subs(x,39。) %計算相對誤差rho=1(*u(1))/(1+*u(1))*lamda39。) %計算級比的范圍x1=cumsum(x0)。Dx+a*x=b39。 %求微分方程的符號解x=subs(x,{39。,39。t39。 %計算級比偏差值，u(1)=a程序結(jié)果：yuce = +004 * Columns 1 through 5 Columns 6 through 10 epsilon = +003 * 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 程序11第三產(chǎn)業(yè)煤炭消費量clc,clearx0=[]。 %累加運算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。,39。a39。x039。,[0:n1])%求已知數(shù)據(jù)的預(yù)測值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運算，還原數(shù)據(jù)epsilon=x