【正文】 進一步優(yōu)化處理。此外筆者認為馬爾科夫鏈模型預(yù)測函數(shù)中取預(yù)測值區(qū)間中點作為最優(yōu)預(yù)測值的做法欠妥，可以考慮對灰色—馬爾科夫鏈模型進一步優(yōu)化處理。這些都有待作進一步深入研究。參考文獻[1]楊榮英，張輝，[J].武漢理 工大學(xué)學(xué)報(交通科學(xué)與工程版)，2001，25(3): 353355．[2]湯中明，[J]. 武漢理工大學(xué)學(xué)報(交通科學(xué)與工程版)，2003，27(5): 728730．[3][D].武漢理工大學(xué)碩士研究生學(xué)位論文， 2005(6).3436．[4][J] .商場現(xiàn)代化，2006(12). 5657．[5]張鳳榮，[J].工業(yè)技術(shù)經(jīng)濟，2003(6): 7273．[6]王曉原，(1,1)模型在區(qū)域物流規(guī)模預(yù)測中的應(yīng)用[J].武漢理工大學(xué)學(xué)報(交通科學(xué)與工程版)，2005,29(3): 415417．[7]陳森，[J].統(tǒng)計與決策，2006(2): 5960．[8]張冬柏，[J]. 黑龍江交通科技，2006(9): 116117．[9](1,1)模型在公路貨運量預(yù)測中的應(yīng)用[J].(6): 2123．[10]蓋春英，—馬爾可夫鏈預(yù)測方法研究[J].中國公路學(xué)報，2003,16(3): 113116．[11]李莉，張建華．物流產(chǎn)業(yè)發(fā)展與國民經(jīng)濟整體水平提升的相關(guān)性分析[J]．中國機械工程．2003(10)：2428．[12]Acevedo J．a(chǎn)mp。Pistikopouios，E．N．Stochastic optimization based algorithms for Process synthesis under Uncertainty[J]． Computers and Chemical Engineering，1998，(2)：6470．[13]Brown，J H．An economic model of airline hub and spoke[J]．Logistics and Transportation Review Volume，1999，(24)：4065．[14]Timo Tera svirta, Diek van Dijk, Marcelo C. Medeiros. Linear models, smooth transition autoregessions, and neural networks for forecasting macroeconomic time series: A reexaurination[J].International Journal of Forecasting, 2005, 21: 75574．[15]孫艷玲．灰色預(yù)測模型在物流需求預(yù)測中的應(yīng)用[J]．物流平臺，2007，(1)：145146．[16]Yi Lin, Sifeng Liu. A Historical Introduction to Grey Systems Theory [C]. IEEE International Conference on Systems Man and Cybemetics,2004:24032408.[17]楊燦，楊敏華，[J].測繪科學(xué),2009, 34(6): 182183.[18]陳永紅，[J].. 30(2): 9598.[19]劉耀林， 劉艷芳， 張玉梅. 基于灰色 馬爾柯夫鏈預(yù)測模型的耕地需求量預(yù)測研究[J] . 武漢大學(xué)學(xué)報: 信息科學(xué)版, 2004(29)。576.致謝從論文選題到搜集資料，從寫稿到反復(fù)修改，在這期間經(jīng)歷了喜悅、痛苦以及彷徨。現(xiàn)在伴隨著本文的最終定稿，我復(fù)雜的心情已經(jīng)煙消云散了，自己甚至有那么點成就感。再就是我想在此抒發(fā)我發(fā)自肺腑的謝意。首先我要感謝我的導(dǎo)師柳虹，她為人隨和熱情，治學(xué)嚴謹細心。在閑聊中她總是像朋友一樣給予鼓勵，而在論文的寫作和措辭等方面她會以專業(yè)標準嚴格要求我。萬事開頭難，當我在選題時陷入迷茫時，她在詢問了我比較感興趣的研究方向后幫助我選題、定題。到最后論文的修改、潤色時，柳老師更是認真地給予我深刻并且細致的指導(dǎo)和熱忱的鼓勵?？傊怯捎诹蠋熓俏覠o私幫助和鼓勵，我的畢業(yè)論文得以順利完成，在此謝謝劉老師。再者我要感謝周圍的同學(xué)們，在我寫論文時，你們給了我許多寶貴的意見和建議。在學(xué)習(xí)生活中也受到你們不少照顧，在此我要向你們表達我真摯的感謝。最后我還要感謝一直在背后默默支持我的家人，沒有你們的付出就沒有我今天的成績。畢業(yè)論文完成后意味著我本科即將畢業(yè)，但是我清楚地明白這僅僅意味著我要帶著老師、同學(xué)、家人的鼓勵踏上新的征程。 附錄%對寧波港總吞吐量預(yù)測的MATLAB代碼clear allclc %清屏z=xlsread(39。C:\Documents and Settings\Administrator\桌面\39。,39。sheet139。,39。A1:H139。)%調(diào)取數(shù)據(jù)format long gxp=mean(z)n=length(z)x(1)=z(1)。for i=2:nx(i)=z(i)+x(i1)。endfor i=1:n1b(i)=(x(i)+x(i+1))/2。 %構(gòu)造矩陣Bendb=[b(1),1。b(2),1。b(3),1。b(4),1。b(5),1。b(6),1。b(7),1]。y=[z(2),z(3),z(4),z(5),z(6),z(7),z(8)]39。 %構(gòu)造向量Yt=inv(b39。*b)*(b39。*y)。a=t(1) %得到發(fā)展系數(shù)u=t(2) %得到調(diào)整系數(shù)for k=2:n+5x(k) = (x(1)u/a)*exp(a*(k1))+u/a。 %時間響應(yīng)函數(shù)endfor i=n+5:1:2x(i)=x(i)x(i1)。 %累加預(yù)測值累減后得到原始數(shù)據(jù)預(yù)測值endp=[0,0,0,1。1/2,0,1/2,0。1/2,0,0,1/2。0,1/2,1/2,0]。 %一步轉(zhuǎn)移概率矩陣vs=[1,0,0,0]。 %開始狀態(tài)向量for i=1:5v=vs*p^i %各預(yù)測年份的狀態(tài)向量for s=4:1:1if v(s) = max(v)e(i) = s。endendende=[2,1,4,2,3,4,3,1,e] %得到各年份所處的狀態(tài)區(qū)間E=[*xp,*xp,*xp] %劃分狀態(tài)區(qū)間的各點for i=1:n+5if e(i)=1 y(i)=x(i)+E(1)。else if e(i)=2 y(i)=x(i)+(E(1)+E(2))/2。else if e(i)=3 y(i)=x(i)+(E(2)+E(3))/2。else y(i)=x(i)+E(3)。endendx %灰色GM(1,1)模型預(yù)測結(jié)果y %灰色馬爾科夫鏈模型預(yù)測 %對寧波港集裝箱吞吐量預(yù)測的MATLAB代碼clear allclc %清屏z=xlsread(39。C:\Documents and Settings\Administrator\桌面\39。,39。sheet139。,39。A1:H139。)%調(diào)取數(shù)據(jù)format long gxp=mean(z)n=length(z)x(1)=z(1)。for i=2:nx(i)=z(i)+x(i1)。endfor i=1:n1b(i)=(x(i)+x(i+1))/2。 %構(gòu)造矩陣Bendb=[b(1),1。b(2),1。b(3),1。b(4),1。b(5),1。b(6),1。b(7),1]。y=[z(2),z(3),z(4),z(5),z(6),z(7),z(8)]39。 %構(gòu)造向量Yt= inv (b39。*b)*(b39。*y)。a=t(1) %得到發(fā)展系數(shù)u=t(2) %得到調(diào)整系數(shù)for k=2:n+5x(k)=(x(1)u/a)*exp(a*(k1))+u/a。 %時間響應(yīng)函數(shù)endfor i=n+5:1:2x(i)=x(i)x(i1)。 %累加預(yù)測值累減后得到原始數(shù)據(jù)預(yù)測值endp=[0,0,1/2,1/2。1,0,0,0。0,1/2,0,1/2。1/2,0,0,1/2]。 %一步轉(zhuǎn)移概率矩陣vs=[0,1,0,0]。 %開始狀態(tài)向量for i=1:5v=vs*p^i %各預(yù)測年份的狀態(tài)向量for s=4:1:1if v(s)=max(v)e(i)=s。endendende=[2,1,4,4,1,3,3,2,e] %得到各年份所處的狀態(tài)區(qū)間E=[*xp,*xp,*xp] %劃分狀態(tài)區(qū)間的各點for i=1:n+5if e(i)=1 y(i)=x(i)+E(1)。else if e(i)=2 y(i)=x(i)+(E(1)+E(2))/2。else if e(i)=3 y(i)=x(i)+(E(2)+E(3))/2。else y(i)=x(i)+E(3)。endendx %灰色GM(1,1)模型預(yù)測結(jié)果y %灰色馬爾科夫鏈模型預(yù)測