数学建模社区-数学中国

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

等级, L3 F/ Z% c7 k; Z& O+ C9 @3 V	相对误差q; O% b- o9 m6 a$ ^/ Y7 a4 \0 @	方差比C0 O1 X* @9 V0 Y7 {. H2 f	小误差概论P ! V! {9 s+ C0 e4 L9 V- I1 [! [' X7 Y
I级- n( ?! F w% F3 i9 O# W	<0.01 7 c0 o- V3 K! M5 p4 \! n5 Q! ^	<0.355 Z* M- Y$ g9 @7 H	>0.95$ M' N3 {3 j5 F% R- F. r
II级2 {" c& w5 _# E2 ]/ i	<0.05 # {1 j7 `; w% z: G	<0.50* q! j0 \& b) B$ Q" V	<0.804 L3 V, \6 C Z. Q; Q( `
III级 + L, ~+ l- W- \" v8 _	<0.10 U4 W2 C% Y" P6 j Q/ M	<0.65/ j) k- K% U2 a2 b Z	<0.70 " A" K' P. m* k& j. y
IV级9 p Q2 ]# c( x9 m4 V	>0.20 X5 a' Q* J4 t/ @9 J	>0.806 I5 L4 \+ u8 Y% k: c6 l0 ^	<0.605 B; W A" h4 x1 a. Y m( P* v

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

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

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