05第5章插值与拟合资料.pdf

第5章 插值与拟合 数学建模算法与应用 3/63 数学 建模 在实际问题中，一个函数 ( )y f x? 往往是从 实验观测得到的，所知道的仅为函数 ( )f x 在某 区间[ , ]a b 上一系列点上的值 ( )i iy f x? ， 0,1, ,i n? ? ， 当需要在这些结点 0 1, , , nx x x? 之间的点 x上的 函数值时，常用较简单的、满足一定条件的函 数 ( )x? 去代替 ( )f x ，插值法是一种常用方法， 其插值函数 ( )x? 满足条件 ( )i ix y? ? ， 0,1, ,i n? ? . 4/63 数学 建模 n nn xaxaxaaxp ++++= ? 2 210)( ),,2,1,0()()( niyxfxp iiin ?=== 5/63 数学 建模 拟合也是已知有限个数据点，求近似函数，不要求 过已知数据点，只要求在某种意义下它在这些点上的总 偏差最小。 6/63 数学 建模 5.1 插值方法 在工程和数学应用中，经常有这样一类数据处理问 题，在平面上给定一组离散点列，要求一条曲线，把这 些点按次序连接起来，称之为插值。 已知 1n? 个点( , )i ix y （ 0,1, ,i n? ? ），下面求各种插值函 数。 7/63 数学 建模 现象 Runge现象 8/63 数学 建模 5.1.1 拉格朗日插值多项式 拉格朗日插值的基函数为 0 1 1 0 1 1 ( ) ( )( ) ( )( ) ( ) ( )( ) ( ) i i n i i i i i i i n x x x x x x x x l x x x x x x x x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? , ( ). 0 0 1 n j j i j j i x x i , , ,n x x? ? ? ? ? ?? ? ( )il x 是n次多项式，满足 0, , ( ) 1, .i j j i l x j i ? ?????? ??? 拉格朗日插值函数 00 0 ( ) ( ) n n n j n i i i ji i i j j i x x L x y l x y x x?? ? ? ? ??? ? ?? ??? ? ?? ??? ???? ? ? ? ? 9/63 数学 建模 5.1.2 分段线性插值 简单地说，将每两个相邻的节点用直线连起来，如此形 成的一条折线就是分段线性插值函数，记作 ( )nI x ，它满足 ( )n i iI x y? ，且 ( )nI x 在每个小区间 1[ , ]i ix x ? 上是线性函数 ( 0,1, , )i n? ? 。 ( )nI x 可以表示为 0 ( ) ( ) n n i i i I x y l x ? ?? ，其中 1 1 1 1 1 1 . , [ , ] ( 0 , ( ) [ , ] , 0 i i i i i i i i i i i x x x x x i x x x x l x x x x i n x x ? ? ? ? ? ? ? ??? ? ?? ????? ??? ? ??? ????????? 其它 ） ，（） ， 10/63 数学 建模 5.1.3 样条插值 许多工程技术中提出的计算问题对插值函数的 光滑性有较高要求，如飞机的机翼外形，内燃机的进、 排气门的凸轮曲线，都要求曲线具有较高的光滑程度 不仅要连续，而且要有连续的曲率，这就导致了样条 插值的产生。 所谓样条（Spline）本来是工程设计中使用的一种 绘图工具，它是富有弹性的细木条或细金属条。绘图员 利用它把一些已知点连接成一条光滑曲线（称为样条曲 线），并使连接点处有连续的曲率。三次样条插值就是 由此抽象出来的。 11/63 数学 建模 数学上将具有一定光滑性的分段多项式称为样条函数。具 体地说，给定区间[ , ]a b 的一个分划 0 1 1: n na x x x x b? ?? ? ? ? ? ?? . 如果函数 ( )S x 满足 （1）在每个小区间 1[ , ]( 0,1, , 1)i ix x i n? ? ?? 上 ( )S x 是m次 多项式； （2） ( )S x 在[ , ]a b 上具有 1m? 阶连续导数。 则称 ( )S x 为关于分划?的m次样条函数，其图形为m次样条曲 线。 显然，折线是一次样条曲线。 12/63 数学 建模 5.1.3.2 三次样条插值 利用样条函数进行插值，即取插值函数为样条函数， 称为样条插值。例如分段线性插值是一次样条插值。这