QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 9043|回复: 18
打印 上一主题 下一主题

BP神经网络预测问题

[复制链接]
字体大小: 正常 放大
jthawking        

2

主题

6

听众

56

积分

升级  53.68%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-9-1 19:44 |只看该作者 |正序浏览
|招呼Ta 关注Ta
以下是我写的一段BP神经网络预测的代码,模型建好后,不知怎样预测下一年的T值,哪位高手来指导下...; f" }# w1 U  K
%X为原始输入数据
( H, J+ ]- l: ~6 u' X2 Z  u0 ?1 wP=[727000 297.35  327.53  2658 6.31  68.3 15948 737  16316 583.89 1.96;
  L4 g1 {- Z9 a: m6 W. V5 X& N30619 325.00  390.51  2763 6.94  69.6 17130 789  18600 707.67 2.28;
1 U1 R, f! n9 G" C33282 423.00  474.63  2801 5.73  69.7 17866 894  19886 861.88 2.32;9 r$ U' O. e0 l; x* b* d) _
33689 467.57  569.55  2887 7.76  69.8 18961 1015 20249 941.99 2.38;0 q7 k/ i2 ^( r& D
41020 538.17  616.25  2920 8.23  69.7 19688 1198 21626 1082.6 2.73;4 @1 I# U+ J5 Q# O. v
43344 832.04  686.37  3495 8.87  68.3 24465 1361 23544 1373.4 2.89;
0 L" S3 Z3 u6 L3 _: l5 m46030 941.94  788.15  3495 9.57  68.2 24814 1710 24941 1756.5 2.82;
* f- u! I( U2 S$ Y, U9 J. A53887 1095.13 969.1   4885 9.65  68.7 36149 2100 25936 2199.5 3.08;
2 v  m. l! c5 ?! J: ?59271 1250.64 1092.6  5376 10.11 69.2 43008 2314 26596 2625.4 3.56;2 A# c6 I4 B$ z: i9 X, l
64507 1437.67 1176.1  6091 15.05 69.2 51946 2500 28494 3229.4 3.98;
: p" V. J  x5 R9 o2 Y70597 1671.29 1273.7  7305 16.81 69.2 65745 2614 29628 3744.7 4.98;5 c5 B, _8 X" ]' v( j
79221 1905.03 1345    8023 18.6  69.1 70346 2798 30063 3792.6 5.1;, k) c2 q/ D: E' a9 p/ C! d
78812 2207.88 927.16  8328 20.52 69.3 71088 2964 30904 3951.64 5.21]';
: y, g4 t% s' M% r0 R) M%为原始目标数据
0 F+ ?3 B2 k" @) l% oT=[3.1,3.6,4.2,4.8,6.7,9.1,13,18.9,29,51.1,78.2,103.5,128];- g5 E$ o# K5 y) N
%原始数据归一化
3 l+ F- x+ K# l7 y[pn,minp,maxp]=premnmx(P);" a. ^2 c- o# L. r
[tn,mint,maxt]=premnmx(T);
: E4 n: {, S  S  o%设置网络隐单元的神经元数%建立相应的BP网络/ d9 {1 O1 {) t5 t) ^; z  p
net=newff(minmax(pn(:,1:13)),[10,1],{'tansig','purelin'},'trainlm');  D" m0 C- g3 g0 Y
%当前输入层权值和阈值
( Q5 ?) J$ O! D- v5 r& BinputWeights=net.IW{1,1}
- Z: n7 [9 x; F5 \1 einputbias=net.b{1}
4 n3 [2 k0 v$ Q%当前网络层权值和阈值 4 w! ^9 ], i8 k6 p
layerWeights=net.LW{2,1} ' V8 s" p$ |1 C. g1 }# n
layerbias=net.b{2}6 c6 w$ Z/ s' E* k9 E& F8 z; @
net=init(net);
5 \' s3 \5 `3 F: _) ?9 D% 网络训练代数
7 M; K* V2 ^( k5 X/ h; l1 Xnet.trainParam.epochs=10000;
1 b! h  O; W  j7 F1 j0 k. |% 网络训练目标误差
. v* A& K/ v7 j$ ]( n) onet.trainParam.goal=0.000001;
6 o8 }! S3 A, G, k6 t/ k; ?% 网络训练动量系数
  J& L) f/ [2 |% m$ ]: r3 ^net.trainParam.mc = 0.9;   , C- |  c; s# p% D7 a. f
% 网络学习速率" k! c/ [% R  w' t
net.trainParam.lr=0.05;# [: j4 a" k# @7 Y  X
%调用TRAINGDM算法训练BP网络
: c% c9 i$ P( L( R' X  ^; c7 ynet=train(net,pn(:,1:13),tn(1:13));. x+ y) [; V; c
t_sim=sim(net,pn);
) T+ F( ^0 Y! s+ o% ot_sim=postmnmx(t_sim,mint,maxt)' x  t; \/ \" @
E=t_sim-T! A$ T4 B* ~  N* \# g# w8 g
M=sse(E)4 T" L3 _9 {( v# E* y# O5 `* A
N=mse(E)- O3 ^) r2 l/ X$ g: |3 x$ z& m
figure(1)8 h7 s1 a8 H* {0 f4 k" x$ o
i=1:1:length(t_sim);
# r. l5 e, n" t4 Y% J6 mj=1:length(T);
: o! g, @+ F# b) f7 N) v# ~# ^plot(i+1995,t_sim,'r',j+1995,T,'bo')" h- L8 {; R0 n5 P3 E' M9 {
legend('模拟数据','原始数据')
7 x. _% h- M9 r. @3 k& bXlabel('年份');Ylabel('数量')+ i! t7 T/ M+ T
title('BP神经网络')
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
杨辰宇        

6

主题

13

听众

71

积分

升级  69.47%

  • TA的每日心情
    开心
    2018-2-7 22:53
  • 签到天数: 52 天

    [LV.5]常住居民I

    自我介绍
    正在参加数学建模比赛

    群组2016数学建模算法集锦

    群组2017国赛备战群组

    回复

    使用道具 举报

    19

    主题

    11

    听众

    272

    积分

    升级  86%

  • TA的每日心情
    郁闷
    2014-10-23 09:29
  • 签到天数: 59 天

    [LV.5]常住居民I

    楼主的问题最后修改完了么?能不能发一份代码学习学习,现在有也遇到同样的问题了,不知道怎么解决了,求帮助。1345134018@qq.com
    回复

    使用道具 举报

    月关        

    0

    主题

    10

    听众

    249

    积分

    升级  74.5%

  • TA的每日心情
    擦汗
    2015-5-23 18:57
  • 签到天数: 84 天

    [LV.6]常住居民II

    自我介绍
    新手,请多多关照!!!

    社区QQ达人

    回复

    使用道具 举报

    月关        

    0

    主题

    10

    听众

    249

    积分

    升级  74.5%

  • TA的每日心情
    擦汗
    2015-5-23 18:57
  • 签到天数: 84 天

    [LV.6]常住居民II

    自我介绍
    新手,请多多关照!!!

    社区QQ达人

    回复

    使用道具 举报

    月关        

    0

    主题

    10

    听众

    249

    积分

    升级  74.5%

  • TA的每日心情
    擦汗
    2015-5-23 18:57
  • 签到天数: 84 天

    [LV.6]常住居民II

    自我介绍
    新手,请多多关照!!!

    社区QQ达人

    回复

    使用道具 举报

    月关        

    0

    主题

    10

    听众

    249

    积分

    升级  74.5%

  • TA的每日心情
    擦汗
    2015-5-23 18:57
  • 签到天数: 84 天

    [LV.6]常住居民II

    自我介绍
    新手,请多多关照!!!

    社区QQ达人

    回复

    使用道具 举报

    19

    主题

    11

    听众

    272

    积分

    升级  86%

  • TA的每日心情
    郁闷
    2014-10-23 09:29
  • 签到天数: 59 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    1

    主题

    5

    听众

    163

    积分

    升级  31.5%

  • TA的每日心情
    无聊
    2013-4-5 21:32
  • 签到天数: 51 天

    [LV.5]常住居民I

    自我介绍
    数模新手
    回复

    使用道具 举报

    3

    主题

    3

    听众

    161

    积分

    升级  30.5%

  • TA的每日心情
    奋斗
    2013-9-20 20:37
  • 签到天数: 33 天

    [LV.5]常住居民I

    新人进步奖

    群组A题讨论群

    群组B题讨论群

    群组C题讨论群

    群组D题讨论群

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-12 14:25 , Processed in 0.486137 second(s), 100 queries .

    回顶部