函数的数值逼近.pptxVIP

会计学1函数的数值逼近 2§1 引言一、函数的工程化表达对于很多实际工程计算问题，函数是通过实验或观测得到的，表达形式上为函数表，无解析表达形式。2. 虽然有些函数存在解析的表达式，但形式过于复杂而不易使用。第1页/共121页 3二、问题的提出设 是R中若干个不同的点，每个点 对应一个数值它们可以是实测得到的，也可以是一个已知函数的值 。如何近似由这组数据 确定的函数？并由此可提出两类问题：作一条曲线，其类型是事先给定的（如：代数多项式），使该曲线经过给定点 。这就是所谓的插值问题。作一条指定的曲线，使该曲线能在“一定意义”下逼近这一组数据。这就是所谓的曲线拟合问题。第2页/共121页 4(1)复杂函数的计算;(2)函数表中非表格点计算(3)光滑曲线的绘制;(4)提高照片分辩率算法(5)定积分的离散化处理;(6)微分方程的离散化处理;(7)积分方程的离散化处理;插值方法的应用:第3页/共121页 5三、插值的定义与存在性求 P( x ) 的方法就是插值法。若存在一简单函数P(x)，使得 P( x )为 f (x) 的插值函数点 x0 , x1, ?, xn 为插值节点 (1)式为插值条件 f ( x ) 为被插函数 [a , b] 为插值区间设 f(x)∈C [a , b], 取点 a ≤x0＜x1＜···＜xn≤b成立，则称定义：第4页/共121页 6若P ( x ) 是次数不超过n 的实系数代数多项式, 即则称P( x)为n 次插值多项式. 相应的插值法称为多项式插值法(代数插值法)。P ( x ) = a0 + a1 x + ? + an x n从几何上看x 0 y y = P(x) a=x0 x1 x2 x3 xn=b ?(xi, yi)y = f (x)曲线 P ( x) 近似 f ( x) 第5页/共121页 7研究问题：（1）满足插值条件的P ( x) 是否存在唯一？（2）若满足插值条件的P ( x) 存在，如何构造P ( x)？（3）如何估计用P ( x)近似替代 f ( x) 产生的误差？第6页/共121页 82、插值多项式的存在唯一性证明：由(1)式(2)定理 若插值结点 x0,x1,…, xn是(n+1)个互异点,则满足插值条件 P(xk)= yk (k = 0,1,…,n) 的n次插值多项式 P(x)=a0 + a1x +……+ anxn 存在且唯一。第7页/共121页 9点是互异的其系数行列式 ：为范德蒙行列式。只要插值节点互不相同，则系数矩阵非奇异。故方程组解存在且唯一。第8页/共121页 10插值多项式的唯一性表明，对同一组节点，它们的 插值多项式是唯一的，可能由不同的方法，会得到不同形式的插值多项式，但它们之间可以相互转化，本质相同，当然误差也一样。 n +1组节点只能确定一个不超过n次的多项式，若n次，如设为?n+1(x)，则有n+2有待定参数a0,a1,…,an, an+1需确定，而n +1个组节点，只构成n +1个插值条 件，即构成n+1个方程，只能确定n+1个变量的方程组。上述证明是构造性的（给出解决问题的方法）即 以通过解线性方程组来确定插值多项式，但这种方法的计算量偏大，计算步骤较多，容易使舍入误差增大。因此实际计算中需要用其它方式进行。说明：第9页/共121页 11x 0 y y = f (x) 的几何意义一、线性插值与抛物线插值1. 线性插值(n =1) 设已知区间[ xk , xk+1]端点处的函数值yk= f (xk)，yk+1 = f (xk+1)，y = L1(x)xk xk+1 求线性插值多项式L 1(x ) ，使其满足 ——过两点 (xk , y