QQ登录

只需要一步,快速开始

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

[代码资源] 请问大神这个代码错在哪里,应该如何修改.以及回答一个简单的题目,谢谢。

[复制链接]
字体大小: 正常 放大
昌辉9        

3

主题

13

听众

72

积分

升级  70.53%

  • TA的每日心情
    开心
    2016-5-23 22:22
  • 签到天数: 54 天

    [LV.5]常住居民I

    自我介绍
    初学者

    群组: Matlab讨论组

    跳转到指定楼层
    1#
    发表于 2014-9-3 16:19 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    20体力
    day=[8.6187,8.3507,8.3142,8.2898,8.2791,8.2796,8.2784,8.2770,8.2770,8.2774,8.2780,8.1013,7.8087,7.3872,6.85,6.81,6.622,6.61,6.25,6.07];
    % U" D) D# T2 W9 U: b9 `dayhistory=day(1:16);%取其中三十天作为历史数据样本5 L6 S# V9 k2 r! ]' N" E
    dayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的
    % V- V: q) n2 ndayday=day(1:12);% 取其中的前25天6 y. a0 ^9 M& W4 M" f! ^' m" o
    daypost=day(5:16);%取其中的随后25天4 R# s1 k" Z" u, ~+ m7 S7 ~" V
    p=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本
    # B+ c* q  F+ M1 `/ _t=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量; e" v  K3 z4 P/ ^# y6 s7 G
    daylast=day(13:16);; s$ q9 F! V0 b- Q$ H  N
    h3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本
    6 V' R1 U5 q' \2 t+ ?r=5:16;
    - X: ^7 M4 M1 j8 v* j- M1 d+ t! [rr=reshape(r,4,3);
    8 ~* e% i5 t8 X! V%%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%
    9 ], ~/ ^' R% D, L! @. {. N  N+ Inet=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');: I+ L# n' A# ]* a$ v* b
    y1=sim(net,p);
    # `, U8 ?) d( N5 ?. d3 g% 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量
    ! ?$ J, W- M. X/ i+ A1 Z% 两层的传递函数均为purelin
    8 ?9 E$ C# G5 u7 H  b0 G$ e% 训练函数为trainlm, i0 o5 p( r# Y8 D1 j3 j
    % 所训练的网络大小为[5,5]
      u6 Z# }9 h) Y! X' L% 仿真训练前的网络  `5 H1 @6 \0 I
    % l6 t- T$ p' i8 o
    %%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%  @% e7 O' V5 u1 h( q
    % network parameters:4 H; `. K8 c1 Y) Y' z0 T$ q1 O* f. S
    %   epochs--epochs of the train
    . J# |  k- @! v' x1 H0 z%   goal--errors goal of the network
    1 q# Q4 H$ K' s' E; M%   lr--learning rate
    : l4 b+ z/ k0 L( ?( W* |%   shows--epochs between the displays
    $ |/ m2 @9 j# |, ~%   time--Maximum time to train in seconds
    $ k! Y$ F& D8 w- I5 u! d- y" Dnet.trainParam.epochs=200000;  % 训练次数
    , E3 R9 _: a, p! m, y# l' pnettrainParam.goal=0.0001;  % 误差期望值" P- f: d5 s, F
    % returns of the train:
    1 @! b( A5 R( y8 o+ p%   net--New network( g1 M' z4 G, o) S0 @, y: q( Z
    %    tr--Training record (epoch and perf).2 B7 _6 L6 J! D1 R+ r
    %     Y--Network outputs.- X- E: z4 E+ P4 n. Q: f
    %     E--Network errors.3 n) D# G8 C6 y% ]
    [net,tr,Y,E]=train(net,p,t);   g0 E! q; v. ?/ Q, L5 K
    %%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%7 H( A9 h2 F' _+ C0 F3 W
    % input the testing points here %3 s0 L( t4 F3 x/ X: R# o. x! r
    title(\'神经网络训练结果\');  A: G' D0 F( g0 ^1 K; b
    xlabel(\'时间(天)\');8 S' ^; F* Z) ^1 G+ f
    ylabel(\'仿真输出结果\');
    * g% P9 e- f; o" P# A, Xlegend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    - p+ D# }# P( d9 B$ `%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%+ p4 B, B9 F% L$ p4 G
    x=1:4;& s7 t' g% i0 i2 |+ d
    y=1:4;, S5 r  o5 Q8 Z+ [% A
    y21=sim(net,p);
    6 [3 k  ?+ m2 Y) i( {1 U8 _y2=reshape(y21,1,12);
    2 `9 e" [. ]0 I' e7 b* yclf
    ' S3 }9 M: h$ a4 F+ pplot(r,y2,\'b-^\')
    - O7 T! C: E$ chold on' p0 F% M- h# p* ~9 ?$ _; f
    plot(1:20,day,\'r-*\')
    9 ]" I' ]! l  ~2 i%%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%) S+ c. G8 y7 e' _. S
    y3=sim(net,h3);
    % F: e5 B6 ?. jplot(17:2:20,y3,\'-*\')' W% |3 E' B! T/ U8 L. V
    hold on/ _0 Z( |/ y  q8 I' C' U6 |
    title(\'神经网络训练结果\');
    5 G" k/ b$ l5 D; Qxlabel(\'时间(天)\');0 w$ k' X% s- p% o/ r
    ylabel(\'仿真输出结果\');
    ! ?5 T  u# n" c8 ~/ E- V0 p8 Jlegend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    8 S. I) Y; a5 K%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    " h/ C+ X3 |: x& M+ wx=1:4;+ m7 i8 d0 ~/ C% p4 A+ D' ?% O5 Q' f( o
    y=1:4;
    : Q1 ^' J6 N5 O' n$ Rplot3(x(1:3),y(1:3),E(x(1:3),y(1:3)))/ \0 V8 E. ?& A  C
    运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢
      u4 w0 |" j* W8 e

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-10-10 04:40 , Processed in 2.238882 second(s), 49 queries .

    回顶部