QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3156|回复: 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];
    ( ~8 m& J7 S' l$ L. s8 B- O0 mdayhistory=day(1:16);%取其中三十天作为历史数据样本
    / Z9 g4 Y0 H* T! D4 xdayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的2 a  I4 X, y/ D  }
    dayday=day(1:12);% 取其中的前25天8 J4 e* ^0 @- s; U& h! K/ m7 D% m
    daypost=day(5:16);%取其中的随后25天. i7 Q- h7 [1 M9 X
    p=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本8 ?" Y3 ^7 h% A* }7 K
    t=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量- x! V/ u! ]1 S
    daylast=day(13:16);" Z0 g+ h. M$ ]+ D5 p4 J
    h3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本
    2 x$ u+ ^; B9 y' b# @. H9 Zr=5:16;8 e! |. n+ T3 S. y, s1 U
    rr=reshape(r,4,3);6 I, A) X: z5 {
    %%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%
    4 {6 ~2 ~9 l  V3 z) w' y, j7 Hnet=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');
    $ A+ @9 U% p4 Ty1=sim(net,p);5 |0 o: g* ?" n* w# d
    % 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量
    4 q, U# n& m. C1 F% 两层的传递函数均为purelin* D* W. G  Q! c  Q/ W7 D
    % 训练函数为trainlm5 F" l4 L( `0 R# r) }* `" G
    % 所训练的网络大小为[5,5]
    ) U. _$ A) y, Q1 l& R' u% 仿真训练前的网络1 H6 `4 I  H- i* k0 a+ k
    & S" f& o8 e$ r- J: T3 k& `3 T$ J
    %%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%
    3 K# q2 N2 L# f% Z9 z- ?+ C7 t% network parameters:6 n% e; E' H, g1 z3 V7 g4 O8 w
    %   epochs--epochs of the train) P/ D  ~' ?; w1 J' C( m
    %   goal--errors goal of the network' F# J* E1 E/ n' F$ d- i+ _5 ~' ?, a
    %   lr--learning rate
    / U4 I# i! S; a8 B" ?%   shows--epochs between the displays
    ) e: {! H2 Z1 i! r5 }: N%   time--Maximum time to train in seconds6 U" |8 _* p1 ?; ~' S
    net.trainParam.epochs=200000;  % 训练次数
    ! Y( |, ~& w) D' P& d! U" qnettrainParam.goal=0.0001;  % 误差期望值
    5 w3 v. Z8 n$ n4 N% returns of the train:
    6 h, a1 c$ z1 {" |) K%   net--New network
    # I9 c; n5 P+ E%    tr--Training record (epoch and perf).- Y7 E8 Z) p( Q
    %     Y--Network outputs.
    1 d# w$ c8 i4 x* V9 \%     E--Network errors.& w+ E; P/ B. [7 n
    [net,tr,Y,E]=train(net,p,t);
    8 B: T* y' r) }3 _" d: d" u, V. Z%%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%
    ; T* B! C, y. y* V: U/ |! y; {% input the testing points here %  e% N8 E- S9 S
    title(\'神经网络训练结果\');
    ' c0 V( D% _9 dxlabel(\'时间(天)\');5 Y' \4 G8 `. C: V2 \( a
    ylabel(\'仿真输出结果\');
    ) f2 s3 s3 o% E) m2 Flegend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    . r! ]9 K& f: J%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%' Y* d2 [& j  E' I5 v
    x=1:4;
    : Z' S/ M5 l; x/ P  J2 |4 v8 ey=1:4;" J! \0 o3 Y" e- T: S  z( N  `9 H
    y21=sim(net,p);
    % ?( M$ i  c: p% g; wy2=reshape(y21,1,12);
    / V' @6 A+ {9 s# wclf- [  s3 x9 d  j( u7 J
    plot(r,y2,\'b-^\')
    " k9 e) P2 V0 _4 Z0 Chold on- X% I/ T- `6 b: v& `" q8 `* a
    plot(1:20,day,\'r-*\')- k4 g& A* O0 K! e2 \
    %%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%9 S7 S. N6 I0 R# i5 a9 d" T* B
    y3=sim(net,h3);) X( i4 u& `% |6 G  _" `& f# M
    plot(17:2:20,y3,\'-*\')! q3 M( y) a$ r3 O9 l8 J
    hold on
    - T* ?  D5 J  q9 f% mtitle(\'神经网络训练结果\');3 v  K" R4 @8 O$ K3 ~
    xlabel(\'时间(天)\');
    2 k: Z! g0 X1 U; f- tylabel(\'仿真输出结果\');
      z) f0 {& m( d- p; d+ i/ {legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    ) B3 b$ f8 H1 o4 q: X- f%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    . M* Y- B: m2 nx=1:4;
    ( h; O+ L; N5 i2 J1 ky=1:4;
    - h. q/ a) {0 |3 @7 Tplot3(x(1:3),y(1:3),E(x(1:3),y(1:3)))6 w; K  `5 f! O0 U
    运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢

    ) x4 L5 Z. O$ w, A! W3 K; P6 y5 q

    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 01:52 , Processed in 0.477012 second(s), 48 queries .

    回顶部