QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3025|回复: 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];. B4 u7 O6 _: F, Q# t
    dayhistory=day(1:16);%取其中三十天作为历史数据样本
    9 |: B, a- B7 t! zdayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的, L( s) M7 `+ G  ]! N6 M9 w2 J
    dayday=day(1:12);% 取其中的前25天
    % S4 N+ l6 G9 i% w1 gdaypost=day(5:16);%取其中的随后25天$ L5 P$ u. J* n2 C  Z
    p=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本
    7 g6 i/ Y  s1 P7 V4 i* Y1 Ft=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量
    & B8 D+ N  C( G4 Q6 q# fdaylast=day(13:16);0 J8 @! r+ Z5 H5 y0 h$ l' j
    h3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本; t+ Y" {) ^9 m! x+ ~
    r=5:16;0 P6 h. I" c" {7 p* R
    rr=reshape(r,4,3);
    , N% v7 T7 G# U/ V7 X( _%%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%
    9 `9 @6 \+ z, H4 Q3 B- znet=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');1 r/ M% z5 p5 ]5 p8 v3 V3 T/ q
    y1=sim(net,p);
    ! S4 v: a: d) {- s5 t$ \/ p% 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量8 {0 |# x- F( ~$ ^& c) K8 P
    % 两层的传递函数均为purelin
    4 N4 }( A8 O) s( |! ]7 e4 i% 训练函数为trainlm
    3 j6 b* o8 d6 D% _5 x% 所训练的网络大小为[5,5]+ V: p& t7 i: V* ?) f$ k8 w
    % 仿真训练前的网络8 }# l& V% f5 I6 B5 H- u
    ! r$ G6 n& b/ }
    %%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%
    # x$ i) y# t3 @; [, p4 m% network parameters:' l" N  n; I3 U+ A# J- P+ B
    %   epochs--epochs of the train
    & F+ e" ^% W8 P/ B; o) f4 J( x%   goal--errors goal of the network( W( j! i5 d+ {' J3 a, _% y- K9 Z
    %   lr--learning rate+ C& U+ @; \5 {( y% P4 p  b
    %   shows--epochs between the displays
    $ z1 g* E8 u# F! H# I" Q- W. T%   time--Maximum time to train in seconds
    & z1 Y$ A- T$ [$ pnet.trainParam.epochs=200000;  % 训练次数* A0 u: r7 |+ c7 T7 @
    nettrainParam.goal=0.0001;  % 误差期望值
    9 \- c7 l: W, {$ q7 w% returns of the train:* M: N5 d( Q! P& X
    %   net--New network
    $ w, a/ l# E  R2 s& _%    tr--Training record (epoch and perf).
    ' J8 j! N, F/ c7 r- Z2 C$ w%     Y--Network outputs.
    . S' ^: l( D5 o+ N%     E--Network errors.
    . O( Q( w' u% i8 F1 c[net,tr,Y,E]=train(net,p,t);
    5 N- w5 ~  I5 _' ?& X%%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%/ J1 Q! J( {$ \1 P: {4 J
    % input the testing points here %
    7 A* A, u/ S  ]( Ttitle(\'神经网络训练结果\');+ }  p) |& S, R- X
    xlabel(\'时间(天)\');
    & w- t, W+ J* A- h0 ^ylabel(\'仿真输出结果\');/ G, J4 ]' T. ?! r+ f5 z8 h/ m
    legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    3 G* l' b6 B3 \# I% n8 ?%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    8 G& H, P9 R, Yx=1:4;
    ( b9 s/ B6 a. T6 [* @: s) |( by=1:4;" a5 m5 W8 O( Q2 y7 u9 ~& a( M
    y21=sim(net,p);
    $ Y) G. U  V2 \" _y2=reshape(y21,1,12);
    2 \4 n8 P2 [5 y4 B6 sclf
    1 ]# H6 {2 D, F( x5 Oplot(r,y2,\'b-^\')
      I3 d& d4 [& M4 t; K2 L$ `hold on
    6 i9 E1 ^6 W; [" b- A% [plot(1:20,day,\'r-*\')3 u( J0 ^; D; M
    %%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%3 t3 t% x& a; N1 e% P
    y3=sim(net,h3);6 p+ i7 t6 m" z2 [
    plot(17:2:20,y3,\'-*\')6 M/ R$ G1 G+ {2 t& `
    hold on
    $ j1 m* i3 u8 }( C, `title(\'神经网络训练结果\');
    5 }7 P! G1 ]* Gxlabel(\'时间(天)\');8 R2 S/ e2 {, v0 d: \8 J- F
    ylabel(\'仿真输出结果\');
    * u( E/ ?" S- Y8 E/ Ilegend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');7 Z) N, R3 n8 V: `, Q
    %%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%- }" M1 T, d( D7 X
    x=1:4;7 N  V. M8 U4 a- B% P( A
    y=1:4;
    7 l4 _9 G. M9 D) M) k$ ]* r  ~/ T2 c1 nplot3(x(1:3),y(1:3),E(x(1:3),y(1:3)))
    . A5 M' k4 o+ F0 `运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢
    0 S. p  r( A% X* Z( w

    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-6-11 09:11 , Processed in 0.337757 second(s), 49 queries .

    回顶部