QQ登录

只需要一步,快速开始

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

BP神经网络预测问题

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

2

主题

6

听众

56

积分

升级  53.68%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-9-1 19:44 |只看该作者 |正序浏览
|招呼Ta 关注Ta
以下是我写的一段BP神经网络预测的代码,模型建好后,不知怎样预测下一年的T值,哪位高手来指导下...
; Z' ^& U6 X& z% o2 S- Q6 |%X为原始输入数据" E1 N% h& Y0 p2 l3 F' w7 g) u
P=[727000 297.35  327.53  2658 6.31  68.3 15948 737  16316 583.89 1.96;
# ?. G3 a+ z! D6 O4 Q30619 325.00  390.51  2763 6.94  69.6 17130 789  18600 707.67 2.28;
1 g5 ~) A  g0 j9 J9 {0 u33282 423.00  474.63  2801 5.73  69.7 17866 894  19886 861.88 2.32;9 V8 I: \* V8 Z0 Y2 h3 G; b
33689 467.57  569.55  2887 7.76  69.8 18961 1015 20249 941.99 2.38;: s: c% L# }% F  v/ a6 W
41020 538.17  616.25  2920 8.23  69.7 19688 1198 21626 1082.6 2.73;: ~) {- Q& [: L, ]; o- ^
43344 832.04  686.37  3495 8.87  68.3 24465 1361 23544 1373.4 2.89;  a. g% |" n/ W
46030 941.94  788.15  3495 9.57  68.2 24814 1710 24941 1756.5 2.82;
2 X. d7 V4 B, ?" |% {% w7 B% j53887 1095.13 969.1   4885 9.65  68.7 36149 2100 25936 2199.5 3.08;
! j& [4 Y6 e+ \$ [59271 1250.64 1092.6  5376 10.11 69.2 43008 2314 26596 2625.4 3.56;
1 j+ v8 R$ D* H" o64507 1437.67 1176.1  6091 15.05 69.2 51946 2500 28494 3229.4 3.98;
3 F9 p8 I" K8 Z$ Q- e2 U70597 1671.29 1273.7  7305 16.81 69.2 65745 2614 29628 3744.7 4.98;
, [; [* T/ B1 k, A( W8 g79221 1905.03 1345    8023 18.6  69.1 70346 2798 30063 3792.6 5.1;
1 K+ _$ W# I* D$ E' s78812 2207.88 927.16  8328 20.52 69.3 71088 2964 30904 3951.64 5.21]';& T* r; N% T5 |) [8 g7 A
%为原始目标数据: m* t3 u7 V9 t
T=[3.1,3.6,4.2,4.8,6.7,9.1,13,18.9,29,51.1,78.2,103.5,128];% P$ E9 q" r/ K
%原始数据归一化
% T! k; B  \* O6 B/ u" p- D[pn,minp,maxp]=premnmx(P);
$ C( h% ^  i" B; H/ ~[tn,mint,maxt]=premnmx(T);
. _1 x. p5 j( `$ F%设置网络隐单元的神经元数%建立相应的BP网络
4 c/ N; u% {, Y# Hnet=newff(minmax(pn(:,1:13)),[10,1],{'tansig','purelin'},'trainlm');
# y/ f# J0 H' \* y8 g* X%当前输入层权值和阈值 0 B3 b0 d6 v9 f. F) y2 z
inputWeights=net.IW{1,1}
* D& I: b% G7 l  ~inputbias=net.b{1}   r% G3 e1 s2 ~- Y! P
%当前网络层权值和阈值 / q/ w7 h! T5 c8 V, r7 B( A
layerWeights=net.LW{2,1}
. E5 e5 J( A' k9 \) M0 wlayerbias=net.b{2}. w9 T, [; [: V5 C8 J- E; v
net=init(net);% D: s* A4 M( ]3 Q) l1 M& K
% 网络训练代数: \( Z' U, p3 [* m. Y* Z
net.trainParam.epochs=10000;3 C2 k# z/ Y1 \6 g+ ?
% 网络训练目标误差
  S  Z; ]  e- H; L+ K7 A* i% bnet.trainParam.goal=0.000001;
; c, I( {6 M. `% Z' V" V% 网络训练动量系数
4 a' H1 B/ ~5 [( unet.trainParam.mc = 0.9;   
1 c7 x, v4 T2 Y6 Y% 网络学习速率) y" a; @8 C' s1 v4 I, l! p$ L  T
net.trainParam.lr=0.05;
+ V+ ]2 P% H) ?9 B2 o( @%调用TRAINGDM算法训练BP网络$ R- E; C# L+ F/ g9 a: [' {* q& z& X
net=train(net,pn(:,1:13),tn(1:13));
+ g* v4 E8 {; D  T4 e9 Tt_sim=sim(net,pn);8 J% _: ^. J, }' @" B& A. R
t_sim=postmnmx(t_sim,mint,maxt)" G8 x. a* @6 ?% @' s
E=t_sim-T/ c8 L0 O' H' g3 I1 {
M=sse(E)
" D$ D% s' T2 ~4 _5 n4 y3 r/ SN=mse(E)7 L' `7 k' t4 P( w# R
figure(1)
) ^% [: b: e8 Y& ci=1:1:length(t_sim);1 Z& N4 e) e: j) F
j=1:length(T);
- q( l3 b/ j4 c; Yplot(i+1995,t_sim,'r',j+1995,T,'bo')7 {/ v6 Y) i4 p0 }
legend('模拟数据','原始数据')
7 S1 D; s- v! r& ^/ u* b9 MXlabel('年份');Ylabel('数量')8 b9 @% A: C9 p( f* C+ o  `
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-4-24 01:22 , Processed in 0.528129 second(s), 100 queries .

    回顶部