第4章插值与曲线拟合讲述.ppt

第4章　插值法与曲线拟合 南京中医药大学信息技术学院 制作:张季 第4章　插值法与曲线拟合 　4.1　Lagrange插值法 　4.2　埃特金算法 　4.3　Newton插值法 　4.4　差分与等距结点插值 　4.5　埃尔米特插值法 4.6　有理分式插值法 4.7　函数逼近 4.8　曲线拟合 内容提要： 概述 设 是在 上有一定光滑性的函数， 是 上 个互不相同的点， 在这 个点上的取值分别 。所谓的插值法就是求一个形式简单的函数 ,使得 i=0,1,2,…,n 通常称给定的点 为插值节点，称函数 为函数 关于插值节点xi的插值函数(插值多项式)，称 为被插函数。 插值多项式的存在唯一性 4.1.1　线性插值 程序框图设计 4.1.2　抛物插值 程序框图设计 4.1　Lagrange插值方法 4.1.3　插值余项 4.2　分段插值 选择插值结点的原则 4.2.1　分段线性插值 4.3.2　分段抛物插值 4.3　Newton插值方法 差商的性质 性质1： 性质2： 性质3： 性质4： 性质5： 差商的计算 牛顿插值多项式 牛顿插值程序框图 例题 已知x=1,4,9 y=1,2,3 求 4.4　等距结点插值 差分的性质 性质1 性质2 性质3 性质4 等距结点插值公式 　　将牛顿插值多项式中各阶差商用相应差分代替，就可以得到相应的等距结点插值公式。 牛顿前插公式： 牛顿后插公式： 向前差分递推表 程序框图 小结 习题 用C语言、matlab写出线性插值、抛物插值、一般拉格朗日插值的源程序代码。 用C语言、matlab写出牛顿插值的源程序代码。 用C语言、matlab写出分段线性插值、分段抛物插值的源程序代码。 3 2 1 9 4 1 x x 1 2 3 1 4 9 二阶差商 一阶差商 x 差分： 一般的，n阶差分递推地定义如下： 即： 并规定　　　　，称之为零阶差分。 设等距结点xi=x0+ih(i=0, ……)相应的函数值为yi=f(xi),并定义　　　　　　　　　　　　　 为其一阶差分。 向后差分： 四阶 三阶 二阶 一阶 函数值 算法设计： 1)升阶。 2)使a[i]=f[i],y=y1。 3)进入循环，使用公式f[i]=a[i]-a[i-1]计算差分值。 4)在K循环中计算连乘积q=q·(t-k+1) 5)在K循环中计算阶乘p=p·k 6)计算插值结果y1=y1+q·f(k)/p 7)判断 |y-y1|E是否成立， 若成立，则程序转去输出并结束， 　　若不成立，则返回步骤1执行。 使用等距结点插值法求插值结点x处的函数值 输入 使用等距结点插值法求插值结点x处的函数值 输出 升阶 使a[i]=f[i],y=y1 在i循环中计算差分 f[i]=a[i]-a[i-1] 计算连乘积q=q·(t-k+1) 计算阶乘值p=p·k 计算插值结果y1=y1+q·f(k)/p 判断 |y-y1|E是否成立， 若成立，则程序转去输出并结束，若不成立，则返回步骤1执行 0.82534 0.87758 0.92106 0.95534 0.98007 0.99500 1.00000 cosx 0.6 0.5 0.4 0.3 0.2 0.1 0.0 x 0.00001 -0.00002 -0.00001 0.00012 0.00010 0.00009 0.00013 0.00025 0.00035 0.00044 -0.00993 -0.00980 -0.00955 -0.00920 -0.00876 -0.00500 -0.01493 -0.02473 -0.03428 -0.04348 -0.05224 1.00000 0.99500 0.98007 0.95534 0.92106 0.87758 0.82534 0 0.1 0.2 0.3 0.4 0.5 0.6 Δ6y Δ5 y Δ4 y Δ3 y Δ2 y Δy cosx x 求cos0.048和cos0.575之值差分表如下：　 计算cos0.048时，用牛顿向前插值公式，取 代入公式得： 计算cos0.575时，用牛顿向后插值公式，取 代入公式得： 　　本章按基本思想，计算公式、算法设计、程序框图设计和源程序的顺序分别介绍了线性插值、抛物插值、一般拉格朗日插值、分段线性和抛物插值、牛顿插值等插值方法。　　 * * 　　插值是数值逼近的重要方法之一。是根据给定的一组自变量和函数值，求取未给出的结点处函数值的近似值。本章主要介绍拉格朗日插值，牛顿插值和分段低阶多项式