QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3071|回复: 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];
    - Y7 S) M4 q. S# Cdayhistory=day(1:16);%取其中三十天作为历史数据样本6 }/ v2 _6 J5 b' A
    dayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的
    ) v7 @, n: Y4 \/ y3 g4 h/ adayday=day(1:12);% 取其中的前25天1 F8 X0 z: Q, K% i. h
    daypost=day(5:16);%取其中的随后25天% q, s* y- j9 [! x
    p=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本" E6 k7 v& U( Q2 O! e
    t=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量
    + H4 u7 I( a. f! S7 E9 F: |9 D- Edaylast=day(13:16);
      [5 l! y0 [3 Ah3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本3 X% ]  l( v: F' A4 U, D
    r=5:16;
    * i. Q8 z& G7 E5 q' Trr=reshape(r,4,3);
    + V" T% C  z/ x/ f4 m0 x%%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%
    ; d: B4 A4 {  {, ~7 enet=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');4 l, X2 n, h% Z% y7 U+ s
    y1=sim(net,p);
    4 i6 E( |/ Z3 h! I% 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量
    7 t- L9 l3 F( ]! t% @# l+ h% 两层的传递函数均为purelin
    - Y, `0 C+ I4 U6 X% 训练函数为trainlm
    % t4 W- d( E) }( e% 所训练的网络大小为[5,5]0 P) K0 M6 z, T! \
    % 仿真训练前的网络
    3 P" I& K% w) a) d2 N$ q/ d# }6 U+ V: Y& B' r
    %%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%
    - e! H! O3 I7 Y5 z- O$ A9 [% network parameters:
    8 z2 @- l, f2 D& \, Q%   epochs--epochs of the train0 A* L7 `' ]4 m
    %   goal--errors goal of the network+ x2 @8 j! t) c5 G9 ~( e
    %   lr--learning rate
    # |# f5 A6 v. W7 K5 {%   shows--epochs between the displays4 `7 U$ l. i! t
    %   time--Maximum time to train in seconds1 a  u: A3 d# Q: I0 [7 P& a4 r- K
    net.trainParam.epochs=200000;  % 训练次数7 i8 a+ J8 `: m0 @1 L6 V0 X/ v
    nettrainParam.goal=0.0001;  % 误差期望值
    8 i7 n  ?' E9 K% returns of the train:
    : q9 x- X8 H# s6 v, Y# }%   net--New network
    & Y: b# _& Y# a%    tr--Training record (epoch and perf).5 _# o' n1 _! k- m/ G8 ~
    %     Y--Network outputs.% q( n$ F! g* [- o0 \1 N
    %     E--Network errors.# P) {. C- c5 D4 O# w& X
    [net,tr,Y,E]=train(net,p,t);
    ' s$ C& R  t6 W0 Q& g( ]%%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%2 |. }- H. N. X
    % input the testing points here %
    * T$ \0 @& `! }title(\'神经网络训练结果\');
    : _9 z9 O  B) q) u7 lxlabel(\'时间(天)\');
    : H) w5 p; i. C% V4 Kylabel(\'仿真输出结果\');
    $ j* |' O0 Q2 `legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');1 N$ ^3 f* @1 j9 ^. t. K
    %%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%# X6 l3 {* Z% f, [! q$ l( ~
    x=1:4;
    ! j- |+ |5 l! b: m/ W# Y8 S3 ^. ey=1:4;; M, S$ O- V3 x6 v  R" k
    y21=sim(net,p);) U8 i4 y' I  C0 l3 ?9 Y
    y2=reshape(y21,1,12);
    " ^+ H8 B5 R) Q, F0 cclf/ P& c/ K/ V0 g& D
    plot(r,y2,\'b-^\')
    ) L7 j! G5 ^* V7 T- a: Whold on  ~* @' E" W/ `: ?& o. g
    plot(1:20,day,\'r-*\')$ a  }0 ?# o  j% D$ P. P5 P5 F
    %%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%
    7 c. ^, _6 y( |# _8 \y3=sim(net,h3);1 H# ]6 |: ~9 W; x& Q6 j1 D( B
    plot(17:2:20,y3,\'-*\')) c8 h$ _1 O  E" m
    hold on
    1 O* }4 _4 z! R! btitle(\'神经网络训练结果\');* g# q/ F  G. o2 |8 F% a
    xlabel(\'时间(天)\');/ b. Q! u! e* w9 C
    ylabel(\'仿真输出结果\');4 X8 }' v% o1 }/ R
    legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    2 D. f; B9 z! ~% l) S& [%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%5 w. w5 g) ]# t1 ?0 r. W0 x6 N
    x=1:4;( s6 B" ^0 F) Y
    y=1:4;
    9 ^/ {6 R( X0 v; x" _. v+ V2 d4 t" ~+ Zplot3(x(1:3),y(1:3),E(x(1:3),y(1:3)))( f6 |; j3 ~: N4 y( M' Q
    运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢
    6 k; n- i$ F: r: S

    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-8-24 19:41 , Processed in 0.768152 second(s), 49 queries .

    回顶部