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

yqxxqy

小弟正在编写GNNM（1,1）的程序，但是就是调试不出来，心里那个急啊~~~
0 ]- G+ ^+ |: F1 J9 x  x  h有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~
有没有现成的程序供参考参考啊~~~
  r+ B4 g* Q. T; w8 x) _5 P* j

zjqylcy

现在的人怎么都那么懒了~~~诶。。。
2 G- X% c* Z. i8 X这个是用来做人口预测的，你参看吧~~~
) ]4 V9 Y$ d4 f; M

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

! e: d# L7 M4 y7 n+ l0 w# Tx1(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;
, y; C% s8 x: W- H* y6 b2 r%计算参数a,u矩阵

i=1:m;
! }5 `( X( [- [/ M/ o) `%计算预测累加数列的值

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);
%计算残差值

yc(1)=ago(1);

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

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

c=std(error)/std(x0);
: t% j& P9 r2 ^! }%计算后验差比

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

[nrow,ncol]=size(relerror);

p=0;

for i=2:ncol

# D! f+ v, d! M9 A. Z7 F, Cif
3 z- G4 B  G+ i/ g  brelerror(1,i)<0.6745*std(x0)

1 D6 }; c8 b0 X( O! Y4 ^( U4 Jp=p+1;

' T+ n3 v9 p) D" B& [1 hend

end

p=p/(n-1);

6 A) \! E9 w, h: v" S

w1=min(abs(error));

w2=max(abs(error));

i=1:n;

! Z" f4 O1 H( @* w0 X%计算关联度

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

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

au
8 _- @2 e) a/ r( \! e%输出参数a,u的值

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

x0;
# C, p2 e, H  ^! X! G* L) G+ e. M%输出原始序列值

f=yc;
; Y! s2 d- k; S; T1 ^- O
%输出预测的值

error;
3 {, h3 \4 K3 B5 Z# ^- \%输出残差的值

c;
4 W1 J: `1 I8 q: g%输出后验差比的值

p;
%输出小误差概率的值

w
%输出关联度

end
2 }/ Z; r; F' Y, p) R

: ]$ v% V6 b- y! i& S( |/ O

diliwang

同求，挺感兴趣

yqxxqy

回复 zjqylcy 的帖子
( ?6 S3 M& P8 m. v* _8 m9 ?4 U
) O, X( b+ W+ o我是想求助关于灰色神经网络GNNM（1,1）的。。。不过还是谢了~~~

gaoshanliu水

linmatsas

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

jshzncd

其实我一直觉得楼主的品味不错！呵呵！
* N' f5 c2 |3 ~. c( M数学中国社区 不走平凡路

海之舟

好高深啊~~~~~~