数学建模社区-数学中国
标题: 紧急求助:关于灰色神经网络GNNM(1,1)的编程问题 [打印本页]
作者: yqxxqy 时间: 2011-3-19 20:56
标题: 紧急求助:关于灰色神经网络GNNM(1,1)的编程问题
小弟正在编写GNNM(1,1)的程序,但是就是调试不出来,心里那个急啊~~~: _6 a# c. e( r0 o4 _' n
有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~# r) x; s! l; ]' L0 d* V( g! v6 U
有没有现成的程序供参考参考啊~~~7 H, e8 Z, g5 S
作者: zjqylcy 时间: 2011-3-19 22:09
现在的人怎么都那么懒了~~~诶。。。
^4 L: I" ~# B: _这个是用来做人口预测的,你参看吧~~~9 F( V9 f9 k, i, m/ r. L
GM( 1, 1) 的matlab程序
gm.m:
function f=gm(x0,m)- K4 [5 Z4 q9 p
%定义为函数gm(x)
n=length(x0);
x1=zeros(1,n);
x1(1)=x0(1);
for i=2:n
% a% }1 j% t y3 f' M# _, `%计算累加序列x1
% Q, X4 s4 |3 K) Ox1(i)=x1(i-1)+x0(i);
end
i=2:n;
E6 b5 e% c" X; E9 p8 g%对原始数列平行移位并负值给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& \1 d0 c: f, m) e% v/ {
%计算参数a,u矩阵
i=1:m;5 E. o' {4 u4 ]1 X8 n' R( G
%计算预测累加数列的值
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;
) R# ^* m6 N* E, O$ g0 _4 g
+ c6 j2 ^+ s7 I) Z: Y%还原数列的值
yc(i+1)=ago(i+1)-ago(i);
i=2:n;
error(i)=yc(i)-x0(i);- {$ d9 S% t% `) y0 Y8 H
%计算残差值
yc(1)=ago(1);
i=1:m-1;
9 @- w: R P2 q- N# J5 r: p; G; \%修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);
c=std(error)/std(x0);
; Z9 O+ c; w2 [! T: `+ I- s: d4 P; ?3 @8 g%计算后验差比
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);
p=0;
for i=2:ncol
! V' H0 Q& Y. W- e0 c* o$ _if
' N5 o& `+ ?5 @relerror(1,i)<0.6745*std(x0)
6 R! Z; c$ Q3 D- T2 x
p=p+1;
/ ^+ }$ Y8 }) |
end
end
p=p/(n-1);& [7 W5 a2 a9 N# X" k; n
0 n) p& s, ^+ k4 S5 f( l0 T
w1=min(abs(error));
w2=max(abs(error));
i=1:n;2 @5 f) K/ U- Y
) T1 V) W! D( W7 w% ~9 G%计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au
7 [ Q( f7 x2 v8 i' h2 Q t8 J%输出参数a,u的值
ago;5 H; z7 }# t. M0 z" G9 b; a# v
%输出累加数列ago的值
x0;* q3 a4 O: _3 ]5 ?
%输出原始序列值
f=yc;% \2 o7 h: u6 W# ~* V F
+ X# n) ]- ^5 e4 a, d# [
%输出预测的值
error;& @4 M( I! N) u# c2 r
%输出残差的值
c;$ n e: p' }' F) s7 V+ R( w
%输出后验差比的值
p;
5 w5 \4 m2 W8 X# |- j* s) D%输出小误差概率的值
w
) @5 H* ~3 h& X0 P) o%输出关联度
end) V& Z4 `$ D' G
! J& u. M. Z1 m U0 P
& g, ?. W3 X) K6 y; r: ?: a4 [
作者: diliwang 时间: 2011-3-19 22:09
同求,挺感兴趣
作者: yqxxqy 时间: 2011-3-20 00:59
回复 zjqylcy 的帖子
3 n b4 J/ _$ o0 d3 @8 K+ ? o6 r; Z: W
我是想求助关于灰色神经网络GNNM(1,1)的。。。不过还是谢了~~~
作者: gaoshanliu水 时间: 2011-3-20 20:07



作者: linmatsas 时间: 2011-3-20 23:25
好高深呀。。。不懂不懂、、、、、
作者: jshzncd 时间: 2011-5-2 02:20
其实我一直觉得楼主的品味不错!呵呵!& E9 i" h% S4 E- _( x! ?6 Z
数学中国社区 不走平凡路
作者: 海之舟 时间: 2011-12-18 17:21
好高深啊~~~~~~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) |
Powered by Discuz! X2.5 |