hongcaisheng
2013-05-03, 12:28
http://d.hiphotos.baidu.com/zhidao/pic/item/b8014a90f603738d7a260b7eb21bb051f819ec7a.jpg
函数：
% 作1-AGO 生成序列x(1)
for i=1:n
x1(i)=sum(x0(1:i));
% 紧邻均生成
for k=2:n % 紧邻均生成 z
z(k)=0.5*x1(k)+0.5*x1(k-1);
end
% 产生数据矩阵B，计算系数a 和u
for i=1:n-1
b(i,1)=-z(i+1);
y(i)=x0(i+1);
end
b(:,2)=1;
y=y'; % 转置为列向量
au=b\y; % 作矩阵除法，计算a u
% 计算GM(1,1)模型的值
yc1(1)=x0(1);
for k=1:n
c=x0(1)-au(2)/au(1);
yc1(k+1)=c*exp( -au(1)*k)+au(2)/au(1);
end
% 显示预测结果
yc0(1)=x0(1);
for k=1:n
yc0(k+1)=yc1(k+1)-yc1(k);
end
disp(uint16(yc0(2:1:n+1)));
函数：
% 作1-AGO 生成序列x(1)
for i=1:n
x1(i)=sum(x0(1:i));
% 紧邻均生成
for k=2:n % 紧邻均生成 z
z(k)=0.5*x1(k)+0.5*x1(k-1);
end
% 产生数据矩阵B，计算系数a 和u
for i=1:n-1
b(i,1)=-z(i+1);
y(i)=x0(i+1);
end
b(:,2)=1;
y=y'; % 转置为列向量
au=b\y; % 作矩阵除法，计算a u
% 计算GM(1,1)模型的值
yc1(1)=x0(1);
for k=1:n
c=x0(1)-au(2)/au(1);
yc1(k+1)=c*exp( -au(1)*k)+au(2)/au(1);
end
% 显示预测结果
yc0(1)=x0(1);
for k=1:n
yc0(k+1)=yc1(k+1)-yc1(k);
end
disp(uint16(yc0(2:1:n+1)));