QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3692|回复: 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 编辑
    ( u1 p; N' w1 H' q& z/ E
    , L. K/ y/ F8 t# s2 }% T各位朋友,有一个问题请教你们。   我用BP神经网络做预测的时候(下面是我的代码)遇到一个问题。  首先我尝试在代码里面不把数据归一化,这样的出来的结果很差(R=0.46).但我把数据作归一化 与 反归一化 之后 确实R提高了很多(R=0.99),这时问题就来了,请看最下面的测试结果代码,那里只有输入(没有输出,因为是预测),premnmx 归一化要有输出的变量T,但我只有输入数据,没有输出数据,因此反归一的结果(最后的预测结果)也无法做到。那么我的预测就不是做不了吗?  很纠结,查了很多资料了,还是不知道怎么解决。& H1 Y$ C9 F0 S# b7 }2 h
    国赛准备到了,希望热心的朋友懂的不吝赐教,我把问题写得很认真,希望各位能帮我解决这个问题,谢谢了。0 c/ ]% c/ D( i0 c3 l( y
    : T: m; d; Q; N! V. e! `
    clear
    . C/ @! D; K% K* ?  |clc
    ) ~9 u4 Z/ e0 s%输入数据--以一列作为一个输入数据" X" B; A, U  a& r. N. N
    x1=[227.5 286.9 236.3 300.6 243.0];  S0 U7 B( a- C1 `; Y9 f$ v/ Q6 e
    x2=[286.9 236.3 300.6 243.0 354.1];6 p' z% c* q$ {3 \8 y/ m" B
    x3=[236.3 300.6 243.0 354.1 294.1];  i# F+ M) d% e# C
    x4=[300.6 243.0 354.1 294.1 361.4];; z" k8 k2 }" y6 B. _* [
    x5=[243.0 354.1 294.1 361.4 333.7];; _6 L5 m) x& O' u  C8 D
    x6=[354.1 294.1 361.4 333.7 426.4];# \% R' R% m  L" Y4 n  z3 E+ P
    x7=[294.1 361.4 333.7 426.4 434.4];9 g' @; t+ D3 y5 Z
    x8=[361.4 333.7 426.4 434.4 449.9];: x2 ]: H+ Y6 j3 K
    x9=[333.7 426.4 434.4 449.9 431.9];
    2 W8 ]# a& W4 z5 {
    / C  A) C' w) ^3 I$ j  E) cP=[x1;x2;x3;x4;x5;x6;x7;x8;x9]';
    # X& t2 ?& A3 Q%输出数据--以一列作为一个输出数据6 R) }2 ]) y+ I( U
    y1=[354.1 294.1 361.4 333.7 426.4];6 F6 {3 P1 `, j5 |
    y2=[294.1 361.4 333.7 426.4 434.4];
      |" }( w& I# H. Gy3=[361.4 333.7 426.4 434.4 449.9];
    2 }; q. |, m/ }" `7 a3 G9 Qy4=[333.7 426.4 434.4 449.9 431.9];0 e9 x. Q8 ]3 b
    y5=[426.4 434.4 449.9 431.9 454.0];0 |; ^6 x7 {6 c8 Y, ^
    y6=[434.4 449.9 431.9 454.0 655.6];
    ' U4 c" e, M  v9 ~8 ]% I' @- sy7=[449.9 431.9 454.0 655.6 907.0];2 J. j* f2 n& r& G
    y8=[431.9 454.0 655.6 907.0 1112.5];
    6 y" q. @8 d; @# f4 Ey9=[454.0 655.6 907.0 1112.5 1401.1];
    / c2 o  t8 L6 l8 j
    / g2 W" M4 X8 OT=[y1;y2;y3;y4;y5;y6;y7;y8;y9]';/ T4 O% S; I* V& E( Y! t0 ^6 j
    , k7 U+ y) c$ m* j$ P5 ?" ]
    %归一化
    3 B) n2 t2 i6 W% [P,minp,maxp,T,mint,maxt] = premnmx(P,T);/ k, V, W% Y% P
    %建立网络1 \) `& [( x) j9 V
    %其中[10,3]代表有一个隐含层,里面有10个神经元;“3”代表有三个输出数据* G' o- u9 P+ S' O5 V# {
    %{'tansig','purelin'}与[10,3]对应,意思为隐含层 与 输出层
    ! }0 ~' R: s9 m! h%'traingdm' 表示调整权重,缩小误差的数学函数+ @1 B* P% }1 x* K& {4 |+ f# ]- ]
    %↓创建一个新的前向神经网络  9 X9 o/ }3 n0 S/ |/ n- E
    net = newff(minmax(P),[6,5],{'tansig','purelin'},'traingdm')7 ]6 S/ n& O8 \; P: J
    ( ], H1 u' B5 L, [% n- {: f0 v
    %  当前输入层权值和阈值
    ' P& [: k3 n- |" f: e7 VinputWeights = net.IW(1,1)7 j! v5 A, F4 c) B" m
    inputbias = net.b(1)7 t+ o  g/ X& _4 r# a2 Y/ I- ^
    5 t' Q/ ]: T; E; f
    %  当前网络层权值和阈值
    7 {/ l# B; B8 i$ k: h/ d5 dlayerWeights = net.LW(2,1)
    ) b& M( P) O# ]! ]layerbias = net.b(2)
    ( l3 i8 m5 v4 x: A5 @* O# R* _) A1 W- {4 W+ f
    %  设置训练参数
    " b6 z, \$ \- _0 h7 a* ~net.trainParam.show = 50;3 h9 m! p6 F+ y7 M' H2 G/ |
    net.trainParam.lr = 0.05;. [  n, y8 Z+ E! V
    net.trainParam.mc = 0.9;% h  E  X0 Q8 C* K3 y5 \
    net.trainParam.epochs = 10000;  %每计算10000次出一次结果, B( j3 A- ^3 m2 @
    net.trainParam.goal = 1e-6;     %期望达到的误差(认为是合适的误差)
    $ k( I  `% T3 f; Z. D: h4 b' T: u4 B5 p9 ^
    %  调用 TRAINGDM 算法训练 BP 网络6 p9 T- H; y+ A6 l* j
    [net,tr] = train(net,P,T);      4 j2 P8 m2 P0 L% I7 d. w  @0 v

    3 n9 f4 W* L8 J%  对 BP 网络进行仿真
    ) _1 K% {7 s: k9 Q0 CA = sim(net,P);                 %A为输出结果4 U5 v- O  }! k  v5 r

    % t5 f- h0 q/ W" M%反归一化8 O, O6 H% F, b9 @# T( n
    %  A = postmnmx(A,mint,maxt);9 {! ~/ H' I- P* X

    1 z* L, ]: T6 T4 r2 [0 G9 u$ f%  计算仿真误差  & P- i( I8 k0 \/ C2 Q8 o$ L, R
    E = T-A;                        %原来的输出 - 训练输出 = 误差
      b: h5 G1 p4 @# }MSE = mse(E)                    %输出误差
    $ M* D2 G+ D3 K: Z% N
    3 W# Y( u1 r! \2 c$ c: f%下面是输入数据,测试结果- b- F; s  Z) @: e1 q8 m- H6 |  Z
    x=[454.0 655.6 907.0 1112.5 1401.1]';- U6 T5 F/ l+ B5 G8 {
    Y = sim(net,x)% D0 Q. F. @4 Y# q: B5 e
    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-6-30 02:24 , Processed in 0.711010 second(s), 71 queries .

    回顶部