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

yqxxqy

小弟正在编写GNNM（1,1）的程序，但是就是调试不出来，心里那个急啊~~~
有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~
0 U' \2 c) V# P( E5 ^& W有没有现成的程序供参考参考啊~~~

海之舟

好高深啊~~~~~~

jshzncd

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

linmatsas

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

gaoshanliu水

yqxxqy

回复 zjqylcy 的帖子

1 L3 O7 U3 ^- L7 K3 D9 r7 g6 V0 }0 |7 q我是想求助关于灰色神经网络GNNM（1,1）的。。。不过还是谢了~~~

diliwang

同求，挺感兴趣

zjqylcy

现在的人怎么都那么懒了~~~诶。。。
9 i) ^: T: {& O: `/ _这个是用来做人口预测的，你参看吧~~~
& ]' z8 X% E6 ~  R: Q7 Q; e. J, p, O

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
; A/ [) W0 I: }1 c6 y0 H5 o1 F%计算累加序列x1

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

end

i=2:n;
; v* s) q4 ~! R8 H%对原始数列平行移位并负值给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;
% g$ J7 P+ O. m, }8 X%计算参数a,u矩阵

i=1:m;
% @0 h3 I  W3 Y8 j$ U* T. Y%计算预测累加数列的值

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;
9 I9 E0 v7 x# d* m$ r2 W9 G  S' }
" N* S7 |; ?3 R9 |1 w& `; K%还原数列的值

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

i=2:n;

error(i)=yc(i)-x0(i);
8 D) s7 }4 {, p- _%计算残差值

yc(1)=ago(1);

i=1:m-1;
7 R, `% i- b- P6 N0 C%修正还原数列的值

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

c=std(error)/std(x0);
2 {# I' Y2 A1 Q/ m6 K9 G: W9 ^%计算后验差比

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)

p=p+1;

end

end

p=p/(n-1);
; v) f* f" m8 k9 o! \  a

w1=min(abs(error));

w2=max(abs(error));

i=1:n;
2 v, h( ?" @! @+ b( L
! c# T, ?. D4 S%计算关联度

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

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

au
4 \% C2 T6 z# T3 @& u- P" x%输出参数a,u的值

ago;
%输出累加数列ago的值

x0;
4 k4 `! y! d% z' }%输出原始序列值

f=yc;

%输出预测的值

error;
%输出残差的值

c;
%输出后验差比的值

p;
, Y* g5 u  q2 H2 j$ e: X8 I' R%输出小误差概率的值

w
%输出关联度

end

7 t$ i! K7 U3 x0 q7 K: i

5 R- @2 ?& _( m2 ]/ Z7 u