BP神经网络预测问题

jthawking

以下是我写的一段BP神经网络预测的代码,模型建好后,不知怎样预测下一年的T值,哪位高手来指导下...
: a: B# b2 O! ~# O%X为原始输入数据
P=[727000 297.35  327.53  2658 6.31  68.3 15948 737  16316 583.89 1.96;
30619 325.00  390.51  2763 6.94  69.6 17130 789  18600 707.67 2.28;
33282 423.00  474.63  2801 5.73  69.7 17866 894  19886 861.88 2.32;
33689 467.57  569.55  2887 7.76  69.8 18961 1015 20249 941.99 2.38;
& T, @8 U% z) x0 n% ?( f/ |41020 538.17  616.25  2920 8.23  69.7 19688 1198 21626 1082.6 2.73;
43344 832.04  686.37  3495 8.87  68.3 24465 1361 23544 1373.4 2.89;
46030 941.94  788.15  3495 9.57  68.2 24814 1710 24941 1756.5 2.82;
7 R  r( Y) V+ @) z4 z53887 1095.13 969.1   4885 9.65  68.7 36149 2100 25936 2199.5 3.08;
: b' Y1 ^1 H" K4 a8 t6 d* g4 [; t8 T$ W3 Y59271 1250.64 1092.6  5376 10.11 69.2 43008 2314 26596 2625.4 3.56;
64507 1437.67 1176.1  6091 15.05 69.2 51946 2500 28494 3229.4 3.98;
1 L9 B8 e" }/ I8 t, D; |5 I0 l! R70597 1671.29 1273.7  7305 16.81 69.2 65745 2614 29628 3744.7 4.98;
79221 1905.03 1345    8023 18.6  69.1 70346 2798 30063 3792.6 5.1;
; p, L' P  u6 w4 {, K7 a0 }3 q3 c78812 2207.88 927.16  8328 20.52 69.3 71088 2964 30904 3951.64 5.21]';
%为原始目标数据
T=[3.1,3.6,4.2,4.8,6.7,9.1,13,18.9,29,51.1,78.2,103.5,128];
8 W* _8 H8 _+ @1 k2 x7 X, N& m%原始数据归一化
8 V( E$ Z5 j: d" o1 @9 j' ?[pn,minp,maxp]=premnmx(P);
[tn,mint,maxt]=premnmx(T);
2 w) x/ z' _5 S  J6 u( N% r% j%设置网络隐单元的神经元数%建立相应的BP网络
% a) m, X# x1 F6 s+ D7 A1 y# ?net=newff(minmax(pn(:,1:13)),[10,1],{'tansig','purelin'},'trainlm');
& q% R! a0 Z9 _& P0 z7 U%当前输入层权值和阈值
& o! a, R$ P' P  Y0 l* qinputWeights=net.IW{1,1}
inputbias=net.b{1}
, r# I2 _4 h' P' z: \%当前网络层权值和阈值
  S) g( y  L6 i: S( ?layerWeights=net.LW{2,1}
& b% j$ i) x" i8 p5 O4 Llayerbias=net.b{2}
: u; D3 e. O! P! D1 r; Tnet=init(net);
" L, G; ]5 u! L; f" b  r) H% 网络训练代数
net.trainParam.epochs=10000;
% 网络训练目标误差
net.trainParam.goal=0.000001;
& R" c2 V; J" O- s/ Z# @# S& W% 网络训练动量系数
6 b. R: Q4 `; ^; P3 N' T4 Dnet.trainParam.mc = 0.9;   
( o' a3 z3 g6 s! `9 D8 ?7 |% 网络学习速率
net.trainParam.lr=0.05;
%调用TRAINGDM算法训练BP网络
net=train(net,pn(:,1:13),tn(1:13));
t_sim=sim(net,pn);
% f' F( b/ o2 U7 Y7 {) rt_sim=postmnmx(t_sim,mint,maxt)
0 X8 l5 `( p: \" p0 U# e+ G( J. lE=t_sim-T
' j9 @: j, v- hM=sse(E)
4 m! G: c3 M( ^# A9 v* k# Z' vN=mse(E)
figure(1)
i=1:1:length(t_sim);
  T( m" B! C* \" r6 s3 P7 R/ t/ sj=1:length(T);
' a+ l- L2 w, P# kplot(i+1995,t_sim,'r',j+1995,T,'bo')
8 t9 j- y. ?8 @5 s# G$ }: Rlegend('模拟数据','原始数据')
Xlabel('年份');Ylabel('数量')
1 R% S8 w3 h" j' p7 a3 O2 Dtitle('BP神经网络')

追梦者

看看O(∩_∩)O~

Rosierl

狂谢Ing……

A123J

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!      eeeeeeeee

泽泽

这个问题我也正想问，你解决了没？

zhuangqg6855

嗯，先看看再说！！！！！！！！！！！！！！！！！！！！

zerodingying

感谢楼主，很不错的说。。。。。

123水无痕

我也想知道……求高手指点

0йly_

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

0.9清1.8清2.7清

厉害啊！楼主！！！！！牛！