QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4116|回复: 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 编辑 6 O' l5 U. [" c: i; \( n" N& x/ c

    7 \. ]0 U% Z1 e% _4 R9 l各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。. u6 d1 z; l* v5 i7 q4 [6 ?
    国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。. x* |" \6 o+ ]5 M1 M
    ( C6 {. d* k, N
    clear
    - s/ m! g+ @! p- Y. Uclc+ R8 h% `" q# ~5 J8 U$ {" l  S
    %输入数据--以一列作为一个输入数据/ K& v. C7 N! J6 p: u$ |1 _7 m8 X
    x1=[227.5 286.9 236.3 300.6 243.0];# L' n2 I' ^: Z
    x2=[286.9 236.3 300.6 243.0 354.1];
    , O. |9 ?: r* O, `+ jx3=[236.3 300.6 243.0 354.1 294.1];
    & y5 t. a" \/ G8 M3 Q' D) O: Gx4=[300.6 243.0 354.1 294.1 361.4];( g$ {' l3 i2 R( H6 B, F0 d
    x5=[243.0 354.1 294.1 361.4 333.7];: B2 X3 I! S. h* X2 Q4 J
    x6=[354.1 294.1 361.4 333.7 426.4];
    $ F. I: j# u" ]- y; sx7=[294.1 361.4 333.7 426.4 434.4];. S, B& o" p  [/ b+ s) m
    x8=[361.4 333.7 426.4 434.4 449.9];
    - {: k/ {$ A/ F2 K2 ix9=[333.7 426.4 434.4 449.9 431.9];
    ' S* _0 q, J' _) p0 R5 {- x2 _- o1 c% V; X( v2 N3 e
    P=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';
    # ^2 D6 \1 K" E, C2 R! V%输出数据--以一列作为一个输出数据$ {% \; t# z6 C% g( y  k$ p
    y1=[354.1 294.1 361.4 333.7 426.4];
    " u. |0 M$ R# J5 Uy2=[294.1 361.4 333.7 426.4 434.4];$ Z/ W3 O& a+ ]1 B# p2 ^3 ~5 ?  {
    y3=[361.4 333.7 426.4 434.4 449.9];
    % c5 e+ N9 z7 u0 Oy4=[333.7 426.4 434.4 449.9 431.9];" |  l- S" g: N* d. S
    y5=[426.4 434.4 449.9 431.9 454.0];
    5 A1 U7 b  S9 O$ Y8 Jy6=[434.4 449.9 431.9 454.0 655.6];
    , @  L0 y0 h) T9 `" Ay7=[449.9 431.9 454.0 655.6 907.0];. H# a# c8 c& e& U8 k( ~0 r, k
    y8=[431.9 454.0 655.6 907.0 1112.5];, m9 J5 u1 F0 `3 g  Y0 c
    y9=[454.0 655.6 907.0 1112.5 1401.1];5 l/ T3 s* g& b. m- W% T# S7 {

    . p2 o: }% p. j" x' Y& |! wT=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
    , X# n! q# K$ A+ d$ c" Y: _& |$ Z! E* Z
    %归一化5 ?, @! _; K$ u) W$ A5 P
    % [P,minp,maxp,T,mint,maxt] = premnmx(P,T);; ?2 G' u' x7 j+ \
    %建立网络
    7 i9 X% U$ U/ k4 ^7 N4 a%其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据- e  H8 h& L' S- K
    %{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层: e+ f7 d, b9 b6 `4 {1 X
    %'traingdm' 表示调整权重,缩小误差的数学函数) }# f4 c  S" X5 ^
    %↓创建一个新的前向神经网络  
    + {/ z- C+ m9 @5 t* t7 ?9 u2 ~net = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')
    1 l" B2 m6 G. t# {4 H
    , ]2 T8 b0 R6 Z& v- ]! V* ^* W%  当前输入层权值和阈值 & j4 E4 A- L6 }+ _8 T7 J& U
    inputWeights = net.IW(1,1)7 Z" ?# B6 I# g) e0 O6 Y1 W
    inputbias = net.b(1)1 C. w* Z9 v7 c/ z# B2 K. d0 ^

    & p- k" r0 k' V. `- _%  当前网络层权值和阈值 4 V) r$ w! T! [, d( b7 N$ `! @- h
    layerWeights = net.LW(2,1)" r+ b# v! W0 a4 G' j  x
    layerbias = net.b(2)
    - T5 i. z' ~8 p5 N( R* E$ u
    # E/ m# ]5 [) m%  设置训练参数( n! G9 A% U; o5 a0 w  a9 l
    net.trainParam.show = 50;
    : a& l3 h  q" y3 P( bnet.trainParam.lr = 0.05;
    * F' |9 [( t1 E: C/ }/ \net.trainParam.mc = 0.9;
    & q9 N$ x% d5 V1 enet.trainParam.epochs = 10000;  %每计算10000次出一次结果
    % v; _% z& x5 U0 Q6 P3 L; W  o" ]net.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)$ d9 ^7 O, y# J& J+ G
    4 a, d0 t% c6 e* i
    %  调用 TRAINGDM 算法训练 BP 网络! g/ |) `0 f+ K( T' B( Y
    [net,tr] = train(net,P,T);      
    . m& |6 c3 a& d& X% ]
    : t+ `# O2 }+ D5 e4 a9 Y) n%  对 BP 网络进行仿真* ^) T. q# i3 q5 f- d5 O; b0 f
    A = sim(net,P);                 %A为输出结果
    * X5 t1 A) I2 \# V9 _) O* R6 k" K' Z6 x& E  I
    %反归一化
    8 `) ?, Q0 o; W2 c4 p6 f) B%  A = postmnmx(A,mint,maxt);
    2 X& c+ t2 H) h; J5 x5 o, x
    * v3 p5 k' _+ A) T- }: J- p%  计算仿真误差  
    " E% `7 l" ?6 Y; Z# @! J$ iE = T-A;                        %原来的输出 - 训练输出 = 误差
    2 }5 t$ t0 Z3 w) w: s. HMSE = mse(E)                    %输出误差; |( A( s# b' c( P4 X

    & I( W7 j& S7 s6 d. \/ y/ k/ C%下面是输入数据,测试结果% o, z% v+ s# Z- K0 ~
    x=[454.0 655.6 907.0 1112.5 1401.1]';* B$ X+ e" n7 {
    Y = sim(net,x)
    1 I- a3 P: r/ q) g* N
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    magic2728 实名认证    中国数模人才认证   

    61

    主题

    478

    听众

    4861

    积分

    升级  95.37%

  • 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-4-16 12:49 , Processed in 0.428151 second(s), 72 queries .

    回顶部