第二部分非线性方程的数值解法8.pptVIP

提示 压缩影响原理的应用 迭代法的局部收敛性 定义2.1：对于方程 定理2.3 迭代法的局部收敛性 第二章 非线性方程的数值解法 2.1 二分法 2.2 非线性方程求解的迭代法 2.3 非线性方程求解的matlab函数 历史背景 代数方程的求根问题是一个古老的数学问题。理论上， 次代数方程在复数域内一定有 个根(考虑重数)。早在16世纪就找到了三次、四次方程的求根公式，但直到19世纪才证明大于等于5次的一般代数方程式不能用代数公式求解，而对于超越方程就复杂的多，如果有解，其解可能是一个或几个，也可能是无穷多个。一般也不存在根的解析表达式。因此需要研究数值方法求得满足一定精度要求的根的近似解。 在科学研究和工程设计中, 经常会遇到的一大类问题是非线性方程: f(x)=0 (2.1) 的求根问题，其中f(x)为非线性函数。 方程f(x)=0的根, 亦称为函数f(x)的零点。 当f(x)不是x的线性函数时，称对应的函数方程为非线性方程。如果f(x)是多项式函数，则称为代数方程，否则称为超越方程（三角方程，指数、对数方程等）。 求解方程 的数值解法大致可以分为三步骤： （1）根的存在性，方程是否有根？如果有，有几个 对于多项式方程，n次方程有n个根。 （2）根的隔离，把区间分为较小的子区间，每个子区 间或者有一个根，或者没有根。这样可以将有根 区间内的任一点都可看成根的一个近似值 （3）根的精确化，对某个近似值设法逐步精确化，使 其满足一定的精确要求 2.1 二分法 二分法又称二分区间法,是求解方程(2.1)的近似根的一种常用的简单方法。 设函数f(x)在闭区间[a,b]上连续,且f(a)f(b)0, 根据连续函数的性质可知, 方程 f(x)= 0 在(a,b)内必有实根,称区间[a,b]为有根区间。 为明确起见,假定方程 f(x)= 0 在区间[a,b]内有惟一实根x*。 二分法的基本思想是: 首先确定有根区间,将区间二等分, 通过判断f(x)的符号, 逐步将有根区间缩小, 直至有根区间足够地小, 便可求出满足精度要求的近似根。 由高等数学知识可知, 设 f (x) 为区间 [a,b] 上的单值连续函数, 如果 f (a)·f (b)0 , 则[a,b]中至少有一个实根。如果 f (x) 在[a,b]上还是单调递增或递减函数，则方程f (x)=0仅有一个实根。 记笔记 由此可大体确定根所在子区间，方法有： (1) 画图法 (2) 逐步搜索法 y=f(x) a b y x ① 取有根区间[a,b]之中点, 将它分为两半,分点 ,这样就可缩小有根区间。 2.1.2 二分法求根 设方程f(x)=0在区间[a,b]内有根,二分法就是逐步收缩有根区间，最后得出所求的根。具体过程如下： ② 对压缩了的有根区间 施行同样的手法, 即取中点 ,将区间 再分为 两半,然后再确定有根区间 , 其长度是 的二分之一。 ③ 如此反复下去,若不出现 ,即可得出 一系列有根区间序列： 上述每个区间都是前一个区间的一半,因此 的长度： 当k→∞时趋于零,这些区间最终收敛于一点 x* 即为所求的根 。 每次二分后,取有根区间 的中点 作为根的近似值，得到一个近似根的序列： 该序列以根x*为极限 　 只要二分足够多次(即k足够大),便有 这里ε为给定精度,由于 ,则： 当给定精度ε＞0后,要想 成立,只 要取k满足 即可，亦即当: 时,做到第 k+1 次二分,计算得到的 就是满足精度要求的近似根 。 在程序中通常用相邻的 与 的差的绝对值或 与 的差的绝对值是否小于ε来决定二分区间的次数。 二分法算法实现 二分法的优点是不管有根区间 [a,b]多大,总能求出满足精度要求的根,且对函数f(x)的要求不高,