灰色預(yù)測MATLAB程序

1、.灰色預(yù)測 作用：求累加數(shù)列、求a b的值 、求預(yù)測方程、求殘差clc %清屏，以使結(jié)果獨立顯示x=71.1 72.4 72.4 72.1 71.4 72.0 71.6;format long; %設(shè)置計算精度if length(x(:,1)=1 %對輸入矩陣進行判斷，如不是一維列矩陣，進行轉(zhuǎn)置變換 x=x;endn=length(x); %取輸入數(shù)據(jù)的樣本量z=0;for i=1:n %計算累加值，并將值賦予矩陣be z=z+x(i,:); be(i,:)=z;endfor i=2:n %對原始數(shù)列平行移位 y(i-1,:)=x(i,:);endfor i=1:n-1 %計算數(shù)據(jù)矩陣B的第一

2、列數(shù)據(jù) c(i,:)=-0.5*(be(i,:)+be(i+1,:);endfor j=1:n-1 %計算數(shù)據(jù)矩陣B的第二列數(shù)據(jù) e(j,:)=1;endfor i=1:n-1 %構(gòu)造數(shù)據(jù)矩陣B B(i,1)=c(i,:); B(i,2)=e(i,:); endalpha=inv(B*B)*B*y; %計算參數(shù) 矩陣 即a b的值for i=1:n+1 %計算數(shù)據(jù)估計值的累加數(shù)列，如改為n+1為n+m可預(yù)測后m-1個值 ago(i,:)=(x(1,:)-alpha(2,:)/alpha(1,:)*exp(-alpha(1,:)*(i-1)+alpha(2,:)/alpha(1,:);%顯示輸出

3、預(yù)測值的累加數(shù)列endvar(1,:)=ago(1,:) %顯示輸出預(yù)測值for i=1:n %如改n為n+m-1，可預(yù)測后m-1個值 var(i+1,:)=ago(i+1,:)-ago(i,:); %估計值的累加數(shù)列的還原，并計算出下一預(yù)測值endfor i=1:n error(i,:)=x(i,:)-var(i,:); %計算殘差endc=std(error)/std(x); %調(diào)用統(tǒng)計工具箱的標(biāo)準(zhǔn)差函數(shù)計算后驗差的比值cago %顯示輸出預(yù)測值的累加數(shù)列alpha %顯示輸出參數(shù) 數(shù)列var %顯示輸出預(yù)測值error %顯示輸出誤差c %顯示后驗差的比值作用：數(shù)據(jù)處理 判斷是否可以用灰

4、色預(yù)測 、求級比、求累加數(shù)列、求a b的值 、求預(yù)測方程clc,clearx0=71.1 72.4 72.4 72.1 71.4 72.0 71.6; %注意這里為列向量n=length(x0);lamda=x0(1:n-1)./x0(2:n) %計算級比range=minmax(lamda) %計算級比的范圍x1=cumsum(x0) %累加運算B=-0.5*(x1(1:n-1)+x1(2:n),ones(n-1,1);Y=x0(2:n);u=BY %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(Dx+a*x=b,x(0)=x0); %求微分方程的符號解x=subs(x,a,b,x0,u(1),u(2),x0(1) %代入估計參數(shù)值和初始值yuce1=subs(x,t,0:n-1); %求已