QQ登录

只需要一步,快速开始

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

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

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

3

主题

13

听众

72

积分

升级  70.53%

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

    [LV.5]常住居民I

    自我介绍
    初学者

    群组Matlab讨论组

    跳转到指定楼层
    #
    发表于 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];  w% M' ~3 c, |6 |% W
    dayhistory=day(1:16);%取其中三十天作为历史数据样本
    1 \6 T! `2 L# N7 K+ J- udayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的5 z! g# F8 `1 V0 _0 p* G
    dayday=day(1:12);% 取其中的前25天( T/ p2 P  o1 b3 u& {8 \
    daypost=day(5:16);%取其中的随后25天
    ( l% t$ O$ b& ]! P, R6 Dp=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本
    , @) T6 ~7 m! G, d: Ft=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量
    7 M. j: I) N) S5 U$ r$ ddaylast=day(13:16);; {) }* ?. P. z* z9 K
    h3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本# m! X8 n; R/ w, H+ Z% r  t
    r=5:16;: H: B7 W/ o8 h" {
    rr=reshape(r,4,3);
    , C; W) ^3 t+ D, H7 @# y5 _%%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%8 C( P* q  q! F4 O1 O9 j2 F
    net=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');
      i8 d  M( M8 ^3 |  S+ l* zy1=sim(net,p);, t1 K; g$ P5 Y8 T9 X
    % 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量7 _  S% z8 m  C5 u% o
    % 两层的传递函数均为purelin. P/ q/ ?% d2 b. F6 @0 {" J6 D  H; p
    % 训练函数为trainlm
    & J, U! R1 ~1 J9 n# |3 r% 所训练的网络大小为[5,5]$ j* H& J: l! {* ~! Y& w% b  q
    % 仿真训练前的网络
    5 E" E" o9 C4 k) D
    ) {/ a: Z8 D+ k9 w%%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%
    0 A6 m9 Y9 J9 |/ w$ _, f: K% network parameters:2 R9 q5 Q& f0 Y" m- s0 {8 f
    %   epochs--epochs of the train
    ; b& s+ V% ]5 _%   goal--errors goal of the network, v% [" P8 C. F6 z3 G: l8 v6 Z; Y$ t
    %   lr--learning rate4 h# C5 X+ A& Q6 X2 c
    %   shows--epochs between the displays
    2 L8 n' J" C% a) _2 X%   time--Maximum time to train in seconds3 |; e$ b9 M) O3 G
    net.trainParam.epochs=200000;  % 训练次数
    5 Z, G2 ?8 ?5 u+ j. Z9 inettrainParam.goal=0.0001;  % 误差期望值
    2 T" \( L% x/ h$ `  v% returns of the train:, V$ w  a4 |# |, z0 N
    %   net--New network6 u$ x6 z5 N2 f( Y. j; h
    %    tr--Training record (epoch and perf).
    7 t6 e9 V% B  A1 k" f+ _! p%     Y--Network outputs.* T% e6 T+ P# u% ?. o0 h+ t1 @
    %     E--Network errors.) F# [6 U9 C  B7 d0 v( D7 r
    [net,tr,Y,E]=train(net,p,t); % w$ C# {* \7 L3 j& r' s& h: `- O" ^# D
    %%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%
    5 f# A, L5 H1 y% |: Y+ I% input the testing points here %
    / ~9 K# C2 d- w; w3 ntitle(\'神经网络训练结果\');+ m7 _9 B9 o1 G1 _! z' e* e  x5 c
    xlabel(\'时间(天)\');# b$ F, O8 h. s3 G
    ylabel(\'仿真输出结果\');; ^( C) ?0 Z+ p4 ]; X
    legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');* \8 u+ y% `- {
    %%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    ) F6 P, l2 ^  \/ \$ Y0 Xx=1:4;! Q2 @8 M8 [  F, d
    y=1:4;
    * c! {  V, P5 e9 j# H" By21=sim(net,p);& m% `7 I, y( k( F6 i) I% h
    y2=reshape(y21,1,12);% Q; y8 a. z& Z2 [* f$ Y& |
    clf) d( I1 O' y& O0 @+ i
    plot(r,y2,\'b-^\')0 Z5 }, ~0 l" j; T/ l
    hold on# C' J* D! }! z& ?3 m
    plot(1:20,day,\'r-*\')) |& t6 h0 q# v4 P
    %%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%
    - A; M/ P, c  _6 m- Ny3=sim(net,h3);
    ' _( A6 b3 I. i) I4 Xplot(17:2:20,y3,\'-*\')# \6 r# d$ B( e' ^8 p5 f4 a( x3 H
    hold on! L# m" v$ j9 n: G' @, Z( c
    title(\'神经网络训练结果\');
    5 {6 C6 _) U  u6 O3 F2 ?xlabel(\'时间(天)\');
    / v2 i- t, v4 q& u5 W2 n! @: zylabel(\'仿真输出结果\');, |3 |, h8 Z: l% L- \2 W6 i
    legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    ! O( G+ ?8 m" I3 k7 @%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    ! ?& B" O! c8 v* d% vx=1:4;
    * E! D8 M, G1 t$ y' \& wy=1:4;. }( `+ R4 G8 g5 S/ ]
    plot3(x(1:3),y(1:3),E(x(1:3),y(1:3)))
    ) a4 Z9 f' J3 C: G& K/ r( J+ k运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢
    ; F) }$ W. P$ g: r+ c+ ]

    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-13 12:28 , Processed in 0.368423 second(s), 50 queries .

    回顶部