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

yqxxqy

小弟正在编写GNNM（1,1）的程序，但是就是调试不出来，心里那个急啊~~~
有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~
有没有现成的程序供参考参考啊~~~
# S. \, w! e, d- X9 d# A  N

zjqylcy

现在的人怎么都那么懒了~~~诶。。。
% a7 ]  ?  Q" x. a3 R3 S5 X* }' Z) D这个是用来做人口预测的，你参看吧~~~
$ J+ q$ [4 O- j5 o2 j0 P

GM( 1, 1) 的matlab程序

gm.m:

function f=gm(x0,m)
" X! `5 l( Y: @, \; A; ~! |%定义为函数gm(x)

n=length(x0);

x1=zeros(1,n);

x1(1)=x0(1);

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

5 f& z6 ^( A  ]! O( Nx1(i)=x1(i-1)+x0(i);

end

i=2:n;
* @4 v0 k( C# z+ M%对原始数列平行移位并负值给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;
%计算参数a,u矩阵

i=1:m;
/ B  H- c: Q+ n8 B$ `$ A" X1 i8 B%计算预测累加数列的值

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;

%还原数列的值

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

i=2:n;

error(i)=yc(i)-x0(i);
3 \+ Y& t1 t6 O. i1 L- z%计算残差值

yc(1)=ago(1);

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

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

c=std(error)/std(x0);
%计算后验差比

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

[nrow,ncol]=size(relerror);

p=0;

for i=2:ncol

9 f9 }; f6 Q; p: Z* Aif
3 `3 D+ [' Z  n& K+ L/ xrelerror(1,i)<0.6745*std(x0)

; K+ O' u' Y6 u/ ep=p+1;

end

end

p=p/(n-1);
& I. W' w2 t7 B! j7 i
+ l+ Z+ u$ g$ q% O

w1=min(abs(error));

w2=max(abs(error));

i=1:n;
& _7 s( C2 g4 b2 O1 d  p' K# k
%计算关联度

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

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

au
% p, ^8 b! {) _  C# z& Y%输出参数a,u的值

ago;
5 _" z8 n0 r/ c9 }% p6 w8 p9 o: I2 P5 E%输出累加数列ago的值

x0;
. t% B  U- [) ^# d' S$ B1 H  H%输出原始序列值

f=yc;
0 x1 j" W5 O# z2 @" u5 J
%输出预测的值

error;
%输出残差的值

c;
1 a/ S8 W' r" {9 r. @) @, A%输出后验差比的值

p;
; |( Z' U' c7 c) F# i* `0 s) R%输出小误差概率的值

w
- F: W* o5 |4 @# {" D" T, ]& _% M%输出关联度

end

' Q0 y. m+ j6 \9 @! F3 y

& t5 Z. Z$ U* v' M0 e% x; c

diliwang

同求，挺感兴趣

yqxxqy

回复 zjqylcy 的帖子

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

gaoshanliu水

linmatsas

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

jshzncd

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

海之舟

好高深啊~~~~~~