第二章数值分析ppt.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 三、误差界与收敛性 作业 P49, 20(1)只列出方程组; 20(2)应得结果. * * * * * * * * * * * * * * * * * * * * * * * * * 上节讨论任意分布节点的插值公式，应用时常碰到等距节点的情形，此时插值公式可简化，为此先介绍差分． 差分及其性质 四、差分形式的牛顿插值多项式 差分的基本性质: 差分表: ?2f0 ?2f1 … ┆ ?2f2 ┆ ┆ ?f0 ?f1 ?f2 ?f3 ┆ f0 f1 f2 f3 f4 ┆ 0 1 2 3 4 ┆ ?2 ?3 … ? fk k 等距节点插值公式 作业 P48, 1，6，8. §4 埃尔米特插值 一、重节点差商与泰勒插值 二、两个典型的埃尔米特插值 §5 分段低次插值 一、高次插值的病态性质 龙格(Runge)现象…. h=1,x=-5:h:5, y=1./(1+x.^2), p=polyfit(x,y,length(x)), xi=linspace(-5,5, 100/h), yi=polyval(p,xi), plot(x,y, o,xi,yi, -) 二、分段线性插值 所谓分段线性插值就是用通过插值点的折线段逼近f(x). 三、分段三次埃尔米特插值 分段线性插值函数导数间断，若已知节点上函数值和导数，可构造一个导数连续的插值函数Ih(x)，满足 作业 P49, 16. §6 样条插值 问题背景… 一、样条插值的概念 二、三次样条插值函数的建立 §1 引 言 第2章 插值法 一、插值问题的提出 应用：例如程控加工机械零件等。 本章求出插值多项式, 分段插值函数, 样条插值函数；讨论P(x)的存在唯一性、收敛性及误差估计. 几何上、发展和实践上. 二、多项式插值 §2 拉格朗日插值 一、线性插值和抛物插值 对给定插值点,求出形如 的插值多项式的方法有多种. 二、拉格朗日插值多项式 练习 给定数据表 xi 0 1 2 3 yi 0 1 5 14 求三次拉格朗日插值多项式L3(x). 三、插值余项与误差估计 例1 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用线性插值计算和抛物插值计算sin0.3367的值, 并估计误差. 例1 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并估计误差. 作业 P48, 2. §3 差商与牛顿插值多项式 一、插值多项式的逐次形成 拉格朗日插值优缺点…. 二、差商及其性质 差商的基本性质 由(3.4)得差商表: k xk f(xk) 一阶差商 二阶差商 三阶差商 … 0 1 2 3 4 ┆ x0 x1 x2 x3 x4 ┆ f(x0) f(x1) f(x2) f(x3) f(x4) ┆ f[x0, x1] f[x1, x2] f[x0,x1,x2] f[x2, x3] f[x1,x2,x3] f[x0,x1,x2,x3] f[x3, x4] f[x2, x3,x4] f[x1,x2,x3,x4] … ┆ ┆ ┆ 三、牛顿插值多项式 k xk f(xk) 一阶差商 二阶差商 三阶差商 四阶差商 0 1 2 3 4 1 2 3 4 5 1 4 7 8 6 3 3 0 1 -1 -1/3 -2 -3/2 -1/6 1/24