插值与拟合模型二题库.pptVIP

作业： 在某一反应工程实验中，测得了如下表所示的实验数据 序号 1 2 3 4 5 6 7 8 温度T 10 20 30 40 50 60 70 80 转化率y 0.1 0.3 0.7 0.94 0.95 0.68 0.34 0.13 现在要确定在其他条件不变的情况下，转化率y和温度T的具体关系，现拟用两种模型去拟合实验，两种模型分别是 试求取上述模型中的参数，并绘图表示拟合效果。 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * 三、拟合的使用及求解 3.1 引言 数据拟合就是从一组数据出发，寻找函数y=f(x)的一个近似表达式 ，要求该函数在某种准则下能尽量反应数据的整体变换趋势，而不一定要经过所有数据点(xi,yi).数据拟合问题也称为曲线拟合问题，而y=f(x)称为拟合曲线。 3.2 拟合模型 直 线 拟 合 问 题 引 例 1 温度t(oC) 20.5 32.7 51.0 73.0 95.7 电阻R(?) 765 826 873 942 1032 已知热敏电阻数据： 求60oC时的电阻R． 在直角坐标系下作图如下(plot) x=[20.5 32.7 51.0 73.0 95.7]; y=[765 826 873 942 1032]; plot(x,y,’+’) xlabel(‘温度’)，ylabel(‘电阻’) 设 R=a1t+a2 a1,a2为待定系数 曲 线 拟 合 问 题 引 例 2 t (h) 0.25 0.5 1 1.5 2 3 4 6 8 c (?g/ml) 19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01 已知一室模型快速静脉注射下的血药浓度数据(t=0注射300mg) 求血药浓度随时间的变化规律c(t). 在直角坐标系下作图如下(plot) 拟合曲线： 方案一 抛物线 a0,a1,a2为待定系数 方案二 c0,k为待定系数 曲 线 拟 合 问 题 的 提 法 已知一组（二维）数据，即平面上 n个点（xi,yi) i=1,…,n, 寻求一个函数（曲线）y=f(x), 使 f(x) 在某种准则下与所有数据点最为接近，即曲线拟合得最好． + + + + + + + + + x y y=f(x) (xi,yi) ?i ?i 为点（xi,yi) 与曲线 y=f(x) 的距离 最小二乘准则： 使n个点（xi,yi) 与曲线 y=f(x) 的距离?i 的平方和最小 ． 记 根据最小二乘准则确定拟合函数的方法称为最小二乘法。 最小二乘法又分为线性最小二乘拟合和非线性最小二乘拟合 线性最小二乘拟合：拟合函数的待定系数a1,a2, …,am+1全部以线性形式出现 例如拟合函数：f(x)=a1xm+ …+amx+am+1 非线性最小二乘拟合：拟合函数的待定系不能全部以线性形式出现 例如拟合函数： 注：最小二乘法中拟合函数的确定 最小二乘法中确定拟合函数是很关键的，通常有以下两种形式， 1：如果没有现成的规则，可以通过散点图，结合曲线的形状变化趋势进行分析，建立经验模型； 2：可以通过机理分析建立数学模型来确定拟合函数，如人口增长的Logistic模型就是通过机理分析法推导出来的。 用MATLAB作线性最小二乘拟合 式f(x)=a1xm+ …+amx+am+1拟合,可利用已有程序: 在MATLAB中线性最小二乘拟合就是做多项 a=polyfit(x,y,m) 输入同长度 的数组x，y 拟合多项 式次数 说明:1.a-输出拟合多项式系数a=[a1,a2, …, am+1] 2.多项式在x处的值y可用以下命令计算： y=polyval（a, x） 温度t(oC) 20.5 32.7 51.0 73.0 95.7 电阻R(?) 765 826 873 942 1032 例：已知热敏电阻数据： 求60oC时的电阻R． x=[20.5 32.7 51.0 73.0 95.7];