QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2982|回复: 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];
    5 p  P+ _! o" y6 V3 ]dayhistory=day(1:16);%取其中三十天作为历史数据样本
    # ?* \- D3 C1 c) y$ Udayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的, ^7 |, B6 T8 x4 A. E2 d
    dayday=day(1:12);% 取其中的前25天
    8 {( w9 W6 B* E6 M- t  v# Mdaypost=day(5:16);%取其中的随后25天
    ! G' H: K" H7 N/ e! b* N6 {3 np=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本: L& t7 t6 c$ i9 p$ G; U
    t=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量
    ) {: P- I6 y, wdaylast=day(13:16);
    - H% x( T% R. j- R. w! V3 L; B# Bh3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本  T: E9 A  i) O4 }
    r=5:16;
    " G8 g8 n1 H' b1 _( C8 T, grr=reshape(r,4,3);' u1 f2 a" n+ S  B+ S" O1 ^1 @
    %%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%/ _' @% ~* `) a6 l  @. q% H7 c
    net=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');
    6 [" K, F6 L9 R# y% u, A( h" {% ey1=sim(net,p);
    * j0 I2 Z. r0 t% ?4 ]% 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量
    - L$ x9 D: @' S& x, X2 f4 `% 两层的传递函数均为purelin
    $ j6 I" ^) [. m# b$ v" B  I% 训练函数为trainlm
    9 E, Q: s3 [2 n- \0 j% 所训练的网络大小为[5,5]6 S. n# q1 Q, i* x
    % 仿真训练前的网络+ S! @# H8 o" O3 A2 G* m

    ( w1 Y0 a1 }9 H! y%%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%* \& `$ a, ?. j% h9 T
    % network parameters:
    , |8 C1 t- z$ m. s: s%   epochs--epochs of the train1 y( |/ a$ n1 T
    %   goal--errors goal of the network
    . y8 d+ j& D/ b1 Q+ E%   lr--learning rate
    9 v/ ~' c9 v% C! N, X: m%   shows--epochs between the displays
    ! h. M6 G7 K: K4 x%   time--Maximum time to train in seconds
    5 B- `1 P7 u% T5 }8 S: P3 Q. snet.trainParam.epochs=200000;  % 训练次数
    + ]4 F9 M' S, g7 x) G5 RnettrainParam.goal=0.0001;  % 误差期望值
    - e- {- D5 @& p: c- ]% returns of the train:
    ( z: @7 r  ^4 |* h/ d  S, @2 k) c%   net--New network0 }9 G; |/ l5 w* w
    %    tr--Training record (epoch and perf).$ z" t  Y( _' }' W% v% `; ]' |+ g. O
    %     Y--Network outputs.- D- R: K, r, M; `- j' Q
    %     E--Network errors.
    , _0 B+ O8 x7 y6 C/ o. o; y. P& c[net,tr,Y,E]=train(net,p,t); 6 U) n) [. |  W" E. T! ]0 D
    %%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%
    3 a5 m3 `, N8 s  J% input the testing points here %
    2 F/ [, n$ `% f" K# M; Etitle(\'神经网络训练结果\');
    3 L4 {* w/ q& @xlabel(\'时间(天)\');  a9 A) q8 e. T' {" Q# D& [6 q
    ylabel(\'仿真输出结果\');0 H8 Q/ {0 g4 O* _1 v  R7 {" m2 \
    legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');) u6 P1 f8 ~5 M
    %%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    $ R4 p( _, K$ s0 K: i! h, Yx=1:4;, G, l8 Z2 {' H$ j1 f% F$ F2 U6 u
    y=1:4;. }/ {5 u6 w7 X$ g
    y21=sim(net,p);- o( l' K$ [3 \: h1 M$ p& R
    y2=reshape(y21,1,12);
    ) W7 Z6 k5 Y3 H/ A0 M9 y# Vclf& T% L' f  e" {0 ~) p
    plot(r,y2,\'b-^\')
    7 {7 Z% A8 P8 P2 c, R7 }5 b& Nhold on8 X- g0 R" G" u$ U
    plot(1:20,day,\'r-*\')
    0 ~. B* O( r1 D4 \6 L1 N% i: |# s%%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%/ h% t5 r; h# h* u1 m
    y3=sim(net,h3);0 v0 M% C$ M3 k/ S
    plot(17:2:20,y3,\'-*\')
    $ K  U3 \# K3 \- U  J( jhold on* }9 r2 O7 ~5 H
    title(\'神经网络训练结果\');: v: K& l3 A5 c8 J! a+ O
    xlabel(\'时间(天)\');2 ?1 i7 P/ [$ M2 |
    ylabel(\'仿真输出结果\');
    9 A* y- ^$ \/ u( E9 H, {# Xlegend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    . @- W+ I* O6 N+ Z/ A, m6 N%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    1 [. h# L& T7 X  f* a# ?$ Nx=1:4;8 `, r, F' N- S/ S& O! a
    y=1:4;
    7 b4 O( b" |$ V* c$ z0 vplot3(x(1:3),y(1:3),E(x(1:3),y(1:3))): F. u  e8 G8 W& y* l! U
    运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢

    & U7 U; K5 t* D7 t9 e$ p" b  i3 n

    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-4-10 19:17 , Processed in 1.831826 second(s), 49 queries .

    回顶部