QQ登录

只需要一步,快速开始

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

[已经解决] BP神经网络预测遇到一个疑问,麻烦过来看看

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

23

主题

8

听众

372

积分

升级  24%

  • TA的每日心情
    开心
    2014-10-4 18:19
  • 签到天数: 102 天

    [LV.6]常住居民II

    自我介绍
    本人是在校学生,即将参加有史以来的第一场国赛,希望能在这里获得帮助。
    跳转到指定楼层
    1#
    发表于 2013-9-8 17:22 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    本帖最后由 数学不烦 于 2013-9-9 21:56 编辑
    0 m2 c9 U' K, R4 I: g  D4 }! m% U6 ^3 j9 V5 L! m
    各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。- b3 e& x% |( W  S1 H9 D
    国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。8 l6 L4 u" ?5 ?4 l

    ! l0 X# r. v) \clear
    8 S- o! c- L' i: S4 qclc+ j& v; y& T" \& N) }
    %输入数据--以一列作为一个输入数据
    % e8 `( D! M8 O  o3 Ax1=[227.5 286.9 236.3 300.6 243.0];
    . `8 l% N4 L$ |! Yx2=[286.9 236.3 300.6 243.0 354.1];
    : ~, M* }$ ?3 u' L( O7 L9 \x3=[236.3 300.6 243.0 354.1 294.1];7 g- m) e9 R* v* V5 ]* X
    x4=[300.6 243.0 354.1 294.1 361.4];6 r/ ]! d2 `& `- `" O8 ~- ~
    x5=[243.0 354.1 294.1 361.4 333.7];
    3 b* t" _0 A& H  Fx6=[354.1 294.1 361.4 333.7 426.4];/ x, U3 H+ h: Z1 H: C; M7 W! J
    x7=[294.1 361.4 333.7 426.4 434.4];7 }( m0 \0 w: O5 I
    x8=[361.4 333.7 426.4 434.4 449.9];( l: }  H( ^9 ?
    x9=[333.7 426.4 434.4 449.9 431.9];/ b4 t. H' @6 W6 O! Z9 R) `  Q6 `
    ! r$ B  i: [; x* G5 |, ~9 W) A
    P=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';+ E/ |. o2 }7 B! S# F: Z* q
    %输出数据--以一列作为一个输出数据& n' @6 g6 Y6 Q; O& b; B# `( y8 ?2 f
    y1=[354.1 294.1 361.4 333.7 426.4];
    # r% t4 Y- B; K+ T6 u+ u, \y2=[294.1 361.4 333.7 426.4 434.4];
    6 m# q5 m  F: q" l. M% E+ Yy3=[361.4 333.7 426.4 434.4 449.9];9 D# H5 r' {1 A- s2 b0 J
    y4=[333.7 426.4 434.4 449.9 431.9];% n9 h  ?/ F3 v0 i4 ~% b
    y5=[426.4 434.4 449.9 431.9 454.0];7 {2 V* k! Y2 |* Z' N" @% m) L
    y6=[434.4 449.9 431.9 454.0 655.6];7 b9 `3 C( U- N
    y7=[449.9 431.9 454.0 655.6 907.0];8 d4 Z9 q$ _! a! t: H! R% K
    y8=[431.9 454.0 655.6 907.0 1112.5];
    ' ~( t+ W- H5 i! ^5 W2 Qy9=[454.0 655.6 907.0 1112.5 1401.1];
    3 B7 f- D2 ~/ b: u' E4 t, E6 _, r3 j8 V. d
    T=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
    / F2 `' O0 E- n9 j
    ! K) |: E* A$ R* Q%归一化
    & C' v" a6 Q' I8 l; p! y% [P,minp,maxp,T,mint,maxt] = premnmx(P,T);, |( V9 r+ j) K. w2 D/ Y; m" v- }" ^
    %建立网络
    ! I8 {" V* ~- g# O1 w2 A2 h7 n/ h+ q%其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据( x: _6 T! I% L( W4 A5 d! x
    %{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
    & ]& \: {2 g- b" z%'traingdm' 表示调整权重,缩小误差的数学函数
    5 C8 [9 B; i7 \5 k) ]9 j. K%↓创建一个新的前向神经网络  
    1 Q& i  D  C8 |% Y* {net = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')
    ' F7 U3 Z5 g' z6 _. n2 i
    ) F$ u0 W2 @7 v, }9 G6 L%  当前输入层权值和阈值 - Y3 z* G% p. I$ a) _, `
    inputWeights = net.IW(1,1)5 ~3 U. T$ v. [! @  a
    inputbias = net.b(1)
    9 h9 m5 \- N  Z+ u3 d. ~* s# ?! U% x) E, ]) b2 u5 j* _' h0 y
    %  当前网络层权值和阈值 " u% E8 m  h. g/ J$ T4 J+ l' x
    layerWeights = net.LW(2,1)1 y% E; ~) Q4 Q) B
    layerbias = net.b(2)
    * c6 Z2 X4 i* H+ s% y4 i1 P* ^- K. p" W; S/ V
    %  设置训练参数
      `4 g) F- j  V8 f5 wnet.trainParam.show = 50;9 ~  t: Y. r2 b, Z# D2 x2 }
    net.trainParam.lr = 0.05;6 \/ D+ P2 @9 r/ F
    net.trainParam.mc = 0.9;
    4 L, l0 x0 S& S/ `3 v- ^net.trainParam.epochs = 10000;  %每计算10000次出一次结果
    $ j. g8 Q' B: V2 Cnet.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)$ A$ s1 V2 |; a2 O. ]- n  Z- s
    1 T+ ~1 ~/ A# v7 l
    %  调用 TRAINGDM 算法训练 BP 网络
    ! m+ B0 i# F7 l9 V  b: L[net,tr] = train(net,P,T);      
    0 l# O4 c" t$ L: P+ _# y1 T3 Y. ^1 x: b
    %  对 BP 网络进行仿真
    + h+ C4 H5 G: {/ N2 b( a6 G. L9 FA = sim(net,P);                 %A为输出结果- t# o" ^9 t3 L" h) a! F

    * b& A" z+ z5 A; `%反归一化6 Z5 |, U! ^8 ?3 U1 ~7 N
    %  A = postmnmx(A,mint,maxt);
    . [0 b- X' ]2 d& h
    . O6 \/ g9 S' j( }5 l%  计算仿真误差  
    8 \+ n5 ~" A" V% W9 ZE = T-A;                        %原来的输出 - 训练输出 = 误差: k0 X. h3 S( N5 D& ~& S. F# r$ _
    MSE = mse(E)                    %输出误差1 G: m* |* O& b) D

    / _/ r: Q: q6 U$ y" p' L) c%下面是输入数据,测试结果
    * o+ H5 w% P. A' P1 ~+ K' Z9 `x=[454.0 655.6 907.0 1112.5 1401.1]';
    ; U0 H( q) s8 ]. S* `Y = sim(net,x)
    2 D7 s9 f; X- Y$ q- u- r4 ?
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

    19

    主题

    11

    听众

    272

    积分

    升级  86%

  • TA的每日心情
    郁闷
    2014-10-23 09:29
  • 签到天数: 59 天

    [LV.5]常住居民I

    不知道你最后是怎么解决的,感觉你最后预测的时候没有把输入数据进行归一化,即要对x进行归一化,再对y进行反归一化,这是我的看法,不知道对不对
    回复

    使用道具 举报

    kaito1412        

    63

    主题

    8

    听众

    3098

    积分

    升级  36.6%

  • TA的每日心情
    开心
    2020-7-21 15:56
  • 签到天数: 829 天

    [LV.10]以坛为家III

    回复

    使用道具 举报

    magic2728 实名认证    中国数模人才认证   

    61

    主题

    478

    听众

    4851

    积分

    升级  95.03%

  • TA的每日心情
    慵懒
    2014-9-29 19:37
  • 签到天数: 409 天

    [LV.9]以坛为家II

    群组数学中国 2015美赛护航

    群组数模专题强化培训

    群组建模思维养成培训

    群组2015美赛护航(强化)

    群组2013年数学建模国赛备

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-9-3 09:43 , Processed in 0.508393 second(s), 75 queries .

    回顶部