灰色系統(tǒng)關(guān)聯(lián)和GM模型

1、灰色系統(tǒng)理論灰色系統(tǒng)理論（GM和關(guān)聯(lián)度分析）和關(guān)聯(lián)度分析） 1.本節(jié)內(nèi)容（1）灰色系統(tǒng)理論的產(chǎn)生和發(fā)展動(dòng)態(tài)；（2）灰色系統(tǒng)的基本概念；（3）灰色系統(tǒng)與模糊數(shù)學(xué)、黑箱方法的區(qū)別；（4）灰色系統(tǒng)GM（1，1）模型，灰色系統(tǒng)模型的檢驗(yàn)，應(yīng)用舉例 (G表示gray（灰色），M表示model（模型），GM（1，1）表示1階的、1個(gè)變量的模型) 。1.1 灰色系統(tǒng)理論的產(chǎn)生和發(fā)展 動(dòng)態(tài) 1982我國(guó)學(xué)者鄧聚龍教授發(fā)表第一篇中文論文灰色控制系統(tǒng)標(biāo)志著灰色系統(tǒng)這一學(xué)科誕生。1985灰色系統(tǒng)研究會(huì)成立，灰色系統(tǒng)相關(guān)研究迅速發(fā)展。1989海洋出版社出版英文版灰色系統(tǒng)論文集，同年，英文版國(guó)際刊物灰色系統(tǒng)雜志正式創(chuàng)刊

2、。目前，國(guó)際、國(guó)內(nèi)200多種期刊發(fā)表灰色系統(tǒng)論文，許多國(guó)際會(huì)議把灰色系統(tǒng)列為討論專題。國(guó)際著名檢索已檢索我國(guó)學(xué)者的灰色系統(tǒng)論著500多次。灰色系統(tǒng)理論應(yīng)用范圍已拓展到工業(yè)、農(nóng)業(yè)、社會(huì)、經(jīng)濟(jì)、能源、地質(zhì)、石油等眾多科學(xué)領(lǐng)域，成功地解決了生產(chǎn)、生活和科學(xué)研究中的大量實(shí)際問題，取得了顯著成果。 1.2灰色系統(tǒng)的基本原理 1.2.1灰色系統(tǒng)的基本概念我們將信息完全明確的系統(tǒng)稱為白色系統(tǒng)，信息未知的系統(tǒng)稱為黑色系統(tǒng)，部分信息明確、部分信息不明確的系統(tǒng)稱為灰色系統(tǒng)灰色系統(tǒng)。系統(tǒng)信息不完全的情況有以下四種：（1）元素信息不完全。（2）結(jié)構(gòu)信息不完全。（3）邊界信息不完全。（4）運(yùn)行行為信息不完全。 1.2.

3、2灰色系統(tǒng)與模糊數(shù)學(xué)、黑箱方法的區(qū)別 區(qū)別主要在于對(duì)系統(tǒng)內(nèi)涵與外延處理態(tài)度不同；研究對(duì)象內(nèi)涵與外延的性質(zhì)不同。研究對(duì)象內(nèi)涵與外延的性質(zhì)不同?；疑到y(tǒng)著重外延明確、內(nèi)涵不明確的對(duì)象，模糊數(shù)學(xué)著重外延不明確、內(nèi)涵明確的對(duì)象。“黑箱”方法著重系統(tǒng)外部行為數(shù)據(jù)的處理方法，是揚(yáng)外延而棄內(nèi)涵的處理方法，而灰色系是揚(yáng)外延而棄內(nèi)涵的處理方法，而灰色系統(tǒng)方法是外延內(nèi)涵均注重的方法。統(tǒng)方法是外延內(nèi)涵均注重的方法。 1.2.3灰色系統(tǒng)的基本原理 公理1：差異信息原理?！安町悺笔切畔ⅲ残畔⒈赜胁町?。公理2：解的非唯一性原理。信息不完全，不明確地解是非唯一的。公理3：最少信息原理?；疑到y(tǒng)理論的特點(diǎn)是充分開發(fā)利用已有

4、的“最少信息”。 公理4：認(rèn)知根據(jù)原理。信息是認(rèn)知的根據(jù)。公理5：新信息優(yōu)先原理。新信息對(duì)認(rèn)知的作用大于老信息。公理6：灰性不滅原理。“信息不完全”是絕對(duì)的。 1.2.4灰色系統(tǒng)理論的主要內(nèi)容 灰色系統(tǒng)理論經(jīng)過10多年的發(fā)展，已基本建立起了一門新興學(xué)科的結(jié)構(gòu)體系，其主要內(nèi)容包括以“灰色朦朧集”為基礎(chǔ)的理論體系、以灰色關(guān)聯(lián)空間為依托的分析體系、以灰色序列生成為基礎(chǔ)的方法體系，以灰色模型（G，M）為核心的模型體系。以系統(tǒng)分析、評(píng)估、建模、預(yù)測(cè)、決策、控制、優(yōu)化為主體的技術(shù)體系。1.2.5灰色系統(tǒng)的應(yīng)用范疇 灰色系統(tǒng)的應(yīng)用范疇大致分為以下幾方面：（1）灰色關(guān)聯(lián)分析。（2）灰色預(yù)測(cè)：人口預(yù)測(cè)；初霜預(yù)測(cè)

5、；災(zāi)變預(yù)測(cè).等等。（3）灰色決策。（4）灰色預(yù)測(cè)控制。1.3 灰色關(guān)聯(lián)分析法灰色關(guān)聯(lián)分析法 灰色關(guān)聯(lián)分析是灰色系統(tǒng)理論的一個(gè)分支應(yīng)用灰色關(guān)聯(lián)分析方法對(duì)受多種因應(yīng)用灰色關(guān)聯(lián)分析方法對(duì)受多種因素影響的事物和現(xiàn)象從整體觀念出發(fā)進(jìn)素影響的事物和現(xiàn)象從整體觀念出發(fā)進(jìn)行綜合評(píng)價(jià)是一個(gè)被廣為接受的方法行綜合評(píng)價(jià)是一個(gè)被廣為接受的方法 1.3.11.3.1灰色關(guān)聯(lián)分析法的建模過灰色關(guān)聯(lián)分析法的建模過程和機(jī)理為程和機(jī)理為 利用灰色關(guān)聯(lián)分析進(jìn)行綜合評(píng)價(jià)的步驟是：1根據(jù)評(píng)價(jià)目的確定評(píng)價(jià)指標(biāo)體系，收集評(píng)價(jià)數(shù)據(jù)。設(shè)n個(gè)數(shù)據(jù)序列形成如下矩陣： mxmxmxxxxxxxXXXnnnn21212121222111,其中 為指

6、標(biāo)的個(gè)數(shù)， m nimxxxXTiiii,2,1,2,12確定參考數(shù)據(jù)列 參考數(shù)據(jù)列應(yīng)該是一個(gè)理想的比較標(biāo)準(zhǔn)，可以以各指標(biāo)的最優(yōu)值 （或最劣值）構(gòu)成參考數(shù)據(jù)列，也可根據(jù)評(píng)價(jià)目的選擇其它參照值記作 mxxxX0000,2,)1(3對(duì)指標(biāo)數(shù)據(jù)進(jìn)行無量綱化 無量綱化后的數(shù)據(jù)序列形成如下矩陣： mxmxmxxxxxxxXXXnnnn10101010222111,常用的無量綱化方法有均值化法（見（123）式）、初值化法（見（124）式）和 變換等sxx .,2,1,1,0)412(1) 312(11mknixkxkxkxmkxkxiiimkiii；或采用內(nèi)插法使各指標(biāo)數(shù)據(jù)取值范圍（或數(shù)量級(jí)）相同例如，某

7、地縣級(jí)醫(yī)院病床使用率最高為90%，最低為60%，我們可以將90%轉(zhuǎn)化10，60%轉(zhuǎn)化為1，其它可以通過內(nèi)插法確定其轉(zhuǎn)化值如80%轉(zhuǎn)化為多少？可進(jìn)行如下計(jì)算： 解之得，即80%轉(zhuǎn)化為7608060901110 x4逐個(gè)計(jì)算每個(gè)被評(píng)價(jià)對(duì)象指標(biāo)序列（比較序列）與參考序列對(duì)應(yīng)元素的絕對(duì)差值 即 （ , 為被評(píng)價(jià)對(duì)象的個(gè)數(shù)）5確定 與)()(0kxkximk, 1ni, 1n)()(minmin011kxkximkni)()(maxmax011kxkximkni6計(jì)算關(guān)聯(lián)系數(shù) 由（125）式，分別計(jì)算每個(gè)比較序列與參考序列對(duì)應(yīng)元素的關(guān)聯(lián)系數(shù)0000minmin( )( )maxmax( )( )( )1

8、25)( )( )maxmax( )( )iiikikiiiikxkx kx kx kkx kx kx kx k（mk, 10.5式中 為分辨系數(shù)，在（0，1）內(nèi)取值，若 越小，關(guān)聯(lián)系數(shù)間差異越大，區(qū)分能力越強(qiáng)。通常 取當(dāng)用各指標(biāo)的最優(yōu)值 （或最劣值），構(gòu)成參考數(shù)據(jù)列計(jì)算關(guān)聯(lián)系數(shù)時(shí)，也可用改進(jìn)的更為簡(jiǎn)便的計(jì)算方法：改進(jìn)后的方法不僅可以省略第三步，使計(jì)算簡(jiǎn)便，而且避免了無量綱化對(duì)指標(biāo)作用的某些負(fù)面影響)()(max)()()()(max)()(min)(0000kxkxkxkxkxkxkxkxkiiiiiiiimk, 17計(jì)算關(guān)聯(lián)序 對(duì)各評(píng)價(jià)對(duì)象（比較序列）分別計(jì)算其個(gè)指標(biāo)與參考序列對(duì)應(yīng)元素的關(guān)

9、聯(lián)系數(shù)的均值，以反映各評(píng)價(jià)對(duì)象與參考序列的關(guān)聯(lián)關(guān)系，并稱其為關(guān)聯(lián)序，記為：00 ( )( ) ( )( )x kkx kkii如果為最優(yōu)值數(shù)據(jù)列，越大，越好；如果為最劣值數(shù)據(jù)列，越大，越不好。011( )miikrkm8如果各指標(biāo)在綜合評(píng)價(jià)中所起的作用不同，可對(duì)關(guān)聯(lián)系數(shù)求加權(quán)平均值即 9依據(jù)各觀察對(duì)象的關(guān)聯(lián)序，得出綜合評(píng)價(jià)結(jié)果011( ),mikikkrWkmmW （k=1,）式中為各指標(biāo)權(quán)重。2.2.灰色關(guān)聯(lián)分析的應(yīng)用舉例灰色關(guān)聯(lián)分析的應(yīng)用舉例 例1：利用灰色關(guān)聯(lián)分析對(duì)6位教師工作狀況進(jìn)行綜合評(píng)價(jià)1評(píng)價(jià)指標(biāo)包括：專業(yè)素質(zhì)、外語(yǔ)水平、教學(xué)工作量、科研成果、論文、著作與出勤2對(duì)原始數(shù)據(jù)經(jīng)處理后得

10、到以下數(shù)值，見下表 編號(hào)編號(hào)專業(yè)專業(yè)外語(yǔ)外語(yǔ)教學(xué)教學(xué)量量科研科研論文論文著作著作出勤出勤1898752927875738397966474688843658669838689576483確定參考數(shù)據(jù)列： 4計(jì)算 ， 見下表09,9,9,9,8,9,9x)()(0kxkxi編號(hào)編號(hào)專業(yè)專業(yè)外語(yǔ)外語(yǔ)教學(xué)教學(xué)量量科研科研論文論文著作著作出勤出勤1101237022124161302032524311146351330061610422515求最值6依據(jù)（125）式， 取計(jì)算，得 011minmin( )( )min(0,1,0,1,0,0)0nmiikx kx k011maxmax( )( )max

11、(7,6,5,6,6,5)7nmiikx kx k111111100.5 700.5 7(1)0.778(2)1.0001 0.5 700.5 7(3)0.778(4)0.636(5)0.467(6)0.333(7)， ，1.000，0.5同理得出其它各值，見下表編號(hào)10.7781.0000.7780.6360.4670.3331.00020.636 0.778 0.636 0.467 0.636 0.368 0.778 31.000 0.636 1.000 0.538 0.538 0.412 0.636 40.538 0.778 0.778 0.778 0.412 0.368 0.538 5

12、0.778 0.538 0.538 1.000 0.778 0.368 0.778 60.778 1.000 0.467 0.636 0.538 0.412 0.778 (1)i(2)i(3)i(4)i(5)i(6)i(7)i7分別計(jì)算每個(gè)人各指標(biāo)關(guān)聯(lián)系數(shù)的均值（關(guān)聯(lián)序）： 8如果不考慮各指標(biāo)權(quán)重（認(rèn)為各指標(biāo)同等重要），六個(gè)被評(píng)價(jià)對(duì)象由好到劣依次為1號(hào)，5號(hào)，3號(hào)，6號(hào)，2號(hào)，4號(hào)即 713. 07000. 1333. 0467. 0636. 0778. 0000. 1778. 001r02030405060.6140.6800.5990.6830.658rrrrr，010503060204r

13、rrrrr灰色系統(tǒng)建模 1 灰色系統(tǒng)理論概述2 灰色GM(1.1)模型3 序列光滑度的理論分析4 灰色GM(1.1)優(yōu)化模型分析 5 灰色模型的應(yīng)用1 1 灰色系統(tǒng)概述灰色系統(tǒng)概述1.1 灰色系統(tǒng)理論的產(chǎn)生及發(fā)展動(dòng)態(tài) 1.2 灰色系統(tǒng)的研究?jī)?nèi)容1.3 灰色系統(tǒng)理論在建模中的應(yīng)用 1.1 灰色系統(tǒng)理論的產(chǎn)生及發(fā)展灰色系統(tǒng)理論的產(chǎn)生及發(fā)展動(dòng)態(tài)動(dòng)態(tài)定義定義1.11.1 系統(tǒng)是客觀世界普遍存在的一種物質(zhì)運(yùn)動(dòng)形式系統(tǒng)是客觀世界普遍存在的一種物質(zhì)運(yùn)動(dòng)形式, ,它和它和運(yùn)動(dòng)性一樣運(yùn)動(dòng)性一樣, ,是物質(zhì)存在的一種根本屬性是物質(zhì)存在的一種根本屬性. .定義定義1.21.2 灰色系統(tǒng)是指灰色系統(tǒng)是指“部分信息已知

14、部分信息已知, ,部分信息未知部分信息未知”的的“小樣本小樣本”,“,“貧信息貧信息”的不確定性系統(tǒng)的不確定性系統(tǒng), ,它通過對(duì)它通過對(duì)“部分部分”已已知信息的生成、開發(fā)去了解、認(rèn)識(shí)現(xiàn)實(shí)世界，實(shí)現(xiàn)對(duì)系統(tǒng)運(yùn)行知信息的生成、開發(fā)去了解、認(rèn)識(shí)現(xiàn)實(shí)世界，實(shí)現(xiàn)對(duì)系統(tǒng)運(yùn)行行為和演化規(guī)律的正確把握和描述行為和演化規(guī)律的正確把握和描述. .灰色系統(tǒng)模型的灰色系統(tǒng)模型的特點(diǎn)：特點(diǎn)：對(duì)試驗(yàn)觀測(cè)數(shù)據(jù)及其分布沒有特殊的要對(duì)試驗(yàn)觀測(cè)數(shù)據(jù)及其分布沒有特殊的要求和限制，是一種十分簡(jiǎn)便的新理論，具有十分寬廣的應(yīng)用領(lǐng)求和限制，是一種十分簡(jiǎn)便的新理論，具有十分寬廣的應(yīng)用領(lǐng)域。域?；疑到y(tǒng)理論，是在一般系統(tǒng)理論的基礎(chǔ)上產(chǎn)生的，灰色

15、系統(tǒng)理論，是在一般系統(tǒng)理論的基礎(chǔ)上產(chǎn)生的，它是系統(tǒng)科學(xué)思想發(fā)展的必然產(chǎn)物，是社會(huì)經(jīng)濟(jì)深它是系統(tǒng)科學(xué)思想發(fā)展的必然產(chǎn)物，是社會(huì)經(jīng)濟(jì)深入發(fā)展對(duì)科學(xué)刺激和需要的產(chǎn)物。當(dāng)我們認(rèn)識(shí)與研入發(fā)展對(duì)科學(xué)刺激和需要的產(chǎn)物。當(dāng)我們認(rèn)識(shí)與研究自然和社會(huì)時(shí)，要從系統(tǒng)的角度出發(fā)，從宏觀上究自然和社會(huì)時(shí)，要從系統(tǒng)的角度出發(fā)，從宏觀上對(duì)其進(jìn)行深入的剖析和整體把握。在實(shí)際中，我們對(duì)其進(jìn)行深入的剖析和整體把握。在實(shí)際中，我們首先要對(duì)事物進(jìn)行系統(tǒng)性認(rèn)識(shí)，進(jìn)而對(duì)已有的系統(tǒng)首先要對(duì)事物進(jìn)行系統(tǒng)性認(rèn)識(shí)，進(jìn)而對(duì)已有的系統(tǒng)進(jìn)行有效控制以及設(shè)計(jì)一些最優(yōu)系統(tǒng)來為人民服務(wù)。進(jìn)行有效控制以及設(shè)計(jì)一些最優(yōu)系統(tǒng)來為人民服務(wù)。對(duì)系統(tǒng)進(jìn)行控制就要通過系統(tǒng)

16、內(nèi)部和外部的信息和對(duì)系統(tǒng)進(jìn)行控制就要通過系統(tǒng)內(nèi)部和外部的信息和信息流來加以實(shí)施，通過對(duì)信息的控制進(jìn)而達(dá)到對(duì)信息流來加以實(shí)施，通過對(duì)信息的控制進(jìn)而達(dá)到對(duì)系統(tǒng)本身的控制。系統(tǒng)本身的控制。但是無論是現(xiàn)代控制理論還是經(jīng)典控制理論，但是無論是現(xiàn)代控制理論還是經(jīng)典控制理論，它們都要依賴正確而精確的數(shù)學(xué)模型，否則，它們都要依賴正確而精確的數(shù)學(xué)模型，否則，一切都很難取得滿意的結(jié)果。然而，在現(xiàn)實(shí)生一切都很難取得滿意的結(jié)果。然而，在現(xiàn)實(shí)生活中，有許多情況不大可能求得精確的數(shù)學(xué)?；钪?，有許多情況不大可能求得精確的數(shù)學(xué)模型，如工業(yè)系統(tǒng)、生物系統(tǒng)、經(jīng)濟(jì)系統(tǒng)、社會(huì)型，如工業(yè)系統(tǒng)、生物系統(tǒng)、經(jīng)濟(jì)系統(tǒng)、社會(huì)系統(tǒng)等。若得不出精

17、確的數(shù)學(xué)模型，現(xiàn)代控制系統(tǒng)等。若得不出精確的數(shù)學(xué)模型，現(xiàn)代控制理論的方法和手段就無法施行，因而，現(xiàn)代控理論的方法和手段就無法施行，因而，現(xiàn)代控制理論對(duì)一些研究對(duì)象也鞭長(zhǎng)莫及。制理論對(duì)一些研究對(duì)象也鞭長(zhǎng)莫及。 當(dāng)人們對(duì)這些問題進(jìn)行潛心研究時(shí)，當(dāng)人們對(duì)這些問題進(jìn)行潛心研究時(shí)，查查德德于于19651965年首創(chuàng)年首創(chuàng)模糊理論模糊理論，第一次用精確的數(shù)，第一次用精確的數(shù)學(xué)方式來分析和研究模糊量，取得了新的突破，學(xué)方式來分析和研究模糊量，取得了新的突破，隨后，模糊集合論迅速應(yīng)用于控制領(lǐng)域，隨后，模糊集合論迅速應(yīng)用于控制領(lǐng)域，收到了良好的效果。模糊控制能夠?qū)σ恍o法構(gòu)造收到了良好的效果。模糊控制能夠?qū)σ恍?/p>

18、無法構(gòu)造數(shù)學(xué)模型的系統(tǒng)進(jìn)行控制，但模糊控制也表現(xiàn)出固數(shù)學(xué)模型的系統(tǒng)進(jìn)行控制，但模糊控制也表現(xiàn)出固有的弱點(diǎn)，即信息利用率低，控制粗糙、精度低等有的弱點(diǎn)，即信息利用率低，控制粗糙、精度低等。因而，在要求高精度的情況下，這種控制難以勝。因而，在要求高精度的情況下，這種控制難以勝任，并且它也未能對(duì)被控對(duì)象的運(yùn)動(dòng)規(guī)律作深刻的任，并且它也未能對(duì)被控對(duì)象的運(yùn)動(dòng)規(guī)律作深刻的闡明，故模糊控制有它的局限性，只適應(yīng)于一些特闡明，故模糊控制有它的局限性，只適應(yīng)于一些特有的模糊系統(tǒng)有的模糊系統(tǒng)。 經(jīng)典控制理論、現(xiàn)代控制理論和模糊控制理論經(jīng)典控制理論、現(xiàn)代控制理論和模糊控制理論都有一個(gè)共同點(diǎn)，那就是它們所研究的對(duì)象系統(tǒng)必

19、都有一個(gè)共同點(diǎn)，那就是它們所研究的對(duì)象系統(tǒng)必須是須是白色系統(tǒng)白色系統(tǒng)（信息完全確知的系統(tǒng)），而事實(shí)上（信息完全確知的系統(tǒng)），而事實(shí)上，無論是自然系統(tǒng)還是社會(huì)系統(tǒng)，宏觀系統(tǒng)還是微，無論是自然系統(tǒng)還是社會(huì)系統(tǒng)，宏觀系統(tǒng)還是微觀系統(tǒng)，無生命系統(tǒng)還是有生命系統(tǒng)，對(duì)我們認(rèn)識(shí)觀系統(tǒng)，無生命系統(tǒng)還是有生命系統(tǒng)，對(duì)我們認(rèn)識(shí)的主體來說，總是信息的主體來說，總是信息不完全的，艱難說明一個(gè)系統(tǒng)的內(nèi)部參數(shù)是完不完全的，艱難說明一個(gè)系統(tǒng)的內(nèi)部參數(shù)是完全的。毫無疑問，內(nèi)部參數(shù)不完全的系統(tǒng)具有全的。毫無疑問，內(nèi)部參數(shù)不完全的系統(tǒng)具有極為普遍的意義。就像模糊理論的誕生一樣，極為普遍的意義。就像模糊理論的誕生一樣，灰色系統(tǒng)理論

20、也應(yīng)運(yùn)而生了?；疑到y(tǒng)理論也應(yīng)運(yùn)而生了。 灰色系統(tǒng)理論是我國(guó)學(xué)者灰色系統(tǒng)理論是我國(guó)學(xué)者鄧聚龍教授于鄧聚龍教授于1919世紀(jì)世紀(jì)8080年代年代初創(chuàng)立并發(fā)展的理論，它把一般系初創(chuàng)立并發(fā)展的理論，它把一般系統(tǒng)論，信息論和控制論的觀點(diǎn)和方法延伸到社統(tǒng)論，信息論和控制論的觀點(diǎn)和方法延伸到社會(huì)，經(jīng)濟(jì)，生態(tài)等抽象系統(tǒng)，結(jié)合運(yùn)用數(shù)學(xué)方會(huì)，經(jīng)濟(jì)，生態(tài)等抽象系統(tǒng)，結(jié)合運(yùn)用數(shù)學(xué)方法發(fā)展的一套解決灰色系統(tǒng)的理論和方法，法發(fā)展的一套解決灰色系統(tǒng)的理論和方法，2020多年來，灰色系統(tǒng)理論引起了國(guó)內(nèi)外學(xué)者的廣多年來，灰色系統(tǒng)理論引起了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注?；疑到y(tǒng)理論已成功應(yīng)用到工業(yè)，農(nóng)泛關(guān)注?；疑到y(tǒng)理論已成功應(yīng)用到工

21、業(yè)，農(nóng)業(yè)，社會(huì)，經(jīng)濟(jì)等眾多領(lǐng)域，解決了生產(chǎn)，生業(yè)，社會(huì)，經(jīng)濟(jì)等眾多領(lǐng)域，解決了生產(chǎn)，生活和科學(xué)研究中的大量實(shí)際問題?；詈涂茖W(xué)研究中的大量實(shí)際問題。1.2 灰色系統(tǒng)理論研究的內(nèi)灰色系統(tǒng)理論研究的內(nèi)容容 灰色系統(tǒng)理論經(jīng)過20年的發(fā)展，已基本建立起一門新興的結(jié)構(gòu)體系，其研究?jī)?nèi)容主要包括：灰色系統(tǒng)建模理論、灰色系統(tǒng)控制理論、灰色關(guān)聯(lián)分析方法、灰色預(yù)測(cè)方法、灰色規(guī)劃方法、灰色決策方法等。 今天我們主要介紹灰色系統(tǒng)建模理論及灰色數(shù)列預(yù)測(cè)?；疑珨?shù)列預(yù)測(cè)是指利用動(dòng)態(tài)GM模型，對(duì)系統(tǒng)的時(shí)間序列進(jìn)行數(shù)量大小的預(yù)測(cè)，即對(duì)系統(tǒng)的主行為特征量或某項(xiàng)指標(biāo)，發(fā)展變化到未來特定時(shí)刻出現(xiàn)的數(shù)值進(jìn)行預(yù)測(cè)。1.3 灰色系統(tǒng)理論在建

22、模中的灰色系統(tǒng)理論在建模中的應(yīng)用應(yīng)用 灰色系統(tǒng)理論在建模中被廣泛用來處理數(shù)據(jù)。與插值擬合相比，利用灰色模型處理數(shù)據(jù)不僅對(duì)數(shù)據(jù)沒有很強(qiáng)的限制，而且精度更高，計(jì)算更簡(jiǎn)便。2 2 灰色灰色GM(1.1)模型模型2.1 灰色生成 2.2 GM（1.1）模型建模機(jī)理2.3 GM(1.1)模型的精度檢驗(yàn) 2.1 灰色生成 將原始數(shù)據(jù)列中的數(shù)據(jù),按某種要求作數(shù)據(jù)處理稱為生成. 客觀世界盡管復(fù)雜,表述其行為的數(shù)據(jù)可能是雜亂無章的,然而它必然是有序的,都存在著某種內(nèi)在規(guī)律,不過這些規(guī)律被紛繁復(fù)雜的現(xiàn)象所掩蓋,人們很難直接從原始數(shù)據(jù)中找到某種內(nèi)在的規(guī)律.對(duì)原始數(shù)據(jù)的生成就是企圖從雜亂無章的現(xiàn)象中去發(fā)現(xiàn)內(nèi)在規(guī)律.

23、常用的灰色系統(tǒng)生成方式有: 累加生成,累減生成,均值生成,級(jí)比生成等,下面對(duì)這幾種生成做簡(jiǎn)單介紹.2.1.1 累加生成 累加生成：即通過數(shù)列間各時(shí)刻數(shù)據(jù)的依個(gè)累加以得到新的數(shù)據(jù)與數(shù)列。累加前的數(shù)列稱原始數(shù)列，累加后的數(shù)列稱為生成數(shù)列。累加生成是使灰色過程由灰變白的一種方法，它在灰色系統(tǒng)理論中占有極其重要地位，通過累加生成可以看出灰量積累過程的發(fā)展態(tài)勢(shì)，使離亂的原始數(shù)據(jù)中蘊(yùn)含的積分特性或規(guī)律加以顯化。累加生成是對(duì)原始數(shù)據(jù)列中各時(shí)刻的數(shù)據(jù)依次累加，從而生成新的序列的一種手段。(0)(0)(0)(0)(0)(1)(1)(1)(1)(1)(1)(0)(1),(2),( ),(1),(2),( ),:x

24、xxxxnxxxxxnxx 令令為為原原始始序序列列, ,記記生生成成數(shù)數(shù)為為如如果果與與之之間間滿滿足足如如下下關(guān)關(guān)系系(1)(0)1( )( );1,2,(21)kixkxikn ,1()AGO AccumulatingGenerationOperator 一一次次累累加加生生成成則則稱稱為為記記為為:r次次累累加加生生成成有有下下述述關(guān)關(guān)系系( )(1)1( )( )(22)krrixkxi (22),1:rr 從從式式 又又有有次次到到 次次的的累累加加為為1( )(1)(1)(1)(1)1( )( )( )(1)( )krrrrrixkxixkxkxk ( )(1)(2)111( )

25、( )( )kkirrriijxkxixj 累累加加生生成成在在灰灰色色系系統(tǒng)統(tǒng)理理論論中中有有著著非非常常重重要要的的地地位位, ,它它能能使使任任意意非非負(fù)負(fù)數(shù)數(shù)列列, ,擺擺動(dòng)動(dòng)的的或或非非擺擺動(dòng)動(dòng)的的, ,轉(zhuǎn)轉(zhuǎn)化化為為非非減減的的的的, ,遞遞增增的的數(shù)數(shù)列列. .2.1.2 累減生成 累減生成，即對(duì)數(shù)列求相鄰兩數(shù)據(jù)的差，累減生成是累加生成的逆運(yùn)算，常簡(jiǎn)記為IAGO (Inverse Accumulated Generating Operation), 累減生成可將累加生成還原為非生成數(shù)列，在建模過程中用來獲得增量信息,其運(yùn)算符號(hào)為. ( )( )( ),:rrixrxi 令令為為 次

26、次生生成成數(shù)數(shù)列列 對(duì)對(duì)作作 次次累累減減生生成成記記為為其其基基本本關(guān)關(guān)系系式式為為(0)( )( )(1)( )(0)( )(0)( )(2)( )(1)( )(1)( )( )( )(1)( )(1)( )( )( )( )( )(1)( )( )(1)(25)( )( )(1)rrrrrrrriririrxkxkxkxkxkxkxkxkxkxkxk (0)(1)( ),(0)0,;(0)1(0)11.ikkikki 式式中中為為 次次累累減減 即即無無累累減減為為1 1次次累累減減, ,即即與與 時(shí)時(shí)刻刻兩兩個(gè)個(gè)零零次次累累減減量量求求差差, ,為為 次次累累減減, ,即即與與 時(shí)時(shí)刻

27、刻兩兩個(gè)個(gè)次次累累減減量量求求差差(25): 從從式式還還可可得得到到以以下下關(guān)關(guān)系系(1)( )(0)( )(0)( )( )( )1(1)(1)11(1)( )( )(1)( )(1)(26)( )( )( )rrrrrkkrriirxkxkxkxkxkxixixk (2)( )(1)( )(1)( )(1)(1)1(2)(2)11(2)( )( )(1)( )(1)(27)( )( )( )rrrrrkkrriirxkxkxkxkxkxixixk :同同理理可可得得( )( )()( )( )(28)irr ixkxk ( )( )(0)( )( )(29)rrxkxk (29),.,.

28、:1,rrr 從從式式可可以以看看出出 對(duì)對(duì) 次次生生成成數(shù)數(shù)列列作作 次次累累減減即即還還原原為為非非生生成成數(shù)數(shù)列列事事實(shí)實(shí)上上 累累加加中中包包含含著著累累減減 累累減減中中包包含含著著累累加加比比如如時(shí)時(shí) 有有1(1)(0)(0)(0)11(1)(0)( )( )( )( )(1)( )(210)kkiixkxixixkxkxk (0)(1)(1)( )( )(1)xkxkxk進(jìn)進(jìn)一一步步有有(1)( )( )( )( )(1)(211)rrrxkxkxk .上上述述關(guān)關(guān)系系式式經(jīng)經(jīng)常常被被用用在在從從生生成成數(shù)數(shù)列列求求還還原原數(shù)數(shù)列列中中2.1.3 均值生成.均均值值生生成成分分為

29、為鄰鄰均均值值生生成成與與非非鄰鄰均均值值生生成成兩兩種種, (1), (2), ( ),( ),( )0.5 ( )0.5 (1),( )Xxxx nkz kz kx kx kz k 所所謂謂就就是是對(duì)對(duì)于于等等時(shí)時(shí)距距的的數(shù)數(shù)列列, ,用用相相鄰鄰數(shù)數(shù)據(jù)據(jù)的的平平均均值值構(gòu)構(gòu)造造新新的的數(shù)數(shù)據(jù)據(jù). .即即若若有有原原始始數(shù)數(shù)列列記記 點(diǎn)點(diǎn)的的生生成成值值為為且且則則稱稱為為鄰鄰均均值值生生成成數(shù)數(shù), ,顯顯然然, ,這這種種生生成成是是相相鄰鄰值值的的等等鄰鄰均均值值生生成成權(quán)權(quán)生生成成. ., (1), (2), ( ), (1), ( ),( ),( ),( )0.5 (1)0.5 (

30、1),( )Xxxkx kx nkkz kz kx kx kz k 所所謂謂就就是是對(duì)對(duì)于于非非等等時(shí)時(shí)距距的的數(shù)數(shù)列列, ,或或雖雖為為等等時(shí)時(shí)距距數(shù)數(shù)列列, ,但但剔剔除除異異常常值值之之后后出出現(xiàn)現(xiàn)空空穴穴的的數(shù)數(shù)列列, ,用用空空穴穴兩兩邊邊的的數(shù)數(shù)據(jù)據(jù)求求平平均均值值構(gòu)構(gòu)造造新新的的數(shù)數(shù)據(jù)據(jù)以以填填補(bǔ)補(bǔ)空空穴穴, ,即即若若有有原原始始數(shù)數(shù)據(jù)據(jù)這這里里為為空空穴穴 記記 點(diǎn)點(diǎn)的的生生成成值值為為且且則則稱稱為為非非鄰鄰均均值值生生成成數(shù)數(shù), ,顯顯然然, ,這這種種生生成成是是空空穴穴前前后后信信息息的的非非鄰鄰均均值值生生成成等等權(quán)權(quán)生生成成. .2.1.4 級(jí)比生成 級(jí)級(jí)比比生生

31、成成是是一一種種常常用用的的填填補(bǔ)補(bǔ)序序列列端端點(diǎn)點(diǎn)空空穴穴的的方方法法. .對(duì)對(duì)數(shù)數(shù)列列端端點(diǎn)點(diǎn)值值的的生生成成, ,我我們們無無法法采采用用均均值值生生成成填填補(bǔ)補(bǔ)空空缺缺, ,只只能能采采用用級(jí)級(jí)比比生生級(jí)級(jí)比比生生成成. .成成是是級(jí)級(jí)比比級(jí)級(jí)比比生生( (k k成成在在建建模模中中可可以以獲獲得得較較好好的的灰灰) )與與光光滑滑比比 ( (k k) )生生成成指指數(shù)數(shù)律律. .的的總總稱稱. .(0)(0)(0)(0)(1),(2),( ),(),( ),XxxxnKk 設(shè)設(shè)序序列列為為原原始始序序列列稱稱為為級(jí)級(jí)比比為為光光滑滑比比 其其表表達(dá)達(dá)式式為為(0)(0)(0)(1)(

32、 )( )/(1)( )( )/(1)(212)kxkxkkxkxk (0)(0)(0)(0)(0)(0)(0) (1),(2),(1), ( ),(1)(1),( )( ),(1)( )Xxxnnxnxnxxn設(shè)設(shè)為為端端點(diǎn)點(diǎn)是是空空穴穴的的序序列列 若若用用右右鄰鄰的的級(jí)級(jí)比比生生成成用用的的左左鄰鄰級(jí)級(jí)比比生生成成則則稱稱和和為為級(jí)級(jí)比比生生成成2.2 GM(1.1)模型建模機(jī)理,(1.1)GM灰灰色色系系統(tǒng)統(tǒng)是是對(duì)對(duì)離離散散序序列列建建立立的的微微分分方方程程是是一一階階微微分分方方程程模模型型, ,其其形形式式為為: :(2(1.1)13)dGMxaxudt :由由導(dǎo)導(dǎo)數(shù)數(shù)定定義義知

33、知0()( )limtdxx ttx tdtt 1,t 當(dāng)當(dāng)很很小小時(shí)時(shí)并并且且取取很很小小的的 單單位位時(shí)時(shí) 則則近近似似地地有有(1)( )xx tx tt 寫寫成成離離散散形形式式為為(1)(1)( )( (1)xx kx kx kt (1),(1)( ),(1)( )(1), ( ).(1), ( ):xxx kttx kx kx kx kx kx kxx kx kt 這這表表示示是是的的一一次次累累減減生生成成 因因此此是是和和二二元元組組合合等等效效值值 則則稱稱與與的的二二元元組組合合為為偶偶對(duì)對(duì), ,記記為為 于于是是我我們們可可以以定定義義一一個(gè)個(gè)從從 到到的的一一個(gè)個(gè)映映射

34、射: (1), ( )(214)dxFx kx kdt ( )(),( ).dxR ttxdtdxR tdtdxaxudt 若若定定義義是是 時(shí)時(shí)刻刻背背景景的的就就是是對(duì)對(duì)應(yīng)應(yīng)的的 的的值值那那么么 每每一一個(gè)個(gè)都都有有一一個(gè)個(gè)偶偶對(duì)對(duì)背背景景值值與與之之對(duì)對(duì)應(yīng)應(yīng)現(xiàn)現(xiàn)在在考考慮慮一一階階微微程程值值分分方方,1,( )(),dxx udtdxdtxx tx tt 它它是是與與的的線線性性組組合合. .那那么么, ,作作這這種種線線性性組組合合時(shí)時(shí), ,所所對(duì)對(duì)應(yīng)應(yīng)的的背背景景值值究究竟竟取取偶偶對(duì)對(duì)是是的的哪哪一一個(gè)個(gè)呢呢? ?如如果果認(rèn)認(rèn)為為在在的的很很短短時(shí)時(shí)間間內(nèi)內(nèi) 變變量量之之間間不

35、不會(huì)會(huì)出出現(xiàn)現(xiàn)突突變變量量 那那么么可可取取偶偶對(duì)對(duì)的的平平均均值值作作為為背背景景值值1( ) ( )(1)(215)2z tx kx k ,(1.1)GM基基于于上上述述機(jī)機(jī)理理 下下面面介介紹紹的的具具體體模模型型及及計(jì)計(jì)算算式式, ,設(shè)設(shè)非非負(fù)負(fù)原原始始序序列列 (0)(0)(0)(0)(1),(2),( )Xxxxn (0),X對(duì)對(duì)作作一一次次累累加加 得得到到生生成成數(shù)數(shù)列列為為 (1)(1)(1)(1)(1),(2),( )Xxxxn (1)0,( )( )kixkx i 其其中中(0)( )(1.1)xkGM于于是是的的白白化化形形式式的的微微分分方方程程為為(1)(1)(21

36、6)dxaxudt , , a u其其中中為為待待定定參參數(shù)數(shù), ,將將( (2 2- -1 16 6) )式式離離散散化化, ,即即得得(1)(1)(1)(1)( (1)(2 17)xkazx ku (1)(1)(1)(1)(1),(1)(1),(1)(1)xkxkdxzkkdt其其中中為為在在時(shí)時(shí)刻刻的的累累減減生生成成序序列列為為在在時(shí)時(shí)刻刻的的背背景景值值. .因因?yàn)闉?1)(1)(1)(1)(0)(1)(1)( )(1)(218)xkxkxkxk (1)(1)(1)1(1)(1)( )(219)2zkxkxk(218),(219) 將將式式代代入入( (2 2- -1 17 7) )

37、式式, ,得得(0)(1)(1)1(1)( )(1)(220)2xkaxkxku(220) 將將式式展展開開得得(1)(1)(0)(1)(1)(0)(0)(1)(1)1(1)(2)12(2)1(2)(3)1(3)(221)2( )1(1)( )12xxxxxxxnxnxn (1)(1)(0)(1)(1)(0)(0)(1)(1)1(1)(2)12(2)1(2)(3)1(3),2( )1(1)( )12xxxxxxYBxnxnxn 令令,(221)Tau 為為待待辨辨識(shí)識(shí)參參數(shù)數(shù)向向量量 則則可可寫寫成成(222)YB 參參數(shù)數(shù)向向量量可可用用最最小小二二乘乘法法求求取取, ,即即1 , ()(2

38、23)TTTa uB BB Y (216), 把把求求取取的的參參數(shù)數(shù)代代入入式式 并并求求出出其其離離散散解解為為(1)(1)(1)(1)(224)akuuxkxeaa 還還原原到到原原始始數(shù)數(shù)據(jù)據(jù)得得(0)(1)(1)(1)(1)(1)( )(1)(1)(225)aakxkxkxkuexea (224),(225)(1.1),(1.1).GMGM 式式稱稱為為模模型型的的時(shí)時(shí)間間相相應(yīng)應(yīng)函函數(shù)數(shù)模模型型 它它是是模模型型灰灰色色預(yù)預(yù)測(cè)測(cè)的的具具體體計(jì)計(jì)算算公公式式2.3 GM(1.1)模型的精度檢驗(yàn)?zāi)ＤＰ托瓦x選定定之之后后, ,一一定定要要經(jīng)經(jīng)過過檢檢驗(yàn)驗(yàn)才才能能判判定定其其是是否否合合理

39、理, ,只只有有通通過過檢檢驗(yàn)驗(yàn)的的模模型型才才能能用用來來作作預(yù)預(yù)測(cè)測(cè), ,灰灰色色模模型型的的精精度度檢檢驗(yàn)驗(yàn)一一般般有有三三種種方方法法: :相相對(duì)對(duì)誤誤差差大大小小檢檢驗(yàn)驗(yàn)法法, ,關(guān)關(guān)聯(lián)聯(lián)度度檢檢驗(yàn)驗(yàn)法法和和后后驗(yàn)驗(yàn)差差檢檢驗(yàn)驗(yàn)法法. .下下面面對(duì)對(duì)這這三三種種方方法法做做個(gè)個(gè)簡(jiǎn)簡(jiǎn)單單介介紹紹. .2.3.1 相對(duì)誤差檢驗(yàn)法(1)(1)(0)(1.1),GMXXX設(shè)設(shè)按按建建模模法法已已求求出出并并將將做做一一次次累累減減轉(zhuǎn)轉(zhuǎn)化化為為即即(0)(0)(0)(0)(1),(2),( )(231)Xxxxn 計(jì)計(jì)算算殘殘差差得得(0)(0) (1), (2), ( )(232)Eeee

40、nXX (0)(0), ( )( )( ),1,2,e kxkxk kn其其中中計(jì)計(jì)算算相相對(duì)對(duì)誤誤差差得得(0)( )( )100%,1,2,(233)( )e krel kknxk計(jì)計(jì)算算平平均均相相對(duì)對(duì)誤誤差差得得11( ),(234)nkrelrel kn 2.3.2 后驗(yàn)差檢驗(yàn)法(0)(0)2212(1.1)(231),(232),GMXXESS 設(shè)設(shè)按按建建模模法法所所求求出出的的如如所所示示 殘殘差差如如所所示示 原原始始序序列列及及殘殘差差序序列列 的的方方差差分分別別為為和和則則2(0)21122211( )1 ( )(235)nknkSxkxnSe ken (0)1111,

41、( ),( )nnkkxxk ee knn 其其中中計(jì)計(jì)算算后后驗(yàn)驗(yàn)差差比比為為21/(236)CSS 計(jì)計(jì)算算小小誤誤差差概概率率為為 1( )0.6745(237)pP e keS 1212,.CpCCSSSSC指指標(biāo)標(biāo) 和和 是是后后驗(yàn)驗(yàn)差差檢檢驗(yàn)驗(yàn)的的兩兩個(gè)個(gè)重重要要指指標(biāo)標(biāo). .指指標(biāo)標(biāo) 越越小小越越好好越越小小表表示示大大而而越越小小大大表表示示原原始始數(shù)數(shù)據(jù)據(jù)方方差差大大, ,即即原原始始數(shù)數(shù)據(jù)據(jù)離離散散程程度度大大. .小小表表示示殘殘方方差差小小, ,即即殘殘差差離離散散程程度度小小. . 小小就就表表明明盡盡管管原原始始數(shù)數(shù)據(jù)據(jù)很很離離散散, ,而而模模型型所所得得計(jì)計(jì)算算

42、值值與與實(shí)實(shí)際際值值之之差差并并不不太太離離散散. .1,0.6745,.ppC p指指標(biāo)標(biāo) 越越大大越越好好越越大大 表表明明殘殘差差與與殘殘差差平平均均值值之之差差小小于于給給定定值值的的點(diǎn)點(diǎn)較較多多 即即擬擬合合值值( (或或預(yù)預(yù)測(cè)測(cè)值值) )分分布布比比較較均均勻勻. .按按兩兩個(gè)個(gè)指指標(biāo)標(biāo) 可可綜綜合合評(píng)評(píng)定定預(yù)預(yù)測(cè)測(cè)模模型型的的精精度度 模模型型的的精精度度由由后后驗(yàn)驗(yàn)差差和和小小誤誤差差概概率率共共同同刻刻劃劃. .一一般般地地, ,將將模模型型的的精精度度分分為為四四級(jí)級(jí), ,見見表表2 2- -1 1模型精度等級(jí) 均方差比值C小誤差概率p1級(jí)（好）C=0.350.95=p2級(jí)

43、（合格）0.35C=0.50.80=p0.953級(jí)（勉強(qiáng)）0.5C=0.650.70=p0.804級(jí)（不合格） 0.65CP0.702 1 表表精精度度檢檢驗(yàn)驗(yàn)等等級(jí)級(jí)參參照照表表 ,Max pC 模模型型的的精精度度級(jí)級(jí)別別的的級(jí)級(jí)別別于于是是的的級(jí)級(jí)別別2.3.3 關(guān)聯(lián)度檢驗(yàn)法灰灰關(guān)關(guān)聯(lián)聯(lián)分分析析實(shí)實(shí)質(zhì)質(zhì)上上就就是是比比較較數(shù)數(shù)據(jù)據(jù)到到曲曲線線幾幾何何形形狀狀的的接接近近程程度度, ,一一般般來來說說, ,幾幾何何形形狀狀越越接接近近, ,變變化化趨趨勢(shì)勢(shì)也也就就越越接接近近, ,關(guān)關(guān)聯(lián)聯(lián)度度就就越越大大. .因因而而在在進(jìn)進(jìn)行行關(guān)關(guān)聯(lián)聯(lián)分分析析時(shí)時(shí), ,必必須須先先確確定定參參考考數(shù)數(shù)列

44、列, ,然然后后比比較較其其它它數(shù)數(shù)列列同同參參考考數(shù)數(shù)列列的的接接近近程程度度, ,這這樣樣才才能能對(duì)對(duì)其其它它數(shù)數(shù)列列進(jìn)進(jìn)行行比比較較, ,進(jìn)進(jìn)而而做做出出判判斷斷. . 000001(1),(2),( ),(1),(2),( ) ,1,2,:iiiiXXXXnXXXX nimXX 設(shè)設(shè)為為參參考考序序列列為為其其它它序序列列 則則與與的的關(guān)關(guān)聯(lián)聯(lián)系系數(shù)數(shù)為為0000min( )( )maxmax( )( )( )( )maxmax( )( )iijijijiiijXjXjXjXjXjXjXjXj ,1,2,jn 其其中中從從關(guān)關(guān)聯(lián)聯(lián)系系數(shù)數(shù)的的計(jì)計(jì)算算來來看看, ,我我們們得得到到比比較

45、較數(shù)數(shù)列列與與參參考考數(shù)數(shù)列列在在各各點(diǎn)點(diǎn)的的關(guān)關(guān)聯(lián)聯(lián)系系數(shù)數(shù)值值, ,結(jié)結(jié)果果較較多多, ,信信息息過過于于分分散散, ,不不便便于于比比較較, ,因因而而有有必必要要將將每每一一比比較較數(shù)數(shù)列列各各個(gè)個(gè)時(shí)時(shí)刻刻的的關(guān)關(guān)聯(lián)聯(lián)系系數(shù)數(shù)集集中中體體現(xiàn)現(xiàn)在在一一個(gè)個(gè)值值上上, ,這這一一數(shù)數(shù)值值就就是是灰灰關(guān)關(guān)鄧鄧氏氏灰灰聯(lián)聯(lián)度度. .色色關(guān)關(guān)聯(lián)聯(lián)度度為為11(239)niijjn (0)(0)(0)(0)0000(0)(0)(0)(0)(1),(2),( ),(1),(2),( ) ,1,2,.(1(,2,),.1,2,),iiiiXxxxnmXxxxnimmimmii i ii i設(shè)設(shè)原原始始

46、數(shù)數(shù)據(jù)據(jù)序序列列為為參參考考序序列列 用用 種種灰灰色色如如果果r r在在所所有有關(guān)關(guān)聯(lián)聯(lián)度度中中建建模模方方法法所所得得模模型型值值分分別別為為求求出出該該個(gè)個(gè)序序列列與與參參考考序序列列的的鄧鄧最最大大 則則第第 種種灰灰色色建建模模方方法法為為所所建建模模型型中中最最好好關(guān)關(guān)聯(lián)聯(lián)度度的的模模型型氏氏3 3 序列光滑度的理論分序列光滑度的理論分析析3.1 序列的光滑性 3.2 提高數(shù)據(jù)序列光滑度的方法3.3 基于函數(shù)sinx變換的改進(jìn)GM(1.1)模型 3.1 序列的光滑性* (1), (2), ( ),: (1), (2), ( ), ( )0.5( )0.5 (1),1,2, ;,(1

47、),:3.1Xxxx nZXZzzz nz kz kz kkn Xx nX 設(shè)設(shè)序序列列是是 的的均均值值生生成成序序列列其其中中是是某某一一可可導(dǎo)導(dǎo)函函數(shù)數(shù)的的代代表表序序列列 我我們們將將刪刪去去后后所所得得的的序序列列仍仍記記為為定定若若滿滿足足義義11*11(1), ( )( );(2)max( )( )max( )( )kik nk nkx kx ix kx kx kz k 當(dāng)當(dāng) 充充分分大大時(shí)時(shí),(1),(2).X則則稱稱 為為光光滑滑序序列列稱稱為為序序列列光光滑滑條條件件(0)(0)(0)(0)(0)(0)1(1)(0)1(0)(1),(2),( ),( )( )( ),2,3

48、,(1)(.)3 2.kiXxxxnxkxkkknxkxiX 定定義義設(shè)設(shè)序序列列稱稱為為序序的的光光滑滑比比列列,( )1( )(,( ).,.ikx kkx ikXXk k-1k-1=1=1光光滑滑比比從從另另一一個(gè)個(gè)側(cè)側(cè)面面反反映映了了序序列列的的光光滑滑性性 即即用用序序列列中中第第 個(gè)個(gè)數(shù)數(shù)據(jù)據(jù)與與前前個(gè)個(gè)數(shù)數(shù)據(jù)據(jù)的的和和的的比比值值來來考考察察序序列列 中中數(shù)數(shù)序序列列中中的的數(shù)數(shù)據(jù)據(jù)變變化化越越平平穩(wěn)穩(wěn) 其其光光滑滑比比據(jù)據(jù)變變化化是是否否平平緩緩 顯顯越越小小然然 (0)1(0)1(0)(0)(1( )0)00( )( )( )( ),( ),1,2,( ),1,2,0,.1)

49、,.,3 3kixkkxixkxk knxk knkkxkk 設(shè)設(shè)為為非非負(fù)負(fù)數(shù)數(shù)據(jù)據(jù)序序列列當(dāng)當(dāng)時(shí)時(shí) 如如果果光光滑滑定定義義光光滑滑為為比比稱稱離離散散序序列列則則(1)( ):(1)1;( )2,3,1;(2) ( )0, ;2,3, ;(33.4)0.5.kX kkknkknX 若若序序列列滿滿足足則則定定義義準(zhǔn)準(zhǔn)光光稱稱滑滑序序列列為為 (0)1(0)1(0)()001()( )( ).( )( ),( ),1,2,(1( ),1,2.1.,),3kixk kxkkkxixkxk knnxk 設(shè)設(shè)為為光光滑滑離離散散序序列列的的充充是是 的的遞遞減減函函數(shù)數(shù)則則稱稱要要條條件件是是定

50、定理理光光滑滑離離為為: :光光滑滑散散序序列列比比3.2 提高數(shù)據(jù)序列光滑度的方法3.2.1 基于函數(shù)lnx變換提高數(shù)據(jù)序列的光滑度(0)(0)(0)(0)(0)(0)(0)11(0)(0)11(1)( )(1)1,ln( )(2)( )(1),ln( )( )ln( )( )kkiixkxxkxkxexkxkxixi 若若為為遞遞增增序序列列, ,且且則則是是光光滑滑離離散散序序列列; ;若若為為遞遞增增序序列列, ,且且則則002( )( )( ),ln( )( )xkxk由由知知 變變換換序序列列具具有有比比原原始始序序列列更更好好的的光光滑滑度度(0)(0)1(0)(0)111(0)

51、(0)11(3)( )(1)1,1,ln( )ln( )ln( )ln( )TkkTiixkxTxkxkxixi 若若為為遞遞增增序序列列, ,且且則則1003( )( )( ),ln( )ln( )Txkxk由由知知 變變換換序序列列 具具有有比比原原始始序序列列更更好好的的光光滑滑度度(0)(0)1(0)(0)(0)1111(0)(0)(0)11111(0)(0)(0)11111(0)(0)(0)111(4)( )(1),1,ln( )ln( )( )ln( )ln( )( )ln( )( )( )ln( )( )( )TkkkTiiiTTkkkTTiiixkxe Txkxkxkxixix

52、ixkxkxkxixixi若若為為遞遞增增序序列列, ,且且則則(0)(0)(0)(0)(0)1(0)(0)(0)(0)1111(0)(0)(0)(0)1111(1)( )(1)1,( ) (0)(2)( )(1)1,1,0( )ln( )ln( )( )( )ln( )ln( )( )aaTkkkkaTiiiixkxxkaxkxTaxkxkxkxkxixixixi 若若為為遞遞增增序序列列, ,且且則則 是是光光滑滑離離散散序序列列; ;若若為為遞遞增增序序列列, ,且且則則對(duì)對(duì)于于有有1 1002( )( )( ),( )ln( )aTxkxk由由知知 變變換換序序列列 具具有有比比原原始

53、始序序列列 更更好好的的光光滑滑度度3 2 30. .()axa基基于于函函數(shù)數(shù)變變換換提提高高光光滑滑度度的的方方法法2( ),:在在應(yīng)應(yīng)用用時(shí)時(shí) 應(yīng)應(yīng)注注意意以以下下兩兩點(diǎn)點(diǎn)000021( )( )( )( )( ),),rrxxxx1 1. .若若的的條條件件( (1 1) ) e e或或( (1 1) )不不成成立立 可可選選擇擇某某個(gè)個(gè)常常數(shù)數(shù)c c, , 使使( (1 1) )( (1 1) )+ +c ce e( (或或最最后后再再還還原原. .02( )( ),( )xka2 2. .對(duì)對(duì)的的a a可可根根據(jù)據(jù)實(shí)實(shí)際際情情況況定定 一一般般情情況況下下越越大大, ,取取得得越越

54、小小. .3 2 3. .()xaa基基于于函函數(shù)數(shù)變變換換提提高高光光滑滑度度的的方方法法(0)(0)(0)(0)(0)( )(0)(0)1( )(0)(0)(0)11111( )(0)(0)(0)1111(1)( )(1)1,(1)(2)( )(1)1,1,1ln( )ln( )( )ln( )ln( )( )xkxkTkkkkxiTiiiixkxaaxkxTaaxkxkxkaxixixi 若若為為遞遞增增序序列列, ,且且則則是是光光滑滑離離散散序序列列; ;若若為為遞遞增增序序列列, ,且且則則對(duì)對(duì)于于有有01021( )( )( )( ),ln( )xkTaaxk由由知知 變變換換序

55、序列列( () )具具有有比比序序列列 更更好好的的光光滑滑度度(0)(0)(0)(0)11( )(0)(0)11111( )(0)(0)111(3)( )(1)1,1,1ln( )ln( )ln( )ln( )xkTTkkkxkTTiiixkxTaaxkxkaxixi 若若為為遞遞增增序序列列, ,且且則則對(duì)對(duì)于于有有011031( )( )( )( ),()ln( )xkTTaaxk由由知知 變變換換序序列列具具有有比比原原始始序序列列 更更好好的的光光滑滑度度3.3 基于函數(shù)sinx變換的改進(jìn)GM(1,1)模型(0)(0)(0)( ),( ),2sin( ),( 1,2, )3.2.xk

56、kxkxkn 若若為為遞遞增增序序列列, ,且且則則是是定定理理光光滑滑離離散散序序列列1(0)(0)(0)1111(0(0)(0)(0)11(0)1)sin( )ln( )ln( )sin( )( ),( ),21ln( )ln).,33(:TkkkTiiixkxkxkxixixixkixiT 若若為為遞遞增增序序列列, ,且且則則有有 式式成成立立定定理理正正111(0)(0)(0)(0)1111111(0)(0(0)(0)(0)(0)1111)sin( )ln( )( )( )sin( )ln( )( )( ),( ),2.413,:)TTTkkkkTTTiiiixkxkxkxkxkxi

57、xixixiixiT 若若為為遞遞增增序序列列, ,且且則則有有正正式式成成立立定定理理1001( )( ), sin( )( ),TaxkxkT由由以以上上定定理理知知( () )能能夠夠2 2有有效效地地提提高高原原始始序序列列的的光光滑滑度度. .基于sinx變換的GM(1,1)的建模過程:設(shè)設(shè)非非負(fù)負(fù)原原始始數(shù)數(shù)據(jù)據(jù)序序列列為為遞遞增增序序列列 (0)(0)(0)(0)(1),(2),( )Yyyyn (0)( )1,1,2,yiin 其其中中 (0)(0)(0)(0)(0)(1),(2),( ),( ), .2rYyyyniyi 將將序序列列進(jìn)進(jìn)行行標(biāo)標(biāo)準(zhǔn)準(zhǔn)化化處處理理 使使得得對(duì)對(duì)

58、任任意意 標(biāo)標(biāo)準(zhǔn)準(zhǔn)化化后后的的均均在在區(qū)區(qū)間間內(nèi)內(nèi)(0)sin( ):rxyk基基于于函函數(shù)數(shù)變變換換對(duì)對(duì)數(shù)數(shù)據(jù)據(jù)序序列列進(jìn)進(jìn)行行處處理理 (0)(0)(0)(0)(1),(2),( )Xxxxn (0)(0)( )sin( ),1,2, .rxkykkn 其其中中, ,(1):AGO 對(duì)對(duì)處處理理后后的的數(shù)數(shù)據(jù)據(jù)序序列列作作一一次次累累加加得得 (1)(1)(1)(1)(1),(2),( )Xxxxn (1)(0)1( )( ),1,2, .kixkxi kn 其其中中, ,(1)( )(1.1),XkGM由由一一階階生生成成模模塊塊建建立立模模型型對(duì)對(duì)應(yīng)應(yīng)的的白白化化微微分分方方程程為為:

59、 :(1)(1)( )( )dxtaxtudt (1.1)GM模模型型對(duì)對(duì)應(yīng)應(yīng)的的差差分分方方程程為為: :(0(1)( )( ),2,3,xkaztu kn 展展開開得得(0)(0)(0)(0)(0)(0)(2)(2)1(3)(3)1( )( )1xzaxzuxnzn (0)(0)(0)(2),(3),( ), , ,.Yxxxna ua u 令令稱稱為為待待辨辨識(shí)識(shí)參參數(shù)數(shù)向向量量, ,其其中中為為待待辨辨識(shí)識(shí)常常數(shù)數(shù)(0)(0)(0)(2)1(3)1( )1zzBzn (1),( ):zk其其中中為為背背景景值值(1)(1)(1)1(1)(1)( ),1,2,1.2zkxkxkkn 待待

60、辨辨識(shí)識(shí)向向量量的的最最小小二二乘乘解解為為: :1()TTB BB Y 白白化化微微分分方方程程的的離離散散解解為為: :(1)(1)(1)(1)akuuxkxeaa 還還原原到到原原始始數(shù)數(shù)據(jù)據(jù)為為: :(0)(1)(1)(1)(1)(1)( )(1)(1)aakxkxkxkuexea (0)(0)( )arcsin( )rykxk (0)(0),( )( ).rykyk最最后后 將將標(biāo)標(biāo)準(zhǔn)準(zhǔn)化化后后的的還還原原為為4 4 灰色灰色GM(1.1)GM(1.1)優(yōu)化模型優(yōu)化模型分析分析4.1 傳統(tǒng)GM(1.1)模型背景值對(duì)預(yù)測(cè)精度的影響 4.2 基于積分背景值構(gòu)造法的改進(jìn)GM(1.1)模型對(duì)