MATLAB 实现灰色预测程序
灰色预测
很好的东西呐,······~~··`~··~~~~~~~~~~~~~~~~~~~~~`````````````
fon [feval,au,ec,C,P]=GM1_1(x, r)
if nrgin<2
myar=0;
end
[mx,nx]=size(x);
if mx==1
x=x';
end
n=length(x);
for i=2:n
z(i-1)=0.5*x1(i)+0.5*x1(i-1);
end
Y=x(2:end);
B(:,1)=-z;
2)/au(1));
yc(1)=x(1);
for k=1:n+myear-1
y1(k+1)=pm*exp(-au*k)+a(2)/au(1);
yc(k+1)=y1(k+1)-y1(k);
end
feval=yc';
ex=ec./x;
r=0;
rou=0.5;
for k=1:n
r=r+rou* s(ec(k))+rou*max(a (ec)));
end
r=r/n;
%原始序列的标准差
s1=std(x);