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

yqxxqy

小弟正在编写GNNM（1,1）的程序，但是就是调试不出来，心里那个急啊~~~
0 X; F: O& {2 @2 U有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~
有没有现成的程序供参考参考啊~~~

zjqylcy

现在的人怎么都那么懒了~~~诶。。。
这个是用来做人口预测的，你参看吧~~~
  l9 \7 S* p3 r! ~3 M7 y: ]4 e

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
' K9 {3 h5 N7 z) U# o3 U%计算累加序列x1

x1(i)=x1(i-1)+x0(i);

end

i=2:n;
2 Y. S$ A8 p; l3 f8 N4 ]. l1 i8 v7 y%对原始数列平行移位并负值给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;
8 h0 g  ~9 ~: D3 x  f( v7 ~%计算参数a,u矩阵

i=1:m;
+ T/ a! Z8 L' o. o0 k, A) A% f%计算预测累加数列的值

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;
( I  M$ m4 R5 b- A  L
5 o* E( X3 R) b9 l%还原数列的值

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

i=2:n;

error(i)=yc(i)-x0(i);
/ R. C/ f1 M  ]! a0 A! k$ p%计算残差值

yc(1)=ago(1);

i=1:m-1;
5 O9 G2 {" g3 i4 b: ^% x6 q%修正还原数列的值

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

c=std(error)/std(x0);
$ v% _+ D. B. c7 x7 t%计算后验差比

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

[nrow,ncol]=size(relerror);

p=0;

for i=2:ncol

if
relerror(1,i)<0.6745*std(x0)

( }% c. m; A  o0 r% s; b$ gp=p+1;

* ^0 f7 ^( z; q( [9 h: Nend

end

p=p/(n-1);
- A2 E0 P+ Y1 f! D8 l
; h  c) O7 s7 _! B' x* j

w1=min(abs(error));

w2=max(abs(error));

i=1:n;
( S* f7 @, K/ ?! |4 q) B
%计算关联度

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

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

au
1 f4 x" F9 s1 d& z1 y%输出参数a,u的值

ago;
# q. r6 j+ u. O8 |0 V: L3 N%输出累加数列ago的值

x0;
%输出原始序列值

f=yc;
: K; H, G4 C) h
%输出预测的值

error;
* E0 s' ?# O- T4 P& I%输出残差的值

c;
%输出后验差比的值

p;
# }2 B) a1 j! u- x%输出小误差概率的值

w
4 n: {9 V% q0 J% {/ q1 k+ h! k%输出关联度

end

1 i, ^& s5 p. }9 X  o0 ]  L; _

3 A) r1 o/ ]/ @% r" U  ^) Q0 v8 K

diliwang

同求，挺感兴趣

yqxxqy

回复 zjqylcy 的帖子

我是想求助关于灰色神经网络GNNM（1,1）的。。。不过还是谢了~~~

gaoshanliu水

linmatsas

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

jshzncd

其实我一直觉得楼主的品味不错！呵呵！
数学中国社区 不走平凡路

海之舟

好高深啊~~~~~~