QQ登录

只需要一步,快速开始

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

BP神经网络预测问题

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

2

主题

6

听众

56

积分

升级  53.68%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-9-1 19:44 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
以下是我写的一段BP神经网络预测的代码,模型建好后,不知怎样预测下一年的T值,哪位高手来指导下...5 Y, \9 y$ f' H
%X为原始输入数据/ J4 l2 o5 r4 w0 d6 A. o
P=[727000 297.35  327.53  2658 6.31  68.3 15948 737  16316 583.89 1.96;) g% f1 P, p; ]1 Q$ `" q2 v0 T
30619 325.00  390.51  2763 6.94  69.6 17130 789  18600 707.67 2.28;
% z! ]7 R* y# M33282 423.00  474.63  2801 5.73  69.7 17866 894  19886 861.88 2.32;
* l7 l- j8 A" W% A% t33689 467.57  569.55  2887 7.76  69.8 18961 1015 20249 941.99 2.38;
- S3 m' T/ G2 |: G4 q; g41020 538.17  616.25  2920 8.23  69.7 19688 1198 21626 1082.6 2.73;0 o+ F% ]1 C* m
43344 832.04  686.37  3495 8.87  68.3 24465 1361 23544 1373.4 2.89;
! G  D& `# t2 o6 p46030 941.94  788.15  3495 9.57  68.2 24814 1710 24941 1756.5 2.82;+ {$ [0 k9 V- D' q* F  [
53887 1095.13 969.1   4885 9.65  68.7 36149 2100 25936 2199.5 3.08;; j( x" z: v$ Q. C7 w
59271 1250.64 1092.6  5376 10.11 69.2 43008 2314 26596 2625.4 3.56;. A* o5 F& |( L4 \1 W" _" i
64507 1437.67 1176.1  6091 15.05 69.2 51946 2500 28494 3229.4 3.98;, ]0 s& b* B* a4 P. ~
70597 1671.29 1273.7  7305 16.81 69.2 65745 2614 29628 3744.7 4.98;
' y1 [# }" z  |9 g% V6 X! F" w- t79221 1905.03 1345    8023 18.6  69.1 70346 2798 30063 3792.6 5.1;, o! P+ Q; Y6 Y. w8 S/ M, k
78812 2207.88 927.16  8328 20.52 69.3 71088 2964 30904 3951.64 5.21]';: O) j# J- i0 ~8 ^+ J" z( J
%为原始目标数据0 @8 ~1 m% a2 [- d9 I/ I
T=[3.1,3.6,4.2,4.8,6.7,9.1,13,18.9,29,51.1,78.2,103.5,128];) O! y+ y+ S6 z- s
%原始数据归一化 0 U6 |* T9 y- ~2 ]/ J
[pn,minp,maxp]=premnmx(P);
8 L, h2 B- f9 ?8 K% d# u) a; A  d0 w/ S[tn,mint,maxt]=premnmx(T);  c! f# l* N% \; b8 V  ]9 ?$ U( Z  I
%设置网络隐单元的神经元数%建立相应的BP网络
, {3 c# b7 h( R: E0 Xnet=newff(minmax(pn(:,1:13)),[10,1],{'tansig','purelin'},'trainlm');
0 |* g+ U# T! d( \%当前输入层权值和阈值 ' \2 n* I* _0 n; u
inputWeights=net.IW{1,1}
- n9 ~0 r) Z4 q3 `inputbias=net.b{1}
2 ]0 y- ~+ y5 F) M6 k" R%当前网络层权值和阈值
5 \* u# q' |4 k: _) ?  klayerWeights=net.LW{2,1}
5 P2 a0 |" p9 A6 [% Wlayerbias=net.b{2}
0 ~5 b4 f$ V  e' x: inet=init(net);
/ N9 W9 t: R, Q1 C" O% 网络训练代数3 P; ?) p: S% {# D3 j- D
net.trainParam.epochs=10000;! y/ g( T: _' M9 E
% 网络训练目标误差
% k' m& r9 g) ?net.trainParam.goal=0.000001;# |  v7 {% x2 S3 G8 D$ k6 P
% 网络训练动量系数
$ m# K% f. c/ V* l: bnet.trainParam.mc = 0.9;   
; r$ P+ I, {+ n% 网络学习速率
* N( O- Q# G' I$ tnet.trainParam.lr=0.05;
4 s; B$ W; n4 [# G8 k3 ~7 _+ X%调用TRAINGDM算法训练BP网络
: L" Y, ^& r; @) K( nnet=train(net,pn(:,1:13),tn(1:13));  U. k- d& \: N3 u/ l  l* I
t_sim=sim(net,pn);
0 Q; [. ~! z' n7 x9 W: Ft_sim=postmnmx(t_sim,mint,maxt)
) X; T+ d  P. WE=t_sim-T; |4 P4 ^" y, T0 H- v% ^2 j
M=sse(E)/ Z% o' T9 _/ B9 @( J' _* b7 T- D
N=mse(E)
+ {' E& ^# X1 G) `0 f8 o8 O  Ffigure(1)
% P1 t3 j% ]+ f7 f$ xi=1:1:length(t_sim);# g  L/ X# e2 K4 N; c2 p" }% G
j=1:length(T);% ^: A  I4 ~( U; P" _
plot(i+1995,t_sim,'r',j+1995,T,'bo')
% @' V" [/ C4 f0 Olegend('模拟数据','原始数据')1 {5 g& `: ]5 v" N6 S
Xlabel('年份');Ylabel('数量')' U& }% l% {0 `$ o0 q) |7 M
title('BP神经网络')
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
追梦者 实名认证       

13

主题

4

听众

594

积分

升级  98%

该用户从未签到

国际赛参赛者

新人进步奖

群组数模应用

群组LINGO

群组四川农业大学数学建模协会

群组重庆交通大学数学建模协会

群组Mathematica研究小组

回复

使用道具 举报

Rosierl        

0

主题

4

听众

69

积分

升级  67.37%

该用户从未签到

回复

使用道具 举报

A123J        

7

主题

4

听众

104

积分

升级  2%

该用户从未签到

回复

使用道具 举报

泽泽 实名认证       

4

主题

4

听众

533

积分

升级  77.67%

  • TA的每日心情
    开心
    2011-6-2 09:51
  • 签到天数: 1 天

    [LV.1]初来乍到

    新人进步奖

    回复

    使用道具 举报

    4

    主题

    4

    听众

    289

    积分

    升级  94.5%

  • TA的每日心情
    奋斗
    2013-1-30 16:57
  • 签到天数: 12 天

    [LV.3]偶尔看看II

    自我介绍
    有梦,便去追!!!!!!!!

    群组加油数模,加油美国赛。

    回复

    使用道具 举报

    1

    主题

    4

    听众

    169

    积分

    升级  34.5%

  • TA的每日心情
    开心
    2012-6-21 21:19
  • 签到天数: 13 天

    [LV.3]偶尔看看II

    群组数学建模培训课堂1

    回复

    使用道具 举报

    0

    主题

    4

    听众

    17

    积分

    升级  12.63%

  • TA的每日心情
    开心
    2012-2-24 09:05
  • 签到天数: 1 天

    [LV.1]初来乍到

    自我介绍
    说什么捏才是对的!
    回复

    使用道具 举报

    0йly_        

    2

    主题

    5

    听众

    20

    积分

    升级  15.79%

    该用户从未签到

    自我介绍
    学生

    群组学术交流A

    如果你要预测的话,你可以先输入一个预测的值,P_test=[];out=sim(net,P_test')~~希望对你有用
    回复

    使用道具 举报

    0

    主题

    7

    听众

    435

    积分

    升级  45%

  • TA的每日心情

    2014-10-29 08:07
  • 签到天数: 96 天

    [LV.6]常住居民II

    群组Matlab讨论组

    群组数学趣味、游戏、IQ等

    群组交朋友

    群组中国科学院软件研究所

    群组学术交流A

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-31 01:50 , Processed in 0.732323 second(s), 104 queries .

    回顶部