function C_net_BP01

u_hebing=textread('u_hebing.txt');

f_hebing=textread('f_hebing.txt');

%数据验证===============================================================

data16test=textread('val_in_03_all_fifth.txt','%s');

data16to2test=hex2dec(data16test)';

Changdu16to2test=length(data16to2test);

lietest=Changdu16to2test/6;

grouptest=zeros(lietest,6);

for itest=1:lietest

grouptest(itest,:)=data16to2test(6*itest-5:6*itest);

end

data_read0test=textread('standard_reading_03_all_fifth.txt');

jiequ_n=96;jiequ_m=121;%截取行数96-121比较匹配

u_jiequ=grouptest(jiequ_n:jiequ_m,:);

f_jiequ=data_read0test(jiequ_n:jiequ_m,1:6);

for iitest=1:(jiequ_m-jiequ_n+1)/2

u_chazhitest(iitest,:)=u_jiequ(2*iitest,:)-u_jiequ(2*iitest-1,:);

f_chazhitest(iitest,:)=f_jiequ(2*iitest,:)-f_jiequ(2*iitest-1,:);

end

% 归一化处理==============================================

[u_guiyi,u_guiyi_ps]=mapminmax(u_hebing.');

[f_guiyi,f_guiyi_ps]=mapminmax(f_hebing.');

u_guiyi_test=mapminmax('apply',u_chazhitest.',u_guiyi_ps);

f_guiyi_test=mapminmax('apply',f_chazhitest.',f_guiyi_ps);

%训练BP神经网络=========================================================

net=feedforwardnet(3,'trainlm');

view(net)

net=train(net,u_hebing.',f_hebing.');

view(net)

f_wucha=net(u_hebing.');

perf=perform(net,f_wucha,f_hebing.')

save 'C_net_BP01' net

%带入BP神经网络

f_out=sim(net,u_chazhitest.').';

%网络权值阈值矩阵=======================================================

b1=net.b{1};

iw=net.iw{1,1};

b2=net.b{2};

lw=net.lw{2,1};

%离线神经网络表达式

f_offline_guiyi=(purelin(lw*tansig(iw*u_guiyi_test+b1)+b2));

f_offline=mapminmax('reverse',f_offline_guiyi,f_guiyi_ps).';

end