查看: 10725|回复: 3

[代码资源] 灰色预测的小误差概率怎么算呢？

字体大小: 正常放大

7 主题	8 听众	169 积分

升级 34.5%

TA的每日心情

	开心 2014-9-1 21:31

签到天数: 38 天

[LV.5]常住居民I

自我介绍: 本科生，新手。

电梯直达

1^#

发表于 2014-7-8 22:38 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

哪位大神来帮帮我吧，那个灰色预测中在检验的时候会有一个小误差概率的计算，就是这个P=P{|ε(i)-ε拔|＜0.674S平方}，怎么用matlab实现呢？

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

平凡之不凡

90 主题	105 听众	9482 积分

TA的每日心情

	怒 2017-1-30 23:35

签到天数: 318 天

[LV.8]以坛为家I

2010挑战赛参赛者

2011挑战赛参赛者

2012挑战赛参赛者

2013挑战赛参赛者

群组: 物联网工程师培训

群组: Matlab讨论组

群组: 2013认证赛C题讨论群组

群组: 数学建模

群组: 2013认证赛A题讨论群组

2^#

发表于 2014-7-8 23:44 |只看该作者 |招呼Ta 关注Ta

很多帖子有，你看一下

1.png (3.99 KB, 下载次数: 906)

2.png (1.27 KB, 下载次数: 821)

点评

红莲业火谢谢！！！十分感谢！！！！我又找到了一份。 %附：预测精度等级 P C % 好 >0.95 0.80 0.70 详情回复发表于 2014-7-17 23:16

我也说一句

发表

使用道具举报

红莲业火

7 主题	8 听众	169 积分

升级 34.5%

TA的每日心情

	开心 2014-9-1 21:31

签到天数: 38 天

[LV.5]常住居民I

自我介绍: 本科生，新手。

3^#

发表于 2014-7-17 23:16 |只看该作者 |招呼Ta 关注Ta

平凡之不凡发表于 2014-7-8 23:44
0 F) I/ x! }7 I% V很多帖子有，你看一下

谢谢！！！十分感谢！！！！我又找到了一份。

%附：预测精度等级 P C
% 好 >0.95 <0.35
% 合格 >0.80 <0.45
% 勉强 >0.70 <0.50
% 不合格 ≤0.70 ≥0.65

disp('请输入初始序列x （例 [1 2 3 4 5] ）');
x=input('x=');
j=length(x);
disp('请输入要预测的年份数k');
k=input('k=');
z(1)=x(1);
for i=2:j
z(i)=z(i-1)+x(i);
end
for i=1

j-1)
b(i,1)=-0.5*(z(i)+z(i+1));b(i,2)=1;
end
for i=1

j-1)
y(i,1)=x(i+1);
end
g=inv(b'*b)*b'*y;
a=g(1,1);u=g(2,1);
for i=1

j+k)
t(i)=(z(1)-u/a)*exp((-a)*(i-1))+u/a;
end
ans(1)=t(1);
for i=2:j
ans(i)=t(i)-t(i-1);
end
for i=(j+1)

j+k)
ansf(i-j)=t(i)-t(i-1);
end
for i=1:j
E(i)=x(i)-ans(i);
end
e=E./x;
xj=sum(x)/j;
S1=sqrt(sum((x-xj).^2)/j);
Ej=sum(E)/j;
S2=sqrt(sum((E-Ej).^2)/j);
C=S2/S1;
Index=abs(E-Ej)<0.6745*S1;
P=length(find(Index>0.5))/length(Index);

disp('原始数列x=');disp(x);
disp('拟合数列ans=');disp(ans);
disp('预测数列ansf=');disp(ansf);
disp('后验差比值C=');disp(C);
disp('小误差概率P=');disp(P);
if C<0.35&&

>0.95
disp('此模型的预测精度为：好！')
else if C<0.45&&

>0.80
disp('此模型的预测精度为：合格')
else if C<0.50&&

>0.70
disp('此模型的预测精度为：勉强')
else disp('此模型的预测精度为：不合格，请尝试其他方法！')
end
end
end