MATLAB线性参数模型估计

一、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

如果疑问，欢迎留言交流，不定期查看谢谢。

                                   分享不易，一分一毛都是激励