算法与例子.pptVIP

* 差商(均差)的概念 算法与例子 牛顿插值公式 取x0, x1, x2,求二次函数 P(x)=a0 + a1(x – x0) + a2 (x – x0)(x – x1) 满足条件 P(x0)=f(x0), P(x1)=f(x1), P(x2)=f(x2) 插值条件引出关于a0, a1, a2方程 解下三角方程组过程中引入符号 a0 = f(x0), a1 = f[x1, x2], a2 = f[x0, x1, x2] P(x)=a0 + a1(x – x0) + a2 (x – x0)(x – x1) 定义5.3 若已知函数 f(x) 在点 x0, x1, ···, xn 处的 值 f(x0), f(x1), ···, f(xn).如果 i ≠ j ,则 ( j = 0,1,…,n-1 ) 一阶差商 二阶差商 ( j = 0,1,…,n-2 ) n阶差商 差商的对称性 由此可以看出，差商的值与节点x0, x1, ···, xn 的排列次序无关，即 x - 2 -1 0 1 3 y -56 -16 -2 -2 4 例 由函数表求各阶差商 x f(x) 一阶差商 二阶差商 三阶差商 -2 -56 -1 -16 40 0 -2 14 -13 1 -2 0 -7 2 3 4 3 1 2 解 按公式计算一阶差商、二阶差商、三阶差商如下 MATLAB程序计算 x=[-2 -1 0 1 3]; y=[-56 -16 -2 -2 4]; f=y -56 -16 -2 -2 4 -56 40 14 0 3 -56 40 -13 -7 1 -56 40 -13 2 2 -56 40 -13 2 0 n=length(x); for k=2:n for j=n:-1:k f(j)=(f(j)-f(j-1))/(x(j)-x(j+1-k)); end f end 牛顿插值公式 f(x)=f(x0)+(x-x0)f[x0, x] n次牛顿基本插值公式 对于插值节点，有 根据代数插值存在唯一性定理，n次牛顿插值公式恒等于n次拉格朗日插值公式，插值余项 根据插值多项式的定义， Nn(x)是f(x)的一个n次插值多项式 *