QQ登录

只需要一步,快速开始

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

紧急求助:关于灰色神经网络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)的程序,但是就是调试不出来,心里那个急啊~~~8 [: J& K& \, {0 J9 @( O
    有没有资深大侠给点帮助啊~~~~~~~~~~~~~~~~~~~~~~~~~
    2 ~! `! ^& k( z' {有没有现成的程序供参考参考啊~~~3 }& w3 C  G$ E( B
    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等

    现在的人怎么都那么懒了~~~诶。。。0 S9 \. C' z9 `) L2 U9 M
    这个是用来做人口预测的,你参看吧~~~
    3 Q- m0 `3 f& p
    GM( 1, 1) 的matlab程序
    gm.m:
    function f=gm(x0,m)" o6 Q$ o0 U6 C7 i5 ^- |
    %定义为函数gm(x)
    n=length(x0);
    x1=zeros(1,n);
    x1(1)=x0(1);
    for i=2:n
    ) z# W: i  f) M6 n# s%计算累加序列x1
    $ T! ~/ n% b% g# y* c
    x1(i)=x1(i-1)+x0(i);
    end
    i=2:n;
    4 @5 ]. b* I5 o8 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;
    - A8 z& m. c* \2 P%计算参数a,u矩阵
    i=1:m;
    % r; a: L+ p) |% 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;
    . A3 z0 I4 Z6 H+ Z
    + G% w+ Z! [9 ~$ P" w7 w7 C%还原数列的值
    yc(i+1)=ago(i+1)-ago(i);
    i=2:n;
    error(i)=yc(i)-x0(i);6 t" c+ C! d# O$ i
    %计算残差值
    yc(1)=ago(1);
    i=1:m-1;
    % G; Z& J! k$ I( k+ v# p. u%修正还原数列的值
    yc(i+1)=ago(i+1)-ago(i);
    c=std(error)/std(x0);/ r  `" M2 e8 J% |4 u+ P8 u
    %计算后验差比
    relerror=abs((error)-mean(error)*ones(size(error)));
    [nrow,ncol]=size(relerror);
    p=0;
    for i=2:ncol
    4 J# c* ^5 f$ C% O' a( S* _
    if
    # p) f" r- R2 u4 j+ J2 xrelerror(1,i)<0.6745*std(x0)
    3 A' @! u7 S, |, H% I" m: s
    p=p+1;

    $ I+ \; T- F  {% vend
    end
    p=p/(n-1);
    . C4 F3 A! Z: ^( j) t
    8 r- j" O5 D! p
    w1=min(abs(error));
    w2=max(abs(error));
    i=1:n;
    - T! i' A# ]) h$ O5 i7 }$ u0 `$ o9 J1 j# h+ `
    %计算关联度
    w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
    w=sum(w)/(n-1);
    au
    ! l7 s! T; h: u: {! X%输出参数a,u的值
    ago;1 a% ]: ]" b3 h  Z
    %输出累加数列ago的值
    x0;( ~" Z! J4 h( H, @3 b1 _8 F% L
    %输出原始序列值
    f=yc;# C' i5 ^/ t8 b8 m
    " u9 T0 L3 O- X. y: v  E$ V0 P
    %输出预测的值
    error;
    , @  A( |" o  r1 a%输出残差的值
    c;
    3 @* ^: D; d+ ^4 ~) v; y; m%输出后验差比的值
    p;
    4 M+ @' C, k# d* d' A%输出小误差概率的值
    w6 O; Q9 W! Q5 E/ [9 ^2 A! L
    %输出关联度
    end) V% G* V: u. {+ K

    5 {5 F# M9 X+ H- _
    , n$ s" U2 J& z' T" o( k
    回复

    使用道具 举报

    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 的帖子
    ' m# A% N  S0 n3 ^
    ( O1 t; f" M  G我是想求助关于灰色神经网络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讨论组

    群组数学建摸协会

    其实我一直觉得楼主的品味不错!呵呵!6 g  K# q3 k( F9 |: b- ^0 X
    数学中国社区 不走平凡路
    回复

    使用道具 举报

    海之舟        

    0

    主题

    3

    听众

    118

    积分

    升级  9%

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

    [LV.5]常住居民I

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-4-17 04:23 , Processed in 0.455914 second(s), 91 queries .

    回顶部