精通MATLAB—综合辅导与指南Chap11 曲线拟合与插值.pdfVIP

第 11 章 曲线拟合与插值 在大量的应用领域中，人们经常面临用一个解析函数描述数据(通常是测量值) 的任务。 对这个问题有两种方法。在插值法里，数据假定是正确的，要求以某种方法描述数据点之 间所发生的情况。这种方法在下一节讨论。这里讨论的方法是曲线拟合或回归。人们设法 找出某条光滑曲线，它最佳地拟合数据，但不必要经过任何数据点。图11.1 说明了这两种 方法。标有o的是数据点；连接数据点的实线描绘了线性内插，虚线是数据的最佳拟合。 11.1 曲线拟合 曲线拟合涉及回答两个基本问题：最佳拟合最佳拟合意味着什么？应该用什么样的曲线？可用 最佳拟合最佳拟合 许多不同的方法定义最佳拟合最佳拟合，并存在无穷数目的曲线。所以，从这里开始，我们走向何 最佳拟合最佳拟合 方？正如它证实的那样，当最佳拟合最佳拟合被解释为在数据点的最小误差平方和，且所用的曲线 最佳拟合最佳拟合 限定为多项式时，那么曲线拟合是相当简捷的。数学上，称为多项式的最小二乘曲线拟合。 如果这种描述使你混淆，再研究图 11.1。虚线和标志的数据点之间的垂直距离是在该点的 误差。对各数据点距离求平方，并把平方距离全加起来，就是误差平方和。这条虚线是使 误差平方和尽可能小的曲线，即是最佳拟合最佳拟合。最小二乘这个术语仅仅是使误差平方和最小 最佳拟合最佳拟合 的省略说法。 Second Order Curve Fitting 12 10 8 ) 6 x ( f = y 4 2 0 -2 0 0.2 0.4 0.6 0.8 1 x 图11.1 2 阶曲线拟合 在MATLAB 中，函数polyfit 求解最小二乘曲线拟合问题。为了阐述这个函数的用法， 让我们以上面图11.1 中的数据开始。 » x=[0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1]; » y=[-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; 为了用 polyfit ，我们必须给函数赋予上面的数据和我们希望最佳拟合数据的多项式的 阶次或度。如果我们选择n=1 作为阶次，得到最简单的线性近似。通常称为线性回归。相 反，如果我们选择n=2 作为阶次，得到一个2 阶多项式。现在，我们选择一个2 阶多项式。 » n=2; % polynomial order » p=polyfit(x, y, n) p = -9.8108 20.1293 -0.0317