数学建模社区-数学中国

标题: 紧急求助:关于灰色神经网络GNNM(1,1)的编程问题 [打印本页]

作者: yqxxqy    时间: 2011-3-19 20:56
标题: 紧急求助:关于灰色神经网络GNNM(1,1)的编程问题
小弟正在编写GNNM(1,1)的程序,但是就是调试不出来,心里那个急啊~~~
0 z9 q) Q# w+ }有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~: t' W: c! _2 z# v5 c
有没有现成的程序供参考参考啊~~~6 q, g! u+ @: b7 H

作者: zjqylcy    时间: 2011-3-19 22:09
现在的人怎么都那么懒了~~~诶。。。
8 `- R5 [6 ^# v- Q( ~/ X( x2 [这个是用来做人口预测的,你参看吧~~~
% d: J0 q( p- u; y
GM( 1, 1) 的matlab程序
gm.m:
function f=gm(x0,m)
3 @7 p* H# |. C: L%定义为函数gm(x)
n=length(x0);
x1=zeros(1,n);
x1(1)=x0(1);
for i=2:n. R$ @" a% {7 V; P( p8 h
%计算累加序列x1
, z9 Q7 W+ d% \" X  L8 u
x1(i)=x1(i-1)+x0(i);
end
i=2:n;
2 D. L4 `% D9 E; Y%对原始数列平行移位并负值给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;) ?9 A+ }3 b8 E- A. _" u5 M" x+ v
%计算参数a,u矩阵
i=1:m;: ~( r; B; \0 j6 R
%计算预测累加数列的值
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;# x! ~* M5 J- [/ T0 U, a0 b, i4 `

) h' s" Y8 w) \# L0 ^  o) Z7 F" M. F%还原数列的值
yc(i+1)=ago(i+1)-ago(i);
i=2:n;
error(i)=yc(i)-x0(i);
/ @( s0 R1 j# [% u( f) E%计算残差值
yc(1)=ago(1);
i=1:m-1;  d! r+ o' G! u- N
%修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);
c=std(error)/std(x0);6 L# f" |0 I9 w  m: _% Q. q* Q
%计算后验差比
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);
p=0;
for i=2:ncol
( c  B( T& G; o
if
, o% s/ l$ x% I) nrelerror(1,i)<0.6745*std(x0)

5 [. }' B, W- e) Zp=p+1;

4 `8 D. ?: {" ^& rend
end
p=p/(n-1);
5 d0 K. M8 w4 L# I" |2 e# o
. C9 A  ~, N4 c/ \0 P
w1=min(abs(error));
w2=max(abs(error));
i=1:n;
; `- J. j8 }9 X+ }
5 e4 V( z+ S  @- B' _2 \%计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au
! B3 ?* Z5 i# n%输出参数a,u的值
ago;6 u6 ]6 u+ y' a$ v2 F( l
%输出累加数列ago的值
x0;
# |3 Y! u6 l, P%输出原始序列值
f=yc;0 `, w8 b8 ?: B; I1 R

; }6 J" ]4 x; Z%输出预测的值
error;
; v: G, m& G$ \%输出残差的值
c;
( _& T3 E2 @  j( X* A%输出后验差比的值
p;& _1 C# ?" D& f) ~" }
%输出小误差概率的值
w
/ Z' |6 Z1 O$ o9 v, w%输出关联度
end
6 D, p/ _8 G. q7 e
1 P1 W' ]+ V7 Y3 O$ L

' j! X0 s% s, Y/ ~9 O; H2 l
作者: diliwang    时间: 2011-3-19 22:09
同求,挺感兴趣
作者: yqxxqy    时间: 2011-3-20 00:59
回复 zjqylcy 的帖子, C" C5 |: B5 L6 S4 O+ ~
  B  @+ ]1 f5 s/ P
我是想求助关于灰色神经网络GNNM(1,1)的。。。不过还是谢了~~~
作者: gaoshanliu水    时间: 2011-3-20 20:07

作者: linmatsas    时间: 2011-3-20 23:25
好高深呀。。。不懂不懂、、、、、
作者: jshzncd    时间: 2011-5-2 02:20
其实我一直觉得楼主的品味不错!呵呵!
5 t5 D0 a1 H3 L2 u0 w9 P数学中国社区 不走平凡路
作者: 海之舟    时间: 2011-12-18 17:21
好高深啊~~~~~~




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5