【正文】 相對(duì)殘差有所減小，因此精度更高。下面再對(duì)該模型進(jìn)行后驗(yàn)差檢驗(yàn)。首先計(jì)算原始數(shù)列的均值 (651)再計(jì)算原始數(shù)列均方差 (652)求殘差均值 (653)求殘差均方差 (654)計(jì)算均方差比較值 (655)計(jì)算小殘差概率 (656) (657)由以上數(shù)據(jù)可得。由表51可以看出，對(duì)于給定的，且，所以改進(jìn)模型可以被評(píng)為“優(yōu)”。但精度要高于GM(1,1)模型，因此用該模型預(yù)測(cè)，結(jié)果將更加接近于真實(shí)值。下面同樣預(yù)測(cè)未來(lái)五年內(nèi)的負(fù)荷，如下:表68 未來(lái)五年電力負(fù)荷預(yù)測(cè)量 （單位：億度）年份20122013201420152016預(yù)測(cè)量 GM(1,N)模型及應(yīng)用實(shí)例GM(1,N)即表示模型是1階的，包含有N個(gè)變量的灰色模型。設(shè)系統(tǒng)有N個(gè)行為因子，即原始數(shù)列為 (658)記為的AGO數(shù)列，即 (659)其中。求的均值數(shù)列 (660)則。于是可得到GM(1,N)的灰微分方程為 (661)其中為灰導(dǎo)數(shù)，為背景值，為參數(shù)。如果對(duì)于一切時(shí)刻，引入向量矩陣記號(hào)：， (662)則GM(1,N)的灰微分方程為 (663)其中為已知數(shù)據(jù)向量，為GM(1,N)已知數(shù)據(jù)矩陣，為參數(shù)向量。用表示的估計(jì)值，如果存在，則有類似GM(1,1)模型，可得預(yù)測(cè)值為 (664)其中取為。累減還原成原始數(shù)列得 (665)使用該模型時(shí)，需要考慮到多種不確定因素，例如天氣、節(jié)假日、大工業(yè)用戶突發(fā)事件、國(guó)家政策、自然災(zāi)害等，都會(huì)對(duì)負(fù)荷預(yù)測(cè)產(chǎn)生一定的影響。由于這些不確定 因素難以數(shù)字化，并且數(shù)據(jù)不易搜集，因此該模型并不常用，充分利用好上述兩種模型足以。7 總結(jié)本文對(duì)灰色預(yù)測(cè)模型進(jìn)行了推導(dǎo)研究，對(duì)GM(1,1)模型及改進(jìn)模型進(jìn)行了殘差、后驗(yàn)差檢驗(yàn)，并推廣了GM(1,1)模型到GM(1,N)模型。作為一種實(shí)用性很高的理論，灰色系統(tǒng)理論具有寬廣的應(yīng)用領(lǐng)域，在電力負(fù)荷預(yù)測(cè)中也得到日益廣泛的應(yīng)用。當(dāng)電力負(fù)荷呈指數(shù)規(guī)律持續(xù)增長(zhǎng)時(shí)，用普通GM(1,1)模型即可進(jìn)行預(yù)測(cè)，該預(yù)測(cè)模型具有精度高、所需樣本數(shù)量少、計(jì)算簡(jiǎn)單、可進(jìn)行多種檢驗(yàn)等優(yōu)點(diǎn)。然而，實(shí)際電網(wǎng)存在數(shù)據(jù)波動(dòng)大、數(shù)據(jù)突變等不確定情況，此時(shí)再用GM(1,1)模型，得到的預(yù)測(cè)值可能會(huì)使得誤差變大，不能準(zhǔn)確代表負(fù)荷消耗趨勢(shì)。因此這個(gè)時(shí)候可應(yīng)用GM(1,1)改進(jìn)模型，可獲得較理想的效果。另外Matlab語(yǔ)言具有良好的運(yùn)行環(huán)境、強(qiáng)大的函數(shù)資源，其編程效率遠(yuǎn)遠(yuǎn)高于其他高級(jí)語(yǔ)言。多變量灰色預(yù)測(cè)模型廣泛的應(yīng)用于許多領(lǐng)域。但該模型參數(shù)估計(jì)以及預(yù)測(cè)都需要經(jīng)過(guò)比較復(fù)雜的計(jì)算，本文灰色預(yù)測(cè)模型通過(guò)Matlab程序能夠方便的解決模型的計(jì)算問(wèn)題。致謝首先我要向張濤老師表示衷心的感謝。在本人做畢業(yè)設(shè)計(jì)期間，一直得到張老師的悉心指導(dǎo)，從論文的選題、研究、審查直至完成，都是在張老師的精心指導(dǎo)下進(jìn)行的，在此謹(jǐn)向張老師表示崇高的敬意和最衷心的感謝。老師淵博的知識(shí)、豐富的實(shí)踐經(jīng)驗(yàn)和嚴(yán)謹(jǐn)?shù)目茖W(xué)作風(fēng)讓我十分敬佩，老師不辭辛勞，對(duì)每一個(gè)細(xì)節(jié)的詳細(xì)講解、指導(dǎo)和答疑解惑讓我深受啟發(fā)。張老師總是在百忙之中抽出寶貴的時(shí)間給我以指導(dǎo)和教誨，其學(xué)者風(fēng)范和嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度令人欽佩。值此學(xué)位論文完成之際，衷心地感謝老師給予的悉心關(guān)懷和諄諄教誨！我要感謝電氣學(xué)院的各位領(lǐng)導(dǎo)、老師幾年來(lái)對(duì)我學(xué)習(xí)上的教誨和支持，感謝他們?cè)趯W(xué)習(xí)和工作上給予我的支持和幫助！此外，在論文撰寫(xiě)過(guò)程中，我得到了許多朋友的大力支持和無(wú)私指導(dǎo)，共同解決了很多問(wèn)題，讓我感到十分幸運(yùn)，在此向許多同學(xué)和朋友，表示衷心的感謝。他們勤勞工作的態(tài)度和樸素善良的本質(zhì)一直是我學(xué)習(xí)的榜樣，祝愿他們永遠(yuǎn)健康、永遠(yuǎn)快樂(lè)！我還要感謝我所列參考文獻(xiàn)的作者，正是他們的許多研究成果給了我很大的幫助，在此表示誠(chéng)摯的謝意。衷心感謝各位專家在百忙之中對(duì)本文的審閱和指導(dǎo)！ 參考文獻(xiàn)[1] （第二版）. 高等教育出版社,2009[2] 司守奎,2012[3] 劉保柱,蘇彥華, . 人民郵電出版社,2010[4] :[碩士學(xué)位論文].上海:上海交通大學(xué),2006[5] :[碩士學(xué)位論文].湖南:湖南大學(xué),2009[6] 孫輝,姜梅,1997年第2期[7] 吳熳紅,2004年第1期[8] ,2012年第14期[9] ,1990年第1期[10] :[本科生學(xué)年論文].杭州師范大學(xué)錢(qián)江學(xué)院,2010[11] ,2002[12] 孟祥斌. 基于小波分析的中長(zhǎng)期負(fù)荷預(yù)測(cè)模型研究與實(shí)現(xiàn):[碩士學(xué)位論文].大連:大連理工大學(xué),2008[13] 陳軍,姜國(guó)麟,2000年第1期[14] :中國(guó)統(tǒng)計(jì)出版社 附錄附件1:GM(1,1)模型程序 exp(1/7)。求級(jí)比ans = 。最大值 exp(1/7) ans =。最小值 x0=[ ]。原始數(shù)列 x1=cumsum(x0)x1 =[ ]。一次累加值 b=length(x1)b = 13。長(zhǎng)度 for i=1:b1ave(i)=1/2*(x1(i)+x1(i+1))。end ave z=ave39。z =[ ]。平均值 a=ones(12,1)。 B=[z,a]。 Y=[ ]。構(gòu)造特征矩陣和特征向量 c=B39。 s=c*B。 d=inv(s)。 f=d*c。 g=f*Y。g = 。求u a= 。 求預(yù)測(cè)表達(dá)式 b= 。 b/aans = m= mansans = for i=1:14ave(i)=*exp(*i)。end。模型表達(dá)式 x0=[ ]。后驗(yàn)差檢驗(yàn) x1=[ ]。殘差數(shù)列 x2=mean(x0)x2 =。平均值 x3=mean(x1)x3 = s=sum((x0x2).^2)/12s = t=sum((x1x3).^2)/12t = s1=sqrt(s)。求均方差s1 = s2=sqrt(t)s2 = c=s2/s1c = s0=*s1s0 = p=abs(x1x3)p = 附件2:殘差模型程序e0=[ 0]。部分殘差 e1=cumsum(e0)e1 =[ ] b=length(e1)。 for i=1:b1ave(i)=1/2*(e1(i)+e1(i+1))。end ave。 z=ave39。z =[ ] a=ones(7,1)。 B=[z,a]。 Y=[ ]。 c=B39。 s=c*B。 d=inv(s)。 f=d*c。 g=f*Yg = x0=[ ]。 殘差模型后驗(yàn)差檢驗(yàn) x1=[0 0 ]。 x2=mean(x0)x2 = x3=mean(x1)x3 = s=sum((x0x2).^2)/12s = t=sum((x1x3).^2)/12t = s1=sqrt(s)s1 = s2=sqrt(t)s2 = c=s2./s1c = s0=*s1s0 = p=abs(x1x3)p = for i=1:12a(i)=*exp(*(i+5))+*exp(*i)。殘差模型總表達(dá)式enda = 畢業(yè)論文（設(shè)計(jì)）原創(chuàng)性聲明本人所呈交的畢業(yè)論文（設(shè)計(jì)）是我在導(dǎo)師的指導(dǎo)下進(jìn)行的研究工作及取得的研究成果。據(jù)我所知，除文中已經(jīng)注明引用的內(nèi)容外，本論文（設(shè)計(jì)）不包含其他個(gè)人已經(jīng)發(fā)表或撰寫(xiě)過(guò)的研究成果。對(duì)本論文（設(shè)計(jì)）的研究做出重要貢獻(xiàn)的個(gè)人和集體，均已在文中作了明確說(shuō)明并表示謝意。 作者簽名： 日期： 畢業(yè)論文（設(shè)計(jì)）授權(quán)使用說(shuō)明本論文（設(shè)計(jì)）作者完全了解**學(xué)院有關(guān)保留、使用畢業(yè)論文（設(shè)計(jì)）的規(guī)定，學(xué)校有權(quán)保留論文（設(shè)計(jì)）并向相關(guān)部門(mén)送交論文（設(shè)計(jì)）的電子版和紙質(zhì)版。有權(quán)將論文（設(shè)計(jì)）用于非贏利目的的少量復(fù)制并允許論文（設(shè)計(jì)）進(jìn)入學(xué)校圖書(shū)館被查閱。學(xué)?？梢怨颊撐模ㄔO(shè)計(jì)）的全部或部分內(nèi)容。保密的論文（設(shè)計(jì)）在解密后適用本規(guī)定。 作者簽名： 指導(dǎo)教師簽名： 日期： 日期： 注 意 事 項(xiàng)（論文）的內(nèi)容包括：1）封面（按教務(wù)處制定的標(biāo)準(zhǔn)封面格式制作）2）原創(chuàng)性聲明3）中文摘要（300字左右）、關(guān)鍵詞4）外文摘要、關(guān)鍵詞 5）目次頁(yè)（附件不統(tǒng)一編入）6）論文主體部分：引言（或緒論）、正文、結(jié)論7）參考文獻(xiàn)8）致謝9）附錄（對(duì)論文支持必要時(shí)）：理工類設(shè)計(jì)（論文）正文字?jǐn)?shù)不少于1萬(wàn)字（不包括圖紙、程序清單等）。：任務(wù)書(shū)、開(kāi)題報(bào)告、外文譯文、譯文原文（復(fù)印件）。、圖表要求：1）文字通順，語(yǔ)言流暢，書(shū)寫(xiě)字跡工整，打印字體及大小符合要求，無(wú)錯(cuò)別字，不準(zhǔn)請(qǐng)他人代寫(xiě)2）工程設(shè)計(jì)類題目的圖紙，要求部分用尺規(guī)繪制，部分用計(jì)算機(jī)繪制，所有圖紙應(yīng)符合國(guó)家技術(shù)標(biāo)準(zhǔn)規(guī)范。圖表整潔，布局合理，文字注釋必須使用工程字書(shū)寫(xiě)，不準(zhǔn)用徒手畫(huà)3）畢業(yè)論文須用A4單面打印，論文50頁(yè)以上的雙面打印4）圖表應(yīng)繪制于無(wú)格子的頁(yè)面上5）軟件工程類課題應(yīng)有程序清單，并提供電子文檔1）設(shè)計(jì)（論文）2）附件：按照任務(wù)書(shū)、開(kāi)題報(bào)告、外文譯文、譯文原文（復(fù)印件）次序裝訂3）其它