《第4章非线性方程(组)的求解》-公开课件.pptVIP

4.1 二分法 4.4 抛物线法 4.5 非线性方程组的求解 非线性方程组可以看作非线性方程的推广，而非线性方程就是非线性方程组的特例。非线性方程组的一般数学描述为： 非线性方程的MATLAB函数求解 1、fzero()函数 MATLAB优化工具箱提供的fzero()函数是专门用于求解单变量非线性方程根的函数，该函数的调用格式为： [x,fval,exitflag,output]=fzero(fun,x0,options,p1,p2,…) 其中，fun表示函数表达式，x0是初始值，可以是标量或长度为2的向量，options是设置的过程参数，它主要包括Display和TolX两个选项，options选项可以通过函数optimset来设定，p1,p2,…是函数表达式中附加的参数，x是返回的根，fval是根x处的目标函数的值，exitflag表明解存在的情况，正数表明解存在，负数表示解不存在(遇到复数、NaN或者无穷大等)。Output包含计算过程中的信息，它是一个结构体，output.algorithm是所选用的算法，output.funcCount是函数赋值次数，output.iterations是迭代次数。 2、fsolve()函数 MATLAB最优化工具箱提供的fsolve()函数是专门用来求解多元方程的实数根的函数。它的调用格式如下： [x,fval,exitflag,output,jacobian]=fsolve(fun,x0,options,p1,p2,…) 其中，大多数参数的意义同fzero()函数，返回的jacobian是函数在x处的Jacobi矩阵， * * 第4章 非线性方程(组)的求解 4.1 二分法 4.2 简单迭代法 4.3 Newton法 4.4 抛物线法 4.5 非线性方程组的求解 4.6 实例解析 本章目标：求 f (x) = 0 的根 原理：若 f ?C[a, b]，且 f (a) · f (b) 0，则 f 在 (a, b) 上必有一根。 a b x1 x2 a b x* bisect.m 误差 分析： 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： 优点：①简单; ② 对f (x) 要求不高(只要连续即可) . 缺点： ①无法求复根及偶重根 ② 收敛慢 注：用二分法求根，最好先给出 f (x) 草图以确定根的大概位置。或用搜索程序，将[a, b]分为若干小区间，对每一个满足 f (ak)·f (bk) 0 的区间调用二分法程序，可找出区间[a, b]内的多个根，且不必要求 f (a)·f (b) 0 。 多用于为其它求根方法提供初始近似值。 试位法 为了加快二分法根的收敛速度，这里再介绍一种方法——试位法，试位法的一般执行过程见下面动画。 a b (a+b)/2 x* (a, f (a)) (b, f (b)) test_bit.m f (x) = 0 x = g (x) 等价变换 f (x) 的根 g (x) 的不动点 思路 从一个初值 x0 出发，计算 x1 = g(x0), x2 = g(x1), …, xk+1 = g(xk), … 若 收敛，即存在 x* 使得 ，且 g 连续，则由 可知 x* = g(x* )，即x* 是 g 的不动点，也就是f 的根。 逐次逼近： 将隐式方程归结为显式计算 4.2 简单迭代法 fixpt.m x y y = x x y y = x x y y = x x y y = x x* x* x* x* y=g(x) y=g(x) y=g(x) y=g(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 原理：将非线性方程线性化 —— Taylor 展开 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间. 将 (x* ? x0)2 看成高阶小量，则有： 线性 /* linear */ x y x* x0 只要 f ?C1，每一步迭代都有f ’( xk ) ? 0， 而且 ，则 x*就是 f 的根。 切线法 4.3 Newton法 newton.m 牛顿下山法 ——Newton’s Method 局部微调： 原理：若由 xk 得到的 xk+1 不能使 | f