紧急求助：关于灰色神经网络GNNM（1,1）的编程问题

yqxxqy

小弟正在编写GNNM（1,1）的程序，但是就是调试不出来，心里那个急啊~~~
9 A) m* t6 ]. a1 {3 |8 D; G! S8 b有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~
有没有现成的程序供参考参考啊~~~
% A+ c( V* V* M; x* h9 O3 d+ A% |

海之舟

好高深啊~~~~~~

jshzncd

其实我一直觉得楼主的品味不错！呵呵！
3 k* i; q9 o* F$ h& x: U+ f' a( u数学中国社区 不走平凡路

linmatsas

好高深呀。。。不懂不懂、、、、、

gaoshanliu水

yqxxqy

回复 zjqylcy 的帖子
/ a/ b" {' F8 k& [0 w
6 o" G/ E. C( w+ I4 I, x我是想求助关于灰色神经网络GNNM（1,1）的。。。不过还是谢了~~~

diliwang

同求，挺感兴趣

zjqylcy

现在的人怎么都那么懒了~~~诶。。。
这个是用来做人口预测的，你参看吧~~~
4 v/ {* ^# C! ~2 O( q1 N

GM( 1, 1) 的matlab程序

gm.m:

function f=gm(x0,m)
%定义为函数gm(x)

n=length(x0);

x1=zeros(1,n);

x1(1)=x0(1);

for i=2:n
%计算累加序列x1

: D' O+ g* t9 R8 n# P1 F0 d, ix1(i)=x1(i-1)+x0(i);

end

i=2:n;
' [! ]; }) `5 u  t7 W, ?%对原始数列平行移位并负值给y

y(i-1)=x0(i);

y=y';

i=1:n-1;

c(i)=-0.5*(x1(i)+x1(i+1));

B=[c' ones(n-1,1)];

au=inv(B'*B)*B'*y;
1 g: k0 \* m5 g& ]" P%计算参数a,u矩阵

i=1:m;
: p: w( ], `2 t; v) @%计算预测累加数列的值

ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1);

yc(1)=ago(1);

i=1:m-1;

. F2 n) |- g0 c6 Q2 I3 n" j" t( L%还原数列的值

yc(i+1)=ago(i+1)-ago(i);

i=2:n;

error(i)=yc(i)-x0(i);
) Z6 B) W. _5 \1 d  q; u8 m% _%计算残差值

yc(1)=ago(1);

i=1:m-1;
4 o% X3 T. B# A, ?%修正还原数列的值

yc(i+1)=ago(i+1)-ago(i);

c=std(error)/std(x0);
! Q% x* b+ P2 h! \( V$ o. k%计算后验差比

relerror=abs((error)-mean(error)*ones(size(error)));

[nrow,ncol]=size(relerror);

p=0;

for i=2:ncol

  a# \. N( I, U# T. E. Iif
relerror(1,i)<0.6745*std(x0)

+ H5 w: ?* ]) w" P8 Gp=p+1;

end

end

p=p/(n-1);

# T2 b% F- S+ t6 z

w1=min(abs(error));

w2=max(abs(error));

i=1:n;
+ l+ C6 M( m  h- V- x1 Q  q
+ L# P: H5 q* W* l! z1 f, w. M' ?%计算关联度

w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);

w=sum(w)/(n-1);

au
- U6 _+ p& f, Q/ V) @7 c%输出参数a,u的值

ago;
8 L" A1 N& p4 J% z+ Z%输出累加数列ago的值

x0;
3 @7 B0 R$ o' p5 l* w%输出原始序列值

f=yc;
+ v9 e- \( H) _5 d% K2 O
  a! q3 X- I2 z; N2 A%输出预测的值

error;
%输出残差的值

c;
%输出后验差比的值

p;
# `+ V$ D! ~$ I, p# n) t9 I%输出小误差概率的值

w
%输出关联度

end
! p' g9 n) X" J/ b

3 k  t% n# J( z5 B) S