模糊聚類(lèi)分析例子1-資料下載頁(yè)
2025-06-26 07:29本頁(yè)面
【正文】 u(2),u(3),x1(1),x1(6)})。yuce=subs(x,39。t39。,0:n1)。digits(6),x=vpa(x) %x為時(shí)間響應(yīng)式x0_hat=[yuce(1),diff(yuce)] %預(yù)測(cè)的數(shù)據(jù)epsilon=x0x0_hat %計(jì)算殘差delta=abs(epsilon./x0) %計(jì)算相對(duì)誤差上海證券交易所綜合指數(shù)的波形預(yù)測(cè)。根據(jù)上海證券交易所綜合指數(shù)的周收盤(pán)指數(shù)數(shù)據(jù)，從1997年2月21日到1998年10月31日的周收盤(pán)指數(shù)曲線(xiàn)如圖所示：解：取9條等間隔的等高線(xiàn)，分別為=1140，=1170，=1200，=1230，=1260，=1290，=1320，=1350，=1380的等高時(shí)刻序列分別為對(duì)應(yīng)于=1140，=（,41,）對(duì)應(yīng)于=1170，=（,23,,76,）對(duì)應(yīng)于=1200，=1230，=1260，=1290，=1320，=1350，=1380分別為=（,,,）=（,,,85）=（7,,）=（,）=（,,）=（,,）=（,69）對(duì)（i=1,2,3,…9）序列，進(jìn)行GM（1,1）預(yù)測(cè)，起響應(yīng)時(shí)分別為y1=+*exp(*t)y2=+*exp(*t)y3=+*exp(.166865*t)y4=+*exp(.159938*t)y5=+*exp(.446077*t)y6=+*exp(.550388*t)y7=+*exp(.191636*t)y8=+*exp(.185059*t)y9=+*exp(.488018*t)對(duì)在1998年11月到2000年3月這5個(gè)月進(jìn)行預(yù)測(cè)，可得等高時(shí)刻的預(yù)測(cè)序列 =（,） =（,） =（,） =（,） =（,） =（） =（,） =（,） =（） 據(jù)此可以畫(huà)出上海證券交易所綜合指數(shù)1998年11月到2000年3月的預(yù)測(cè)波形圖 如下 Matlab程序如下： %GM（1,1）函數(shù) function gm11(x0,t)n=length(x0)。x1=cumsum(x0)。lamda=x0(1:n1)./x0(2:n)。for i=2:nz(i)=*(x1(i)+x1(i1))。endB=[z(2:n)39。,ones(n1,1)]。Y=x0(2:n)39。u=B\Y。x=dsolve(39。Dx+a*x=b39。,39。x(0)=x039。)。x=subs(x,{39。a39。,39。b39。,39。x039。},{u(1),u(2),x1(1)})。yuce1=subs(x,39。t39。,[0:n1])。digits(6),y=vpa(x) yuce1=subs(x,39。t39。,[0:n1+t])。yuce=[x0(1),diff(yuce1)]%預(yù)測(cè)數(shù)據(jù)clc,clearx1=[,41,]。x2=[,23,,76,]。x3=[,,,]。x4=[,,,85]。x5=[7,,]。x6=[,]。x7=[,,]。x8=[,,]。x9=[,69]。gm11(x1,4)gm11(x2,4)gm11(x3,4)gm11(x4,4)gm11(x5,4)gm11(x6,4)gm11(x7,4)gm11(x8,4)gm11(x9,4)
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1