第七节 方程求根.pptVIP

* 第七节 方程求根 求初始近似根的“逐步扫描法”。 一、迭代法：xk+1=g(xk) 二、牛顿法： 三、弦截法 牛顿法的突出优点：收敛速度快。 缺点：需要计算导数 。若计算 较困难，则用牛顿法不方便。 为避免计算 ，可用差商 近似代替 ，于是得 弦截法的 迭代格式： 上式有明显的几何意义：设曲线y=f(x)上横坐标为x0和xk的点分别为P0 和Pk ，则差商 表示弦 的斜 率,弦 的方程为 O x* xk+1 xk Pk P0 x0 x y y=f(x) 可见，前述迭代格式求得的xk+1实际上是弦 与x轴交点的横坐标（令y=0解出x即可） 因而此法可形象地称为弦截法。 弦截法的收敛速度比牛顿法慢得多。为了 加快收敛速度，我们改用差商 来 代替牛顿迭代格式中的 。于是得到 快速弦截法 迭代格式： 快速弦截法也是迭代法。不过，前面讨论的迭代法在计算xk+1时只用到上一步的结果xk（一步迭代）。而快速弦截法则是一种多步迭代（两步迭代：在计算xk+1时要用到前两步的结果xk和xk-1，因此在使用快速弦截法时，必须给出两个初始近似根x0和x1）。 例 用快速弦截法求方程 的根。设方程的两个初始近似根为x0=0.5 , x1=0.6 表2—5 k xk xk-xk-1 0 0.5 1 0.6 0.1 2 0.56532 -0.03468 3 0.56709 0.00177 4 0.56714 0.00005 与例1（P81）中牛顿法的计算结果相比较，可以看出快速弦截法的收敛速度也是相当快的，迭代到第4步就得到精度 的结果。 第八节 函数插值 在生产实践和科学实验中，经常要遇到研究变量之间的关系的问题，但在许多情况下，往往只能通过实验和观测得到若干个点上的函数值 f(xi)=yi(i=0,1,2, … ,n) 如何根据这些数据估计函数f(x)在给定点x的函数值呢？ 我们设法寻找一个简单的函数p(x)来估计（近似代替）f(x)，使p(xi)=yi(i=0,1,2, … ,n)。 这样的问题称为函数插值问题，p(x)称为插值函数，x0 ,x1, … ,xn称为插值节点。 一、拉格朗日插值公式 因为多项式便于求值、求导等，所以常用多项式作为插值函数，称为代数插值。即：已知函数y=f(x)在n+1个节点x0x1 … xn的函数值y0,y1, … ,yn，求一个次数不超过n次的多项式pn(x)，使函数值pn(xi)=yi(i=0,1,2,…,n)。 几何意义：通过给定的n+1个点(xi,yi)(i=0, 1,2, … ,n)作一条次数不超过n次的代数曲线y= pn(x)，来近似地表示曲线y=f(x)。 可以证明（从略）：代数插值问题的解是唯一的。 1. 线性插值 先讨论最简单的情况。设已知函数y=f(x)在两个点x0,x1的函数值y0和y1,求一个一次多项式p1(x)，满足 p1(x0)=y0 , p1(x1)=y1 从几何上看，y=p1(x)表示一条通过两点(x0, y0),(x1,y1)的直线。此直线的点斜式方程为 若把它按y0,y1整理，则得 （线性插值公式，两点插值公式） O x x0 x1 y A B y=p1(x) y=f(x) y0 y1 （线性插值几何意义图） 为了导出一般的n次插值多项式pn(x)，我们研究一次插值多项式p1(x)的表达式的结构。令 则有 A0(x0)=1, A0(x1)=0, A1(x0)=0, A1(x1)=1 称A0(x), A1(x)是以x0 ,x1为节点的基本插值多项式。于是上述公式可写成 此式