《数值分析》-- 拉格朗日插值

文章目录

  • 问题
  • 一、拉格朗日插值基函数
  • 二、拉格朗日插值多项式
  • 三、n次Lagrange插值多项式余项
  • 习题
  • 总结


问题

《数值分析》-- 拉格朗日插值_第1张图片


一、拉格朗日插值基函数

《数值分析》-- 拉格朗日插值_第2张图片
《数值分析》-- 拉格朗日插值_第3张图片
《数值分析》-- 拉格朗日插值_第4张图片

  • n=1时一次基函数
    《数值分析》-- 拉格朗日插值_第5张图片
  • 两点线性插值问题
    《数值分析》-- 拉格朗日插值_第6张图片
  • 问题:即已知函数 f(x)在点 x 0 x_0 x0 x 1 x_1 x1点的函数值
    y 0 y_0 y0=f( x 0 x_0 x0), y 1 y_1 y1=f( x 1 x_1 x1).
    求线性函数 L 1 L_1 L1(x)= a 0 a_0 a0+ a 1 a_1 a1x
    使满足条件: L 1 L_1 L1( x 0 x_0 x0)= y 0 y_0 y0, L 1 L_1 L1( x 1 x_1 x1)= y 1 y_1 y1

在这里插入图片描述
《数值分析》-- 拉格朗日插值_第7张图片

①:把x= x 0 x_0 x0带入 L 1 L_1 L1(x)中, L 1 ( x 0 ) L_1(x_0) L1(x0) = y 0 y_0 y0 + 0 = y 0 y_0 y0
②:把x= x 1 x_1 x1带入 L 1 L_1 L1(x)中, L 1 L_1 L1( x 1 x_1 x1) = y 0 y_0 y0 + ( y 1 y_1 y1- y 0 y_0 y0)( x 1 x_1 x1- x 0 x_0 x0) / ( x 1 x_1 x1- x 0 x_0 x0) = y 0 y_0 y0 + y 1 y_1 y1 - y 0 y_0 y0 = y 1 y_1 y1

《数值分析》-- 拉格朗日插值_第8张图片
在这里插入图片描述

则称 l 0 l_0 l0(x)叫做点 x 0 x_0 x0的一次插值基函数, l 1 l_1 l1(x)叫
做点 x 1 x_1 x1的一次插值基函数

《数值分析》-- 拉格朗日插值_第9张图片

当x= x 0 x_0 x0时, l 0 l_0 l0=1; x= x 1 x_1 x1时, l 0 l_0 l0=0;
当x= x 0 x_0 x0时, l 1 l_1 l1=0; x= x 1 x_1 x1时, l 1 l_1 l1=1;

  • n=2时二次基函数
    《数值分析》-- 拉格朗日插值_第10张图片

二、拉格朗日插值多项式

《数值分析》-- 拉格朗日插值_第11张图片

对上述举例:

begin---------------------------------------------

  • n=1时
    在这里插入图片描述

①:把x= x 0 x_0 x0带入 L 1 L_1 L1(x)中,L1(x0) = y 0 y_0 y0 + 0 = y 0 y_0 y0
②:把x= x 1 x_1 x1带入 L 1 L_1 L1(x)中,L1(x1) = y 0 y_0 y0 + ( y 1 y_1 y1- y 0 y_0 y0)( x 1 x_1 x1- x 0 x_0 x0) / ( x 1 x_1 x1- x 0 x_0 x0) = y 0 y_0 y0 + y 1 y_1 y1 - y 0 y_0 y0 = y 1 y_1 y1
③:由此推广到 L n L_n Ln(x)的多项式, L n L_n Ln( x j x_j xj) = y j y_j yj j=0,1,2,…,n

end-----------------------------------------------

再由插值多项式的唯一性: P n P_n Pn( x) = L n L_n Ln( x)

特别地:
n =1时又叫线性插值,其几何意义为过两点的直线.
n =2时又叫抛物(线)插值, 其几何意义为过三点的抛物线.

  • 注意
    《数值分析》-- 拉格朗日插值_第12张图片
    习题

三、n次Lagrange插值多项式余项

截断误差 R n ( x ) R_n(x) Rn(x)=f(x) - L n L_n Ln(x)也称为n次Lagrange插值多项式的余项

  • 拉格朗日余项定理
    《数值分析》-- 拉格朗日插值_第13张图片

ω n \omega_n ωn + _+ + 1 _1 1( x x x) = ( x x x- x 0 x_0 x0)( x x x- x 1 x_1 x1)( x x x- x 2 x_2 x2)…( x x x- x n x_n xn)
由给定条件可知 R n ( x ) R_n(x) Rn(x)在节点 x k x_k xk( k k k=0,1,…,n)上为0,即 R n ( x ) R_n(x) Rn(x) = 0 ( k k k=0,1,…,n)
插值节点 x i x_i xi 上误差等于零

《数值分析》-- 拉格朗日插值_第14张图片
《数值分析》-- 拉格朗日插值_第15张图片

  • 截断误差限⭐(求误差会用到该公式)
    《数值分析》-- 拉格朗日插值_第16张图片
  • 结论例题⭐
    《数值分析》-- 拉格朗日插值_第17张图片
    习题

习题

  • 例题
  1. 在这里插入图片描述
    《数值分析》-- 拉格朗日插值_第18张图片
  2. 在这里插入图片描述
    《数值分析》-- 拉格朗日插值_第19张图片
    《数值分析》-- 拉格朗日插值_第20张图片
  • 例题
  1. 《数值分析》-- 拉格朗日插值_第21张图片
    《数值分析》-- 拉格朗日插值_第22张图片
2.

在这里插入图片描述
《数值分析》-- 拉格朗日插值_第23张图片
《数值分析》-- 拉格朗日插值_第24张图片

  1. 《数值分析》-- 拉格朗日插值_第25张图片
    《数值分析》-- 拉格朗日插值_第26张图片
    《数值分析》-- 拉格朗日插值_第27张图片
    《数值分析》-- 拉格朗日插值_第28张图片
    《数值分析》-- 拉格朗日插值_第29张图片
    《数值分析》-- 拉格朗日插值_第30张图片
    《数值分析》-- 拉格朗日插值_第31张图片

总结

《数值分析》-- 拉格朗日插值_第32张图片
《数值分析》-- 拉格朗日插值_第33张图片
《数值分析》-- 拉格朗日插值_第34张图片

  • 拉格朗日插值的优缺点
    《数值分析》-- 拉格朗日插值_第35张图片

你可能感兴趣的:(数值分析,数值分析,插值方法,拉格朗日插值)