待解决的问题:
软测量作业第二题:
2. 炼钢厂出钢时盛钢水的钢包在使用过程中受钢水和炉渣对包衬耐火材料的浸蚀,其容积不断增大。表2所示为使用不同次数x时钢包容积y(由于容积不便测量,故以钢包盛满钢水重量来表示)的某组实测数据。试求两者之间的定量关系。
表2 钢包容积数据表
使用次数x 钢包容积y 使用次数x 钢包容积y
2 106.42 11 110.59
3 108.20 14 110.60
4 109.53 15 110.90
5 109.50 16 110.76
7 110.00 18 111.00
8 109.93 19 111.20
10 110.49
① 画出散点图;
② 构造一个非线性函数;
③ 利用样本数据回归出模型参数;
④ 对实测值和模型估计值进行比较。
matlab代码:
我首先根据几个常见的非线性函数的图像,觉得a+bIn(x)最符合散点图的样子,所以应用在了下面的内联函数inline中。
%软测量课程作业
%第二题,绘制散点图
x = [2,3,4,5,7,8,10,11,14,15,16,18,19];
y= [106.42,108.20,109.53,109.50,110.00,109.93,110.49,110.59,110.60,110.90,110.76,111.00,111.20];
scatter(x,y)
hold on
myfunc=inline('beta(1)+beta(2)*log(x)','beta','x');
beta0=[0.2,0.2];
beta = nlinfit(x,y,myfunc,beta0)
x1=0:1:20;
y1=106.2987+1.7197*log(x1);
plot(x1,y1,'-r')