QQ登录

只需要一步,快速开始

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

BP神经网络预测问题

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

2

主题

6

听众

56

积分

升级  53.68%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-9-1 19:44 |只看该作者 |正序浏览
|招呼Ta 关注Ta
以下是我写的一段BP神经网络预测的代码,模型建好后,不知怎样预测下一年的T值,哪位高手来指导下...
" Z4 v% @$ Y& L. L+ r9 w& F! X%X为原始输入数据* G# @7 C2 h, l% {4 O
P=[727000 297.35  327.53  2658 6.31  68.3 15948 737  16316 583.89 1.96;: {2 s8 r, [+ @+ D
30619 325.00  390.51  2763 6.94  69.6 17130 789  18600 707.67 2.28;
1 k) C7 b. N5 z! W: o5 ^4 Z. J8 K1 C33282 423.00  474.63  2801 5.73  69.7 17866 894  19886 861.88 2.32;+ A) s2 a. o  H" j( g8 z4 r
33689 467.57  569.55  2887 7.76  69.8 18961 1015 20249 941.99 2.38;
# t. c$ P! U/ ~% i+ y7 C41020 538.17  616.25  2920 8.23  69.7 19688 1198 21626 1082.6 2.73;
8 q: x- h. P. b4 K43344 832.04  686.37  3495 8.87  68.3 24465 1361 23544 1373.4 2.89;
3 X) Q7 x! L! J, I1 H$ _) G6 }46030 941.94  788.15  3495 9.57  68.2 24814 1710 24941 1756.5 2.82;
3 ^: ^- i$ p9 Y# `53887 1095.13 969.1   4885 9.65  68.7 36149 2100 25936 2199.5 3.08;
& Q1 f6 l1 ]+ r! _' A3 q59271 1250.64 1092.6  5376 10.11 69.2 43008 2314 26596 2625.4 3.56;
5 R# o! H# |6 y8 Z# R64507 1437.67 1176.1  6091 15.05 69.2 51946 2500 28494 3229.4 3.98;% o( y* T5 b9 J7 s0 t
70597 1671.29 1273.7  7305 16.81 69.2 65745 2614 29628 3744.7 4.98;
: }; a- ]3 g) a6 r79221 1905.03 1345    8023 18.6  69.1 70346 2798 30063 3792.6 5.1;. q! `  G, t# A1 x: d7 Y2 O
78812 2207.88 927.16  8328 20.52 69.3 71088 2964 30904 3951.64 5.21]';
# u8 W1 g5 P( _- t" o7 C7 N%为原始目标数据1 Q% Q1 o9 Q* c7 {% y6 e. P
T=[3.1,3.6,4.2,4.8,6.7,9.1,13,18.9,29,51.1,78.2,103.5,128];6 e: C: }4 m$ A* W" P
%原始数据归一化
6 b8 t  e) f) ~: o* b" S% i7 P[pn,minp,maxp]=premnmx(P);
0 Z- G9 T" N. D. e' n5 [" D[tn,mint,maxt]=premnmx(T);
' e. T$ W+ F6 @6 j1 z% s, X. {%设置网络隐单元的神经元数%建立相应的BP网络. p( T) d3 A* N  d* o
net=newff(minmax(pn(:,1:13)),[10,1],{'tansig','purelin'},'trainlm');8 w1 z( E" U6 ]( \/ w3 a
%当前输入层权值和阈值 + P" R* o) Y$ \: P
inputWeights=net.IW{1,1}
2 J- S1 Y1 X! z, i7 Vinputbias=net.b{1}
/ M+ M( F3 K) ~6 D, ^2 G# C# d%当前网络层权值和阈值 ! V7 Q; I& `: c; O' @7 O1 s8 E
layerWeights=net.LW{2,1} 9 d, O, K, s* ^7 _0 Y2 u
layerbias=net.b{2}
0 w' @8 f+ I% _* k& {7 Gnet=init(net);1 F3 Y2 I. x- U! J
% 网络训练代数
4 b, W8 F2 D! W# _' l3 E. d+ Nnet.trainParam.epochs=10000;  Z6 n& O0 p! r
% 网络训练目标误差' d2 L+ ^1 [0 V; @3 b3 g9 V& ^; s- J# x8 j
net.trainParam.goal=0.000001;1 u( e# j" q% b4 O) y3 B
% 网络训练动量系数
7 j. h& Y, H4 \  dnet.trainParam.mc = 0.9;   
: h0 ^0 v9 G9 E) E% 网络学习速率( {0 u  k2 _' e+ ~& d% D: v
net.trainParam.lr=0.05;
7 x3 \& Y- r- F. U; q2 l. E7 X%调用TRAINGDM算法训练BP网络
/ D8 \$ e) N1 B$ b' Cnet=train(net,pn(:,1:13),tn(1:13));
# \# K" p! E9 p& wt_sim=sim(net,pn);
! n" v) H8 [- J7 z& E1 v4 }t_sim=postmnmx(t_sim,mint,maxt)
. J7 g8 b6 }2 U6 sE=t_sim-T% G9 S- @* r) K) H7 t% W
M=sse(E)
4 b6 n( o' p2 Z1 I( AN=mse(E)
# |1 E: b7 S' i+ G) Xfigure(1)$ W& C: \' t4 e  C4 _; n# v
i=1:1:length(t_sim);( h9 S4 [4 C4 j5 F8 s7 u, s
j=1:length(T);5 t! I" L+ r, A! {
plot(i+1995,t_sim,'r',j+1995,T,'bo')" U( _5 b- |6 h+ R3 V" Z% m) `3 u
legend('模拟数据','原始数据')* ]4 `' w* x1 z' j" u5 ?
Xlabel('年份');Ylabel('数量')
9 D+ r& W, n* m+ J: B7 p2 [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-7-31 12:38 , Processed in 0.582215 second(s), 100 queries .

    回顶部