灰色預(yù)測(cè)MATLAB程序

作用：求累加數(shù)列、求ab的值、求展望方程、求殘差clc%清屏，以使結(jié)果獨(dú)立顯示x=[];formatlong;%設(shè)置計(jì)算精度iflength(x(:,1))==1%對(duì)輸入矩陣進(jìn)行判斷，如不是一維列矩陣，進(jìn)行轉(zhuǎn)置變換x=x';endn=length(x);%取輸入數(shù)據(jù)的樣本量z=0;fori=1:n%計(jì)算累加值，并將值給予矩陣bez=z+x(i,:);be(i,:)=z;endfori=2:n%對(duì)原始數(shù)列平行移位y(i-1,:)=x(i,:);endfori=1:n-1%計(jì)算數(shù)據(jù)矩陣B的第一列數(shù)據(jù)c(i,:)=*(be(i,:)+be(i+1,:));endforj=1:n-1%計(jì)算數(shù)據(jù)矩陣B的第二列數(shù)據(jù)e(j,:)=1;endfori=1:n-1%結(jié)構(gòu)數(shù)據(jù)矩陣BB(i,1)=c(i,:);B(i,2)=e(i,:);endalpha=inv(B'*B)*B'*y;%計(jì)算參數(shù)矩陣即ab的值fori=1:n+1%計(jì)算數(shù)據(jù)預(yù)計(jì)值的累加數(shù)列，如改為

n+1為

n+m可展望后

m-1個(gè)值ago(i,:)=(x(1,:)-alpha(2,:)/alpha(1,:))*exp(-alpha(1,:)*(i-1))+alpha(2,:)/alpha(1,:);%顯示輸出展望值的累加數(shù)列endvar(1,:)=ago(1,:)%顯示輸出展望值fori=1:n%如改n為n+m-1，可展望后m-1個(gè)值var(i+1,:)=ago(i+1,:)-ago(i,:);%預(yù)計(jì)值的累加數(shù)列的復(fù)原，并計(jì)算出下一展望值endfori=1:nerror(i,:)=x(i,:)-var(i,:);%計(jì)算殘差endc=std(error)/std(x);%調(diào)用統(tǒng)計(jì)工具箱的標(biāo)準(zhǔn)差函數(shù)計(jì)算后驗(yàn)差的比值ago%顯示輸出展望值的累加數(shù)列alpha%顯示輸出參數(shù)數(shù)列var%顯示輸出展望值error%顯示輸出偏差c%顯示后驗(yàn)差的比值

c作用：數(shù)據(jù)辦理

判斷能否能夠用灰色展望

、求級(jí)比、求累加數(shù)列、求

ab

的值、求展望方程clc,clearx0=[]';%注意這里為列向量n=length(x0);lamda=x0(1:n-1)./x0(2:n)%計(jì)算級(jí)比range=minmax(lamda')%計(jì)算級(jí)比的范圍x1=cumsum(x0)%累加運(yùn)算B=[*(x1(1:n-1)+x1(2:n)),ones(n-1,1)];Y=x0(2:n);u=B\Y%擬合參數(shù)u(1)=a,u(2)=bx=dsolve('Dx+a*x=b','x(0)=x0');%求微分方程的符號(hào)解x=subs(x,{

'a','b','x0'},{u(1),u(2),x0(1)})

%代入預(yù)計(jì)參數(shù)值和初始值yuce1=subs(x,'t'

,[0:n-1]);

%求已知數(shù)據(jù)的展望值