QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3887|回复: 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 编辑
    , b# L; U! @+ J5 i4 V2 _5 K' F5 M9 v7 \' s
    各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。: k+ k; @$ S8 a* M6 V8 n1 G
    国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。3 D( B' l2 o# h2 I+ x

    4 v& c# ]1 P8 c8 f; B+ Nclear
    ; ]+ K3 G5 l1 Oclc
    9 |* p* o% d; ~9 \3 @( S( ?' u* i%输入数据--以一列作为一个输入数据
    4 @! v* _8 G0 z- W. n" [% E/ R8 D; o$ `x1=[227.5 286.9 236.3 300.6 243.0];: K; r' E& n3 \
    x2=[286.9 236.3 300.6 243.0 354.1];- A# y1 x; l4 ]' R. E. J3 B( O0 [
    x3=[236.3 300.6 243.0 354.1 294.1];
    # ?$ b3 Q( C, z: j7 C: }+ v6 x8 jx4=[300.6 243.0 354.1 294.1 361.4];, D7 h  k$ f4 C6 F5 e0 m
    x5=[243.0 354.1 294.1 361.4 333.7];
    , G7 N% K- b& S9 j2 H# N" G% G( `4 wx6=[354.1 294.1 361.4 333.7 426.4];
    , I# s! C8 b2 |+ B: A3 r9 vx7=[294.1 361.4 333.7 426.4 434.4];  X( t8 j: b8 e% }8 ]0 v
    x8=[361.4 333.7 426.4 434.4 449.9];
    * x* {8 A9 n" ~& e  ?2 kx9=[333.7 426.4 434.4 449.9 431.9];( q1 R0 Q( m: s3 J) F

    3 e& q- |- `6 I9 q& H7 X5 v+ NP=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';
    6 l. |) q4 Z( H1 }% w% R6 F1 p3 A%输出数据--以一列作为一个输出数据9 n( g2 K. z$ ~6 v9 R' |* s# J# O
    y1=[354.1 294.1 361.4 333.7 426.4];
    ) j/ w6 l+ N  U3 Cy2=[294.1 361.4 333.7 426.4 434.4];
    " |9 v/ y0 j' M$ o% f) q  {) F5 f8 py3=[361.4 333.7 426.4 434.4 449.9];
    : X" }# Z0 k$ J1 F& fy4=[333.7 426.4 434.4 449.9 431.9];* R# E% L5 L6 w- G- T
    y5=[426.4 434.4 449.9 431.9 454.0];
    * h8 W* `3 _; H* z7 {( ^y6=[434.4 449.9 431.9 454.0 655.6];$ j; |8 P& M6 `% g) i# D7 M
    y7=[449.9 431.9 454.0 655.6 907.0];
    3 c( ]$ U3 j2 F" O$ v& My8=[431.9 454.0 655.6 907.0 1112.5];
    5 v* h5 e, q6 V& C8 R5 ^y9=[454.0 655.6 907.0 1112.5 1401.1];: J8 a' _9 J4 w/ s3 K1 ~. r+ _
      n, Y, S2 ~7 j
    T=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
    $ w4 Y3 }: m  ?6 W! s3 a: `# V# x8 L* W2 P& Q% o3 T
    %归一化
    ) X, C+ z' @1 K2 {6 [/ t% k% [P,minp,maxp,T,mint,maxt] = premnmx(P,T);
    . E# B; I( k4 X! G/ W" o%建立网络, r: U0 Q, Q% A" s% i
    %其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据
    9 T# d8 G; @7 b9 T; d%{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
    $ ]; _/ V1 w# M# u# X4 X%'traingdm' 表示调整权重,缩小误差的数学函数, w7 d& K: S8 E2 [6 N: W; g6 T- v
    %↓创建一个新的前向神经网络  
    4 F  }4 @- b5 H( C, @, m. n3 Q, D; Jnet = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')6 w) P& W: J) ]3 Z) O+ j

    0 G  ?6 y1 `2 Y7 b# e%  当前输入层权值和阈值 ! B' C3 E/ m, m/ q# B$ h
    inputWeights = net.IW(1,1)8 M" \3 f- `  K" o0 D4 |3 J- d
    inputbias = net.b(1)
    4 V  T& f0 d7 W+ N  R. r9 ~6 D( v; W  i8 J
    %  当前网络层权值和阈值
    3 g2 V  K1 n( a% alayerWeights = net.LW(2,1)
    1 A7 K/ {( a8 c1 u. f/ Dlayerbias = net.b(2)
    , B, y: k8 w- k" g
    # |) h- ?, l! L( m' B%  设置训练参数
    2 B3 w% H1 V# mnet.trainParam.show = 50;# G! I" M5 q- |0 ~
    net.trainParam.lr = 0.05;
    8 ]+ r2 ^0 s% k% ^$ u& |0 N8 Knet.trainParam.mc = 0.9;! J9 _; {: n8 \7 K( _( G
    net.trainParam.epochs = 10000;  %每计算10000次出一次结果5 ?7 @9 B" ?' L: F
    net.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)
    / k& ^0 a# O( i- {& O& B5 p9 K: H
    7 k+ O- [, {' G' t6 K- Z# n%  调用 TRAINGDM 算法训练 BP 网络
    3 c4 V8 {# L# C[net,tr] = train(net,P,T);      # F1 N( R& K4 o9 r  d
    / R& ^; B. c# s; c. X* R& m
    %  对 BP 网络进行仿真. J9 N3 c* v% R* q# C9 e. y1 L. v+ `
    A = sim(net,P);                 %A为输出结果
    , [, j! a0 I7 Q0 @) U4 i* H) c' x" N: H5 W& H  B
    %反归一化
    " q3 f+ O) T" T# H6 S- p/ w( T%  A = postmnmx(A,mint,maxt);. U) ?# r3 d5 |7 B
    ; Z# w, M$ U8 Q: l/ E5 U5 D* k
    %  计算仿真误差  & S4 x- F, u1 Y
    E = T-A;                        %原来的输出 - 训练输出 = 误差
    7 C! `$ J1 N  F1 e/ EMSE = mse(E)                    %输出误差% C" o; x+ o5 M: W3 Y
    * u( M0 d  T5 F
    %下面是输入数据,测试结果' ]! M/ q0 U# G8 t1 U9 Q8 j
    x=[454.0 655.6 907.0 1112.5 1401.1]';
    2 p( D$ Q& z5 |, `Y = sim(net,x)& A8 L& y; _( |2 R& K( [: B
    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, 2025-10-11 01:24 , Processed in 1.049602 second(s), 72 queries .

    回顶部