QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4106|回复: 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 编辑
    8 s3 X9 k" o! d% _0 Q. |
    , b3 U. Z9 P& p& F各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。
    ; p; h$ ^2 n7 S9 F( ?+ Z国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。. v" `' y0 g# j5 B
    # s. d( \/ [" R( R+ w
    clear* G5 @4 O/ F* R$ o) e5 x* o
    clc
    $ a' m$ x. L' \- s$ h, _8 I, D%输入数据--以一列作为一个输入数据
    5 b: l( D* [) V, Bx1=[227.5 286.9 236.3 300.6 243.0];
    ) T: E) l* t. t+ o$ J% ux2=[286.9 236.3 300.6 243.0 354.1];
    5 F. K  @# R5 T! \# u/ ?1 l1 Ux3=[236.3 300.6 243.0 354.1 294.1];
    ' p1 F' L3 t7 E5 A- K) |x4=[300.6 243.0 354.1 294.1 361.4];4 K, k/ ^7 t; S& _, o
    x5=[243.0 354.1 294.1 361.4 333.7];
    5 f3 i' \" A  @x6=[354.1 294.1 361.4 333.7 426.4];
    # s# _5 Y) g. h9 i3 Ux7=[294.1 361.4 333.7 426.4 434.4];
    ) F' i& c7 J8 V4 n% `x8=[361.4 333.7 426.4 434.4 449.9];
    + j7 T( e8 N, ~9 u$ ?x9=[333.7 426.4 434.4 449.9 431.9];" H" X; Y' ~) b

    . @% B- I1 c, G$ j% _9 kP=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';
    - r4 c$ {( F. c%输出数据--以一列作为一个输出数据0 ]7 D) ^# O! F* ]* V
    y1=[354.1 294.1 361.4 333.7 426.4];; Y1 D1 Z5 p+ x  `. [
    y2=[294.1 361.4 333.7 426.4 434.4];
    # `6 ?, ?. ^( W; [5 ~) L7 Ly3=[361.4 333.7 426.4 434.4 449.9];0 R' I% [- \( F# \+ I
    y4=[333.7 426.4 434.4 449.9 431.9];( _: o/ S7 a% f# j! s. s7 z
    y5=[426.4 434.4 449.9 431.9 454.0];
    9 p9 V( \& a0 \y6=[434.4 449.9 431.9 454.0 655.6];
    & f6 Y0 Y) \* G& z& hy7=[449.9 431.9 454.0 655.6 907.0];' N& y, X2 y2 r. f9 X
    y8=[431.9 454.0 655.6 907.0 1112.5];3 h6 T+ b( J9 a/ `- i
    y9=[454.0 655.6 907.0 1112.5 1401.1];
    : S% B7 f9 p* H5 t, {9 o3 y; v) G
    / V0 S; T) b; U1 R8 t# V% l3 l6 ET=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';/ g$ t! A, L! B) r
    + U, w$ q- o/ U. P1 X
    %归一化
    ! q; l: u1 O* @# S2 l- h% [P,minp,maxp,T,mint,maxt] = premnmx(P,T);
    % I& X9 z: i* g6 U% _9 h%建立网络
    , h1 B. k' B' L# j; a%其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据
      r$ ^  b* k# H0 a%{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
    / d7 y& ~# G* P%'traingdm' 表示调整权重,缩小误差的数学函数
    - p7 o( Y; I* N+ s%↓创建一个新的前向神经网络  
    5 ]  m6 \3 h7 g+ Unet = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')
    . F8 @  V) \0 X) {+ ?* @- ~
    & I. L/ V& H. S, {5 N%  当前输入层权值和阈值
    + N6 c2 n4 u& n& P  o6 binputWeights = net.IW(1,1)
    / G5 ~5 x/ |% |. [& xinputbias = net.b(1)
    / K# d: V3 n5 u& |" l7 o( r5 j) S/ J& z
    %  当前网络层权值和阈值
      W4 z3 S# O$ c/ F/ Y6 M1 hlayerWeights = net.LW(2,1), y2 K, P, E1 Y
    layerbias = net.b(2); V0 Q7 G0 V1 \8 n' L

    $ h5 k& u! P. Y% r. Q( i%  设置训练参数
    ( N7 T9 Z# [' i9 Pnet.trainParam.show = 50;/ e4 ?+ T% w4 j' j" N2 R% S, k8 z
    net.trainParam.lr = 0.05;+ ?/ g( X" E: n! n! _
    net.trainParam.mc = 0.9;! T) t* g+ W; f2 O4 t" x, k5 I
    net.trainParam.epochs = 10000;  %每计算10000次出一次结果6 M/ o9 G" j& n. `* r6 S% f7 K9 O
    net.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)9 D( }7 j" p+ u7 H# H

    $ A% E! U3 W; U$ k3 c%  调用 TRAINGDM 算法训练 BP 网络
    - l( q, v8 C( E/ ^4 e- o0 R5 m! y[net,tr] = train(net,P,T);      
      p1 y% L4 ?& ~7 y2 P3 P$ _8 c+ x7 C" T- g& O* p6 y5 A7 }
    %  对 BP 网络进行仿真+ Y% B. m- U! P4 q: V
    A = sim(net,P);                 %A为输出结果$ X7 a/ m4 i1 @! y9 k- S
    ; ^& ]4 `9 N9 d3 c4 A+ Q5 l
    %反归一化* V# A' O( G: j0 J  ]" O1 t3 Y
    %  A = postmnmx(A,mint,maxt);
    9 [0 {/ g: w7 G& a
    5 ]! O  z5 q' o# h$ T# f%  计算仿真误差  
    6 q/ ~* q5 y4 C/ s$ M' KE = T-A;                        %原来的输出 - 训练输出 = 误差; c/ Y" n8 n: m6 P; w  D% t6 C$ }
    MSE = mse(E)                    %输出误差
    9 Q& V/ p) o% V$ I8 _
    . _; g: I2 W' I, |+ Q  @9 a4 U%下面是输入数据,测试结果
    9 q) \' ~) R* f6 X, b" D" Q1 xx=[454.0 655.6 907.0 1112.5 1401.1]';) G( i; l" |6 |
    Y = sim(net,x)
    ! x- E2 w' y' c, G; Y+ `8 U! Q
    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-15 08:24 , Processed in 0.821852 second(s), 71 queries .

    回顶部