第2部分+插值方法.ppt

明德新民 止于至善 2.4 分段低次插值公式 明德新民 止于至善 (2) 将区间[-5，5]分割成十等份，得到节点和节点处的函数值为 节点处的插值基函数记为 那么10次插值多项式函数为 其图形如下： 明德新民 止于至善 2.4 分段低次插值公式 明德新民 止于至善 (3) 将区间[-5，5]分割成十六等份，得到节点和节点处的函数值为 节点处的插值基函数记为 那么16次插值多项式函数为 其图形如下： 明德新民 止于至善 从图中，可看见，在0附近插值效果是好的，即余项较小，另一种现象是插值多项式随节点增多而振动更多。 这种插值多项式当节点增加时反而不能更好地接近被插值函数的现象，称为龙格现象 2.4 分段低次插值公式 （1）I(x)在每个小区间[xi，xi+1]上是个次数不超过一 的多项式; y0 y1 y2 … yn y x0 x1 x2 … xn x 求一个多项式I(x), 使其满足如下条件： （3） I(xi)=yi =f(xi) (i=0,1, ..., n). 已知数据表格： 设在[a, b]上取n+1 个节点，且 a= x0 x1x2…xn-1xn=b, f(x)的函数值为yi =f(xi) (i=0,1,2,…,n), 即 （2） I(x）∈C[a, b]; 明德新民 止于至善 称之为f(x)在区间[a, b] 上关于数据(xi , yi ) (i=0,1,2,…,n)的分段线性插值函数. 说明： I(x）的特点（优点和缺点）？ 失去了原函数的光滑性。 (1) 插值公式 明德新民 止于至善 在每个子区间[xi，xi+1]， 或 (2) 插值余项 明德新民 止于至善 因此,在插值区间[a,b]上有 (2) 插值余项 将区间[-5，5]分成10等份，做分段线性插值函数，并 做出图形观察逼近程度。 课下练习： 明德新民 止于至善 2.4 分段低次插值公式 明德新民 止于至善 解 将区间[-5，5]分割成十等份，得到节点和节点处的函数值为 节点处的分段线性插值基函数记为 那么分段线性插值多项式函数为 其图形如下： 明德新民 止于至善 1、两点带导数的三次埃尔米特插值多项式 求多项式 H3(x), 使其满足如下条件： H3(x)称为两点带导数的三次埃尔米特插值多项式。 法1822 -1901 2.5 埃尔米特(Hermite）插值 给定如下数据表： y0 y1 y m0 m1 y′ x0 x1 x H3(x)是一个次数不超过3 的多项式; （2） 明德新民 止于至善 x1 x0 x1 x0 1 0 0 0 μ1(x) 0 1 0 0 μ0(x) 0 0 1 0 λ1(x) 0 0 0 1 λ0(x) 导数值 函数值 节 点 基 函 数 1、两点带导数的三次埃尔米特插值多项式 (1) 埃尔米特插值公式 设节点xi 处的插值基函数分别是λi(x) 和μi (x)(i=0,1) 令 由 得 由 得 于是 2.5 埃尔米特(Hermite）插值 明德新民 止于至善 同理 于是 令 由 得 同理 (1) 埃尔米特插值公式 2.5 埃尔米特(Hermite）插值 x1 x0 x1 x0 1 0 0 0 μ1(x) 0 1 0 0 μ0(x) 0 0 1 0 λ1(x) 0 0 0 1 λ0(x) 导数值 函数值 节 点 基 函 数 ? 设f(x)在包含x0、x1的区间[a,b]内存在4阶导数，则当x∈[a,b]时, 有 设 则当x∈(x0 , x1)时, 余项有如下估计式 且与x有关) (2) 埃尔米特插值余项 1、两点带导数的三次Hermite插值多项式 2.5 埃尔米特(Hermite）插值 明德新民 止于至善 ? 例 已知f(x)=x1/2及其一阶导数的数据见下表,用埃尔米特插值公式计算1251/2的近似值,并估计其截断误差. 1/24 1/22 f (x) 12 11 f(x) 144 121 x 解 得 由 可得 2、n+1个节点带导数的埃尔米特插值多项式 求一个多项式 H2n+1(x), 使其满足如下条件： H2n+1(x)称为n+1