QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4097|回复: 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 编辑 4 A2 Z7 o1 s6 q/ B1 g- }) w* U

    0 L  |. @# V$ s& w" l: V1 o各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。
    " s6 s! u: U+ ?国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。" Q% Z: _" v% O2 D9 M# j: u1 F& l6 U% P

    + O# ^. A/ E  g' ]/ Q, nclear: K' s; D* [# S
    clc$ i$ o' m2 k' k  l6 r1 N
    %输入数据--以一列作为一个输入数据
    2 Z. ^+ E1 K/ e: M# Px1=[227.5 286.9 236.3 300.6 243.0];
    ; |, s8 W) {) D2 E  y& Ox2=[286.9 236.3 300.6 243.0 354.1];5 f$ i& Q8 J( R. V3 H/ l
    x3=[236.3 300.6 243.0 354.1 294.1];8 ]; }7 Z- K  g8 ^
    x4=[300.6 243.0 354.1 294.1 361.4];5 T9 H5 m8 n" x5 F: P- d- m  w
    x5=[243.0 354.1 294.1 361.4 333.7];
    $ I! r) o8 g+ M' Rx6=[354.1 294.1 361.4 333.7 426.4];  Z; {# K: H. z
    x7=[294.1 361.4 333.7 426.4 434.4];5 i1 p, G; a0 n8 D7 a6 ]9 h* h
    x8=[361.4 333.7 426.4 434.4 449.9];
    ; @9 v0 o0 |' u/ Vx9=[333.7 426.4 434.4 449.9 431.9];% z3 M1 D6 p' z# J
      R& D0 K5 f: i
    P=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';, Z& w9 H/ A' |6 O# O
    %输出数据--以一列作为一个输出数据6 d' L) I0 o1 V/ K$ ]* y
    y1=[354.1 294.1 361.4 333.7 426.4];
    * q7 k2 E, h& Jy2=[294.1 361.4 333.7 426.4 434.4];
    - f6 B* r4 y: Ry3=[361.4 333.7 426.4 434.4 449.9];6 I7 _/ C8 ?/ o2 t, U& ]
    y4=[333.7 426.4 434.4 449.9 431.9];+ W4 f( C3 I, P
    y5=[426.4 434.4 449.9 431.9 454.0];2 p# b# L  U1 c& @& Z
    y6=[434.4 449.9 431.9 454.0 655.6];
    0 ]3 S7 o0 t0 Y, T7 \y7=[449.9 431.9 454.0 655.6 907.0];9 ^* K1 \0 @: H( x) W' [
    y8=[431.9 454.0 655.6 907.0 1112.5];
    # s' J9 @: G' c* Wy9=[454.0 655.6 907.0 1112.5 1401.1];; m1 F: X; a. J# V3 a& |

    1 M$ ]9 |' F& B8 I/ Z/ XT=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';
    ! F& c. N* A6 y* J, z9 i' G
    . C, s) V2 L2 f%归一化  K# f: [# H7 i" H. ~7 `/ b  o7 Y
    % [P,minp,maxp,T,mint,maxt] = premnmx(P,T);0 R/ Q" }' a) |. q0 ^! c% z1 p, {8 e
    %建立网络
    . v4 ^# A) w" ]# p%其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据
    % n5 D+ F3 A0 _5 n( a* A6 K$ I%{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
    & W. w1 i" ^- _: N, X! x%'traingdm' 表示调整权重,缩小误差的数学函数! t$ A: D) `$ z. q, ?
    %↓创建一个新的前向神经网络  8 B; w3 S- Y* h  E% J
    net = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')
    % J: H" M. `& O% r& L2 k2 u3 r3 f, P6 W+ ^/ A: |4 b: t
    %  当前输入层权值和阈值
      {) R# L: L. @6 }inputWeights = net.IW(1,1)
    / j6 ^% C1 M% \& c# \" ainputbias = net.b(1)
    4 ~" l3 i2 i. S3 H! s; ~; ]
    - C# j( d- K$ Q%  当前网络层权值和阈值 6 d% O. {6 x% o$ ~0 b) O
    layerWeights = net.LW(2,1)
    - h5 }. _: v& M! [, s& @6 z) [layerbias = net.b(2)
    # S0 @0 O4 D, B: E# p3 P4 C1 g$ Q- N- m: h+ g
    %  设置训练参数& K  v. w4 E3 V0 \
    net.trainParam.show = 50;
    $ }, R3 ]* o7 s; snet.trainParam.lr = 0.05;4 C  U& D+ i% o( r! Q2 H
    net.trainParam.mc = 0.9;
    + @0 G0 q, x; h3 @- x" p8 fnet.trainParam.epochs = 10000;  %每计算10000次出一次结果
    * O" S' o& Y; q& u. ?, Wnet.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)
    & d' L% V3 k4 F* ]2 _8 h( L! l& w# J- n% V4 z; ]3 i, j3 C- U
    %  调用 TRAINGDM 算法训练 BP 网络* h1 h, \7 E+ w( W
    [net,tr] = train(net,P,T);      3 }* K" u$ Q" w% V" [3 K

    , f, M/ r3 P4 S6 w, s% \4 S%  对 BP 网络进行仿真/ m" \1 d' B; x" O
    A = sim(net,P);                 %A为输出结果
    6 {  r. n# n1 ^; _
    4 D: R6 Z8 s8 v. u5 F+ |%反归一化8 O8 J: m' p1 S/ o  X. v% a- A* A
    %  A = postmnmx(A,mint,maxt);0 |9 T. w! C7 S3 t9 C( ]5 h
    5 W9 e" o- C2 g
    %  计算仿真误差  ; M) J# k8 m1 V4 E/ I
    E = T-A;                        %原来的输出 - 训练输出 = 误差
    0 N6 I' J+ _- E$ f; P! j+ q! hMSE = mse(E)                    %输出误差8 C+ B/ D. h9 b! G3 C6 [/ c

    2 S# a, \2 {% W( y# S' I4 h# i%下面是输入数据,测试结果1 r: W  M; Z1 p2 Z
    x=[454.0 655.6 907.0 1112.5 1401.1]';
    / ?1 X! U' B/ J( o3 T9 g9 sY = sim(net,x)5 ^) g' S3 m; \, 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, 2026-1-15 07:13 , Processed in 0.584887 second(s), 70 queries .

    回顶部