%% 清空环境变量
warning off %关闭报警信息
close all %关闭开启的图窗
clear %清空变量
clc %清空命令行
%% 导入数据
res = xlsread('C:\新建文件夹\新建 Microsoft Excel 工作表.xlsx');
%% 划分训练集和数据集
temp = randperm(78);%打乱数据
p_train = res(temp(1:70),2:1543)';%行代表一组数据,最后一列是类,其余列是特征
T_train = res(temp(1:70),1)';
M = size(p_train,2);
p_test = res(temp(71:end),2:1543)';
T_test = res(temp(71:end),1)';
N=size(p_test,2);
%% 数据归一化
[p_train,ps_input] = mapminmax(p_train,0,1);
p_test = mapminmax('apply',p_test,ps_input);
[T_train,ps_output] = mapminmax(T_train,0,1);
T_test = mapminmax('apply',T_test,ps_output);
%% 创建网络
net = newff(p_train,T_train,6);
%% 设置训练参数
net.trainParam.epochs=10000;
net.trainParam.goal=1e-6;
net.trainParam.lr=0.01;
%% 训练网络
net = train(net,p_train,T_train);
%% 仿真训练
t_sim1 = sim(net,p_train);
t_sim2 = sim(net,p_test);
%% 数据返归一化
t_sim1=mapminmax('reverse',t_sim1,ps_output);
t_sim2=mapminmax('reverse',t_sim2,ps_output);
%% 均方根误差
error1 = sqrt(sum((t_sim1-T_train).^2)'/M);
error2 = sqrt(sum((t_sim2-T_test).^2)'/N);
%% 绘图
figure
plot(1:M,T_train,'r-*',1:M,t_sim1,'b-o','LineWidth',1)
legend('真实值','预测值')
xlabel('预测样本')
ylabel('预测结果')
%string = {'训练集预测结果对比';['RMSE='num2str(error1)]};
%title(string)
xlim([1,M])
grid
figure
plot(1:N,T_test,'r-*',1:N,t_sim2,'b-o','LineWidth',1)
legend('真实值','预测值')
xlabel('预测样本')
ylabel('预测结果')
%string = {'训练集预测结果对比';['RMSE='num2str(error2)]};
%title(string)
xlim([1,N])
grid
01-21
10-05
“相关推荐”对你有帮助么?
-
非常没帮助
-
没帮助
-
一般
-
有帮助
-
非常有帮助
提交