- 处理concrete_data
temp=randperm(size(attributes,2));
input_train=attributes(:,temp(1:80));
output_train=strength(:,temp(1:80));
input_test=attributes(:,temp(81:end));
output_test=strength(:,temp(81:end));
N=size(output_test,2)
[fin_train,stru]=mapminmax(input_train,0,1);
fin_test=mapminmax('apply',input_test,stru);
[fou_train,stru1]=mapminmax(output_train,0,1);
net=newff(fin_train,fou_train,9);
net.trainParam.epochs=1000;
net.trainParam.goal=1e-4;
net.trainParam.lr=0.01;
net=train(net,fin_train,fou_train);
test_result=sim(net,fin_test);
output_bp=mapminmax('reverse',test_result,stru1);
err=abs((output_bp-output_test)./output_test)
figure
plot(1:N,output_bp,'r',1:N,output_test,'b')
legend('预测值','实际值')