本帖最后由 鸿蒙123 于 2017-4-24 12:42 编辑
%% 清空环境变量
clc
%输入输出数据
input=data(:,1:2219);
output=[ones(1000,1);zeros(1000,1);ones(509,1);zeros(638,1)];
%随机提取2000个样本为训练样本,1147个样本为预测样本
input_train=input(1:2000,1:2219);
output_train=output(1:2000,1);
input_test=input(2001:3147,1:2219);
output_test=output(2001:3147,1);
%输入数据归一化
[inputn,inputps]=mapminmax(input_train);
%% BP网络训练
% %初始化网络结构
net=newff(inputn,output_train,10);
net.trainParam.epochs=1000;
net.trainParam.lr=0.1;
net.trainParam.goal=0.0000004;
%% 网络训练
net=train(net,input_train,output_train);
%% BP网络预测
%预测数据归一化
input_test=mapminmax('apply',input_test,inputps);
%网络预测输出
BPoutput=sim(net,input_test);
%% 结果分析