数值分析常用的插值方法.docx

1、数值分析常用的插值方法数值分析报告班级: 专业: 流水号: 学号: 姓名:常用的插值方法序言在离散数据的基础上补插连续函数，使得这条连续曲线通过全部给定的离散数 据点。插值是离散函数逼近的重要方法，利用它可通过函数在有限个点处的取值状 况，估算出函数在其他点处的近似值。早在6世纪，中国的刘焯已将等距二次插值用于天文计算。17世纪之后，牛顿、 拉格朗日分别讨论了等距和非等距的一般插值公式。在近代，插值法仍然是数据处 理和编制函数表的常用工具，又是数值积分、数值微分、非线性方程求根和微分方 程数值解法的重要基础，许多求解计算公式都是以插值为基础导出的。插值问题的提法是：假定区间a, b上的实值函数

2、f (x)在该区间上n + 1个 互不相同点Xo，X1 Xn处的值是f ( Xo)， f ( Xn)，要求估算f (x)在a，b中某点的值。其做法是：在事先选定的一个由简单函数构成的有 n + 1个参数Co，Ci,Cn的函数类0(Co, Ci，Cn)中求出满足条件P(Xi) = f (Xi) (i = 0，1, n)的函数P(x),并以P(x)(乍为f(x)的估值。此处f (x)称为被插值函数，xo，X1，Xn 称为插值结(节)点，(Co, C1,Cn)称为插值函数类，上面等式称为插值条件， (Co,Cn)中满足上式的函数称为插值函数，R (x)= f (x)- P ( X)称为插 值余项。求

3、解这类问题，它有很多种插值法，其中以拉格朗日 (Lagra nge)插值和牛顿(Newt on)插值为代表的多项式插值最有特点，常用的插值还有 Hermit插值，分段插值和样条插值。一.拉格朗日插值1问题提出：已知函数y f x在n+1个点Xo*丄x上的函数值y,yi丄，y*，求任意一点X的函数值f X。说明：函数y f x可能是未知的；也可能是已知的，但它比较复杂，很难计算其函数值f X2.解决方法:3.f X，则构造一个n次代数多项式函数Fn x来替代未知(或复杂)函数yao , ai , a2 ,L , an。3.构造Pn X的依据:X。根据这个条件，可以写出非齐次线性方程组:aoax2

4、a2XoLnanXoyoaoa1x12a2xLnanX1y1LLaoax2a?XnLnanXnynR x逼近于原来的函数其系数矩阵的行列式D为范德萌行列式:1 Xo Xo2 L1 x1 x12 LM M M1 Xn Xn2 L故当n+1个点的横坐标X。丄，Xn各不相同时，方程组系数矩阵的行列式D不等于 零，故方程组有唯一解。即有以下结论。结论：当已知的n+1个点的横坐标Xo,Xi,L ,xn各不相同时，则总能够构造唯一的 n 次多项式函数Pn X，使Pn X也过这n + 1个点。4.几何意义5.举例：已知函数f x 二，求f 115。分析：本题理解为，已知“复杂”函数f x ，当x=81,10

5、0,121,144时，其对 应的函数值为：y=9,10,11,12，当x=115时，求函数值f 115 。解:（1）线性插值：过已知的（100,10 ）和（121,11）两个点，构造1次多项式函数P x，于是有（2）抛物插值：构造2次多项式函数巳x，使得它过已知的（100,10）（121,11）和（144,12）三个点。于是有2次拉格朗日插值多项式：则有其表达式为:lk xX Xjj 0 Xk Xjj k拉格朗日公式特点:2.每一项中的分子是关于x的n次多项式，分母是一个常数;3.每一项的分子和分母的形式非常相似，不同的是:分子是x L ，而分母是Xk7.误差分析（拉格朗日余项定理）n 1 n

6、X Xkk 0其中 在X0,X1,L ,Xn, X所界定的范围内针对以上例题的线性插值，有R 115 f 115 115 100 115 1212!函数f x在100,115区间绝对值的极大值为f 100 2.5 10 4,则有：R 115 f 115 0.01125 0.05于是近似值f 115 R 115 10.71428571428572有三位有效数字。针对以上例题的抛物线插值，有fF2 115 f 115 115 100 115 121 115 1443!函数f x在100,115区间绝对值的极大值为f 100 3.75 10 6，则有F2 115 f 115 0.00163125V0

7、.005于是近似值f 115 R2 115 10.72275550536420有四位有效数字。8.拉格朗日插值公式的优点公式有较强的规律性，容易编写程序利用计算机进行数值计算。9.拉格朗日插值通用程序程序流程图如下：y输入nxi,yi, (i=0,1,n )t (即插值点x)开始p=0, k = 0k=ny计算l(k)p=p+l*ykk=k+1输出P结束l = 1,=0 njkyl = l*(t-xj)/(xk-xj)j=j+1j=k+1nj yl = l*(t-xj)/(xk-xj)j=j+1文件lagrange.m 女口下：格朗日插值close alln=input(已知的坐标点数n=?)

8、;x=i nput(x1,x2,.,x n=?);y=i nput(y1,y2,.,y n=?);%定义t为符号量xx=input(插值点=?);syms tp=0;for k=1: nl=1;for j=1:k-1l=l*(t-x(j)/(x(k)-x(j);endfor j=k+1: nI=l*(t-x(j)/(x(k)-x(j);endp=p+l*y(k);endp=inline(p); %把符号算式p变为函数形式fplot(p,mi n(mi n(x),xx)-1,max(max(x),xx)+1); %画多项式函数hold onp(xx) %显示插值点plot(x,y,o,xx,p(

9、xx),*); %画已知点和插值点在MATLAB命令窗口输入：lagra nge然后有以下对话过程和结果，已知的坐标点数n=?6x1,x2,.,x n=?1,3,5,7,9,11y1,y2,.,y n=?-1,20,0,-1,12,3插值点=? 8ans =5.67187500000000有以下图形：二.牛顿插值拉格朗日插值的缺点：无承袭性（继承性）若算出3点的抛物插值精度不够，再进行 4点的3次多项式插值时，必须从头 算起，前面算出的3点抛物插值的计算结果不能利用。而泰勒插值却是具有承袭性的，如线性插值的结果不精确，那么再加上一项，就变成了泰勒抛物插值，如：泰勒1次插值：R X f Xo f

10、 Xo x Xo泰勒2次插值：P2 X f Xo f Xo X Xo而牛顿插值就是具有承袭性的插值公式1.差商的概念设n+1个点Xo, 凡互不相等，则定义：Xi和Xj i j两点的一阶差商为：f X, Xj、亠、 fXi , Xj , Xk 二点的二阶差商为： f Xi, Xj , Xk Xi ,Xj,Xk,X|四点的三阶差商为：f Xi,Xj,Xk,X|Xo 2X Xo O2!Xif XjXi XjXi.Xj f Xj. XkXi Xkn+1个点Xo,捲丄，Xn的n阶差商为:Xo-XL ,Xnf X0,X1,L ,Xn ! f N，X2，L ,XnX Xn差商具有对称性：f Xj,Xj f

11、Xj,Xj ； f Xj,Xj,Xk f Xj-XXk2.牛顿插值解决的问题与拉格朗日插值解决的问题相同只是表述n次多项式pn X的公式不同。3.牛顿插公式的推导根据差商的概念，有：f X f Xo f X, Xo X Xo f X,Xo 是 X, Xo 两点的一阶差商;把以上各式从后向前逐次代入，可以得到:Fn X Rn Xf X, Xo, X,丄，Xn X X X Xj L X Xn以上Pn X的表达式称为牛顿插值公式，可以证明， n次牛顿插值多项式与 n次拉格朗日插值多项式完全相同，只是表达形式不同。故，拉格朗日余项定理与牛顿余项定理相同：Rn X Pn XX, Xo,Xi 丄XnX X

12、kTTkoX Xk，其中 在Xo,Xi,L ,Xn, X所界定的范围内。、 f n 1则有公式：f X, Xo, Xi,L , Xnn 1 !4.牛顿插值差商表Xiyi一阶差商二阶差商n阶差商*x0yo1x1y1fx0,x1(x-x0)x2y2fx1,x2fx0,x1,x2(x-x0)(x-x1)x3y3fx2,x3fx1,x2,x3(x-x0)(x-x2)xn-1yn-1xnynfxn -1,x nfxn-2,x n-1,x nfx0,xn(x-x0)(x-xn-1)5举例已知函数f(x)当x=-2,-1,0,1,2时，其对应函数值为f(x)=13,-8,-1,4,1。求f(0.5)的值。解

13、：该题目与例1相比，就是多了一个点，所以和例1的差商表相比，只需多一列, 多一行：xiyi一阶差商二阶差商三阶差商四阶差商*-2131-1-8-21(x+2)0-1714(x+2)(x+1)145-1-5(x+2)(x+1)x21-3-4-11(x+2)(x+1)x(x-1)而5个点的4次牛顿插值多项式P4 x是在P3 x的基础上多增加1项:F4 x 13 21x2 14x2x1 5x2x1x x2x1xx1则f 0.5 P4 0.5 13 21 0.5 2 14 0.5 2 0.5 1 5 0.5 2 0.5 1 0.5 0.5 2 0.5 1 0.5 0.52.6875可以在MATLAB下

14、运行程序newton02.m :p4=i nli ne(13-21*(x+2)+14*(x+2)*(x+1)-5*(x+2)*(x+1)*x+(x+2)*(x+1)*x*(x-1);fplot(p4,-2.5,2.5,r);hold onxi=-2,-1,0,1,2;yi=13,-8,-1,4,1;plot(xi,yi, *);plot(0.5,p4(0.5),o);可以得到以下图形：6.牛顿插值的优点（1） 具有承袭性质（2） 利用差商表，计算多点插值，比拉格朗日公式计算方便。7.牛顿插值算法的通用程序以下是程序流程图:MATLAB的通用程序newton.m为：%牛顿插值close alln

15、=input(已知的坐标点数n=?);x=i nput(x1,x2,.,x n=?);y=i nput(y1,y2,.,y n=?);xx=input(插值点=?);% 计算差商：fx1,x2,fx1,x2,x3,.,fx1,x2,.,x nf=y;for i=1: n-1%计算第i阶差商for k=n :-1:i+1f(k)=(f(k)-f(k-1)/(x(k)-x(k-i);endendsyms t %定义t为符号量P=f(1)；for k=2:n1=1;for j=1:k-1l=l*(t-x(j);end p=p+l*f(k);endp=i nlin e(p);%把符号算式p变为函数形式

16、fplot(p,mi n(min( x),xx)-1,max(max(x),xx)+1);hold on%画多项式函数p(xx)%显示插值点plot(x,y,o,xx,p(xx),*);%画已知点和插值点在MATLAB命令窗口输入：n ewt on然后有以下对话过程和结果,已知的坐标点数n=?6x1,x2,.,x n=?1,3,5,7,9,11y1,y2,.,y n=?-1,20,0,-1,12,3插值点=? 8 ans =5.67187500000000有以下图形：三.总结和展望插值与逼近都是指用某个简单的函数在满足一定条件下在某个范围内近似代替 另一个较复杂的函数或解析表达式未能给出的函数

17、，以便于简化对后者的各种计算 或揭示后者的某些性质。插值方法理论是近似计算和逼近函数的有效方法。此外，它也是数值微积分， 微分方程数值解等数值分析的基础。在图形处理等很多需要优化的实际中，也有着 很广泛的应用。我们期望在以后的生活中会更加熟练和更好的运用插值方法。参考文献1李庆扬，王能超,易大义数值分析M.:华中科技大学出版社，1982.2 吴才斌.插值方法J.湖北大学成人教育学院学报，1999,(5).3徐萃薇，孙绳武.计算方法引论M.:高等教育出版社,2002.4 林鹭.拉格朗日插值多项式的一种并行算法J.厦门大学学报：自然科学版，2004,43(5):592-595.5吴筑筑.计算方法M.:清华大学出版社，2004:61-84.6 杨士俊，王兴华.Hermite插值多项式的差商表示及其应用J.高校应用数学学 报 A 辑,2006,21(1):70-78.7 齐东旭，李华山.数据逼近的多结点样条技术J.中国科学(E辑),1999,29：334-387.8王仁宏.数值逼近M.:高等教育出版社,1999.9 孙亮.数值分析方法课程的特点与思想J.工科数学,2002,18(1):84-86.