tbb5.txt

神经网络预测bp算法程序

filename=input('请输入训练BP神经元的数据文件名：','s')
%filename='e:\data\bj420to623.txt';
a=csvread(filename,0,0);
ptemp=a';
%div=[1,0,0,0;0,1/1755,0,0;0,0,1/300,0;0,0,0,1/1450]';
%p=div*ptemp(:,41:61);
p1=ptemp([2,3,4],:);
p=p1([1],:)-p1([2],:)-p1([3],:);
pp=p;
p5=p;
ptag=zeros(14,65);
for  j=7:65
    for i=1:7
    tag(i,j)=p5(i+j-7);
end;
end;
p=tag(:,7:2:24)/2000;
tagfilename=input('请输入目标矢量的数据文件名：','s')
%tagfilename='e:\data\bj420to623add.txt';
b=csvread(tagfilename,0,0);
temp=b'/150;
temp=temp(3,:);
for  j=7:65
    for i=8:14
    tag(i,j)=temp(i+j-14);
end;
end;
t1=temp(8:2:25);
%t is 目标矢量
t=t1/150
NNTWARN OFF
%定义BP网络神经元
s1=5;
%初始BP网络神经元
[w1,b1,w2,b2]=initff(p,s1,'logsig',t,'purelin');
%训练BP网络神经元
df=1000;%step number per one disp 
me=90000;%total step
eg=0.0025;%err
lr=2.0;%study rat
[w1,b1,w2,b2,ep,tr]=trainbpx(w1,b1,'logsig',w2,b2,'purelin',p,t,[df me eg lr]);
%绘制误差曲线
%figure
%ploterr(tr,eg);
%利用训练完的BP网络神经元进行预测
p2=tag(:,7:24)/2000;
c=simuff(p2,w1,b1,'logsig',w2,b2,'purelin');
c=double(int16(c*150+0.5))
outfilename=input('请输入保存结果的文件名:','s');
fid=fopen(outfilename,'w');
fprintf(fid,'%f\n',c);
st=fclose(fid);