QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4223|回复: 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 |- u3 V" G* R/ A

    1 z: _$ i0 p, u# s$ u6 v各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。( n' d) _! U. p) o3 t& q2 u
    国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。% G6 \# k  O; X: A, y5 `+ A
    + q5 p% v4 C( |% s' e5 {
    clear
    & P+ e6 J; U& b; v$ D3 _, j/ d& lclc3 ~* ]( {& ~$ k! d
    %输入数据--以一列作为一个输入数据0 [$ k1 J+ {" a
    x1=[227.5 286.9 236.3 300.6 243.0];
    3 B0 A" x/ u5 Hx2=[286.9 236.3 300.6 243.0 354.1];
    / ]$ ~( k& H7 m( D! I6 xx3=[236.3 300.6 243.0 354.1 294.1];, G, s. n  C. Q! _! D7 k, Q
    x4=[300.6 243.0 354.1 294.1 361.4];
    . ~8 X0 y9 L. ]' f4 b6 i% ?x5=[243.0 354.1 294.1 361.4 333.7];0 Z: h# ~9 i. w6 P
    x6=[354.1 294.1 361.4 333.7 426.4];. n" j# e" w0 F2 j# P3 o1 ?
    x7=[294.1 361.4 333.7 426.4 434.4];
    + M& J+ _. y2 g- E* |$ D; hx8=[361.4 333.7 426.4 434.4 449.9];  ~7 V3 k# T$ U' F
    x9=[333.7 426.4 434.4 449.9 431.9];+ F4 u" l5 }2 }8 t/ h3 R4 ?8 J& y
    6 g6 K: v. g3 U2 p
    P=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';0 {7 L, _/ k: J% K1 r6 M9 Z  i
    %输出数据--以一列作为一个输出数据( K+ n0 t& [3 j# d8 H
    y1=[354.1 294.1 361.4 333.7 426.4];
    " E, o. @' a$ ?y2=[294.1 361.4 333.7 426.4 434.4];
    & Y6 r$ d2 @6 [( M2 `1 A' Fy3=[361.4 333.7 426.4 434.4 449.9];, u. x& l( Y' F! p
    y4=[333.7 426.4 434.4 449.9 431.9];) h: p% `5 J7 c# }# X
    y5=[426.4 434.4 449.9 431.9 454.0];
    9 H" M1 ?( ?/ V$ f' w7 Y2 [4 Ey6=[434.4 449.9 431.9 454.0 655.6];9 E% ]3 L7 {6 l$ a& ^6 `, b6 K; z
    y7=[449.9 431.9 454.0 655.6 907.0];* S) ]9 i, m/ E( p! k: ]
    y8=[431.9 454.0 655.6 907.0 1112.5];
    ( T( b$ z5 |' D* f0 T) Iy9=[454.0 655.6 907.0 1112.5 1401.1];! ?- P( S) ~/ K; `- h

    % j/ Q$ M6 w# B. IT=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
    ' a4 l& `1 A* ~' _9 `
    0 h( d/ n5 r$ b' D' H; V# z%归一化. ^$ C$ [) \# G- j& |+ q
    % [P,minp,maxp,T,mint,maxt] = premnmx(P,T);% g3 v6 \/ r' r) k" J. F* Z
    %建立网络
    " \1 c* s: s# h7 f%其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据4 |8 {3 Z: m0 t. g" b
    %{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层; |9 b* i  _9 |# M
    %'traingdm' 表示调整权重,缩小误差的数学函数4 c1 ]: w: B+ v& L/ C7 P
    %↓创建一个新的前向神经网络  
    # ]5 e! _8 h/ Qnet = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')
    0 S, `9 I( k% S- y9 l
    : L) r! I) b; [0 j2 g%  当前输入层权值和阈值 . Y" M+ {0 D9 v$ m3 l
    inputWeights = net.IW(1,1)
    3 m) ~. d4 T8 j/ B6 n* _  einputbias = net.b(1)
    ) e+ f4 l. q  [3 X/ p( r' o2 @) p8 h( ]+ A1 |6 \* n
    %  当前网络层权值和阈值 ' F' f. G, f/ x1 H4 a
    layerWeights = net.LW(2,1)! [5 y) y0 t! g' t% a, K
    layerbias = net.b(2)! _1 w2 H( i% H" s2 N

    6 X0 f  s) O$ y7 f%  设置训练参数
    7 A! K8 L3 u# d" u: B+ e* Znet.trainParam.show = 50;) \9 [  d1 m1 e
    net.trainParam.lr = 0.05;
    7 Z( a1 X. q, pnet.trainParam.mc = 0.9;2 U4 K8 G6 n1 o- L1 |
    net.trainParam.epochs = 10000;  %每计算10000次出一次结果
    ! O% o4 Q) X) R8 g  T5 w& c- hnet.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)
    , D6 a" x6 K0 J3 ]; A8 @* q/ l6 @0 i( I* d% F0 g- R9 W5 V
    %  调用 TRAINGDM 算法训练 BP 网络( u9 y. A0 A7 X. Z& [
    [net,tr] = train(net,P,T);      
    8 Q' B9 s, L5 r
    3 T" S- ~6 L4 e* T%  对 BP 网络进行仿真5 z" j* j+ f' ?8 Q6 M
    A = sim(net,P);                 %A为输出结果
    + o. s/ x& T# P6 D8 R0 T, w  U+ P7 G2 s5 g' U0 g
    %反归一化+ H' [# O9 @6 b0 q! z6 @9 }$ f
    %  A = postmnmx(A,mint,maxt);5 ^; q7 \; x0 Q
    6 Z1 p9 t9 Y& O1 f* ]
    %  计算仿真误差  
    ' \/ f  H( {7 ]) p2 kE = T-A;                        %原来的输出 - 训练输出 = 误差
    ! o7 Q7 [: E3 o$ }MSE = mse(E)                    %输出误差8 Q& v) ^. h0 H) L$ L. x1 M

    9 G# l2 Z( W$ z) u1 k%下面是输入数据,测试结果; h& Z  @& }) G0 D
    x=[454.0 655.6 907.0 1112.5 1401.1]';
      k/ u3 C, K' E; SY = sim(net,x)
    * g5 {; B+ T, P2 H7 G0 F# D$ Z
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    magic2728 实名认证    中国数模人才认证   

    61

    主题

    478

    听众

    4851

    积分

    升级  95.03%

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

    [LV.9]以坛为家II

    群组数学中国 2015美赛护航

    群组数模专题强化培训

    群组建模思维养成培训

    群组2015美赛护航(强化)

    群组2013年数学建模国赛备

    回复

    使用道具 举报

    kaito1412        

    63

    主题

    8

    听众

    3098

    积分

    升级  36.6%

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

    [LV.10]以坛为家III

    回复

    使用道具 举报

    19

    主题

    11

    听众

    272

    积分

    升级  86%

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

    [LV.5]常住居民I

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

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-9-3 08:38 , Processed in 0.341204 second(s), 72 queries .

    回顶部