mgm.m

基于泰勒级数的灰度预测

function [gg,emspp]=mgm(t,x,n1)
cc=size(x);%定义数据的个数.
x0=x;%产生原始数据
nn=cc(1);n=cc(2);
%[a b]'=(A'*A)^(-1)*A'*X
%准备相关数据
B=[];y=[];x1=[];A=[];X=[];
for j=1:
x1(1,1)=x0(1,1);%定义累加第一个数据
%----------计算矩阵A,x1----------------------------------
for i=1:n-1
    x1(1,i+1)=x1(1,i)+x0(1,i+1);   
    B(i,1)=-(x1(1,i+1)+x1(1,i))/2;
    X(i,1)=x0(1,i+1);
end
A=[B,ones(n-1,1)];
auz=(A'*A)^(-1)*A'*X;
%------ 预测值计算---------
g0=[]; ems=[];
g0(1,1)=x(1,1);ems(1,1)=0;a=auz(1);b=auz(2);
for i=2:n1
  g0(1,i)=(x(1,1)-b/a)*(1-exp(a))*exp(-a*i);
end
%for i=1:n
 %   ems(1,i)=abs(g0(1,i)-x(1,i));
%end
%gg=g0';emspp=ems';