QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2987|回复: 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];
    * b9 q) p$ _( q. @* ]$ @dayhistory=day(1:16);%取其中三十天作为历史数据样本4 ]5 a+ l2 F. F( P
    dayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的
    % ~, Q; P" @' p8 Cdayday=day(1:12);% 取其中的前25天$ B. h' Z. w- _' I7 c& w
    daypost=day(5:16);%取其中的随后25天8 m" U6 @% \* ^# `0 G% h" d, _; U0 l
    p=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本/ _# e# k) ~, u
    t=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量% P* `2 n  ?. O' G. R% t
    daylast=day(13:16);
    % f. g% [' l) U4 Z- V' F) e- b+ nh3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本
    4 n5 O: n. e# vr=5:16;+ v# I, F) e% Z0 i0 H8 b3 S
    rr=reshape(r,4,3);# m5 N6 c# o1 `! G0 j
    %%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%
    / u, f- v3 P( R5 b# V6 t; Q0 r) lnet=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');& H* D7 A! H2 w3 ~. z, ~7 I4 U
    y1=sim(net,p);% O* B2 y1 d+ `" u9 o  z* l
    % 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量
    + s! F9 s$ |, q+ D/ i( w4 ^7 o% 两层的传递函数均为purelin4 w* J5 ]1 ]* V+ l- F1 U/ c. L
    % 训练函数为trainlm
    + G! Q* }4 e" s3 h% 所训练的网络大小为[5,5]2 L7 n0 E6 h7 @9 D' i, x( C. y' D2 [
    % 仿真训练前的网络( Z5 q% [' l( o1 {# O+ @  r: E

    1 r1 z; {2 f' ^( p%%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%7 ^# G# y2 R9 J( X
    % network parameters:
    7 E! [. p# J+ O/ `/ |3 a%   epochs--epochs of the train3 _0 z( b8 l% A# b, g4 Z6 _) l
    %   goal--errors goal of the network2 O8 u5 f0 o8 y) d  c
    %   lr--learning rate' _: m9 w+ W$ _8 S
    %   shows--epochs between the displays
    $ l  f( R# O: z+ e7 g%   time--Maximum time to train in seconds
    7 h4 [4 H8 N0 w! C* unet.trainParam.epochs=200000;  % 训练次数0 U! v+ l5 u' B( e4 U5 w) j% h% {
    nettrainParam.goal=0.0001;  % 误差期望值4 ~/ G5 X( @( d
    % returns of the train:  a3 X" U# o; `0 }" l9 e. d
    %   net--New network* d/ }5 ^7 v) u% D
    %    tr--Training record (epoch and perf).8 h8 I) |( A0 u7 y5 z: r
    %     Y--Network outputs.
    4 ?1 x  n& w4 }%     E--Network errors.$ p( H4 z8 C4 e& q
    [net,tr,Y,E]=train(net,p,t);
    $ R: ^* J- ]2 W# @7 G%%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%
    9 h& m9 Q$ z' E% `% input the testing points here %
    ' \' a8 v/ }& B0 \" G' {& k) ~title(\'神经网络训练结果\');, E/ A; b2 b4 G/ Q
    xlabel(\'时间(天)\');
    9 w$ h4 m- w: I# Lylabel(\'仿真输出结果\');
    0 {" q: Y. Z! n+ `$ y8 hlegend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');6 G0 B" L" T- `, ?& b4 p
    %%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%9 N, c, ~3 A) e1 }* G: O5 d) d
    x=1:4;
    2 W" V- `' u: y/ [y=1:4;
    % i) K% f# X) ]5 e" |# j( Ny21=sim(net,p);: e9 I- @0 Q% {* `
    y2=reshape(y21,1,12);
    . Y3 g, P) ^' ~& @2 iclf
    ! l( l/ Z3 E- f+ m- c3 |3 [( Eplot(r,y2,\'b-^\')
    4 S& }# {  r( x0 F- I* e+ Ahold on5 [5 k2 a) U, W6 R
    plot(1:20,day,\'r-*\')3 }8 X) f5 o. U8 Y# Y! B( D! u
    %%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%
    0 N" [4 F* C1 C" N; H8 X3 v' {y3=sim(net,h3);
    . t: |3 m( i1 p1 [, N4 |plot(17:2:20,y3,\'-*\')
    $ Q4 N8 y$ L4 k" T5 U+ V. fhold on; [- f; k: h2 I+ F2 r+ J
    title(\'神经网络训练结果\');
    ' U& V  @1 s1 I. g0 ^5 G2 f9 Kxlabel(\'时间(天)\');) {' j" c9 ^* }/ [3 ]
    ylabel(\'仿真输出结果\');
    7 o7 r* ], j0 x. q; Vlegend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');4 m+ O) l/ M% r3 ?5 H! x
    %%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    ) {4 ]4 p$ c7 H1 }5 Y& ~3 v* \) \x=1:4;
    ) Q: S) ^' M' M. W; Py=1:4;
    7 [8 f8 u. z" a/ v5 \plot3(x(1:3),y(1:3),E(x(1:3),y(1:3)))
    ) N7 P$ y" g6 `6 _4 z7 F运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢
    7 i% c# J! D2 ^. k% D6 e' D6 f

    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-14 04:52 , Processed in 0.413003 second(s), 49 queries .

    回顶部