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

yqxxqy

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

zjqylcy

现在的人怎么都那么懒了~~~诶。。。
6 G2 R8 ~5 D( U* N) ]' q8 C, z( g- r这个是用来做人口预测的，你参看吧~~~
" }7 v6 d) N6 _6 E

GM( 1, 1) 的matlab程序

gm.m:

function f=gm(x0,m)
, R' R3 l) r. |%定义为函数gm(x)

n=length(x0);

x1=zeros(1,n);

x1(1)=x0(1);

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

5 d$ B/ P6 t) {& {8 b* n5 C  Px1(i)=x1(i-1)+x0(i);

end

i=2:n;
%对原始数列平行移位并负值给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;
( V- Y8 X7 S9 w+ D' d%计算参数a,u矩阵

i=1:m;
%计算预测累加数列的值

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;
  x- l; E- i$ O9 i  t, X8 Q# Y
%还原数列的值

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

i=2:n;

error(i)=yc(i)-x0(i);
/ P! n8 D: ]. ]+ T" R' r& ~%计算残差值

yc(1)=ago(1);

i=1:m-1;
%修正还原数列的值

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

c=std(error)/std(x0);
6 c7 p- ^! B0 g  A%计算后验差比

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

[nrow,ncol]=size(relerror);

p=0;

for i=2:ncol

if
# b$ r' a; o9 i: y& c4 Jrelerror(1,i)<0.6745*std(x0)

p=p+1;

3 R$ p% k, l7 D8 yend

end

p=p/(n-1);

w1=min(abs(error));

w2=max(abs(error));

i=1:n;
+ C( S8 G  A2 u/ o$ F: ?; C
%计算关联度

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

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

au
, [: \% _0 s0 j) c3 N%输出参数a,u的值

ago;
: g4 c7 n- k$ K. O%输出累加数列ago的值

x0;
3 Y# X" I5 g1 I$ ]0 A%输出原始序列值

f=yc;
( H! A, p# S0 h4 a9 O1 ~
3 ~7 @: H. a/ ^%输出预测的值

error;
& }6 L( s; O6 G5 x%输出残差的值

c;
7 T* i# }3 C; }- H# E7 n%输出后验差比的值

p;
0 z# i- r  q* n) i0 V& I& E%输出小误差概率的值

w
%输出关联度

end

9 l1 n( W  ^" k8 i! k

7 \0 |" n1 e' f4 g

diliwang

同求，挺感兴趣

yqxxqy

回复 zjqylcy 的帖子
, c- p: E& I4 I4 `: y
3 N' P4 u5 Y3 M( l我是想求助关于灰色神经网络GNNM（1,1）的。。。不过还是谢了~~~

gaoshanliu水

linmatsas

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

jshzncd

其实我一直觉得楼主的品味不错！呵呵！
. H3 ]1 `# X) x* F: u: {数学中国社区 不走平凡路

海之舟

好高深啊~~~~~~