bp神经网络matlab代码82083

BP神经网络算法的matlab代码，本程序根据训练好的网络文件ANN.mat预测新的数据文件，得到均方误差，并画出预测数据和原数据的对比图。

%%
%对程序的说明%
%做生成数（使用0.1做因子）之后13年的数据按统计方法标准化后输入训练仿真对比结果。这个过程相当于插值计算%
%输出结果请最后五句r1rlplotyy（ktokt）plotyy（kakf）%

%这样做的缺点是不能检验预测效果的好坏还有一点我对网络初始化的过程没有把握%

clear
% 数据输入
inv=[189.6  244.2  328.5  405.3  487.5  568.6  669.9  801.6  785.9  859.2  979.7  1086.3  1241];
d1=[15.3  16.6  17.2  22.8  37.6  30.8  38.9  50.8  58.2  41.2  55.8  71.7  82.3];
d2=[117.7  140.7  177.8  209.8  247.7  290.8  343.2  337.1  308.1  353.6  350.1  367.2  415.7];
d3=[56.6  86.9  133.5  172.6  202.2  247  287.8  413.7  419.6  464.4  573.8  647.4  742.1];
f=[822.3  959.7  1198.4  1607.5  1997.6  2380.9  2683.8  2798.9  2897.4  3253  3561  3882 4433];
k=1991:1:2003;
% 使用13年数据训练得到网络仿真计算相对误差
%对原始数据初始化的过程为t（1）=f（1）;  t（i）=t（i-1）+0.1*f（i）（i=2:13）
% 累加 p1（1）=d1（1）;for i=2:13 p1（i）=0.1*d1（i）+p1（i-1）;end %
% 累减 a（1）=v（1）;for i=2:13 a（i）=（v（i）-v（i-1））*10;end %
% 先赋予0值
p1=zeros（113）;p2=zeros（113）;p3=zeros（113）;t=zeros（113）;
p1（1）=d1（1）;for i=2:13 p1（i）=0.*d1（i）+p1（i-1）;end
p2（1）=d2（1）;for i=2:13 p2（i）=0.1*d2（i）+p2（i-1）;end
p3（1）=d3（1）;for i=2:13 p3（i）=0.1*d3（i）+p3（i-1）;end
t（1）=f（1）;for i=2:13 t（i）=0.1*f（i）+t（i-1）;end
b=[p1;p2;p3;t]‘;
p=[（p1-mean（p1））./std（p1）;（p2-mean（p2））./std（p2）;（p3-mean（p3））./std（p3）];
tt=（t- mean（t））./std（t）;%数据标准化处理标准化为网络输入p期望输出tt

%建立网络并训练
net=newff（minmax（p）[41]{‘tansig‘‘purelin‘}‘trainlm‘）;
net.iw{11}=zeros（size（net.iw{11}））+0.5;
net.lw{21}=zeros（size（net.lw{21}））+0.75;
net.b{11}=zeros（size（net.b{11}））+0.5;
net.b{21}=zeros（size（net.b{21}））;
net.trainParam.epochs=3000;
net.trainParam.goal =0.000005;
net=train（netptt）;
%网络仿真得到网络输出，并计算误差
tt1=sim（netp）;
to=tt1.*std（t）+mean（t）;%利用标准化的逆变换得到t1的近似值to
a=zeros（113）;a（1）=to（1）;
for i=2:13 a（i）=（to（i）-to（i-1））*10; end %累减得到近似的实际产值a
r_net=（tt1-tt）./tt;%网络相对误差
r_add=（to-t）./t;%总产值累加数相对误差
r_real=（a-f）./f;%实际总产值相对误差

%输出结果%
r1=[r_net;r_add;r_real]‘%网络相对误差/累加数相对误差/实际值相对误差
r=r_real./r_add %误差放大倍数
l=[to;t;to-t;a;f;a-f]‘%累加数计算值/累加数/绝对误差//计算的实际值/实际值/绝对误差
plotyy（ktokt）;%总产值的仿真累加值（to）随时间变化曲线与实际累加值（t）随时间变化曲线
%在这里设置断点以查看这两条曲线%
plotyy（kakf）;%总产值的仿真值（a）随时间变化曲线与实际值（f）随时间变化曲线

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2474  2005-05-22 15:06  bp源程序\bp.m

     文件       1988  2005-06-13 21:54  bp源程序\bp2.m

     文件      28792  2004-08-11 20:11  bp源程序\sourse\ANN.mat

     文件       8442  2004-08-12 10:20  bp源程序\sourse\AnnData.txt

     文件       1388  2004-08-12 10:32  bp源程序\sourse\Predict.asv

     文件       1371  2004-08-12 10:27  bp源程序\sourse\Predict.m

     文件        273  2004-08-12 10:34  bp源程序\sourse\说明.txt

    ..AD...         0  2006-03-18 19:39  bp源程序\sourse

     目录          0  2006-03-18 19:38  bp源程序

----------- ---------  ---------- -----  ----

                44946                    10