42~数值分析.pptVIP

重点插商 Newton插值计算 插商表1 插商表2 求Nn x 插商表1计算简单，好实现，但数值不稳定。 插商表2在计算机上稳定性好，但算法复杂。 计算Nn x 常采用秦九韶程序（取n 4） 例题 在实际应用中 ，常是等距节点情况，即 这里h 0为常数，称为步长，这时Newton插值公式就可以简化，为此我们引入差分概念。 等距节点Newton插值公式 插商与差分的关系 （1）用前插表示N x 在等距节点条件下有： （2）用后插表示N x 例题 Lagrange插值公式所求得L x 保证了节点处的函数值相等，也就是保证了函数的连续性，但不少实际问题还需要插值得光滑度，也就是还要求它在节点处的导数值也相等，导数的阶数越高则光滑度越高。现代的仿生学就是一个典型的例子。在设计交通具的外形，就是参照海豚的标本上已知点及已知点的导数，做插值在计算机上模拟海豚的外形制成飞机、汽车等外形。 Hermite插值多项式 构造H x 算法实现 Hermite插值余项 特例（n 1 例题 * * F n ……… …… … … … F 3 F 2 F 1 F 0 单元号 三阶插商 二阶插商 一阶插商 f x sin x 关于节点的差商表 一阶差商 二阶差商 三阶差商 四阶差商 0 0 0.954 929 658 0.791 089 691 -0.208 607 6 0.607024424 -0.351 538 65 -0.136 489 09 1 0.255 872 63 -0.447 100 35 -0.091 254 70 0.028 797 106 （Newton插值法） 1 输入： 2 （i 0，1，2，…，n） 3 计算差商 对i 1，2，…，n做 1 对j i，i+1，…，n做; 2 对j i，i+1，…，n做; x f x 0 1 -0.181 269 0.2 0.818 731 0.032 585 -0.148 411 -0.005 955 0.4 0.670 320 0.026 903 0.001 077 -0.121 508 -0.004 878 -0.000 191 0.6 0.548 812 0.022 025 0.000886 -0.099 483 -0.003 992 0.8 0.449 329 0.018 033 -0.018 033 1.0 0.367 879 所谓Hermite 插值是插值函数与被插值函数不仅在节点处函数值相等,而且导数值也相等.设函数f x 在区间[a,b]上具有一阶连续导数,如果有不超过2n+1次多项式满足 则称为f x 关于节点的Hermite插值多项式。 其中,为基本插值多项式 。 由条件可得 将A、B代入上式可得式 ,而式 完全类似可证。 设函数y f x 在区间[a,b]上的n+1个互异节点 处给定了函数值及一阶导数（i 0,1,2…n）,则一定存在唯一的不超过2n+1次的多项式满足 （i 0,1,2…n） 且有 其中为基本插值多项式。 所以有 与真值相比已有三位有小数字。 由，可有，于是可设 由 和得 将A、B代入式得 由此可解得 Newton插值公式 由差商定义 1 一阶差商 二阶差商 三阶差商 n阶差商 单元号 F 0 F 1 F 2 F 3 F n 试用Newton插值公式计算sinx在x π/12处的近似值。 解 先列差商，所以得 等距节点Newton插值公式 设函数f x 在等距节点 i 0,1,2,…，n 上值为 ，则 1 称 i 0,1,2,…，n 为函数f x 在点上的一阶向前差分（简称差分）；又称 k 1,2,…,n;i 0,1,…,n-k 为函数f x 在点上的k阶向前差分，这里约定; 2 称 i n,n-1,…,1 为函数f x 在点上的后差分；又称 k 1,2,…,n; i n-k+1,…,2,1 为函数f x 在点上的k阶向后差分，同样约定。 设函数y f x 在各节点的取值如下表所示，试计算各阶差分值。 x 0 0.2 0.4 0.6 0.8 1.0 f x 1 0.818 731 0.670 320 0.54