现代检测技术课程实验编程:最小二乘法应用编程
现代检测技术课程实验编程:最小二乘法应用编程
- 一、最小二乘法编程题目描述
- 二、最小二乘法编程题目要求
- 三、什么是最小二乘法
- 四、最小二乘法编程步骤
- 4.1、 界面的设计
- 4.2、 程序的编写
- 4.2.1、程序在计算按钮如下的回调函数中编写
- 4.2.2、编辑的文本框输入的数据转换成数字类型的数据
- 4.2.3、将Xi、Yi数据存放与数组中
- 4.2.4、计算最小二乘法直线拟合的回归参数a、b
- 4.2.5、计算五条直线的残差平方和
- 4.2.6、判断最佳的最小二乘法直线的拟合
- 4.2.7、数据和图像的显示
- 五、最小二乘法编程总结
一、最小二乘法编程题目描述
最小二乘法编程题目描述如下所示
在对量程为10MPa的压力传感器进行标定时,传感器输出电压值与压力值之间的关系如下表所示,请简述最小二乘法准则的意义,并分析下列电压-压力直线中哪一条最符合最小二乘法准则?(使用计算机辅助进行计算)
(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,和备选直线方程的参数。
所编程序,要能够直接显示哪条直线为最佳直线,不能人为进行判断。
对所编程序的原理和运行结果进行介绍和分析。
三、什么是最小二乘法
最小二乘法定义如下所示
最小二乘法(又称最小平方法)是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合,其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达 。最小二乘法是解决曲线拟合问题最常用的方法。
四、最小二乘法编程步骤
最小二乘法编程
4.1、 界面的设计
使用MATLAB2014b软件进行GUI页面的设计,如下图所示
- 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: 5if minSub > subArrays(i)minSub = subArrays(i);subJudge = i;end
endstrJudge = '最符合最小二乘法准侧的直线是: 第';
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 ony1 = a1 + b1 * Xi;
axes(handles.axes1);
plot(Xi, y1, 'm');
hold ony2 = a2 + b2 * Xi;
axes(handles.axes1);
plot(Xi, y2, 'r');
hold ony3 = a3 + b3 * Xi;
axes(handles.axes1);
plot(Xi, y3, 'y');
hold ony4 = a4 + b4 * Xi;
axes(handles.axes1);
plot(Xi, y4, 'k');
hold ony5 = a5 + b5 * Xi;
axes(handles.axes1);
plot(Xi, y5, 'g');
hold on
三、 程序的运行结果
- 程序可以自动进行最佳拟合直线的判断。从运行的结果可以确定是第五条直线是最佳最小二乘法直线的拟合。
输入用最小二乘法计算得到的直线拟合回归参数,运行得到如下结果
五、最小二乘法编程总结
- 用MATLAB所编写的GUI页面程序实现了计算最小二乘法直线拟合的回归参数的计算。
- 自动判断哪条直线最符合最小二乘法直线拟合的准侧。
- 最小二乘法直线拟合的图像显示。
现代检测技术课程实验编程:最小二乘法应用编程相关推荐
- 现代检测技术课程实验:电桥测量位移
现代检测技术课程:电桥测量位移 一.内容:电桥测量位移 二.电路连接 三.电路仿真 四.x与Uo的关系式 一.内容:电桥测量位移 电桥测量位移描述如下 使用2个线性电位器(总电阻值为R0=5000Ω, ...
- Web前端开发技术课程实验报告实验3:Vue路由实验
实验代码:实验3第2题实验参考.rar-互联网文档类资源-CSDN下载 Web前端开发技术课程实验报告 实验3:Vue路由实验 姓名:_ __ _ ___ ___ 班级:_ _ _ ___ _ _ ...
- python扫雷游戏实验分析_高级编程技术课程实验报告-扫雷游戏
班级 通 信 一 班 姓名 陈远春 学号 201100120084 课程实验内容 扫雷游戏,是 Windows 操作系统自带的一款经典游戏.其规则简单,上手 容易, 不论男女老少皆可娱乐. 扫雷的目的 ...
- 【接口技术】实验五:可编程并行接口8255
文章目录 一. 实验目的 二. 实验内容 实验5-1 8255并行输入/输出实验 实验5-2 并行接口键盘反转扫描实验 一. 实验目的 了解8255芯片结构及编程方法. 了解8255输入/输出实验方法 ...
- 计算机仿真实验结果的评价与分析,基于问卷调查法的计算机仿真技术课程综合性实验效果分析评价...
摘要:综合性实验对于提高学生综合素质,培养学生实践能力与创新精神具有极其重要的作用,为了科学评价计算机仿真技术课程综合性实验效果,文章基于问卷调查法,阐述如何设计情况调查表,如何采用多种评价方法对学生 ...
- 软件测试文档结课论文,现代检测技术结课论文.doc
结 课 论 文 题 目 名 称 动态分析仪 课 程 名 称 现代检测技术 专 业 学 号 学 生 姓 名 年 月 日 一. 研究内容 (1):设计一个典型系统的动态响应的过度过程的分析仪. (2):输 ...
- 杨月江计算机导论答案,计算机导论课程实验教学大纲-计算机学院-华北科技学院.doc...
计算机导论课程实验教学大纲-计算机学院-华北科技学院.doc 华 北 科 技 学 院 2010级 实验教学大纲 网络工程系 目 录 TOC \o "1-1" \h \z \u HY ...
- 【技术】DTEmpower核心功能技术揭秘(2) - AIOD智能异常点检测技术
一.概述 在上一篇<DTEmpower核心功能技术揭秘(1)--HierarchicalStratify分层分类技术>中,我们提到工业设计方法中往往大量采用了数据驱动的方式构建代理模型.但 ...
- 计算机微机原理及接口技术实训室,《微机原理与接口技术》课程实验报告.doc...
<微机原理与接口技术>课程实验报告 实验一基本操作 一.实验目的 ? 1.熟悉Netlab MCS-51单片机实验系统的设置和使用方法 2.掌握基本MCS-51软件的编写.修改和编译的方法 ...
最新文章
- (C++)1015 德才论
- 为什么现在的人越来越不幸福
- 【计算机网络】手动配置hosts文件解决使用GitHub和Coursera网站加载慢/卡的问题
- QT5日志功能(qDebug、qWarnng、qCritical、qFatal)
- 小程序picker下拉框选择时间,只需要年和月
- 解决Zabbix网页端Get value error: cannot connect to [[192.168.238.139]:10050]: [113] No route to host问题...
- Java为何能一次编写,到处运行
- android 替代map,Android为什么推荐使用SparseArray来替代HashMap?
- 组装台式计算机配置清单,diy之家 - 2017组装电脑配置清单_电脑diy主机配置推荐...
- 8086CPU工作原理
- linux系统查看ssh是否开启,linux查看与开启ssh
- 时空序列预测之PredRNN++(Casual LSTM和GHU解决时空预测学习中的深度困境)
- 世界上最著名也最危险的APT恶意软件清单
- K3 CLOUD计划管理之计划方案(MPS/MRP)
- word如何设置每一章节的页眉都不同
- buck电路 dac stm32_嵌入式设计-基于STM32F334的BUCK同步降压数字电源设计-电路城论坛 - 电子工程师学习交流园地...
- 使用JDBC连接mysql遇到的若干问题及解决办法
- 2021年六级翻译:海南岛
- m.soudashi.cn 地图_搜索引擎网站推广优化有什么技巧?
- setlocale(LC_ALL,) 中的 LC_* 是什么意思,以及 LANGUAGE/LANG 都是些啥?