- 在线时间
- 3 小时
- 最后登录
- 2013-5-18
- 注册时间
- 2008-12-16
- 听众数
- 3
- 收听数
- 0
- 能力
- 0 分
- 体力
- 40 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 16
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 16
- 主题
- 3
- 精华
- 0
- 分享
- 0
- 好友
- 2
升级   11.58% TA的每日心情 | 开心 2013-5-18 14:52 |
---|
签到天数: 1 天 [LV.1]初来乍到
- 自我介绍
- 开朗 活泼
 |
本帖最后由 hanyunxiang 于 2010-5-25 17:03 编辑
' }) a& G7 O) ~+ k7 {* ^8 Y4 U. K9 P$ A; q7 z! X/ B
关于BP神经网络预测
1 d5 W' t0 H+ r1 W. E3 [* {
! `( a$ |4 N* z% f& }+ Y) t8 u注:里面的笑脸代表冒号(不知道会出现笑脸)
- {- A' W( t0 B8 ?3 @
v* }* i" Z E" l3 `* E1 [2 j( A. P. C* N& V
关于BP时间序列预测
, D5 V- @7 m( h% ]8 {# M程序中x是已知的运输量数值(2000年-2009年),以4个数据为一组,前三个数据作为输入,第四个作为输出,共分7组,前六组训练,隐层节点为7,最后一组检验:
E' v. m( W- S5 d0 T+ f: [3 F. ex=[13369.2 14873.7 17137.3 17432.5 24193.5 28435.1....* d8 m, q6 |& N4 ?9 s
33197.3 38758.6 40576.2 48606.3];2 M S* k E, T, P" [4 _' h% B
X(1, =(x(1, -min(x(1, ))/(max(x(1, )-min(x(1, ));(归一化)
4 O& i6 [6 @2 V' w& \/ Wfor i=1:71 n# ~8 \* t' y! ]4 I& Q' k
X1(i, =[X(i) X(i+1) X(i+2) X(i+3)];- b5 _# M3 y: [( P( {& D
end
8 C* n- [. v' f0 d% S0 HX2=(X1(:,1:3))';' r0 Y9 Y& L% x: C
T2=(X1(:,4))';
; c8 y, X. }4 |% ZX3=X2(:,1:6);# l3 m$ Y" ?* b9 A
T3=T2(:,1:6);
4 L2 d* @8 R5 G/ k8 Z8 knet=newff(minmax(X3),[7 1],{'logsig','tansig'},'trainlm');9 o, ?4 Y7 A& v1 o
net.trainParam.epochs=1000;
' f/ ~, t. w5 _9 R' g+ E& p6 @net.trainParam.goal=0.000001;3 ]! {3 c8 Y0 k8 h5 W |$ H
LP.lr=0.1;
C. j3 ]& m3 D3 Anet=train(net,X3,T3);8 C& Y& }$ `3 e+ r3 _7 b
y1=sim(net,X3);
+ {" [8 D. e7 A
& ]3 W# F4 s. X5 g4 M6 d2 e2 H以下是检验:
; m s! t9 I" c& h: j
5 @+ v) \$ S. L$ q) ~0 Sx_test=X2(:,7);
( _2 {: H0 H5 \; Yy_test=sim(net,x_test)
* ~6 r) h9 U* R' Z- _1 u# `+ iY_test=y_test*(max(x(1, )-min(x(1, ))+min(x(1, )(反归一化)
0 w+ K! L/ X3 L. v/ H检验效果还不错; T% @: P6 j8 O- Z% [
6 g! @" ~8 @; R; _( T
但运用以下命令在预测2010年运输量时,向量yc1为2007、2008、2009年运输量归一化的数值;; C) ]7 g9 P7 Q) i( Z* ?* n
为什么反归一化输出和2009年数值一样,怎么回事呢???初学者盼解答??
/ T( C- w+ k4 E* A1 a. F y1 P
4 _6 b" ^- |+ _1 Q( nyc1=[0.7205 0.7721 1.0000]';' \0 e' L t8 l2 O& t
jg1_test=sim(net,yc1(:,1))" S, X3 j7 j+ ?4 Y/ G/ V$ Y
Jg1_test=jg1_test*(max(x(1, )-min(x(1, ))+min(x(1, )(反归一化)
: G9 J; \ r; i0 N* S. K8 T R |
zan
|