九层妖塔 起于垒土


Matlab 最小二乘法拟合一阶线性拟合&传感器实验

  • 一、代码
  • 二、数据处理结果
  • 三、Notes

一、代码

%电容传感器位移实验数据 最小二乘法一阶线性拟合
x = [13.080,12.580,12.080,11.580,11.080,10.580,10.080,9.580,9.080,8.580,8.080,7.580,7.080];  %13
y = [-440,-360,-283,-205,-134,-61,0,67,123,184,235,290,337];
p = polyfit(x,y,1);   %最小二乘法一阶线性拟合
a = p(1)     %截距
b = p(2)     %斜率
%计算在 x 中的点处拟合的多项式 p。用这些数据绘制得到的线性回归模型。
f = polyval(p,x);  %拟合直线
plot(x,y,'o',x,f,'-');
grid on;grid minor; %显示网格
legend('实验数据','拟合直线');
str={'-440','-360','-283','-205','-134','-61','0','67','123','184','235','290','337'};
text(x,y+15,str,'HorizontalAlignment','center');
text(9,200,['y= ',num2str(p(1)),'x+',num2str(p(2))],'FontSize',15);
xlabel('X (mm)'); ylabel('V(mV)');
title('电容传感器位移实验数据  最小二乘法一阶线性拟合');m = max(f-y)    %输出值(多次测量时为平均值)与拟合直线的最大偏差
no_line = m / 800     %非线性误差

二、数据处理结果


三、Notes

 ●函数polyfit—多项式曲线拟合
 ●函数polyval—多项式计算
 ●函数text—向数据点添加文本说明


Matlab 隐函数方程求解&最小二乘法拟合一阶线性拟合&二阶拟合&传感器实验

  • 一、代码
  • 二、数据处理结果
  • 三、Notes
%实验三十 Pt100 铂电阻测温特性实验
%实验数据 最小二乘法一阶线性拟合
t = linspace(40,70,7);
Vo = [-112.5,-119.4,-126.2,-133.6,-140.8,-148.2,-155.5];figure(1)  %电压与温度的分析
Vo_p = polyfit(t,Vo,1);   %最小二乘法一阶线性拟合
Vo_a = Vo_p(1)     %斜率
Vo_b = Vo_p(2)     %截距
%计算在 x 中的点处拟合的多项式 p。用这些数据绘制得到的线性回归模型。
Vo_f = polyval(Vo_p,t);  %拟合直线
plot(t,Vo,'o',t,Vo_f,'-');
grid on;grid minor; %显示网格
legend('实验数据','拟合直线');
str={'-112.5','-119.4','-126.2','-133.6','-140.8','-148.2','-155.5'};
text(t,Vo+1,str,'HorizontalAlignment','center');
text(52,-140,['Vo= ',num2str(p(1)),'t+',num2str(p(2))],'FontSize',15);
xlabel('温度t'); ylabel('电压Vo');
title('实验三十 Pt100 铂电阻测温特性实验  最小二乘法一阶线性拟合');
Vo_m = max(Vo_f-Vo)    %输出值(多次测量时为平均值)与拟合直线的最大偏差
Vo_no_line = Vo_m / 50     %非线性误差figure(2)  %阻值与温度的分析
syms Rt K R2 R3 R4 Vc syms_Vo
K=10;
R2=5000;
R3=5000;
R4=100;
Vc=4;
RT = zeros(1,length(Vo));
for n = 1:length(Vo)syms_Vo=Vo(n) * 0.001;eqn =  K*((R4/(R4+R2)) - (Rt/(R3+Rt)))*Vc == syms_Vo ;solx = solve(eqn,Rt);  %求解Rt表达式RT(n) = vpa(solx);     %转化为小数
endRT
p = polyfit(t,RT,1);   %最小二乘法一阶线性拟合
a = p(1)     %截距
b = p(2)     %斜率
%计算在 x 中的点处拟合的多项式 p。用这些数据绘制得到的线性回归模型。
f = polyval(p,t);  %拟合直线
plot(t,RT,'o',t,f,'-');
grid on;grid minor; %显示网格
legend('实验数据','拟合直线');
str={'114.6727'  '115.5754'  '116.4653'  '117.4341'  '118.3770'  '119.3465'  '120.3033'};
text(t,RT+0.5,str,'HorizontalAlignment','center');
text(55,120,['Rt= ',num2str(p(1)),'t+',num2str(p(2))],'FontSize',15);
xlabel('温度t'); ylabel('阻值Rt');
title('实验三十 Pt100 铂电阻测温特性实验  最小二乘法一阶线性拟合');m = max(f-RT)    %输出值(多次测量时为平均值)与拟合直线的最大偏差
no_line = m / 50    %非线性误差、```

Matlab 隐函数方程求解最小二乘法拟合一阶线性拟合二阶拟合传感器实验相关推荐

  1. matlab最小二乘法拟合图旋转,【Matlab】—{最小二乘法拟合一阶线性拟合传感器实验}...

    [Matlab]-{最小二乘法拟合一阶线性拟合传感器实验} [Matlab]-{最小二乘法拟合一阶线性拟合&传感器实验} ???九层妖塔?起于垒土 [Matlab]-{最小二乘法拟合一阶线性拟 ...

  2. matlab将方程的解代入原方程,matlab代入方程求解

    已知F1=0.3320X1+0.2741X2+0.3294X3+0.3309X4+0.0835X5+0.3321X6+0.3284X7+0.2561X8+0.3311X9+0.3042X10+0.32 ...

  3. matlab求方程实根,matlab怎么求方程的根

    MATLAB解方程_IT/计算机_专业资料.一般的代数方程函数solve用于求解一般代数方程的根,假定S为符 号表达式,命令solve (S)求解表达式等于0的根,也 可以再输入一个...... MA ...

  4. 最小二乘法的一阶、二阶辨识系统

    最小二乘法辨识在控制系统里面运用广泛.下面我分别讲下最小二乘法进行一阶辨识.二阶辨识的公式推导以及利用最小二乘进行参数的估计.为了便于做实验,我设置的输入都为冲击响应,因为冲击响应的s变换就为1. 下 ...

  5. 利用MATLAB求解一阶线性常系数非齐次微分方程组

    用矩阵函数求解一阶线性常系数非齐次微分方程组 主要步骤 1.问题形式 2.求矩阵函数 3.代入矩阵A的指数函数得最终解 主要步骤 本来想用在矩阵论期中开卷考试验证计算结果的,结果一个解方程组的题也没考 ...

  6. matlab 线性最小二乘法,matlab_最小二乘法线性和非线性拟合.ppt

    3.血液容积v, t=0注射剂量d, 血药浓度立即为d/v. 2.药物排除速率与血药浓度成正比,比例系数 k(>0) 模型假设 1. 机体看作一个房室,室内血药浓度均匀--一室模型 模型建立 在 ...

  7. matlab曲面方程,如何用Matlab拟合数据求解曲面方程 - 仿真模拟 - 小木虫 - 学术 科研 互动社区...

    我不知道你的方程,我粗略得到一个R2 = 0.99以上的方程,代码及结果如下 高次多项式拟合代码 : load shuju.mat B=fit([Dotx,Dotx],Dotz,'poly43') 结 ...

  8. matlab 椭圆方程拟合,matlab中如何插值拟合求椭圆方程

    [g_fitting.rar] 使用正交多项式完成数据拟合.程序对读入的gps采样点完成曲线拟合. (2007-08-01, matlab, 1KB, 26次) [曲面拟合.rar] 这是利用matl ...

  9. matlab龙格库塔法求通解,基于matlab及龙格库塔法求解布拉修斯方程.doc

    基于matlab及龙格库塔法求解布拉修斯方程 Runge-Kutta法求解布拉修斯解 摘要 薄剪切层方程主要有三种解法,即相似解,非相似条件下对偏微分方程组的数值解和近似解.布拉修斯解是布拉修斯于19 ...

最新文章

  1. Metadata Service 架构详解 - 每天5分钟玩转 OpenStack(165)
  2. MySQL数据库优化总结
  3. GDI文字旋转90度
  4. xml mysql 模糊查询_mybatis+Spring mysql的模糊查询问题
  5. 分布式系统中,权限设计实践
  6. 【英语学习】【Level 07】U03 Amazing wonders L5 Modern structures in China
  7. 苹果iPhone XI新爆料:用了被小米当噱头的TOF技术
  8. 以太坊区块链浏览器搭建
  9. 微信小程序登陆踩坑记录
  10. Reducer类——hadoop
  11. php免费问答源码,whatsns问答系统PHP免费源码 v4.1
  12. [HTML]解决html5中设置的颜色和浏览器显示的颜色不一致的问题
  13. 志强处理器结尾的ES、QS、正式版的区别
  14. 如何运用python画名字_[原创]如何使用Python在好友画我上画一个标准的汉字震惊朋友圈...
  15. ubuntu工作站配置
  16. haproxy stats enable功能详解
  17. Java是剑客,.NET是刀客
  18. Elasticsearch专栏-1.为什么引入es
  19. 软件测试标准升级|新版25000标准解读
  20. 【成功 debug】Assertion idx_dim >= 0 idx_dim < index_size “index out of bounds“

热门文章

  1. 微信公众平台系列 一
  2. JAVA8用哪个版本的MYSQL_MySQL用哪个版本,5.7还是8.0?
  3. uni-app触发点击事件
  4. lucene java 庖丁解牛_Lucene分词器之庖丁解牛
  5. Dremel made simple with Parquet(CN)
  6. java-net-php-python-ssm巴音学院餐饮安全与卫生防御管理系统计算机毕业设计程序
  7. hdmi怎么支持2k分辨率_选择2K显示器有没有必要?
  8. 疑问代词which/what/who的用法
  9. linux who 时间,linux命令之who、w、whoami
  10. 深入浅出keepalived+nginx实现网关主备高可用