灰色预测工具箱 各种灰色模型

function [px0,ab,rel]=gm11(x0,number) %[px0,ab,rel]=gm11(x0,number) %px0为预测数列，rel为平均相对误差,rel为平均相对误差(为百分比) %默认的number参数为原数组大小 if nargin==1 number=max(size(x0)); end n=max(size(x0)); x1=zeros(size(x0)); for k=1:n for i=1:k x1(k)=x1(k)+x0(i); end end z=zeros(size(x0)); for k=2:n z(k)=0.5*(x1(k)+x1(k-1)); end y=x0'; y(1)=[]; b(:,1)=-z'; b(:,2)=1; b(1,:)=[]; ab=inv(b'*b)*b'*y; a=ab(1); b=ab(2); px0(1)=x0(1); %求还原值系列 for k=1:number-1 px0(k+1)=(1-exp(a)) * ( x0(1)-b/a ) * exp(-a*k); end temp=px0(1:n); x0; temp=(temp-x0)./x0; %相对误差 temp(1)=[]; %删除第一个为零的误差 temp=abs(temp); rel=sum(temp)/(n-1)*100;