数学分析 三次样条插值..pptVIP

三次样条插值 分段低次插值 分段线性插值 问题：已知 及 ，i=0,1,...,n 为[x0, xn]上的分段线性函数且满足 计算公式 余项公式 收敛性 分段3次Hermite插值 已知 及 ， ， 为[x0, xn]上的分段3次多项式且满足 分段3次Hermite插值 三次样条函数 三次样条函数的确定 边界条件 边界条件 三次样条插值求解公式 利用分段3次Hermite插值 设 ，利用分段3次Hermite插值公式(含未知mi)，使得插值条件和连接条件的连续性和一阶光滑性满足； 由n-1个二阶光滑性约束条件和边界条件来求个待定参数 ，转化为解三对角方程组(一阶导数边界条件 ) 三次样条插值求解步骤 由计算机数值求解的步骤 由边界条件和插值条件列三对角方程组; 用追赶法解方程组求得 ; 判断插值点 x 在第i 个小区间; 用第i 个小区间的三次样条插值多项式求插值. 例4.6 求满足下列数据的三次样条插值函数 S(x) 解法一（分析推导, 称为承袭法）; 解法二（待定系数法）; 解法三（用计算公式）. * * 问题 抛物线插值的误差比线性插值要小，是不是插值多项式的次数越高，精度就越好？ NO! 例：在 [?5, 5] 上考察 的Ln(x)。取 ? Ln(x) f (x) - 5 - 4 - 3 - 2 - 1 0 1 2 3 4 5 - 0.5 0 0.5 1 1.5 2 2.5 n 越大，端点附近的抖动越大，称为 Runge现象。 分段低次插值 基本思想：用分段低次多项式来代替单个多项式。 具体作法：(1) 把整个插值区间分割成多个小区间； (2) 在每个小区间上作低次插值多项式； 优点：公式简单、 运算量小、稳定性好、收敛性 … (3) 将所有插值多项式拼接整一个多项式。 缺点：节点处的导数可能不连续，失去原函数的光滑性。 在处理实际问题时，总是希望将所得到的数据点用得越多越好。最简单的方法是用直线段将函数值点直接连接。 但不光滑！ 计算公式 收敛性 余项公式 一阶光滑！而且实际工程中一般不知道 的值！ 样条函数 由一些按照某种光滑条件分段拼接起来的多项式组成的函数。 最常用的样条函数为三次样条函数，即由三次多项式组成，具有二阶连续导数。 设节点 a = x0 x1 … xn-1 xn = b ，若函数s(x)在[a,b]上有二阶连续导数，在每个小区间[ xi , xi+1 ]上是三次多项式，则称其为三次样条函数。如果同时满足 s(xi) = f (xi) (i = 0, 1, 2, …, n)，则称 s (x) 为 f (x) 在 [ a , b ]上的三次样条函数。 定义 三次样条插值 问题：已知 及 ，S(x)为[xi-1, xi]上的不超过3次多项式，且满足 (1) 插值条件： (2) 连接条件: 节点： x0 x1 … xn-1 xn 函数值：yi = f (xi) (i = 0, 1, 2, …, n) 由定义可设： s(x) 满足： (i = 0, 1, 2, …, n) 其中 为[ xk-1 , xk ]上的三次多项式，且满足 (k = 1, 2, …, n) 每个 sk(x) 均为三次多项式，有4个待定系数，所以共有 4n 个待定系数，需 4n 个方程才能确定。前面已经得到 2n +2 (n -1) = 4n –2 个方程，还缺 2 个方程！ (k = 1, 2, …, n-1) 实际问题通常对样条函数在端点处的状态有要求，即所谓的边界条件。 第一类边界条件：给定函数在端点处的一阶导数，即 第二类边界条件：给定函数在端点处的二阶导数，即 当 时，称为自然边界条件，此时的样条函数称为自然样条函数。 第三类边界条件：设 f (x) 是周期函数，并设 xn – x0 是一个周