非线性方程(组)解法.pptVIP

第3讲 非线性方程(方程组)的数值解法 * 介绍求解非线性方程 的几种常见和有效的 数值方法，同时也对非线性方程组 求解，简单介绍一些最基本的解法。无论在理论上， 还是在实际应用中，许多问题的求解，在解析方法无 能为力时，数值方法则可以借助于计算机出色完成。 * 2.1二分法 概念： 有根区间：存在根 隔根区间：唯一根 * 根的存在性 首先确定有限区间：依据零点定理。 设 ，且 ，则方程 在区间 上至少有一个根。如果 在 上恒正或恒负，则此根唯一。 * 二分法 用二分法(将区间对平分)求解。 令 若 ，则 为有根区间，否则 为有根区间 记新的有根区间为 ， 则 且 * 二分法 对 重复上述做法得 且 * 二分法 设所求的根为 ， 则 即 取 为 的近似解。 * 二分法 迭代终止准则 * 即 2.2一般迭代法 2.2.1 迭代法及收敛性 对于 有时可以写成 形式 如： * 迭代法及收敛性 考虑方程 。 这种方程是隐式方程，因而不能直接求出它的根。 但如果给出根的某个猜测值 ， 代入 中的右端得到 ，再以 为一个猜测值， 代入 的右端得 ，反复迭代 得 * 迭代法及收敛性 若 收敛，即 则得 是 的一个根 * 简单迭代法 将 变为另一种等价形式 。 选取 的某一近似值 ，则按递推 关系 产生的迭代序列 。这种方法算为简单迭代法。 * 例题 例2.2.1 试用迭代法求方程 在区间(1，2)内的实根。 解：由 建立迭代关系 k=0,1,2,3…… 计算结果如下: * 例题 精确到小数点后五位 * 例题 但如果由 建立迭代公式 仍取 ，则有 ， 显然结果越来越大， 是发散序列 * 2.3 Newton迭代法 设x*是方程f (x) = 0的根， 又x0 为x* 附近的一个值， 将f (x) 在x0 附近做泰勒展式： * Newton迭代法 有： 即 以 代替 重复以上的过程，继续下去得： * Newton迭代法 * 迭代终止的标准 或 3.非线性方程组的迭代解法 或 * 的Jacobi矩阵 3.非线性方程组的迭代解法 * 将 在 处进行泰勒展开 一元函数 为向量 建立 3.非线性方程组的迭代解法 * 迭代初始向量 迭代终止的标准 或 Matlab中范数：norm(x) 4.Matlab函数 * 非线性方程 Z=fzero(filename,x0,tol) filename: f(x)函数文件 X0: 迭代初值 tol：精度参数 非线性方程组 Z=fsolve(filename,x0,option) option: 最优化工具箱的优化参数 * *