数学建模社区-数学中国

标题: 灰色预测模型的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

等级4 W" c0 v) V9 `; Y! b1 t* r4 m	相对误差q * I R4 ^" o. B0 H" C	方差比C) J6 ]! ?& q- R' i# j	小误差概论P n9 _4 N+ v0 I5 v$ ^( y
I级0 V% m& Z% O2 P! ?& q3 R	<0.01 , h* w5 M5 C2 R$ Y# R2 m! E! L	<0.352 x+ r; F( v6 o5 D: D. _& p	>0.959 M3 t; D" [$ z7 b- V6 D0 l
II级 / ?7 H5 r; B* Z3 ]7 k' ?	<0.05' e1 ]0 T# g5 `& t	<0.50 : P$ h% K9 E1 F( c	<0.800 d; u; P, C9 o) w
III级 " m2 r! ?# T' z. {5 i	<0.10 # W# `! p6 B! i% F, L0 N- ^	<0.65# \- `' K6 v6 q	<0.70' g! p: M! \|0 G2 F
IV级: G% }& \( t' F# `3 h( J/ P	>0.20 3 s/ J! i4 Y' M& y	>0.80- a: F D& Q9 p9 u F( ]	<0.60 : \|) O, [2 E6 r0 W+ P

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

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

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