建立多个(.m)文件,有一个主函数,调用其他函数即可。如果想直接使用的话,要注意代码里面的路径,每个人运行的路径不一样。数据我就不上传了,大家想用,用自己的数据跑一下吧,注意矩阵的行和列,很多新手跑的时候认为跑不出来以为代码有问题,其实往往不是的。有什么好的想法大家可以一起研究一下。

clc
clear all;
filename = 'E:\\研一\\智能传感器\\多元回归\\多元回归\\数据.xlsx';
[data,txt,raw] = xlsread(filename);%读取全部内容
filename = 'E:\\研一\\智能传感器\\多元回归\\多元回归\\融合数据处理.xlsx';
[data_chuli,txt,raw] = xlsread(filename);%读取全部内容
[A,H,P] = get_A(data);%计算出系数矩阵的6个数据
[P_ronghe,P,Wu_P] = get_P(A,H,P);
data_ronghe = get_data_ronghe(P_ronghe,P,Wu_P);
[B_K,data_fit] = get_B_K(data);
[ronghe_shuju_xainxing] = get_shujuronghe(data_chuli);

以下是实现某一个功能的函数,可以放入一个(.m)文件也可以放入单独的文件

function [ronghe_shuju_xainxing] = get_shujuronghe(data_chuli)
ronghe_shuju_xainxing = zeros(3,6,'double');%表5-4for i = 1:1:6      ronghe_shuju_xainxing(1,i) = data_chuli((i*3)-1,9);ronghe_shuju_xainxing(2,i) = data_chuli(i*3,9);ronghe_shuju_xainxing(3,i) = data_chuli(i*3,9) - data_chuli((i*3)-1,9);endcsvwrite('E:\研一\智能传感器\多元回归\多元回归\ronghe_shuju_xainxing.csv',ronghe_shuju_xainxing);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [P_ronghe,P,Wu_P] = get_P(A,H,P)%求P_ronghe,P,Wu_P矩阵
P_ronghe = zeros(1,36);%定义空间
Wu_P = zeros(1,36);%定义空间for j = 1:36for i = 1:6 P_ronghe(1,j) = P_ronghe(1,j) + A(1,i)*H(i,j);endWu_P(1,j) = P_ronghe(1,j) - P(1,j);end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function  data_ronghe = get_data_ronghe(P_ronghe,P,Wu_P)
data_ronghe = zeros(18,6);%与data数据进行对比
P = reshape(P,6,6);               P = P';%P为6*6;
P_ronghe = reshape(P_ronghe,6,6); P_ronghe = P_ronghe';%P_ronghe为6*6
Wu_P = reshape(Wu_P,6,6);         Wu_P = Wu_P';%Wu_P为6*6for i = 1:6for j = 1:6data_ronghe(1:3:16,j) = P(1:1:6,j); data_ronghe(2:3:17,j) = P_ronghe(1:1:6,j);data_ronghe(3:3:18,j) = Wu_P(1:1:6,j);endend
csvwrite('E:\研一\智能传感器\多元回归\多元回归\data_ronghe.csv',data_ronghe);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%最小二乘法,以21.5度为静态,建立最小二乘法模型
function  [B_K,data_fit] = get_B_K(data)%求P矩阵
B_K = zeros(1,2,'double');
data_fit = zeros(3,6);
sun_x = 0;  sun_y = 0;  sun_xy = 0; sun_x_squ = 0;
Num = 6;for i = 1:1:Numsun_x = sun_x + data(2,i+1);sun_y = sun_y + data(1,i+1);sun_xy = sun_xy + data(2,i+1)*data(1,i+1);sun_x_squ = sun_x_squ  + data(2,i+1)^2;endB_K(1,1) = (sun_x_squ*sun_y-sun_x*sun_xy)/(Num*sun_x_squ-sun_x*sun_x);B_K(1,2) = (Num*sun_xy-sun_x*sun_y)/(Num*sun_x_squ-sun_x^2);for i = 1:1:Numdata_fit(1,i) = data(1,i+1);data_fit(2,i) = (B_K(1,1)+B_K(1,2)*data(2,i+1));data_fit(3,i) = (B_K(1,1)+B_K(1,2)*data(2,i+1))-data(1,i+1);end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [A,H,P]=get_A(data)%求A,H,P的系数矩阵
Up = data(2:2:12,2:7);Up = Up(:)';%Up为1行36列的矩阵
Ut = data(3:2:13,2:7);Ut = Ut(:)';%Ut为1行36列的矩阵
%P为1*36的矩阵
P = [zeros(1,6,'double'),ones(1,6,'double'),2*ones(1,6,'double'),3*ones(1,6,'double'),4*ones(1,6,'double'),5*ones(1,6,'double')];
H = zeros(6,36);%H为6*36的矩阵
for i=1:1:36H(1,i) = 1;H(2,i) = Up(1,i);H(3,i) = Ut(1,i);H(4,i) = Up(1,i)^2;H(5,i) = Up(1,i)*Ut(1,i);H(6,i) = Ut(1,i)^2;
end
A = zeros(1,6,'double');% A矩阵为1*6
A = P*H'*inv(H*H');%求A的系数矩阵
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

matlab多元回归相关推荐

  1. matlab多元回归程序,多元回归程序MATLAB程序

    <多元回归程序MATLAB程序>由会员分享,可在线阅读,更多相关<多元回归程序MATLAB程序(45页珍藏版)>请在人人文库网上搜索. 1.程序MATLAB多元回归程序matl ...

  2. matlab多元回归模型分析,matlab多元回归工具箱 Excel数据分析工具进行多元回归分析.doc...

    matlab多元回归工具箱 Excel数据分析工具进行多元回归分析.doc matlab多元回归工具箱 Excel数据分析工具进行多元回归分析 导读:就爱阅读网友为您分享以下"Excel数据 ...

  3. matlab多元回归程序,多元回归程序MATLAB程序.doc

    标准实用 文案大全 matlab回归(拟合)总结 前言 1.学三条命令 polyfit(x,y,n)---拟合成一元幂函数(一元多次) regress(y,x)----可以多元, nlinfit(x, ...

  4. 多元线性回归matlab代码例题_matlab多元线性回归

    1.matlab多元回归示例如下: 解决问题:油价预测 方法:多元线性回归 实现:matlab regress()函数 技巧:通过增加X1^2,X2^2,或者X1*X2等构造的特征项,可以提高回归模型 ...

  5. 【预测模型】基于灰度模型预测房价matlab源码

    1 简介 自2006年开始,房价保持较长时间持续攀升,房价过高已成为一个全社会关注的民生问题.从市场价格内在发展趋势和外部因素影响两方面建立房价构成模型,将居民消费物价指数,房地产开发综合景气指数,土 ...

  6. 毕设论文数据分析记录-part2:相关性分析

    20220225-相关性分析 polyfit线性拟合后进行t检验判断线性关系显著水平-一阶差分去趋势 a=xlsread('gai-maizey.xls','gai-maizey'); a1=a(:, ...

  7. matlab某种水泥在凝固,Matlab实现多元回归实例

    Matlab 实现多元回归实例 假设已有数据X 和Y ,在Matlab 软件包中,使用stepwise 命令进行逐步回归,得到回归方程n n Y a X a X a X ε=++???++1122,其 ...

  8. 多元回归——matlab,R与spss

    多元回归(Multiple Regression) 具体是指在相关变量中指定其中一个为因变量,其余一个或多个变量视为自变量,建立多个变量之间的线性或非线性数学模型. 其中,称为回归参数,我们的目的就是 ...

  9. MATLAB实现一元回归和多元回归

    在回归分析中,如果有两个或两个以上的自变量,就称为多元回归.事实上,一种现象常常是与多个因素相联系的,由多个自变量的最优组合共同来预测或估计因变量,比只用一个自变量进行预测或估计更有效,更符合实际. ...

最新文章

  1. PLM(产品生命周期管理)的前世今生
  2. mysql查询_MySQL基础,查询语句详解
  3. RTC-IC-PCF2129
  4. gdb info frame信息解释
  5. 操作系统中的进程与线程和java中的线程
  6. BZOJ.4516.[SDOI2016]生成魔咒(后缀自动机 map)
  7. 关于Python的一些学习笔记(小白式笔记,持续更新)
  8. 用函数模板实现选择排序算法_干货|STL容器和算法
  9. 程序员出差是去干什么_为什么女生不适合做程序员?原来互联网公司女生少,是由于这几点...
  10. 安全验证框架使用笔记001---Shiro简介
  11. 大学计算机四级报名,2016下半年安徽理工大学计算机四级报名
  12. MTK平台修改分区大小之后,通过fastboot工具烧录的说明
  13. pix4d空三导入cc_CC(Context Capture)软件安装及空三过程中的十大常见报错与解决方法...
  14. 扣费克星 1.72 更新说明
  15. fortran 语法学习之二
  16. ​ 众至科技数据防泄露系统,保护企业办公核心数据
  17. 有软件负载均衡,也有硬件负载均衡,选择哪个?
  18. NSGA-II算法阅读笔记
  19. towhee+elasticsearch实现本地以图搜图
  20. 生理期是女孩变美的黄金期

热门文章

  1. win10磁盘分区解决方案
  2. 基于雷达的安防视频监控系统软件设计
  3. PHOTOSHOP应用技巧系列文章之投影(一)
  4. month(字段)、year(字段)
  5. 读书笔记-只有偏执狂才能生存
  6. 遥感识别4——语义分割样本数据集制作
  7. 大智慧新一代文件下载
  8. 如何开机进入grub界面_开机显示grub怎么办,开机进入grub界面
  9. BENNING46393电源维修SMPS-MC2贝朗北宁电源维修
  10. 使用JavaScript实现GPA计算器(学科实践任务 一)