QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3826|回复: 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 编辑 2 A% [7 d( q5 ?# x4 C0 d

      O6 I" E9 h/ I各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。; t9 \5 V4 i" P4 W8 I
    国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。( }0 q4 y! t( m5 I- [' M5 n/ O, i

    ; T6 `6 [& v& v  `8 z+ |, Dclear
    7 m0 y6 F) z2 t' G2 b! Fclc2 W! ~3 _' A, Y4 E
    %输入数据--以一列作为一个输入数据. H" U* V  j! d
    x1=[227.5 286.9 236.3 300.6 243.0];
    9 W) P( A: k3 l! x) Qx2=[286.9 236.3 300.6 243.0 354.1];
    : t! \) v, _2 J. ~! a+ cx3=[236.3 300.6 243.0 354.1 294.1];. c: L+ V9 ?  J1 b
    x4=[300.6 243.0 354.1 294.1 361.4];, ^6 `- w6 c  X9 |& Q
    x5=[243.0 354.1 294.1 361.4 333.7];. {: ~- j4 ]6 G( c3 e$ j0 z- B, M
    x6=[354.1 294.1 361.4 333.7 426.4];5 G  R  [- H/ o5 j7 E2 I* c& k: u
    x7=[294.1 361.4 333.7 426.4 434.4];
    ; A3 W, v( D; r( x0 sx8=[361.4 333.7 426.4 434.4 449.9];
    6 V1 w/ p# f% ~1 J6 [8 e1 a& bx9=[333.7 426.4 434.4 449.9 431.9];
    2 d  P5 D5 e3 U3 a6 K- m' ?7 H/ ]2 d- u
    P=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';+ X) L( f) e3 l2 g0 {
    %输出数据--以一列作为一个输出数据  t0 F6 B3 t7 [) v1 w
    y1=[354.1 294.1 361.4 333.7 426.4];
    # |+ t: f- }; {& ey2=[294.1 361.4 333.7 426.4 434.4];
    ) F# U- X9 r  \9 dy3=[361.4 333.7 426.4 434.4 449.9];
    ! S3 F1 v1 o4 K9 Z/ W7 k! U$ \y4=[333.7 426.4 434.4 449.9 431.9];
    ( m: m$ ^! J% s: i- _y5=[426.4 434.4 449.9 431.9 454.0];
    5 I) F4 O% A" s0 ]- hy6=[434.4 449.9 431.9 454.0 655.6];
    7 Z2 G3 N0 z' O. L& T7 _2 By7=[449.9 431.9 454.0 655.6 907.0];. v9 z7 ^4 o2 \! s
    y8=[431.9 454.0 655.6 907.0 1112.5];  q- X5 j* z) l4 x
    y9=[454.0 655.6 907.0 1112.5 1401.1];
    ; y! U1 u1 i$ G1 t; C2 s( B5 ]  Z8 q6 E* C' X& I" U9 A
    T=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
    : \. D# R; }, y) b* N' x  r$ H
    ) Z+ ^$ ~2 N+ T* b$ I%归一化
    . d+ X4 e+ ?4 U# z% [P,minp,maxp,T,mint,maxt] = premnmx(P,T);/ I3 }8 N" W8 }! h# f
    %建立网络% r' Z  C6 b! _
    %其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据! j) ~' O# r( X; `' O. ^
    %{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
    + c5 V0 ~7 `4 m0 N1 W%'traingdm' 表示调整权重,缩小误差的数学函数
    6 U! N' {1 B( O2 Z%↓创建一个新的前向神经网络  
    3 M8 I2 B1 q" N  ?+ Ynet = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')
    4 J7 q& L* f8 N- }) `5 U# S8 o
    5 b2 J$ g7 [$ z; F4 s%  当前输入层权值和阈值 7 }7 u& `$ f; Q8 m
    inputWeights = net.IW(1,1)5 ^) s2 Y* l* D8 P* u
    inputbias = net.b(1)0 f1 T$ M1 S7 g( W  x4 Z7 Z! }
    8 ?6 I6 K) g6 Z# U! X
    %  当前网络层权值和阈值
    5 C( N, j: M" k8 E4 AlayerWeights = net.LW(2,1)
    ( Q0 L: A7 @% x' y3 Y' L, O/ elayerbias = net.b(2)
    $ w) }" E# X& }- B, L1 _+ u" ~4 T1 V( B# u; J
    %  设置训练参数5 O( t' W. A8 I& L
    net.trainParam.show = 50;  y' Z& s6 p; G. ]+ k, p
    net.trainParam.lr = 0.05;
    , ?: N3 w: C7 r% _net.trainParam.mc = 0.9;, k8 C+ B+ x" q" D
    net.trainParam.epochs = 10000;  %每计算10000次出一次结果
    ' U2 G' T0 ?9 t$ l5 F( fnet.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)
    ' _7 V0 z7 g$ C: O# A  i
    ' t- W8 K) N& K- |$ t- S) i%  调用 TRAINGDM 算法训练 BP 网络3 w; s) c  Q  X$ W+ m) v% |
    [net,tr] = train(net,P,T);      
    5 ?' R0 G, H7 _' z6 j+ L0 J# l8 l
    1 A4 ]6 q' S/ J: b3 S5 c3 f9 O%  对 BP 网络进行仿真. R# [; q: A7 z* T6 h! x! b
    A = sim(net,P);                 %A为输出结果
    6 l2 B' ^' ^7 T. z. X+ b# r& G: [" u0 ], Y) d/ ^' X
    %反归一化# W6 i+ ^6 }7 I% {# @
    %  A = postmnmx(A,mint,maxt);) a, c( A1 w$ ]! Z( P5 p: f1 `

    7 P" U9 z) C) m& \1 E) h%  计算仿真误差  
      D$ k  Z/ b) X, P/ K( hE = T-A;                        %原来的输出 - 训练输出 = 误差! z* c$ n" [" d9 X) k: v1 A7 P  Y
    MSE = mse(E)                    %输出误差7 h+ w; m6 s: m( W# K
    ( a! \) Y7 }! I
    %下面是输入数据,测试结果$ H) [' I" q- e# U8 _0 F+ c8 ~3 A, [' @
    x=[454.0 655.6 907.0 1112.5 1401.1]';$ L+ p+ B% ], m  }) _% f: J: i
    Y = sim(net,x)
    , b( }, w4 _" U$ {
    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-8-16 16:00 , Processed in 0.410359 second(s), 71 queries .

    回顶部