李庆扬数值分析——插值法.pptVIP

描述事物之间的数量关系：函数。 有两种情况： 一是表格形式——一组离散的数据来表示函数关系；另一种是函数虽然有明显的表达式，但很复杂，不便于研究和使用。 从实际需要出发：对于计算结果允许有一定的误差，可以把函数关系用一个简单的便于计算和处理的近似表达式来代替，从而使问题得到简化。;当精确函数 y = f(x) 非常复杂或未知时，在一系列节点 x0 … xn 处测得函数值 y0 = f(x0), … yn = f(xn)，由此构造一个简单易算的近似函数 g(x) ? f(x)，满足条件g(xi) = f(xi) (i = 0, … n)。这里的 g(x) 称为f(x) 的插值函数。; 根据实际需要，可以用各种不同的函数来近似原来的函数。;§1 拉格朗日多项式 ;1.2 二次插值; 对于线性插值的两种形式解进行适当的分析, 从中寻求规律而得到启发,就有了所谓的拉格朗日插值法(公式)和牛顿插值(公式).; 首先, 线性插值的两点式可看作是两个特殊的一次式的一种线性组合.; 由此启发，我们希望二次插值也能由一些二次插值基函数来线性组合:;同理可得 ;n ? 1;定理;设节点;注：? 通常不能确定 ?x , 而是估计 , ?x?(a,b) 将 作为误差估计上限。;例1 求经过A(0,1)，B(1,2)，C(2,3)三个插值点的插值多项式.;例2：已知;n = 2; 例3 考虑下述的插值法问题：求二次多项式P(x)，满足P(x0) = y0, 其中 是已给的数据并给出使这一问题的解存在且唯一的条件.; 回顾 拉格朗日插值公式 ;n ? 1;Lab 1. Lagrange Polynomial Given a set of sample points of a certain function f , use Lagrange polynomial to approximate function values at some given points . Input There are several sets of inputs. For each set: The 1st line contains an integer 20 ? n ? 0 which is the degree of Lagrange polynomial. n = ?1 signals the end of file. The 2nd line contains n+1 distinct real numbers . The 3rd line contains n+1 real numbers . The last line of a test case consists of an integer m 0 and m real numbers . The numbers are separated by spaces and new lines. ;Output (? represents a space) For each ai , you are supposed to print the function value at this point in the following format: fprintf(outfile, f(%6.3f)?=?%12.8e\n, a, f ); The outputs of two test cases must be seperated by a blank line. ; Lagrange插值公式(利用插值基函数很容易得到): 含义直观,结构紧凑,在理论分析中非常方便; 计算机上实现也很容易. ;差 商;);注: 1. 事实上,从上述可看出二次牛顿插值公式是用待定系数法求得的;