文档介绍：MATLAB中的曲线拟合与插值
在大量的应用领域中,人们经常面临用一个解析函数描述数据(通常是测量值)的任务。对这个问题有两种方法。在插值法里,数据假定是正确的,要求以某种方法描述数据点之间所发生的情况。这种方法在下一节讨论。这里讨论的方法是曲线拟合或回归。人们设法找出某条光滑曲线,它最佳地拟合数据,但不必要经过任何数据点。。标有'o'的是数据点;连接数据点的实线描绘了线性内插,虚线是数据的最佳拟合。
曲线拟合
曲线拟合涉及回答两个基本问题:最佳拟合意味着什么?应该用什么样的曲线?可用许多不同的方法定义最佳拟合,并存在无穷数目的曲线。所以,从这里开始,我们走向何方?正如它证实的那样,当最佳拟合被解释为在数据点的最小误差平方和,且所用的曲线限定为多项式时,那么曲线拟合是相当简捷的。数学上,称为多项式的最小二乘曲线拟合。如果这种描述使你混淆,。虚线和标志的数据点之间的垂直距离是在该点的误差。对各数据点距离求平方,并把平方距离全加起来,就是误差平方和。这条虚线是使误差平方和尽可能小的曲线,即是最佳拟合。最小二乘这个术语仅仅是使误差平方和最小的省略说法。
2阶曲线拟合
在MATLAB中,函数polyfit求解最小二乘曲线拟合问题。为了阐述这个函数的用法,让我们以上面图
。
» x=[0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1];
» y=[-.447 ];
为了用polyfit,我们必须给函数赋予上面的数据和我们希望最佳拟合数据的多项式的阶次或度。如果我们选择n=1作为阶次,得到最简单的线性近似。通常称为线性回归。相反,如果我们选择n=2作为阶次,得到一个2阶多项式。现在,我们选择一个2阶多项式。
» n=2; % polynomial order
» p=polyfit(x, y, n)
p =
- -
polyfit 的输出是一个多项式系数的行向量。其解是y = - +-。为了将曲线拟合解与数据点比较,让我们把二者都绘成图。
» xi=linspace(0, 1, 100); % x-axis data for plotting
» z=polyval(p, xi);
为了计算在xi数据点的多项式值,调用MATLAB的函数polyval。
» plot(x, y, ' o ' , x, y, xi, z, ' : ' )
画出了原始数据x和y,用'o'标出该数据点,在数据点之间,再用直线重画原始数据,并用点' : '线,画出多项式数据xi和z。
» xlabel(' x '), ylabel(' y=f(x) '), title(' Second Order Curve Fitting ')
将图作标志。。
多项式阶次的选择是有点任意的。两点决定一直线或一阶多项式。三点决定一个平方或2阶多项式。按此进行,n+1数据点唯一地确定n阶多项式。于是,在上面的情况下,有11个数据点,我们可选一个高达10阶的多项式。然而,高阶多项式给出很差的数值特性,人们不应选择比所需的阶次高的多项式。此外,随着多项式阶次的提高,近似变得不够光滑,因为较高阶次多项式在变零前,可多次求导。例如,选一个10阶多项式
» pp=polyfit(x, y, 10) ;
» format short e % change display format
» pp.' % display polynomial coefficients as a column
ans =
-+005
+006
-+006
+006
-+006
+006
-+005
+005
-+004
+002
--001
要注意在现在情况下,多项式系数的规模与前面的2阶拟合的比较。还要注意在最小(--001)和最大(+006)系数之间有7个数量级的幅度差。将这个解作图,并把此图与原始数据及2阶曲线拟合相比较,结果如何呢?
» zz=polyval(pp, xi); % evaluate 10th order polynomial
» plot(x, y,