第5章_插值法精品.ppt

第五章 插值法 实践中有些函数解析式未知，或虽有明确解析式，但计算复杂，这时需要用比较简单且易于计算的函数p(x)去近似代替它，使得 p(xi)= yi (i=0,1,2,…,n) 这类问题称为插值问题。函数p(xi)称为插值函数。 x0,x1,… xn称为插值节点或简称节点。插值节点所在的区间称为插值区间。p(xi)= yi 称为插值条件。 多项式的插值问题 构造n次多项式 Pn(x)= a0 + a1x + a2x2+…+ anxn 使满足Pn(xi)= yi (i=0,1,2,…,n)，及利用多项式Pn(x)进行插值计算的问题。 §1.1 差 商 差商的定义 f(x)在xi点的零阶差商为 f[xi]= f(xi) (i=0,1,2,…,n) f(x)在[xi,xj]上的一阶差商为 §1.1 差商 区间[xi,xi+1,…,xi+n]上的n阶差商为： §1.1 差商 差商表 §1.1 差 商 例5.1 试列出f(xi)=x3在节点x=0,2,3,5,6上的各阶差商值。 §1.1 差 商 如以x代表时间t,f(x)代表路程s,则一阶差商为?si/ ?ti=Vi，它相当于在[ti, ti+1]范围内的一种平均速度，二阶差商则为上述平均速度的平均变化率，即平均加速度,…,所以差商表的数值可以直接反映出函数值的变化情况。 差商的重要特性——对称性，即差商的值与同组节点排列的次序无关。 §1.1 差 商 §1.1 差 商 §1.2 牛顿基本差商公式 设x为插值区间内的一个节点，按照差商定义，有如下关系式 §1.2 牛顿基本差商公式 由上式逐次解出f(x), f[x0 ,x], f[x1 ,x0 , x], f[x2 , x1 ,x0 , x],…,并代入f(x)得： §1.2 牛顿基本差商公式 其中Pn(x)称为牛顿基本差商公式， Rn(x)称为牛顿基本差商公式的余式。 若用Pn(x)近似f(x),则误差为Rn(x)。 当x=xi(i=0,1, 2,…,n)时， Rn(xi )=0, Pn(xi)= yi 当x?xi(i=0,1, 2,…,n)时， Rn(xi ) ? 0, Pn(x)? f(x) §1.2 牛顿基本差商公式 §1.2 牛顿基本差商公式 §1.3 牛顿基本差商公式的余式估计 f(x)=Pn(x)+ Rn(x) §1.3 牛顿基本差商公式的余式估计 §1.3 牛顿基本差商公式的余式估计 §1.3 牛顿基本差商公式的余式估计 §1.3 牛顿基本差商公式的余式估计 Pn (x): §1.3 牛顿基本差商公式的余式估计 §1.3 牛顿基本差商公式的余式估计 例5.9 用插值法求 的值。 §1.3 牛顿基本差商公式的余式估计 §1.3 牛顿基本差商公式的余式估计 采用事后估计误差方法： §1.3 牛顿基本差商公式的余式估计 事后估计误差公式： §2.1 差分 函数在等距节点上的值为y0 y1… yn ，称 ?yi-1= yi - yi-1 为函数f(x) 在[xi-1, xi]上的一阶差分。称 ?2yi-1= ? yi - ? yi-1 为函数f(x) 在[xi-1, xi+1]上的二阶差分。称 ?kyi-1= ?k-1yi - ?k-1yi-1 为函数f(x) 在[xi-1, xi+k-1]上的k阶差分。 §2.1 差分 §2.1 差分 ?y0= y1 – y0 ?y1= y2 – y1 ?y2= y3 – y2 §2.1 差分 §2.1 差分 等距节点情况下xi= x0+ih ，用差分表示差商： §2.1 差分 §2.2 牛顿前向插值公式 建立等距节点的牛顿基本差商公式: §2.2 牛顿前向插值公式 §2.2 牛顿前向插值公式 §2.2 牛顿前向插值公式 §2.2 牛顿前向插值公式 §2.3 牛顿后向插值公式 §2.3 牛顿后向插值公式 在等距节点情况下，以xn xn-1… x0顺序建立牛顿基本差商公式 §2.3 牛顿后向插值公式 §2.3 牛顿后向插值公式 §2.3 牛顿后向插值公式 §2.3 牛顿后向插值公式 §2.3 牛顿后向插值公式 §2.3 牛顿后向插值公式 §2.4 斯梯林插值公式（1） §2.4 斯梯林插值公式（2） §2.4 斯梯林插值公式（3） 令P2(x)= §2.4 斯梯林插值公式（3） 令P2(x) §2.4 斯梯林插值公式 §2.4 斯梯林插值公式 §2.4 斯梯林插值公式 §2.4 斯梯林插值公式 §2.4 斯梯林插值公式 §2