数值计算课件第二章.ppt

浙江大学研究生学位课程 《实用数值计算方法》 化学工程学系 第二章 插值和逼近 2.1 多项式插值 2.2 密切多项式插值 2.3 分段插值和样条函数 2.4 有理函数 2.5 多维插值 2.6 函数的逼近 插值多项式 Lagrange 多项式 Aitken-Neville 逐步加精格式 Newton 多项式 Newton-Gregory 前差分插值 Gregory 后差分插值 Gauss 差分插值 Stirling 差分插值 2.2 密切多项式插值 2.3 分段插值和样条函数 2.4 有理函数插值  Rational Function Interpolation 2.5 多维插值 一般化的方法示例 x取3次，y取2次，共12个待定系数 需要12个节点上的型值 解线性代数方程组可以求得 第二章 习题 2-1. Show that the divided difference is a symmetric function. That is, f [x1, x2,……,xk]=f [y1,y2,……,yk] if the set {yi} is a permutation of the set xi . 2-2. Write a FORTRAN program to compute and print the Table of 1st- to 4th-order divided differences. Use following sets of data to test your program. 2.4.1 有理函数的一般情况 有理函数 度型 degree type (n,m) 度数 number of degree d=n+m 标准形式 系数个数 d+1 型值点数 d+1 2.4.1 使有理函数通过d+1个型值点 得到线性方程组 并非所有情况下（2-84）一定有解 例：给出三个型 值点，试确 定有理函数 系数。 2.4.1 解 构造有理函数 满足 得到线性方程组 这是因为数据给出的型值点中存在着 不可达点 inaccessible point 2.4.2 逆差商递推算法 Recursive Procedures Using Inverse Differences 设计一种递推算法，使度型（n,m）中的n和m 的数值交替地逐个增加如下： 2.4.2 设 通过 而 通过 诸节点上之型值点。则 2.4.2 2.4.2 2.4.2 以上公式可由下列计算步骤完成 表2-9 型值和各阶逆差商 2.4.3 反差商算法 Reciprocal Difference Algorithm 为了避免逆差商不对称性的麻烦 定义另一种差商 它具有对称性，称之为反差商 由以上定义，有 又根据逆差商的定义 2.4.3 2.4.3 2.4.3 具体计算可按以下步骤 2.5.1 多维线性插值 以二维为例，最简单的考虑是 由 求 的值 可认为是x,y两个方向线性插值的叠加 2.5.1