- 在线时间
- 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 编辑 ( e/ Z( ~% f% A" M5 a: ]) p
: t; B& \$ ~6 |) x5 n. k
关于BP神经网络预测5 L+ V% {% y; R% r
( s3 R8 g8 Z0 m( e7 w1 z& l注:里面的笑脸代表冒号(不知道会出现笑脸)! [/ X3 b% l" R1 n
7 O6 J+ C6 R" O- J% J: S9 `* b/ E1 x. n& h5 F- @. U
关于BP时间序列预测
& ^5 A" ^! K9 _% z# j5 h程序中x是已知的运输量数值(2000年-2009年),以4个数据为一组,前三个数据作为输入,第四个作为输出,共分7组,前六组训练,隐层节点为7,最后一组检验:1 J7 S8 ]% c' r# f1 U1 P
x=[13369.2 14873.7 17137.3 17432.5 24193.5 28435.1....
. D( q( @( r. m# p* P 33197.3 38758.6 40576.2 48606.3];
; T- p8 Z) a# e- K7 @ gX(1, =(x(1, -min(x(1, ))/(max(x(1, )-min(x(1, ));(归一化)$ j9 R+ Y3 P, F6 ]2 [ E
for i=1:7% S& l5 {1 ?1 g" M8 h3 Z, K; a- \' E
X1(i, =[X(i) X(i+1) X(i+2) X(i+3)];& n& T* ?+ [& l
end! u7 D8 t# A2 T7 l# ]" s
X2=(X1(:,1:3))';
% D( O# [$ W( L' Z+ ^ k( u3 ZT2=(X1(:,4))';
7 M$ F2 {2 ~1 GX3=X2(:,1:6);
) {1 C- D, ?0 z5 b2 |0 a& W: F2 TT3=T2(:,1:6);
) i! I ]1 a' W' r& Onet=newff(minmax(X3),[7 1],{'logsig','tansig'},'trainlm');- n5 p5 ^4 e9 [" Y9 Z
net.trainParam.epochs=1000;
; L$ j' M2 |: wnet.trainParam.goal=0.000001;/ c& Q0 ^" Z' Q3 v0 h% O
LP.lr=0.1;
0 O( j% O# a/ |; Y) C' Znet=train(net,X3,T3);% Y. h* e ?7 @7 Y' q. @. u3 _
y1=sim(net,X3);% s. F/ F% U1 R( b( q. R) `. E
2 M7 `6 N9 l" N7 K3 E; `8 s, z
以下是检验:& U3 d9 X& k" ]( @
( ~' }: L6 U( O) p# E1 B0 @x_test=X2(:,7); E- _2 _' p' s& [$ c& E
y_test=sim(net,x_test)
) B) ?" ] C. w8 j+ }8 q4 m9 QY_test=y_test*(max(x(1, )-min(x(1, ))+min(x(1, )(反归一化)
/ s& r# }7 Y5 M; h* a, c检验效果还不错;6 T6 _( q$ e. h' P$ F. Y: O5 T1 x! Q
, }' l% M! L1 N
但运用以下命令在预测2010年运输量时,向量yc1为2007、2008、2009年运输量归一化的数值;: H+ t7 ?/ N" U- _, `. O" _/ m
为什么反归一化输出和2009年数值一样,怎么回事呢???初学者盼解答??
9 T: C0 m9 d1 i l1 R0 s
6 w+ s; q& i! n" m" Zyc1=[0.7205 0.7721 1.0000]';
8 a8 q% \6 {# u( C- V* O% Wjg1_test=sim(net,yc1(:,1))! {' `: b+ t5 z/ G9 r
Jg1_test=jg1_test*(max(x(1, )-min(x(1, ))+min(x(1, )(反归一化)
$ V e. {+ u1 v6 j |
zan
|