现代检测技术课程实验编程:最小二乘法应用编程

简介: 现代检测技术课程实验编程:最小二乘法应用编程

e2bd1989831d48ff8cc0d78728be70d1.jpg

一、最小二乘法编程题目描述


ea32577bc3da49069a46609727645aa7.jpg



最小二乘法编程题目描述如下所示

在对量程为10MPa的压力传感器进行标定时,传感器输出电压值与压力值之间的关系如下表所示,请简述最小二乘法准则的意义,并分析下列电压-压力直线中哪一条最符合最小二乘法准则?(使用计算机辅助进行计算)




bc66bb826d694026a70fd17fd7c57395.jpg

(1) y=5.02x+0.08;(2)y=5.05x+0.07;

(3) y=4.95x+0.12;(4)y=4.95x+0.15;

(5) y=5.00x+0.07


二、最小二乘法编程题目要求


最小二乘法编程题目要求如下所示


使用计算机软件(VB、VC、JAVA、LabVIEW、Matlab、Python均可)编程完成本次编程题目;


所编程序要有较为美观的GUI界面,可以通过人机界面输入校准数据xi/yi,和备选直线方程的参数。


所编程序,要能够直接显示哪条直线为最佳直线,不能人为进行判断。


对所编程序的原理和运行结果进行介绍和分析。


三、什么是最小二乘法


最小二乘法定义如下所示


最小二乘法(又称最小平方法)是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合,其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达 。最小二乘法是解决曲线拟合问题最常用的方法。


四、最小二乘法编程步骤


9991ecb3904a46db95b0058f1616b140.jpg


最小二乘法编程


4.1、 界面的设计


使用MATLAB2014b软件进行GUI页面的设计,如下图所示


4e9c0940de7c44ce9383a2e002f971ea.png



  • 1 界面可以进行测量数据的压力和电压的输入、五条拟合直线的回归参数a、b。
  • 2 点击计算按钮,可以完成最小二乘法直线拟合的回归参数a、b。五条直线的残差平方和。自动判断哪条直线最符合最小二乘法准则。
  • 3 页面最右边的轴可以显示最小二乘法直线拟合的图像。


4.2、 程序的编写


4.2.1、程序在计算按钮如下的回调函数中编写


1 程序在计算按钮如下的回调函数中编写,代码如下


function pushbutton1_Callback(hObject, eventdata, handles)


4.2.2、编辑的文本框输入的数据转换成数字类型的数据


2 将可以编辑的文本框输入的数据转换成数字类型的数据,代码如下

Xi1 =  get(handles.edit_Xi1, 'String'); Yi1 =  get(handles.edit_Yi1, 'String');
Xi2 =  get(handles.edit_Xi2, 'String'); Yi2 =  get(handles.edit_Yi2, 'String');
Xi3 =  get(handles.edit_Xi3, 'String'); Yi3 =  get(handles.edit_Yi3, 'String');
Xi4 =  get(handles.edit_Xi4, 'String'); Yi4 =  get(handles.edit_Yi4, 'String');
Xi5 =  get(handles.edit_Xi5, 'String'); Yi5 =  get(handles.edit_Yi5, 'String');
Xi1 = str2num(Xi1); Yi1 = str2num(Yi1);
Xi2 = str2num(Xi2); Yi2 = str2num(Yi2);
Xi3 = str2num(Xi3); Yi3 = str2num(Yi3);
Xi4 = str2num(Xi4); Yi4 = str2num(Yi4);
Xi5 = str2num(Xi5); Yi5 = str2num(Yi5);
a1 =  get(handles.edit_a1, 'String'); b1 =  get(handles.edit_b1, 'String');
a2 =  get(handles.edit_a2, 'String'); b2 =  get(handles.edit_b2, 'String');
a3 =  get(handles.edit_a3, 'String'); b3 =  get(handles.edit_b3, 'String');
a4 =  get(handles.edit_a4, 'String'); b4 =  get(handles.edit_b4, 'String');
a5 =  get(handles.edit_a5, 'String'); b5 =  get(handles.edit_b5, 'String');
a1 = str2num(a1); b1 = str2num(b1);
a2 = str2num(a2); b2 = str2num(b2);
a3 = str2num(a3); b3 = str2num(b3);
a4 = str2num(a4); b4 = str2num(b4);
a5 = str2num(a5); b5 = str2num(b5);


4.2.3、将Xi、Yi数据存放与数组中

3 将Xi、Yi数据存放与数组中,代码如下

Xi = [Xi1 Xi2 Xi3 Xi4 Xi5];
Yi = [Yi1 Yi2 Yi3 Yi4 Yi5];


4.2.4、计算最小二乘法直线拟合的回归参数a、b

4 计算最小二乘法直线拟合的回归参数a、b,代码如下

squareXi = Xi .* Xi;
squareYi = Yi .* Yi;
mulXiYi = Xi .* Yi;
sumXi = sum(Xi);
sumYi = sum(Yi);
sumSquareXi = sum(squareXi);
sumSquareYi = sum(squareYi);
sumMulXiYi = sum(mulXiYi);
Lxx = sumSquareXi - sumXi * sumXi / 5;
Lxy = sumMulXiYi - sumXi * sumYi / 5;
b = Lxy / Lxx;
a = sumYi / 5 - b * sumXi / 5;

4.2.5、计算五条直线的残差平方和

5 计算五条直线的残差平方和,代码如下


SubYiXi1 = Yi - (a1 + b1 * Xi);
squareSubYiXi1 = SubYiXi1 .* SubYiXi1;
sumSub1 = sum(squareSubYiXi1);
SubYiXi2 = Yi - (a2 + b2 * Xi);
squareSubYiXi2 = SubYiXi2 .* SubYiXi2;
sumSub2 = sum(squareSubYiXi2);
SubYiXi3 = Yi - (a3 + b3 * Xi);
squareSubYiXi3 = SubYiXi3 .* SubYiXi3;
sumSub3 = sum(squareSubYiXi3);
SubYiXi4 = Yi - (a4 + b4 * Xi);
squareSubYiXi4 = SubYiXi4 .* SubYiXi4;
sumSub4 = sum(squareSubYiXi4);
SubYiXi5 = Yi - (a5 + b5 * Xi);
squareSubYiXi5 = SubYiXi5 .* SubYiXi5;
sumSub5 = sum(squareSubYiXi5);


4.2.6、判断最佳的最小二乘法直线的拟合

6 判断最佳的最小二乘法直线的拟合,代码如下

subArrays = [sumSub1 sumSub2 sumSub3 sumSub4 sumSub5];
minSub = subArrays(1);
subJudge = 1;
for i = 2: 5
    if minSub > subArrays(i)
        minSub = subArrays(i);
        subJudge = i;
    end
end
strJudge = '最符合最小二乘法准侧的直线是: 第';
subJudge = num2str(subJudge);
strJudge1 = '条直线';
strJudge = strcat(strJudge, subJudge, strJudge1);


4.2.7、数据和图像的显示

7 数据和图像的显示,代码如下


set(handles.text_judge, 'String', num2str(strJudge));
set(handles.edit1, 'String', num2str(sumSub1));
set(handles.edit2, 'String', num2str(sumSub2));
set(handles.edit3, 'String', num2str(sumSub3));
set(handles.edit4, 'String', num2str(sumSub4));
set(handles.edit5, 'String', num2str(sumSub5));
set(handles.edit_result_a, 'String', num2str(a));
set(handles.edit_result_b, 'String', num2str(b));
plot(Xi, Yi, '*');
hold on
y1 = a1 + b1 * Xi;
axes(handles.axes1);
plot(Xi, y1, 'm');
hold on
y2 = a2 + b2 * Xi;
axes(handles.axes1);
plot(Xi, y2, 'r');
hold on
y3 = a3 + b3 * Xi;
axes(handles.axes1);
plot(Xi, y3, 'y');
hold on
y4 = a4 + b4 * Xi;
axes(handles.axes1);
plot(Xi, y4, 'k');
hold on
y5 = a5 + b5 * Xi;
axes(handles.axes1);
plot(Xi, y5, 'g');
hold on

三、 程序的运行结果


5bfa975927fb4657b118f601ac4dd058.png


  • 程序可以自动进行最佳拟合直线的判断。从运行的结果可以确定是第五条直线是最佳最小二乘法直线的拟合。

输入用最小二乘法计算得到的直线拟合回归参数,运行得到如下结果


30144cb731f74be6b0d6e127ac62d545.png


五、最小二乘法编程总结


909a93cc45e846509825ae2831a38087.jpg



  • 用MATLAB所编写的GUI页面程序实现了计算最小二乘法直线拟合的回归参数的计算。
  • 自动判断哪条直线最符合最小二乘法直线拟合的准侧。
  • 最小二乘法直线拟合的图像显示。

4cd359e17aac4b98a577da5b946712a0.jpg

相关文章
|
8月前
|
算法
计算机算法设计与分析 第1章 算法概述 (笔记)
计算机算法设计与分析 第1章 算法概述 (笔记)
现代检测技术课程实验编程:波特图分析仪原理仿真:一阶检测系统编程仿真
现代检测技术课程实验编程:波特图分析仪原理仿真:一阶检测系统编程仿真
现代检测技术课程实验编程:波特图分析仪原理仿真:一阶检测系统编程仿真
|
机器学习/深度学习 算法 安全
机器学习的统计方法 贝叶斯决策理论入门
机器学习的统计方法 贝叶斯决策理论入门
213 0
机器学习的统计方法 贝叶斯决策理论入门
|
索引
【仿真建模】第二课:AnyLogic入门基础课程 - 行人仿真空间逻辑讲解
每个智能体都有个属性index,从0到n,所以我们创建多层建筑没必要一层一层建立,只需要把一层建筑封装成一个类,然后拖动出来,使用类似for循环的机制,去复制即可。指定初始位置,index*20,代表每层高度为20,通过index自增的索引进行高度自增。这里要注意,一定要将矩形墙的中心对准原点,因为在复制封装好的类时,它会以原点为参考原点。设置墙的颜色(大家自己选一个颜色即可),透明度设置为100(主要是为了能看清建筑内部)在工程面板找到刚新建的层,然后修改dZ为40(高度为40)
597 0
【仿真建模】第二课:AnyLogic入门基础课程 - 行人仿真空间逻辑讲解
|
网络协议 程序员 测试技术
接口编程的经典案例|学习笔记
快速学习接口编程的经典案例
|
机器学习/深度学习 数据采集 存储
不谈高级原理,只用简单的语言来聊聊机器学习
不谈高级原理,只用简单的语言来聊聊机器学习
379 0
不谈高级原理,只用简单的语言来聊聊机器学习
|
存储 机器学习/深度学习 算法
K近邻算法哪家强?KDTree、Annoy、HNSW原理和使用方法介绍
K近邻算法哪家强?KDTree、Annoy、HNSW原理和使用方法介绍
|
算法 数据可视化
ML之MaL: 流形学习MaL的概念认知、算法分类、案例应用、代码实现之详细攻略
ML之MaL: 流形学习MaL的概念认知、算法分类、案例应用、代码实现之详细攻略
|
算法 机器人 数据挖掘
带你读《数据挖掘导论(原书第2版)》之三:分类:基本概念和技术
本书所涵盖的主题包括:数据预处理、预测建模、关联分析、聚类分析、异常检测和避免错误发现。通过介绍每个主题的基本概念和算法,为读者提供将数据挖掘应用于实际问题所需的必要背景。
带你读《计算思维导论实验 与习题指导》之二:计算基础
本书围绕《计算思维导论》主教材,设计了13个实验,并针对前8章内容设计了习题,包括单选题、多选题、填空题、判断题等。通过实验和习题,能帮助学生:了解计算思维的概念和计算机发展简史;理解进制转换、字符编码和中文编码等相关知识,掌握数制转换的方法和口诀;了解计算机硬件并学会配置与组装计算机,同时能够对简单故障进行判断和排除;掌握上网浏览、查询资料、收发电子邮件等信息时代的必备知识,同时学会局域网的搭建、WWW和FTP服务器的构建;掌握利用Access创建数据库的方法,并能初步设计与管理数据库;掌握命题符号化方法,以及基本的推理理论,并能利用真值表、等值演算等方法进行简单的逻辑推理等能力。

热门文章

最新文章