数学建模社区-数学中国»论坛 › 【科学软件论坛】数学建模 › MATLAB论坛 › Matlab 资料库视频教程讲义代码 › 灰色预测模型的Matlab 程序及检验程序

查看: 3802|回复: 2

[代码资源] 灰色预测模型的Matlab 程序及检验程序

[复制链接]

字体大小: 正常放大

wangkai1

9 主题	8 听众	449 积分

升级 49.67%

TA的每日心情

	擦汗 2016-4-9 09:36

签到天数: 131 天

[LV.7]常住居民III

群组: Matlab讨论组

群组: 国赛讨论

群组: 科技写作基础培训

电梯直达

1^#

发表于 2014-8-31 09:33 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

本帖最后由 wangkai1 于 2014-8-31 09:35 编辑

灰色预测模型的Matlab 程序及检验程序%灰色预测模型程序

clear

syms a b;

c=[a b]';

A=[46.2 32.6 26.7 23.0 20.0 18.9 17.5 16.3];% 原始序列

B=cumsum(A);%累加

n=length(A);

for i=1: (n-1)

C(i)=(B(i)+B(i+1))/2;

end

%计算待定参数

D=A;

D(1)=[];

D=D';

E=[-C; ones(1,n-1)];

c=inv(E*E')*E*D;

c=c';

a=c(1);

b=c(2);

%预测往后预测5个数据

F=[];F(1)=A(1);

for i=2:(n+5)

F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;

end

G=[];G(1)=A(1);

for i=2: (n+5)

G(i)=F(i)-F(i-1);

end

t1=2002:2009;

t2=2002:2014;

plot(t1,A,'o',t2,G)

%灰色预测模型检验程序

function [ q,c,p ] = checkgm( x0,x1 )

%GM 检验函数

%x0 原始序列

%x1 预测序列

%·返回值

% q –- 相对误差

% c -- ·方差比

% p -- 小误差概率

e0=x0-x1;

q=e0/x0;

s1=var(x0);

%qpa=mean(e0);

s2=var(e0);

c=s2/s1;

len=length(e0);

p=0;

for i=1:len

if(abs(e0(i)) < 0.6745*s1)

p=p+1;

end

p=p/len;

end

等级. I9 p5 ]$ c; u. [* M( L Q	相对误差q 9 ?. N0 Y; P4 c5 f. P1 A	方差比C) `, T: I* ^' M	小误差概论P Q0 H- M+ {4 U
I级 * v3 o" W% l: w	<0.01 & A7 M3 Y! l/ h- z: Y: {	<0.35 ' G( K) e. g0 n ]	>0.95" P) J. f, D2 }1 R# }/ c4 E
II级 ' f8 F4 f) {9 p4 i+ P/ t0 A) J% }3 g	<0.05, v# P0 `' x1 Q. w	<0.50 2 p& r% z4 _4 e/ e# R6 z5 h	<0.80/ g3 u$ r* K* D: D: m' m: ^! y* D
III级- ^+ Q: p, y2 J7 E4 ]6 `	<0.10 ) C5 G: o' T6 B6 }% P4 F# N5 _	<0.65 " Q- W6 I/ B' }4 g0 w	<0.70/ y U$ q. e- i+ P8 z6 p
IV级2 k9 f# P3 i9 o5 O/ E0 s' @+ S6 w	>0.209 p1 F# F; k$ Q8 [0 e7 W" d! l( v9 d, k	>0.80 8 w$ w4 [0 w' c0 q. h7 [	<0.60" O2 R6 Z* L: m' S: m6 G+ S