高等数值分析-插值法报告

1、南京理工大学课程考核论文课程名称：高等数值分析论文题 目： 基于 matlab 的函数插值方法性能比较姓名：xxx学号：xxxxxxxxxx成绩：任课教师评语:签名:2016年5月5日摘要函数插值是指已知某函数的在若干离散点上的函数值或者导数信息，通过求解该函数中待定形式的插值函数以及待定系数，使得该函数在给定离散点上满足约束。本文首先介绍了五种插值方法：线性插值、lagrange插值、newdun插值、三次样条插值和三次B样条插值，并从五种插值法的基本思想和具体实例仿真入手，探讨了五种插值 法的优缺点。通过对五种插值法的对比研究及实际应用的总结，从而使我们在以后的应用中能够更好、更快的解决问

2、题 。关键字插值法对比matlab目录 TOC o 1-5 h z HYPERLINK l bookmark2 o Current Document 摘要2. HYPERLINK l bookmark10 o Current Document 0引言4. HYPERLINK l bookmark12 o Current Document 1插值问题的提出、发展史及简单应用 4 HYPERLINK l bookmark14 o Current Document 插值问题的提出4. HYPERLINK l bookmark16 o Current Document 插值法的发展史4. HYPERL

3、INK l bookmark18 o Current Document 插值法的简单应用 4. HYPERLINK l bookmark6 o Current Document 2五种插值法的定义5. HYPERLINK l bookmark8 o Current Document 线性插值法5. HYPERLINK l bookmark20 o Current Document Lagrange 插值法 5. HYPERLINK l bookmark30 o Current Document Newton 插值法6. HYPERLINK l bookmark32 o Current Docu

4、ment 三次样条插值法6. HYPERLINK l bookmark36 o Current Document B样条插值6. HYPERLINK l bookmark43 o Current Document 五种插值法的 matlab仿真实现 .8. HYPERLINK l bookmark45 o Current Document 五种插值方法性能对比 1.1 HYPERLINK l bookmark47 o Current Document 结束语12 HYPERLINK l bookmark49 o Current Document 参考文献120引言近半世纪由于计算机的广泛使用和

5、造船、航空、精密机械加工等世纪问题的需要，使插值法在理论上和实践上得到进一步发展，尤其是20世纪40年代后期发展起来的样条插值等，更获得广泛应用，称为计算机图形学的基础。插值是指已知某函数的在若干离散点上的函数值或者导数信息，通过求解该函数中待定形式的插值函数以及待定系数， 使得该函数在给定离散点上满足约束。 插值函数又叫作基函 数，如果该基函数定义在整个定义域上，叫作全域基，否则叫作分域基。1插值问题的提出、发展史及简单应用插值问题的提出许多实际问题都用函数来表示某种内在规律的数量关系，其中相当一部分函数是通过实验或观测得到的。虽然 f(x应某个区间la,b】上是存在的，有的还是连续的，但却

6、只能给出 b,b】上一系列点Xj的函数值y = f(x (i =0,1,2这只是一张函数表.有的函数虽有解析 表达式，但由于计算复杂，使用不方便，通常也造一个函数表，如大家熟悉的三角函数表、对数表、平方根和立方根表 .为了研究函数的变化规律，往往需要求出不在表中的函数值.因此，我们希望根据给定的函数表做一个既能反映函数f(x)的特性，又便于计算简单函数p(x ),用p(x M似f(x )。通常选一类较简单的函数(如代数多项式或分段代数多项式) 作为f(x ),并使p(x )= f (xi对(i =0,1,2)成立.这样确定的p(x)就是我们希望得到的 插值函数。插值法的发展史插值法是一种古老的

7、数学方法，它来自生产实践.早在一千多年前的隋唐时期制定历法时就应用了二次插值，隋唐刘焯将等距点二次插值应用于天文计算.但插值理论都是在17世纪微积分产生以后才逐渐发展的，牛顿的等距节点插值公式及均差插值公式都是当时的重要成果。插值法的简单应用在现代机械工业中用计算机程序控制加工机械零件，根据设计可给出零件外形曲线的某些型值点仅1=% li =0,1,2,n ),加工时为控制每步走刀方向及步数，就要算出零件外形曲线 其他店的函数值，才能加工出外表光滑的零件。2五种插值法的定义本章比较了五种插值方法，首先介绍了各种插值法的定义。2.1线性插值法线性插值法是指使用连接两个已知量的直线来确定在这两个已

8、知量之间的一个未知量的值的方法。假设我们已知坐标 根据图中所示，我们得到(x0,y2五种插值法的定义本章比较了五种插值方法，首先介绍了各种插值法的定义。2.1线性插值法线性插值法是指使用连接两个已知量的直线来确定在这两个已知量之间的一个未知量的值的方法。假设我们已知坐标 根据图中所示，我们得到(x0,y0)与(xi,yi),要得到x0,xi区间内某一位置x在直线上的值。则直线上，两点间的任意一点y- yox- xoa =yi - yo xi - xox对应的y值可以求出：y = i - a yo + a yi-(2i)-(22)这样通过a就可以直接得到V。Lagrange 插值法Lagrang

9、e插值法是在节点上给出节点基函数，然后做基函数的线性组合，组合系数为节点函数值，这种插值多项式称为拉格朗日插值公式。下面讨论n+i个节点x0 cxixn的n次插值多项式Ln(x)，假定它满足条件Ln xj = yjj =0,i, ,n (2 - 3)为了构造Ln(x ),先定义n次插值基函数.定义：若n次多项式lj(x)(j =0,i，n n+i个节点x0 xixn上满足条件ik = i .lj(xk)=，(j,k = 0,i，n、2 - 4)0,krj就称这n+i个n次多项式lo(x)li(x)，ln(x)为节点xo,xi，xn上的n次插值基函数.推导可得n次插值基函数为lk(x)= (x-

10、x0 -(x-XkJx-Xkdx - Xn)(k = 0,i，,n). (2- 5)xk-xoxk-x-xk-xki xk-xn满足上式的插值多项式可表示nLn x =、ydk x (2- 6)k=0n由 lk(x 2定义知 Ln(Xj)=Z y&(Xj)=yj(j=0,1，n).k=0形如(2 - 6)式的插值多项式 Ln(x)称为Lagrange插值多项式。Newton插值法利用插值基函数很容易得到拉格朗日插值多项式，公式结构紧凑，在理论分析中甚为方便，但当插值节点增减时全部插值基函数均要随之变化，整个公式也将发生变化，这在实际计算中是很不方便的，为了克服这一缺点，提出了牛顿插值。牛顿插值

11、通过求各阶差商，递推得到的一个公式：fn X = f X0 + f X0,X1 X - x0 + fx0,x1,x2 X - X0 X- X1 + ?+ fX0,Xn X- X0 X - Xn- 1-(2)我们称fn x为牛顿Newton插值多项式。E x = f x - N n x = f Lxo , x1x n | . x - x j , x 三 a, b I (2 - 8)为插值余项.三次样条插值法设在区间Nb】上取n+1个节点a = x0 : x1 : xn = b给定这些点函数值 yi = f (xi )若函数Sx瞒足条件：S(X )=yi,i =0,1,2,n;在每个区间 及，xi

12、由( =0,1,2,，n )上是3次多项式；S Xi )三C2 a,b 1;取下列边界条件之一：(i )第一边界条件：S(X0)= f&0)S(Xn)= f M )；(ii )第二边界条件：S(X0)= f&0)S (Xn)= f 8 )或 S(X0)= S (Xn)= 0 ；(iii)周期边界条件：Sk X0 =Sk Xn ,k =1,2,称S仅为3次样条插值函数。B样条插值B样条基函数的定义：给定参数u轴上的一个分割U =Ui(Ui Ui书),(i =0,1,2，m),由下列递推关系定义的Ni,p(u)称为U的p次(p+1阶)B样条基函数。1,ui u ui 1Ni3 0,U-Uiu i

13、_n 1 _ uNi,p(u) = -Ni,p_i(u)p-Ni i,p i(u)(2- 10)ui 1 - uiui-p 1 - ui 1规定0 =00其中，p表示B样条的次数(即为p+1阶)，ui为节点，U为节点矢量。B样条曲线的定义：设Po, p1,Pn为给定空间的n+1个控制顶点，U=“0,5,Um是m+1个节点矢量，称下列参数曲线nC (u) - pi Ni,p(u) a u b(2 - 11) i =0为p次的B样条曲线，折线 p0 , R,.,以为B样条曲线的控制多边形。其中次数p,控制顶点个数n+1,节点个数m+1具有如下关系：m=n+p+1取n=3,则有三次B样条曲线的基函数

14、如下:三次B样条曲线段P0,3取n=3,则有三次B样条曲线的基函数如下:三次B样条曲线段P0,3(t)为:(2 - 12)13)3五种插值法的matlab仿真实现程序简介：将五种插值法写在同一个程序中，并将结果显示在同一个 figure中，可以更好的对比每一种插值法的性能。最后将五种插值法结果显示在同一个分图中，进行插值性能比较。程序中插值原函数：y = x2 ?sin x ,线性插值调用linear函数，三次样条插值调用spline函数，三次 B样条调用spapi函数。%-程序包括五种插值方法：线性插值、三次样条插值、三次 B样条插值、Lagrange插 值、Newton插值clcclear

15、 allclose allx=0:0.4:9;n=length(x);IK=randperm(n);x=x(IK);y=x.A2.*cos(x); m=length(y);IK=randperm(n);x=x(IK);y=x.A2.*cos(x); m=length(y);xi=0:0.02:9;method=linear,spline for k=1:1:2 subplot(3,2,k) yi=interp1(x,y,xi,char(method(k);yg=xi.A2.*cos(xi);plot(x,y,ko,xi,yi,r:,xi,yg,b);text(2,5,char(method(k

16、),Fontsize,12); xlim(min(x),max(x);endsubplot(3,2,3)sp2=spapi(3,x,y);fnplt(sp2,r:); hold on plot(x,y,ko,xi,yg,b);text(2,5,B-spline,Fontsize,12); xlim(min(x),max(x);% set(gcf,Color,w); % x0=xi;% s=0;fori=1:n%线性插值和三次样条插值%计算理论值%画数据点及绘制曲线%3次B样条插值lagrangez=ones(1,length(x0);for j=1:nif j=iz=z.*(x0-x(j)./

17、(x(i)-x(j);%lagrange 插值基函数end ends=s+z*y(i);%lagrange 插值函数end y1=s; % newdun%计算均差表 YY=zeros(n,n); Y(:,1)=y.; for j=2:n for k=j:n Y(k,j)=(Y(k,j-1)-Y(k-1,j-1)/(x(k)-x(k-j+1);endend%计算牛顿插值公式M=length(x0);yn=zeros(1,M);for j=1:Mfori=1:nz=1;for k=1:i-1z=z*(x0(j)-x(k);endyn(j)=yn(j)+Y(i,i)*z;endend% 画图subp

18、lot(3,2,4)plot (x ,y ,ko,x0 ,y1 ,r: ,x0 ,yg ,b);text(2,5,Lagrange,Fontsize,12);xlim(min(x),max(x);subplot(3,2,5)plot (x ,y ,ko,x0 ,yn ,r: ,x0 ,yg ,b);text(2,5,Newton,Fontsize,12);xlim(min(x),max(x);subplot(3,2,6)plot (x ,y ,ko,x0 ,y1 ,y: ,x0 ,yn ,g,x0,yg,b);hold onsp1=csapi(x,y);%三次样条函数fnplt(sp1,k); hold onsp2=spapi(3,x,y);%三次 B 样条插值fnplt(sp2,r); legend (data ,lagrange,newdun,linear,三次样条插值,三次 B 样条插值); 程序运行结果如下：linear-100Lagrange-100-100uaia lagrange newdun-100 linearlinear-100Lagrange-100-100uaia lagrange newdun-100三次样条插值三次 B