第九章-插值与拟合
18页1、第九章 插值与拟合插值:求过已知有限个数据点的近似函数。拟合:已知有限个数据点,求近似函数,不要求过已知数据点,只要求在某种意义下它在这些点上的总偏差最小。插值和拟合都是要根据一组数据构造一个函数作为近似,由于近似的要求不同,二者的数学方法上是完全不同的。而面对一个实际问题,究竟应该用插值还是拟合,有时容易确定,有时则并不明显。1 插值方法下面介绍几种基本的、常用的插值:拉格朗日多项式插值、牛顿插值、分段线性插值、Hermite插值和三次样条插值。1.1 拉格朗日多项式插值 1.1.1 插值多项式用多项式作为研究插值的工具,称为代数插值。其基本问题是:已知函数在区间上个不同点处的函数值,求一个至多次多项式 (1)使其在给定点处与同值,即满足插值条件 (2)称为插值多项式,称为插值节点,简称节点,称为插值区间。从几何上看,次多项式插值就是过个点,作一条多项式曲线近似曲线。次多项式(1)有个待定系数,由插值条件(2)恰好给出个方程 (3)记此方程组的系数矩阵为,则 是范德蒙特(Vandermonde)行列式。当互不相同时,此行列式值不为零。因此方程组(3)有唯一解。这表明,只要个节点互不相
2、同,满足插值要求(2)的插值多项式(1)是唯一的。插值多项式与被插函数之间的差称为截断误差,又称为插值余项。当充分光滑时,其中。1.1.2 拉格朗日插值多项式实际上比较方便的作法不是解方程(3)求待定系数,而是先构造一组基函数 是次多项式,满足令 (4)上式称为次Lagrange插值多项式,由方程(3)解的唯一性,个节点的次Lagrange插值多项式存在唯一。 1.1.3 用Matlab作Lagrange插值Matlab中没有现成的Lagrange插值函数,必须编写一个M文件实现Lagrange插值。设个节点数据以数组输入(注意Matlat的数组下标从1开始),个插值点以数组输入,输出数组为个插值。编写一个名为lagrange.m的M文件:function y=lagrange(x0,y0,x);n=length(x0);m=length(x);for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j=k p=p*(z-x0(j)/(x0(k)-x0(j); end end s=p*y0(k)+s; end y(i)=s;end
3、1.2 牛顿(Newton)插值在导出Newton公式前,先介绍公式表示中所需要用到的差商、差分的概念及性质。1.2.1 差商定义 设有函数为一系列互不相等的点,称为关于点一阶差商(也称均差)记为,即 称一阶差商的差商 为关于点的二阶差商,记为。一般地,称 为关于点的阶差商,记为 容易证明,差商具有下述性质: 1.2.2 Newton插值公式线性插值公式可表成 称为一次Newton插值多项式。一般地,由各阶差商的定义,依次可得 将以上各式分别乘以,然后相加并消去两边相等的部分,即得 记显然,是至多次的多项式,且满足插值条件,因而它是的次插值多项式。这种形式的插值多项式称为Newton插值多项式。称为Newton插值余项。Newton插值的优点是:每增加一个节点,插值多项式只增加一项,即因而便于递推运算。而且Newton插值的计算量小于Lagrange插值。由插值多项式的唯一性可知,Newton插值余项与Lagrange余项也是相等的,即由此可得差商与导数的关系其中。1.2.3 差分当节点等距时,即相邻两个节点之差(称为步长)为常数,Newton插值公式的形式会更简单。此时关于节点间函数
4、的平均变化率(差商)可用函数值之差(差分)来表示。定义 设有等距节点,步长为常数,。称相邻两个节点处的函数值的增量为函数在点处以为步长的一阶差分,记为,即 类似地,定义差分的差分为高阶差分。如二阶差分为 一般地,阶差分为 ,上面定义的各阶差分又称为向前差分。常用的差分还有两种: 称为在处以为步长的向后差分; 称为在处以为步长的中心差分。一般地,阶向后差分与阶中心差分公式为 差分具有以下性质: (i)各阶差分均可表成函数值的线性组合,例如 (ii)各种差分之间可以互化。向后差分与中心差分化成向前差分的公式如下: 1.2.4 等距节点插值公式如果插值节点是等距的,则插值公式可用差分表示。设已知节点,则有 若令,则上式又可变形为上式称为Newton向前插值公式。1.3 分段线性插值1.3.1 插值多项式的振荡用Lagrange插值多项式近似,虽然随着节点个数的增加,的次数变大,多数情况下误差会变小。但是增大时,的光滑性变坏,有时会出现很大的振荡。理论上,当,在内并不能保证处处收敛于。Runge给出了一个有名的例子:对于较大的,随着的增大,振荡越来越大,事实上可以证明,仅当时,才有,而在此区间
《第九章-插值与拟合》由会员F****n分享,可在线阅读,更多相关《第九章-插值与拟合》请在金锄头文库上搜索。
2023-03-09 286页
2023-10-27 32页
2023-11-25 51页
2022-10-02 190页
2023-06-24 8页
2024-01-06 6页
2023-02-14 5页
2023-12-30 5页
2023-11-04 7页
2023-06-09 17页