• 回归预测 | MATLAB实现PCA-BP主成分降维结合BP神经网络多输入单输出回归预测


    回归预测 | MATLAB实现PCA-BP主成分降维结合BP神经网络多输入单输出回归预测

    效果一览

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

    3

    基本介绍

    MATLAB实现PCA-BP主成分降维算法结合BP神经网络多输入单输出回归预测(完整源码和数据)
    1.输入多个特征,输出单个变量;
    2.多输入单输出回归预测;
    3.多指标评价,代码质量极高;excel数据,方便替换,运行环境2018及以上;
    4.评价指标MAE、MBE、R2,代码质量极高,方便学习和替换数据;
    5.要求2018版本及以上。


    程序设计

    • 完整源码和数据获取方式:私信回复MATLAB实现PCA-BP主成分降维结合BP神经网络多输入单输出回归预测
    %%  清空环境变量
    warning off             % 关闭报警信息
    close all               % 关闭开启的图窗
    clear                   % 清空变量
    clc                     % 清空命令行
    
    %%  导入数据
    res = xlsread('data.xlsx');
    
    %%  划分训练集和测试集
    temp = randperm(103);
    
    P_train = res(temp(1: 80), 1: 7)';
    T_train = res(temp(1: 80), 8)';
    M = size(P_train, 2);
    
    P_test = res(temp(81: end), 1: 7)';
    T_test = res(temp(81: end), 8)';
    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);
    
    
    
    %%  仿真测试
    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);
    
    
    
    %%  相关指标计算
    % 决定系数 R2
    R1 = 1 - norm(T_train - T_sim1)^2 / norm(T_train - mean(T_train))^2;
    R2 = 1 - norm(T_test -  T_sim2)^2 / norm(T_test -  mean(T_test ))^2;
    
    disp(['训练集数据的R2为:', num2str(R1)])
    disp(['测试集数据的R2为:', num2str(R2)])
    
    % 平均绝对误差 MAE
    mae1 = sum(abs(T_sim1 - T_train)) ./ M ;
    mae2 = sum(abs(T_sim2 - T_test )) ./ N ;
    
    disp(['训练集数据的MAE为:', num2str(mae1)])
    disp(['测试集数据的MAE为:', num2str(mae2)])
    
    % 平均相对误差 MBE
    mbe1 = sum(T_sim1 - T_train) ./ M ;
    mbe2 = sum(T_sim2 - T_test ) ./ N ;
    
    disp(['训练集数据的MBE为:', num2str(mbe1)])
    disp(['测试集数据的MBE为:', num2str(mbe2)])
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
    • 10
    • 11
    • 12
    • 13
    • 14
    • 15
    • 16
    • 17
    • 18
    • 19
    • 20
    • 21
    • 22
    • 23
    • 24
    • 25
    • 26
    • 27
    • 28
    • 29
    • 30
    • 31
    • 32
    • 33
    • 34
    • 35
    • 36
    • 37
    • 38
    • 39
    • 40
    • 41
    • 42
    • 43
    • 44
    • 45
    • 46
    • 47
    • 48
    • 49
    • 50
    • 51
    • 52
    • 53
    • 54
    • 55
    • 56
    • 57
    • 58
    • 59
    • 60
    • 61
    • 62
    • 63
    • 64

    参考资料

    [1] https://blog.csdn.net/kjm13182345320/article/details/129215161
    [2] https://blog.csdn.net/kjm13182345320/article/details/128105718

  • 相关阅读:
    个人网站制作 Part 25 添加实时聊天功能 | Web开发项目添加页面缓存
    园子的现代化建设-复活:沉睡2年多的新闻评论功能重新开放
    实战:实现一个LRU
    字符与数字的相互转换
    计算机网络基础(二):物理层、数据链路层及网络层
    树莓派安装使用全记录
    Java剑指Offer青蛙跳台阶问题
    R实现地图相关图形绘制
    【动态规划】区间动态规划
    Worthington真菌分泌酶研究 丨半乳糖氧化酶方案
  • 原文地址:https://blog.csdn.net/kjm13182345320/article/details/132744276