QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4221|回复: 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 编辑 5 l4 |% J+ a+ j' _
    7 A0 I0 l+ }- D3 o4 j7 a
    各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。/ Y- r% W' u+ S* B5 x% I; H
    国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。
    1 J* j. d; P* T& w' p8 X4 g) r: [; N" r- e
    clear
    , G4 l5 V: e8 j+ sclc
    9 {4 p. }2 s0 x' }+ {2 M%输入数据--以一列作为一个输入数据
    , w4 ]) N. Z6 z0 |2 Y% V6 n. r0 Qx1=[227.5 286.9 236.3 300.6 243.0];
    3 W/ A% b. P# W4 V1 ]2 w/ e- lx2=[286.9 236.3 300.6 243.0 354.1];7 G' A: b9 Y3 [" p4 c" N- c( y
    x3=[236.3 300.6 243.0 354.1 294.1];
    $ z" L3 o( Y1 @x4=[300.6 243.0 354.1 294.1 361.4];
    7 j8 G3 G* z+ h/ |: Zx5=[243.0 354.1 294.1 361.4 333.7];
    & ^% h8 o/ E# xx6=[354.1 294.1 361.4 333.7 426.4];
    . C# g2 H) |  cx7=[294.1 361.4 333.7 426.4 434.4];; R$ R  L2 ^5 g; J0 U" U
    x8=[361.4 333.7 426.4 434.4 449.9];  P  W# t3 C# M6 [% W9 N- y
    x9=[333.7 426.4 434.4 449.9 431.9];
    # T3 o  z* i* V1 [' a
    0 V" j& V' a8 E7 ~P=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';( u; w/ V, a& ?2 K) _
    %输出数据--以一列作为一个输出数据
    2 Q/ @+ j: a2 V0 J0 ]8 s  v1 s: Gy1=[354.1 294.1 361.4 333.7 426.4];
    ; V$ |0 `/ @4 h+ [. }& dy2=[294.1 361.4 333.7 426.4 434.4];& ^. h: m" \* O
    y3=[361.4 333.7 426.4 434.4 449.9];
    8 }9 P! ~6 ?8 M- K0 Yy4=[333.7 426.4 434.4 449.9 431.9];
    4 @3 ]5 H% q& ey5=[426.4 434.4 449.9 431.9 454.0];
    " o- A9 V1 d) K% Z$ i4 J. Zy6=[434.4 449.9 431.9 454.0 655.6];
    ' O4 u$ |, x, Ny7=[449.9 431.9 454.0 655.6 907.0];' Q  Q( m' J. b% n
    y8=[431.9 454.0 655.6 907.0 1112.5];/ c6 J8 u5 s* t9 |' P  K" h  U
    y9=[454.0 655.6 907.0 1112.5 1401.1];  r/ A, H" _+ ^# A
    $ ]. N+ U9 G6 u1 s
    T=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
    ; M' D6 d+ n" C5 _" E
    , g3 x  t- s# I8 x" s2 Z/ P; N%归一化
    9 y" \) L' U- R  ~% O% [P,minp,maxp,T,mint,maxt] = premnmx(P,T);
    ) l& L, x4 u1 ^# F' D%建立网络* r; D, `! }4 @6 R+ K( L  s1 o
    %其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据1 z  U) I6 L2 S6 v3 `$ O* p- D
    %{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
    ' r! j4 p, k# w%'traingdm' 表示调整权重,缩小误差的数学函数
    ' w8 A* s# |) o: O%↓创建一个新的前向神经网络  2 {: `+ q# |5 h' I# W
    net = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')$ }- {) O3 g9 @, V" ?; s
    0 N/ U. ~* `1 }7 U& ~1 C
    %  当前输入层权值和阈值 - e4 X2 I# n' ~" I( y
    inputWeights = net.IW(1,1)8 I, O" H" D$ D
    inputbias = net.b(1)% h. \  Z8 j7 H/ ?) Z# b6 r
    4 O% V% I3 a2 s: L/ z! b3 U8 o
    %  当前网络层权值和阈值
    ( k: a' m! p2 tlayerWeights = net.LW(2,1)
    ) J8 e! M4 S/ Xlayerbias = net.b(2)
    4 T: e8 a; ~4 E4 j$ j) X5 L
    % ^$ s2 J. P" n/ [9 k8 ~% e+ u7 E% W%  设置训练参数& u/ ^+ j/ w  y* z3 U5 P
    net.trainParam.show = 50;
    2 l) e" |$ ~2 F2 N) rnet.trainParam.lr = 0.05;  k7 i2 R3 C2 x7 C9 E7 p
    net.trainParam.mc = 0.9;) D8 F2 h3 x+ A
    net.trainParam.epochs = 10000;  %每计算10000次出一次结果
    / F; n9 n( m1 I0 |net.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)
    . ?& G! L& m" J+ a' W7 t* q4 y! Z( z% K
    %  调用 TRAINGDM 算法训练 BP 网络( D& O6 b% g7 f! I3 m
    [net,tr] = train(net,P,T);      
    5 N9 r0 W  X2 a! `1 v: [. z8 r, B& b* w5 k4 O
    %  对 BP 网络进行仿真
    $ s8 x( q# O9 }6 L. I. YA = sim(net,P);                 %A为输出结果
    ' k2 {* T# j) |& M- A* U' g4 n" {% ^" S0 D, T3 G+ A# j
    %反归一化
    9 k3 r5 a: U6 `4 D6 |6 Q%  A = postmnmx(A,mint,maxt);- r3 o- l1 v" @
    # _) f, Y# J. r% N- t5 D$ h3 M
    %  计算仿真误差  7 b" z& |  a6 i0 t
    E = T-A;                        %原来的输出 - 训练输出 = 误差4 P; j% {8 Y/ G' M- U  O# k
    MSE = mse(E)                    %输出误差
    ( }+ `) C/ A( a5 `2 ~8 w# e4 R) K5 A) k) c4 ^- Q# l5 r
    %下面是输入数据,测试结果, H& i/ W$ C! P6 U( j+ z, `
    x=[454.0 655.6 907.0 1112.5 1401.1]';
    . W, D& X3 u5 r1 zY = sim(net,x)
    ( Z( G$ y# [! @$ e8 Z$ C. I
    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-1 20:05 , Processed in 0.492797 second(s), 72 queries .

    回顶部