041722sjwbp.m

matlab实现BP神经网络

S1=xlsread('e:data/Kyangben.xls');%困境样本
S2=xlsread('e:data/KTest.xls');%困境测试
S3=xlsread('e:data/Zyangben.xls');%正常样本
S4=xlsread('e:data/ZTest.xls');%正常测试
%输入
P1=S1(1:429,1:27);%困境
P3=S3(1:4195,1:27);%正常
P=[P1',P3'];
%目标输出
A=2*ones(429,429);%困境
A=A(1,:);
B=ones(4195,4195);%正常
B=B(1,:);
T=[A,B];
%设计模型
net=newff(minmax(P),[27 1],{'tansig','purelin'},'trainrp');
%训练
net.trainParam.show=50;
net.trainParam.lr=0.01;
net.trainParam.epochs=600;
net.trainParam.goal=1e-10;
[net,tr]=train(net,P,T);
%仿真
P=S2(1:100,1:27);
Y=sim(net,P')