- 在线时间
- 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 编辑 + i* q2 e+ n6 a% v0 S0 s: N
0 S! X! {; I+ Z \+ W6 ?+ t: i- _
关于BP神经网络预测
' `* t+ @/ N, D2 T6 T, R, S. O3 [6 }! S2 Y$ M0 U
注:里面的笑脸代表冒号(不知道会出现笑脸)
4 g1 Y# x0 ~$ g3 i
]. D e- _* B- M) r. x6 Q+ w; p# s$ k" u5 q
关于BP时间序列预测1 S" j) M- e) v/ _5 f
程序中x是已知的运输量数值(2000年-2009年),以4个数据为一组,前三个数据作为输入,第四个作为输出,共分7组,前六组训练,隐层节点为7,最后一组检验:# D5 m& t# y& |3 ~% h- `
x=[13369.2 14873.7 17137.3 17432.5 24193.5 28435.1....
& j5 T0 N1 e' `6 Q& ^8 | 33197.3 38758.6 40576.2 48606.3];
3 E& k0 j2 a5 {. Y( U8 oX(1, =(x(1, -min(x(1, ))/(max(x(1, )-min(x(1, ));(归一化); V9 v# I* E! F! y+ l
for i=1:7
7 ]" n& m! I: J8 A9 g X1(i, =[X(i) X(i+1) X(i+2) X(i+3)];/ D% a. J/ M4 [8 Z1 d" y
end
8 p& q+ W9 B$ D" \) ZX2=(X1(:,1:3))';2 h! V% O% _2 T6 t
T2=(X1(:,4))';
$ H# R. v$ ~; c6 uX3=X2(:,1:6);
' \$ N) o0 u& nT3=T2(:,1:6);
# w/ R- @" \/ c0 V* Tnet=newff(minmax(X3),[7 1],{'logsig','tansig'},'trainlm');( D% {9 Y9 Q% X8 W) O
net.trainParam.epochs=1000;0 v5 `7 v y2 ?! B& d
net.trainParam.goal=0.000001;* B& k3 P6 n. L" K- M; y3 O
LP.lr=0.1;
* ]% C6 \4 U! W# hnet=train(net,X3,T3);! R0 c0 y) c% U/ {6 h# s/ Y6 \/ U
y1=sim(net,X3);% l' d( D2 l4 p. X! ]7 n; ?
* z$ e, S0 ?9 u9 p% o9 b以下是检验: i, z5 L. w$ i
8 `9 e# D' `8 j! Vx_test=X2(:,7);" e5 g+ h/ p, `6 U5 f9 |( m
y_test=sim(net,x_test)' x0 J1 _+ [- |1 e t
Y_test=y_test*(max(x(1, )-min(x(1, ))+min(x(1, )(反归一化)
5 F0 A: Y- Y+ m0 V检验效果还不错;
0 a: r( b- n4 W: F r* G# f, C' k; j. H- l% R
但运用以下命令在预测2010年运输量时,向量yc1为2007、2008、2009年运输量归一化的数值;$ U4 R" J6 y4 d+ ~0 W0 m( b+ F* G+ k
为什么反归一化输出和2009年数值一样,怎么回事呢???初学者盼解答??! G/ j% Z( A1 r8 l; }$ y0 M& T
% E* |) J' Q. T" z; v0 N% `yc1=[0.7205 0.7721 1.0000]';! U2 |5 N+ l& n: K, G/ B0 F
jg1_test=sim(net,yc1(:,1))* a0 }" W5 V* R1 C
Jg1_test=jg1_test*(max(x(1, )-min(x(1, ))+min(x(1, )(反归一化)
8 ?0 s- U. ^ ~1 N* M |
zan
|