| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1972 人关注过本帖
标题:求助:灰色模型符合预测M程序
只看楼主 加入收藏
hexcimal
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2008-6-16
收藏
 问题点数:0 回复次数:2 
求助:灰色模型符合预测M程序
有人对灰色模型在符合预测中的应用M程序熟悉吗,需要20个历史数据(如20个一点数据) 预测某一天的1点的符合数据,哪位朋友帮帮忙 ,万分感谢!
搜索更多相关主题的帖子: 灰色模型 预测 
2008-06-16 11:56
mervynabc
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2009-8-21
收藏
得分:0 
function []=greymodel(y)
% 本程序主要用来计算根据灰色理论建立的模型的预测值。
% 应用的数学模型是 GM(1,1)。
% 原始数据的处理方法是一次累加法。
% y=input('请输入数据 ');
n=length(y);
yy=ones(n,1);
yy(1)=y(1);
for i=2:n
    yy(i)=yy(i-1)+y(i);
end
B=ones(n-1,2);
for i=1:(n-1)
    B(i,1)=-(yy(i)+yy(i+1))/2;
    B(i,2)=1;
end
BT=B';
for j=1:n-1
    YN(j)=y(j+1);
end
YN=YN';
A=inv(BT*B)*BT*YN;
a=A(1);
u=A(2);
t=u/a;
t_test=input('请输入需要预测个数:');
i=1:t_test+n;
yys(i+1)=(y(1)-t).*exp(-a.*i)+t;
yys(1)=y(1);
for j=n+t_test:-1:2
    ys(j)=yys(j)-yys(j-1);
end
x=1:n;
xs=2:n+t_test;
yn=ys(2:n+t_test);
plot(x,y,'^r',xs,yn,'*-b');
det=0;
for i=2:n
    det=det+abs(yn(i)-y(i));
end
det=det/(n-1);
disp(['百分绝对误差为:',num2str(det),'%']);
disp(['预测值为: ',num2str(ys(n+1:n+t_test))]);





灰度预测.rar (952 Bytes) 可以试试

2009-08-21 22:59
kiku
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2009-10-16
收藏
得分:0 
非常感谢,我去试试看
2009-10-16 16:45
快速回复:求助:灰色模型符合预测M程序
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.025252 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved