多项式插值与样条插值ppt.pptx

现代数值计算第3章 多项式插值与样条插值第3章 多项式插值与样条插值 实践中常有这样的问题： (1) 由实验得到某一函数f(x)在一系列点x0, x1, …, xn处的值y0, y1, …, yn, 其函数的解析表达式是未知的 (2) 或者f(x)虽有解析式, 但计算复杂, 不便于使用 需要构造一个简单函数p(x)近似地代替f(x) —— 这就是函数逼近问题§3.0 基本概念1. 插值与拟合 设已知被逼近函数f (x)在离散点xi ? [a, b]上的值f (xi) = yi, i = 0, 1, 2, …, n (1) 要求p(x)满足（甚至 ）的问题称为函数插值。 (2) 要求p(x)满足 为最小的问题称为数据拟合（曲线拟合）§3.0 基本概念2. 简单函数类 设?0, ?1, …, ?n线性无关, 令? = span{?0, ?1, …, ?n}为简单函数类, 其中?0, ?1, …, ?n称为?的基函数。 逼近问题即用p(x) = a0?0(x) + a1?1(x) +…+ an?n(x)来做逼近, 问题归结为求其中的待定系数a0, a1, …, an。§3.1 多项式插值3.1.1 多项式插值问题的定义 【定义1】设区间[a, b]上定义的函数f(x)在n+1个互异点a ? x0 ? x1 ? ... ? xn ? b的函数值为y0, y1, ..., yn, 若存在一个次数不超过n的多项式pn(x)满足条件：pn(xi) = yi i = 0, 1, 2, …, n(3.1)则称pn(x)为函数f(x)的n次插值多项式 按条件(3.1)求函数f(x)近似表达式pn(x)的方法称为插值法, 条件(3.1)称为插值条件 点xi (i = 0, 1, 2, …, n)称为插值节点 区间[a, b]称为插值区间§3.1 多项式插值3.1.1 多项式插值问题的定义 【定义1】设区间[a, b]上定义的函数f(x)在n+1个互异点a ? x0 ? x1 ? ... ? xn ? b的函数值为y0, y1, ..., yn, 若存在一个次数不超过n的多项式pn(x)满足条件：pn(xi) = yi i = 0, 1, 2, …, n(3.1)则称pn(x)为函数f(x)的n次插值多项式 几何意义：曲线y = pn(x)与曲线y = f(x)在平面上相交于n+1个互异的点(xi, yi) (i = 0, 1, 2, …, n).§3.1 多项式插值3.1.2 插值多项式的存在唯一性 【定理1】存在唯一pn(x) ? Pn[x]满足插值条件(3.1)证明：设pn(x) = a0 + a1x +…+ anxn 满足插值条件pn(xi) = yi i = 0, 1, 2, …, n 从而只需求解线性方程组pn(xi) = a0 + a1xi +…+ anxin = yi (i = 0, 1, 2, …, n) (3.2)(3.2)的系数行列式为范德蒙行列式：§3.1 多项式插值3.1.2 插值多项式的存在唯一性 【定理1】存在唯一pn(x) ? Pn[x]满足插值条件(3.1)证明：只需求解线性方程组pn(xi) = a0 + a1xi +…+ anxin = yi (i = 0, 1, 2, …, n) (3.2)(3.2)的系数行列式为范德蒙行列式：因为x0, x1, …, xn互异, 所以Vn ≠ 0。即(3.2)存在唯一解, 从而存在唯一的pn(x) ? Pn[x] 满足插值条件(3.1)。§3.1 多项式插值3.1.2 插值多项式的存在唯一性 【例1 】给定数据求次数不大于3的插值多项式p3(x)解：设p3(x) = a0 + a1x + a2x2 + a3x3依题意有解之得：a0 = 10, a1 = 5, a2 = – 5, a3 = 2即有 p3(x) = 10 + 5x – 5x2 +2x3xi-1125yi-77-435§3.1 多项式插值3.1.2 插值多项式的存在唯一性 【例1 】给定数据求次数不大于3的插值多项式p3(x)注：(1) 范德蒙矩阵的条件数很大 —— 误差大计算量大(2) 选择适当基函数使插值多项式具特殊形式xi-1125yi-77-435§3.1 多项式插值3.1.3 插值基函数 所求的插值多项式pn(x)是线性空间Pn[x]中的一个点, 而Pn[x]的基是不唯一的。因此, n次代数插值多项式pn(x)可以写成多种形式. 由线性空间的不同基底出发, 构造满足插值条件的多项式的方法称为基函数法。§3.1 多项式插值3.1.