一、MATLAB程序Linear_Model.m
%Linear_Model.m
%@New
clear;clf;
%输入温度t实验数据
t = [-19.75 -10.44 -10.03 -5.00 -0.45 0.04 5.03 9.51 20.51 30.34];
%输入Remington电阻值实验数据
R = [4466.8 3165.3 3211.1 2641.1 2226.4 2266.9 1828.0 1570.6 1115.0 800.5];
a = ones(size(t));
T = 273.16*a+t;%将摄氏温度变为绝对温度
fR=1./log(R);%求电阻值的自然对数的倒数
p=polyfit(fR,T,4);%最小二乘法拟合
R0=800:100:4500;
fR0=1./log(R0);%求电阻值的自然对数的倒数
T0=polyval(p,fR0);
%figure;
plot(R0,T0,R,T,'o');
%plot(fR,T,'o');%绘制热敏电阻温度特性曲线
hold on;
%plot(p(5)+p(4)*fR+p(3)*fR.^2+p(2)*fR*fR*fR+p(1)..,T,'-');%?
xlabel('R/Ommiga ')
ylabel('T/K)')
二、参数估计值:
A4= -3.3877
A3= 1.8276
A2= -0.3697
A1= 0.0334
A0= -0.0011
三、拟合曲线及实验数据点
untitled.bmp
如果疑问,欢迎留言交流,不定期查看谢谢。
分享不易,一分一毛都是激励