QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5454|回复: 9
打印 上一主题 下一主题

紧急求助:关于灰色神经网络GNNM(1,1)的编程问题

[复制链接]
字体大小: 正常 放大
yqxxqy        

21

主题

4

听众

321

积分

升级  7%

  • TA的每日心情
    难过
    2016-4-16 21:09
  • 签到天数: 1 天

    [LV.1]初来乍到

    新人进步奖

    跳转到指定楼层
    1#
    发表于 2011-3-19 20:56 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    小弟正在编写GNNM(1,1)的程序,但是就是调试不出来,心里那个急啊~~~
    ( U. [" S0 W" J! \; a/ c; d' x有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~- t+ \8 k+ \. l
    有没有现成的程序供参考参考啊~~~
    7 C1 W+ n  b0 y1 ], }
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    zjqylcy 实名认证       

    7

    主题

    5

    听众

    1205

    积分

    宇宙无敌霹雳大头衔

    升级  20.5%

  • TA的每日心情
    奋斗
    2013-5-18 15:14
  • 签到天数: 101 天

    [LV.6]常住居民II

    自我介绍
    男~~

    新人进步奖

    群组杭州联盟

    群组小草的客厅

    群组C 语言讨论组

    群组数学建模

    群组数学趣味、游戏、IQ等

    现在的人怎么都那么懒了~~~诶。。。
    # l4 b* p2 [5 c7 @  Y) }这个是用来做人口预测的,你参看吧~~~
    5 F6 Y! [7 |6 v7 H* f& ^1 Z4 q
    GM( 1, 1) 的matlab程序
    gm.m:
    function f=gm(x0,m)& R$ _0 \5 A7 h" A& d, {4 `
    %定义为函数gm(x)
    n=length(x0);
    x1=zeros(1,n);
    x1(1)=x0(1);
    for i=2:n" l- e8 S$ E9 K- H2 K  x/ U# K' |
    %计算累加序列x1

    / i$ Q5 I) i9 Z: E$ T( o$ ux1(i)=x1(i-1)+x0(i);
    end
    i=2:n;9 u( I" d8 r: \" r
    %对原始数列平行移位并负值给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;
    $ e; m+ ^8 G4 A/ o% h+ A%计算参数a,u矩阵
    i=1:m;5 |- a) [* N$ L; D; t
    %计算预测累加数列的值
    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;$ u6 u# M) }- R/ |! f1 h( @1 U4 h0 d

    3 m  x# R0 }" ~9 z7 j/ X# t- e: q. `%还原数列的值
    yc(i+1)=ago(i+1)-ago(i);
    i=2:n;
    error(i)=yc(i)-x0(i);
    " ~. y6 I7 {: _%计算残差值
    yc(1)=ago(1);
    i=1:m-1;. I+ e. f. ?, L8 m
    %修正还原数列的值
    yc(i+1)=ago(i+1)-ago(i);
    c=std(error)/std(x0);3 i4 j. A/ r6 o3 ?
    %计算后验差比
    relerror=abs((error)-mean(error)*ones(size(error)));
    [nrow,ncol]=size(relerror);
    p=0;
    for i=2:ncol

    + W& X, U/ M# j! |% Lif( K) r9 h, O2 ]
    relerror(1,i)<0.6745*std(x0)

    2 l% h2 H3 r4 O& C/ }6 l( t) u$ tp=p+1;
    % T1 r7 h) c1 V6 v* b8 Q5 Q' D
    end
    end
    p=p/(n-1);
    % e& r1 b: C  h; T/ J" |. U% k  H* K5 }
    w1=min(abs(error));
    w2=max(abs(error));
    i=1:n;$ j' Z8 g1 |* M' k7 _/ Z4 g, r0 Z

    : L9 J7 y- `$ K% \: [%计算关联度
    w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
    w=sum(w)/(n-1);
    au
    1 s2 M" B! U, h%输出参数a,u的值
    ago;
    8 H3 n1 K6 c: @* v# T/ h" T%输出累加数列ago的值
    x0;* w3 s  A; N, o
    %输出原始序列值
    f=yc;
    7 t! h' q% i/ R- r: ?
    3 S3 s& v. g* z# r%输出预测的值
    error;7 Z4 |' @1 P$ ~0 T& }9 E
    %输出残差的值
    c;0 e8 |; H$ f; u+ u
    %输出后验差比的值
    p;
    5 R% i1 a* M8 V5 }, C  `%输出小误差概率的值
    w* l' F/ J6 ~# [$ z
    %输出关联度
    end
    ) W5 Y+ U2 X2 g9 V6 {; I7 m2 u; K4 b/ R. z% D) U) C' R
    ; H# ]6 B% q; y5 F1 @. G: D
    回复

    使用道具 举报

    diliwang 实名认证       

    0

    主题

    4

    听众

    64

    积分

    升级  62.11%

  • TA的每日心情
    开心
    2012-2-5 18:47
  • 签到天数: 10 天

    [LV.3]偶尔看看II

    回复

    使用道具 举报

    yqxxqy        

    21

    主题

    4

    听众

    321

    积分

    升级  7%

  • TA的每日心情
    难过
    2016-4-16 21:09
  • 签到天数: 1 天

    [LV.1]初来乍到

    新人进步奖

    回复 zjqylcy 的帖子( T9 b- j+ b& J0 w' |) o7 P: o
    , U2 P3 i3 V. I4 F8 I
    我是想求助关于灰色神经网络GNNM(1,1)的。。。不过还是谢了~~~
    回复

    使用道具 举报

    17

    主题

    3

    听众

    2216

    积分

  • TA的每日心情
    开心
    2012-1-30 23:29
  • 签到天数: 39 天

    [LV.5]常住居民I

    群组小草的客厅

    群组数学建模

    群组Matlab讨论组

    群组LINGO

    群组中南民族大学

    回复

    使用道具 举报

    linmatsas 实名认证       

    53

    主题

    13

    听众

    3591

    积分

    逍遥游

  • TA的每日心情
    奋斗
    2014-12-2 09:53
  • 签到天数: 54 天

    [LV.5]常住居民I

    自我介绍
    额。。。。世界上最讨厌的事情就是自我介绍。。。

    邮箱绑定达人 新人进步奖 发帖功臣 最具活力勋章

    群组Matlab讨论组

    群组数学建模

    群组小草的客厅

    群组2012数学一考研交流

    群组C 语言讨论组

    回复

    使用道具 举报

    jshzncd 实名认证       

    1

    主题

    3

    听众

    27

    积分

    升级  23.16%

  • TA的每日心情
    奋斗
    2014-8-23 05:25
  • 签到天数: 1 天

    [LV.1]初来乍到

    群组Matlab讨论组

    群组数学建摸协会

    其实我一直觉得楼主的品味不错!呵呵!- S7 G. E4 m" h0 F  a, v, S
    数学中国社区 不走平凡路
    回复

    使用道具 举报

    海之舟        

    0

    主题

    3

    听众

    118

    积分

    升级  9%

  • TA的每日心情
    开心
    2012-2-17 00:39
  • 签到天数: 38 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    9#
    无效楼层,该帖已经被删除
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-17 02:20 , Processed in 0.421885 second(s), 97 queries .

    回顶部