数学建模社区-数学中国

标题: 灰色预测模型的Matlab 程序及检验程序 [打印本页]

作者: wangkai1 时间: 2014-8-31 09:33
标题: 灰色预测模型的Matlab 程序及检验程序
本帖最后由 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

等级$ k* q& g6 j) o l! D Q% z: y	相对误差q" K. ]! r3 W& }3 x1 p	方差比C3 X6 ^' Q+ [+ U' m	小误差概论P) y7 M! U V" `- F0 l* G0 h; V% m4 \5 I
I级 % j) p g; h+ J4 y/ D	<0.018 v8 A( W8 q" G' v# R! w	<0.35 ; q1 P" ~8 F& `# V9 u! O3 b	>0.95 2 T1 Y; {2 a, F0 F" u1 c
II级 y: v% y* r$ B9 b- B1 p	<0.05 6 H- G0 ^0 E, C	<0.50 ; Q' B+ n. I) \|* U/ J	<0.801 l; t+ ?6 \|( O+ l
III级7 h. ^( T" R K, E2 D. U	<0.107 C& B* l- I# b* f& `! H	<0.65 + z2 L0 H& ~' E5 e	<0.70$ G: f7 X0 L# A1 Y
IV级( d4 h$ f: \|5 n9 c3 B( [8 m	>0.20 , E7 M7 _$ q6 Y# P6 v( I& ^3 U	>0.80 ) e$ y B. `" x5 Z ~ _	<0.60 5 d" \|: t1 }' u9 s1 r# R) u

作者: 流雨星月 时间: 2014-8-31 14:35
顶一个。。。。

作者: 梦@di?~ 时间: 2014-8-31 18:09
不错，楼主辛苦了。。。

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)