实验三插值.docVIP

试验三Matlab插值与拟合 实验目的 了解插值的基本内容和原理；掌握用matlab求解插值问题,包括一维插值和二维插值的各种常用方法； 实验原理和方法 插值法是实用的数值方法，是函数逼近的重要方法。在生产和科学实验中，自变量x与因变量y的函数y = f(x)的关系式有时不能直接写出表达式，而只能得到函数在若干个点的函数值或导数值。当要求知道观测点之外的函数值时，需要估计函数值在该点的值。 如何根据观测点的值，构造一个比较简单的函数y=φ(x)，使函数在观测点的值等于已知的数值或导数值，进而用简单函数y=φ(x)在点x处的值来估计未知函数y=f(x)在x点的值。寻找这样的函数φ(x)，办法是很多的。φ(x)可以是一个代数多项式，或是三角多项式，也可以是有理分式；φ(x)可以是任意光滑（任意阶导数连续）的函数或是分段函数；函数类的不同，自然地有不同的逼近效果。 一维插值 已知个数据节点: 构造一个(相对简单)函数(称为插值函数),通过全部结点即 （j＝0,1,…n）再用计算插值，即 数学上插值方法非常多，只介绍几种常用方法： （1）拉格朗日(Lagrange)插值 已知函数f(x)（称为被插值函数）在n+1个点x0,x1,…,xn处的函数值为 y0,y1,…,yn 。求一n次多项式函数Pn(x)（称为插值函数），使其满足： Pn(xi)=yi,i=0,1,…,n.. 解决此问题的拉格朗日插值多项式公式如下： 其中为次多项式： 称为拉格朗日插值基函数,可以验证该多项式通过所给数据点。 特别地: 两点一次(线性)插值多项式为： 三点二次(抛物)插值多项式为: matlab没有Lagrange插值的函数，需要自己编写。 Lagrange插值插值在头尾附近会出现一些不好的振荡现象（龙格现象）Lagrange插值如上缺点的原因在于：多项式次数过高不一定是好事。而克服这一缺点的办法之一就是采用分段低次插值，其中最简单的是分段线性插值。 （2）分段线性插值 分段线性插值就是在每个子区间上做通过两个端点的直线。 该方法具有良好的收敛性，即，其中为被插值函数，可以通过增加插值节点的个数来控制插值误差，是一种较实用的插值方法，Matlab有分段线性插值函数，但总体光滑程度不够。 （3）分段三次埃尔米特（Hermite）插值 对插值函数的要求可以根据需要进行改变，由此产生了不同的插值方法，分段三次埃尔米特插值就是其一。在插值问题中，如果除了指定在插值节点的函数值以外，同时指定插值节点处的导数值，这时插值问题变为： 给定一组观察数据及，要求一个分段多项式函数，满足： 这就相当于所求分段多项式在每个小段上满足四个条件： 可以确定四个待定参数，三次多项式恰好有四个系数，所以可以考虑用三次多项式作为插值函数，这就是分段三次Hermite插值。 （4）分段三次样条插值 在数学上，光滑程度的定量描述是：函数(曲线)的k阶导数存在且连续，则称该曲线具有k阶光滑性。 光滑性的阶次越高，则曲线光滑程度越好。分段线性插值具有零阶光滑性，也就是不光滑；分段三次埃尔米特插值具有一阶光滑性，但是这样的光滑程度在工程设计和机械加工等实际应用中是远远不够的。是否存在较低次的分段多项式达到较高阶光滑性的方法？三次样条插值就是一个很好的例子。 什么是样条？所谓样条，就是一个细的、可弯曲的木制或塑料条，在飞机或轮船等的设计制造过程中为描绘出光滑的外形曲线(放样)所用的工具。1946年,Schoenberg将样条引入数学,即所谓的样条函数。三次样条函数本质上是一段一段的三次多项式拼合而成的曲线，在拼接处,不仅函数是连续的,且一阶和二阶导数也是连续的。 满足以上条件的分段函数称为三次样条插值函数。 可以看出，三次样条插值是分段三次多项式插值，但是不仅要求过所给数据节点，而且要求在所给节点处的二阶导数也是连续的，从而构造出来的插值函数比较光滑。 在n个小区间构造S(x)，共有n个三次多项式，每个多项式有4个系数需要确定，总共需确定4n个参数才能给出S(x) （1）在n+1个节点上要求S(x)过点：， 该条件给出n+1个方程。 （2）在n-1个内部节点上要求连续，一阶导数和二阶导数连续： （），该条件给出3（n-1）个方程。 由（1）和（2）总共得到4n-2个方程，还差两个。为此常用的方法是在两个边界节点和附加要求，这就是所谓的边界条件。根据实际问题的不同，三次样条插值常用到下列三类边界条件。 （1）边界条件：。即指定两个边界节点的一阶导数值为给定值 （2）边界条件：。即指定两个边界节点的二阶导数值为给定值 特别地，当和都为零时，称为自然边界条件。 (3)周期性边界条件： 总而言之，理论上三次样条插值是确定的。 matlab 插值命令 命令1 interp1 功能 一维数据插值。该命令