西工大计算方法4重点分析.ppt

要确定4n个系数，还需附加2个约束条件. 常用的约束条件有以下三类： 此时一般有 成立. ③.周期性边界条件 , ②.弯矩边界条件 特别的称 为自然边界条件. ①.转角边界条件 2. 三弯矩构造法 记 , 基本步骤如下： ①.取 为待定参数，并用S(x)的插值条件写出 的表达式。 ④.代入S(x)的表达式，得各个区间上的表达式。 ②.用 在内节点 的连续条件及边界条件导出关于 的方程组。 ③.求解后得到 。 式中 。 对 积分两次, 并利用插值条件 , 确定两个积分常数, 得到 计算 类似可以得到 两边同乘以 , 得 令 , 有 式中 , . 若附加弯矩约束条件, 得 系数矩阵严格对角占优, 故系数矩阵非奇异, 上述线性方程组有唯一解，可用追赶法求解。将解带回到子区间上的表达式中（用二阶导表示），即有s(x)在每个区间上的表达式。 若附加转角边界条件, 得 线性方程组为 对于周期性边界条件, 得： 即 式中 , . 线性方程组为: 将 当作已知参数, 从后 个方程中求解出用 表示的后 个参数, 然后将它们代入第一个方程解得 , 最终得到其它参数. Remark： 1.类似地，可以使用节点处的一阶导数来表示三次样条插值函数。 2.对三次样条插值函数来说，当插值节点逐渐加密时，可以证明，不但样条插值函数收敛于函数本身，而且其导数也收敛于函数的导数。 差商型误差估计 导数和差商的关系 差商型误差估计特点：对被插值函数光滑性要求不高；但不适用于实际计算。 四、例题 解 1）建立差商表 1.0 1.5 2.0 0.8415 0.9975 0.9093 0.312 -0.1764 -0.4884 2）插值 Newton插值多项式适用于节点任意分布的情形。但当节点等距分布时，可以简化Newton插值公式。 §4.4 等距节点插值 设a=x0x1… xn=b，yi=f(xi)为等距节点xi=x0+h(i=0,1,…,n)上的函数值，其中h=(b-a)/n称为步长。 在此基础上我们先定义差分，用差分表示Newton插值多项式，从而得到等距节点的插值公式。 一、差分的定义与性质 定义：称 ?yi=yi+1-yi (i=0,1,…,n-1) 为f(x)在xi处以h为步长的一阶向前差分。 ?2yi＝?yi＋1-?yi =yi+2-2yi+1+yi (i=0,1,…,n-2) 称为f(x)在xi处以h为步长的二阶向前差分。 一般地，?myi＝?m-1yi＋1-?m-1yi (i=0,1,…,n-m) 称为f(x)在xi处以h为步长的m阶向前差分。 差分的性质 性质1：各阶差分可用函数值线性表示，其计算公式为： 其中 性质2：差分与差商满足下述关系： 证明：利用数学归纳法 当k＝1时，有 即结论成立。 设k＝m-1时结论成立，即 则当k＝m时，有 由数学归纳法知，结论成立。 证毕 Remark：类似地可以定义向后差分与中心差分： 性质3：差分与导数满足关系： 证明：利用差商与导数、差分的关系，有： 证毕 二、Newton向前插值公式 令x=x0+th，由xi=x0+ih(i=0,1,…,n)得： x-xi=(t-i)h，则有: 将差商与差分的关系式 带入Newton插值多项式，得: 从而可得Newton向前插值多项式及其余项为： 三、差分表 Newton向前插值公式，又称表初公式，它利用差分表的最上面一个斜行的数值进行计算。 四、例题 解 # 五、Newton向后插值公式 类似于向前差分，也可以得到差商与向后差分的关系： 将插值节点从大到小排列，即 类似于向前插值公式，可得到Newton向后插值公式，又称表末公式，它利用差分表的