QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2983|回复: 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];+ O7 e9 S4 U6 s) p5 b
    dayhistory=day(1:16);%取其中三十天作为历史数据样本
    8 `, o8 {  q( _( m8 A, J* T5 ydayhismod=reshape(dayhistory,4,4);% 将历史数据分为6个样本,每个大小为5,其中reshape是以列排序的9 B  P3 J" Z, G0 Q* [2 J# Q* W
    dayday=day(1:12);% 取其中的前25天$ v" h: D" @  K. h+ X
    daypost=day(5:16);%取其中的随后25天5 N3 [; H/ Z! t2 ~. p. x
    p=reshape(dayday,3,4);% 将前25天数据分为5行5列矩阵作为网络的训练输入样本$ g; L; r. _& l
    t=reshape(daypost,3,4); %将随后的25天分为5行5列矩阵作为网络的目标输出向量
    * N+ j" y) a; a$ t( v; \daylast=day(13:16);
    & x, z9 F& O/ e& L6 Ih3=reshape(daylast,4,1);% 将倒数第二个样本作为网络测试时的输入样本% ^( O- D: J1 U; c
    r=5:16;! q" G2 j# G# |
    rr=reshape(r,4,3);
    1 D6 B/ c  Z6 o  R1 L+ v2 D%%%%%%%%%%%%%% 新建网络bp %%%%%%%%%%%%%%%%/ z+ Q+ H4 \) i" V9 i
    net=newff(minmax(p),[4,3],{\'purelin\' \'purelin\'},\'trainlm\');
    % u2 c; t; f) N9 My1=sim(net,p);1 M( M  ^( @# s' [) K% z  C
    % 新建网络,其中minmax(p)为p的没一次输入的最大最小值向量
    $ X4 c& J8 {0 I$ u% 两层的传递函数均为purelin
    ! }7 c* j/ j# R  M) e6 |2 Z% 训练函数为trainlm
    ! M) n9 `" C) v' ^9 \0 P% 所训练的网络大小为[5,5]
    6 x9 }8 O7 q+ v* C% 仿真训练前的网络: o8 _0 ~8 x8 P- W8 G0 W
    % _8 N& v7 `- l; b
    %%%%%%%%%%%  进行网络训练  %%%%%%%%%%%%%%$ }8 G, h0 t6 ?, s1 p- r  x
    % network parameters:# D& Q" C9 j" U% O
    %   epochs--epochs of the train6 G# U1 T  A8 G, z8 X3 n* x$ F
    %   goal--errors goal of the network
    2 z% K7 A. F" R: B' p3 d%   lr--learning rate+ ^6 ?/ V1 I/ c. |6 ]) ]$ J* o3 x0 w
    %   shows--epochs between the displays
    / m- a9 w! p! B- s6 g/ b" i%   time--Maximum time to train in seconds
    " f0 L( o$ B5 G) Cnet.trainParam.epochs=200000;  % 训练次数( H4 ?% q3 `7 ?1 f: Z2 P
    nettrainParam.goal=0.0001;  % 误差期望值
    6 s+ o1 U+ O/ n7 N2 t4 \  L% returns of the train:
    ) v* s8 x1 ~7 \3 S%   net--New network. [! s' D, c% W* h6 A
    %    tr--Training record (epoch and perf).
    0 W" m9 [  E$ y6 O* Z" U%     Y--Network outputs.
    1 t0 E; c* y/ k0 R  E$ X%     E--Network errors.
    $ _6 p/ T2 E: o- l/ G0 K9 d; y[net,tr,Y,E]=train(net,p,t); , k% z& a* R1 W3 n: F
    %%%%%%%%%%%  网络测试 %%%%%%%%%%%%%%%%
    * x/ i' k7 _) b% R9 G% input the testing points here %4 e1 u' Q) R$ W7 K* e0 v
    title(\'神经网络训练结果\');
    7 P! v4 g: e% n. sxlabel(\'时间(天)\');
    * J/ c9 S1 u  Y& L( z9 o% D' |ylabel(\'仿真输出结果\');2 P2 g7 p3 e8 ]8 X8 p: j  `
    legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');' @2 V! l% x0 t" r) ?* G) K; \, L
    %%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    3 W+ \1 @+ w- ^" qx=1:4;  y, {  E$ S/ z. t
    y=1:4;  _# w9 h$ S% W) F& H
    y21=sim(net,p);
    0 C) C$ h! p1 A0 I+ D( J, B9 c+ T- k" Hy2=reshape(y21,1,12);
    0 B3 w# z0 ]: M: Lclf/ d. j! ]% A) z" f
    plot(r,y2,\'b-^\')& F( \% S# r0 F, G+ \4 L  O& ~( y
    hold on
    - O( x$ o( x5 k8 k' F6 k* O) xplot(1:20,day,\'r-*\')
    ( z; u' _: O/ D; p7 i/ }5 l- v( p- S  P%%%%%%%%%%%%% 预测 %%%%%%%%%%%%%%%
    " c5 c8 m, G& L3 b- B' Y  Dy3=sim(net,h3);$ f5 l9 c  k1 U+ k1 T) {
    plot(17:2:20,y3,\'-*\')
    ( c# C, j9 L$ Rhold on
    - G$ a4 c8 H! w4 o7 Stitle(\'神经网络训练结果\');
    ! u8 A6 L$ C0 ]1 p  [. f3 yxlabel(\'时间(天)\');* f. ~5 O, ~3 b; {
    ylabel(\'仿真输出结果\');2 P" x( J# w& e! Y* J
    legend(\'仿真模拟值\',\'实际值\',\'神经网络预测值\');
    0 g' T9 @% ?- J: r8 A%%%%%%%%%%%%%%%%%%  绘制误差曲面 %%%%%%%%%%%%%
    / l) \& v2 [. D( M6 ox=1:4;
    5 V* |" O/ d, ?y=1:4;) C* {; c6 x) K5 T: j. Y9 `
    plot3(x(1:3),y(1:3),E(x(1:3),y(1:3)))% Z' g1 t' ]1 x! h, V- k
    运用简单的神经网络,怎样处理已经有的数据,一般输出的怎样的结果,有什么用,谢谢
    , e% B& V; V5 N& \( e! h

    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 01:02 , Processed in 0.488302 second(s), 49 queries .

    回顶部