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

yqxxqy

小弟正在编写GNNM（1,1）的程序，但是就是调试不出来，心里那个急啊~~~
有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~
" [! d  N" ?) K: v3 B有没有现成的程序供参考参考啊~~~
, q3 U1 H, a. b3 p! Z4 l4 j* ?5 M# n

zjqylcy

现在的人怎么都那么懒了~~~诶。。。
0 G" |. Z0 o' K/ n3 B1 ?这个是用来做人口预测的，你参看吧~~~

GM( 1, 1) 的matlab程序

gm.m:

function f=gm(x0,m)
) d7 v! D& O* |) z%定义为函数gm(x)

n=length(x0);

x1=zeros(1,n);

x1(1)=x0(1);

for i=2:n
4 T* x' Z5 ?" v%计算累加序列x1

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

end

i=2:n;
2 d( W8 p) @+ z5 ~/ T- b! c%对原始数列平行移位并负值给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;
, h+ W6 ]- G0 L$ _4 \%计算参数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;
0 C" D: L: p8 c! |2 D) k( s1 E
6 Z' y5 h" o8 T2 o%还原数列的值

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;
3 Q% z( q8 \+ }! l$ C  Q6 A6 ~; I. o( e%修正还原数列的值

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

) W) M. _. m+ T/ Oif
* ?4 \7 C$ K# l# `0 Q1 Srelerror(1,i)<0.6745*std(x0)

+ q& b- g+ K+ O1 h5 X$ w* V1 np=p+1;

end

end

p=p/(n-1);
% o$ k- W/ W& q' L" |4 \; M: e4 ?5 Y
; z, R# N! `; N% k7 W

w1=min(abs(error));

w2=max(abs(error));

i=1:n;

5 r' X5 F7 |( s) y( N- S%计算关联度

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

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

au
/ c* [2 V( u' e%输出参数a,u的值

ago;
. {+ l) b3 T: A% P5 H& e0 v%输出累加数列ago的值

x0;
%输出原始序列值

f=yc;

%输出预测的值

error;
%输出残差的值

c;
1 ?" L8 j3 h+ d%输出后验差比的值

p;
5 f0 B% d. c/ C8 }( ?% V) Y3 ~%输出小误差概率的值

w
% U7 {' s9 w& L6 n, i+ Z9 c, y8 \%输出关联度

end

; k- Y& b2 L; o: H: @  T5 c: g

% j0 i+ f* d$ k8 T0 b

diliwang

同求，挺感兴趣

yqxxqy

回复 zjqylcy 的帖子

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

gaoshanliu水

linmatsas

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

jshzncd

其实我一直觉得楼主的品味不错！呵呵！
1 N+ {; @5 n' L: V数学中国社区 不走平凡路

海之舟

好高深啊~~~~~~