BP神经网络预测问题

jthawking

以下是我写的一段BP神经网络预测的代码,模型建好后,不知怎样预测下一年的T值,哪位高手来指导下...
%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;
; N- u* W% p0 t& _: t3 b- N- B33689 467.57  569.55  2887 7.76  69.8 18961 1015 20249 941.99 2.38;
* U4 R8 {, k0 q& p41020 538.17  616.25  2920 8.23  69.7 19688 1198 21626 1082.6 2.73;
! R6 w7 E& C! l+ a! |- W+ j# [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;
53887 1095.13 969.1   4885 9.65  68.7 36149 2100 25936 2199.5 3.08;
$ y9 J$ l3 w2 W% X59271 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;
70597 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;
2 w$ {9 V7 n. `% ?* q+ t1 i7 D78812 2207.88 927.16  8328 20.52 69.3 71088 2964 30904 3951.64 5.21]';
& y9 ]" E3 Q# l3 y' n, [; _* n8 T%为原始目标数据
! @# G/ h4 m( Q" R4 KT=[3.1,3.6,4.2,4.8,6.7,9.1,13,18.9,29,51.1,78.2,103.5,128];
1 ?! c# @( r- o# N2 M4 D0 D%原始数据归一化
[pn,minp,maxp]=premnmx(P);
[tn,mint,maxt]=premnmx(T);
%设置网络隐单元的神经元数%建立相应的BP网络
net=newff(minmax(pn(:,1:13)),[10,1],{'tansig','purelin'},'trainlm');
%当前输入层权值和阈值
, k5 A; X+ ]* t2 ^% R, O# kinputWeights=net.IW{1,1}
inputbias=net.b{1}
%当前网络层权值和阈值
layerWeights=net.LW{2,1}
layerbias=net.b{2}
net=init(net);
% 网络训练代数
net.trainParam.epochs=10000;
3 e2 b, |# I) q# W% 网络训练目标误差
7 a9 T* Q4 f& W  rnet.trainParam.goal=0.000001;
% 网络训练动量系数
net.trainParam.mc = 0.9;   
5 _4 p  h7 M1 Y1 i. I0 A% 网络学习速率
, g! e* U5 {- Z+ p  r: }net.trainParam.lr=0.05;
8 M2 y0 s' K% Y1 u  O3 o%调用TRAINGDM算法训练BP网络
net=train(net,pn(:,1:13),tn(1:13));
3 z% Q5 J1 h! a) ?/ R- kt_sim=sim(net,pn);
t_sim=postmnmx(t_sim,mint,maxt)
E=t_sim-T
M=sse(E)
N=mse(E)
9 D# k9 r0 G7 f8 H6 wfigure(1)
9 d' s" x% T! {$ x' |6 ?i=1:1:length(t_sim);
j=1:length(T);
3 l: \0 X1 K7 Pplot(i+1995,t_sim,'r',j+1995,T,'bo')
legend('模拟数据','原始数据')
Xlabel('年份');Ylabel('数量')
title('BP神经网络')

杨辰宇

X和Y大写了
& K% ]8 m% V" V: ^

狼之魂汪洋

楼主的问题最后修改完了么？能不能发一份代码学习学习，现在有也遇到同样的问题了，不知道怎么解决了，求帮助。1345134018@qq.com

月关

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

月关

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

月关

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

月关

看不懂。。。。。

狼之魂汪洋

楼主威武。。。。。

左手边623

//////////////////////////////////

melissa513

看不是很懂啊，感觉自己还没有学到这里呢