数值方法 第4章　非线性方程.doc

第4章　 非线性方程数值解法 1　引言 其中特别 此称为代数方程。当时，其根不能用和开方的有限次运算公式来表示，因此非线性方程的数值解法是迭代法 ，若有 为正整数，满足，那么称为的重零点，或称为方程的重根，那么称为方程的单根 ，那么有 ，如果，那么方程在内至少有一个根，称为方程的有限区间。 求方程的有根区间，其中 由此看出，的三个有根区间为 ，可以用二分法来缩小有根区 间，设取中点； 若，则令新的有根区间为 ；若，则令； … 设的有根区间，即 ，取的中点，若 ，令，若 ，那么取， 如此可以找到小的有根区间。 二分法可以求方程的根的近似 值，但收敛速度相当慢，一般不采用。 §2　迭代法及其收敛性 (I)　不动点和不动点迭代 是一元连续函数，要解的方程是 　　　　 　　　　　　　　　（1） 为进行迭代，把(1)变成等价形式 　　　　 　　　　　　　　　（2） 其中是连续函数。 由（1）变成等价形式(2)，不是唯一的。 例2.1　，　 令　　 ， 令　　，　　 显然，，但，都 与是等价的。 由可以构造迭代公式 　　　（3） 称为迭代函数。 如果有　，由的连续性得 　　　　　 满足方程（2），称是的一个不动点， 也是方程的根，迭代方法（3）称 为不动点迭代法。 例2.2求方程在附近的根 ，其中 解　将方程写成等价形式 取，迭代有 　　 可以看出，迭代收敛。 的近似值。 再将方程写成另一等价形式 　，　 , 　 发散，不收敛 （II）整体收敛性 定理2.1　设 如果有 那么在上一定存在不动点 如果满足的条件，并且存在常数 使得 那么在上的不动点是唯一的。 证明 由于，所以有 如果上两式中有一等号成立，那么在 上有不动点。若两式中 都为严格不等号，， 由的连续性，必存在??使得 即为的不动点。 （唯一性）设有两个不动点 ，并且，那么 得出矛盾不动点唯一。 推论　　设，并满足 　　　　； 还存在常数，使得 那么在上存在唯一的不动点 定理2.2 设　，并满足 　 　存在常数使得 那么对任，由 产生的序列收敛到在上的不动 点，并对于整数有 证明 由题设条件，在上有唯一的 不动点。由条件保证， ；并有 　 递推有 令有， 即 令　有 。 推论　设　， 并， 那么定理结论成立。 例子2.3　解方程，其中 ，由于， 所以方程在（1,2）内有根（单根！） （i）　 　　　 取，如　 这样不能进行迭代，故不会收敛 （ii）　 　　　　 ，　，　　 从而有 　　 ，； 迭代收敛 （II）局部收敛性与收敛阶 前面讨论了迭代法在区间上的收敛 性，此收敛称为全局收敛性（整体收敛）。 很多情况我们不能这样考虑。实际应用时 只在不动点附近来考虑其收敛性，此称 为局部收敛性。 定义2.1　设在区间I上有不动点， 如果存在的某一邻域 ，对任，由迭代法产生 的序列且收敛到，则称迭代 法局部收敛。 定理2.3　设为的不动点，在 的某一邻域连续，并且。则迭代 法局部收敛。 证明 由于在的邻域连续，所以任 给，使得当时有 取使，从而有 其中　　； 注意到取，那么 根据定理2.2，推论知迭代法 上收敛，即局部收敛。 定理2.3中，导数条件可以改为 。 对于收敛快慢，引入收敛阶的概念 定义2.2　设序列收敛到，记误差 。若存在实数，及非零常数 c，使 则称为阶收敛，c为渐近误差常数。 当时分别称为线性收敛，平方收 敛。 定理2.4　设为的不动点，整数 ，在的邻域上连续，且满足 则由产生的序列在的邻域 是阶收敛的，且有 证：因定理2.3保证 局部收敛性，取充分按近于 的，设，有? 由Taylor展开 在之间，利用条件 利用的连续性，两边取极限 例子2.4 讨论 在例2.3 解法(ii)的收敛阶 解 方法一阶 §3 Newton迭代法 （I）Newton迭代法的计算公式 为求解的根，设有一个近似 值，用Taylor公式 在之间 其中假定，若，可得 若上式中右边最后一项忽略不计，作为 的新近似值，有 这称Newton迭代法，也称Newton- Raphson迭代法 Newton迭代为不动点迭代 定理3.1 设满足，若 （表示为单根），并在的闭邻域上 连续，那么Newton法　　 至少二阶收敛而且 证： 为其不动点 , 收敛至少是二阶的。 Taylor展开有 即 相减有 例3.1 用Newton法解方程，其中 解 有根区间，仍取 ，计算 例3.2 用Newton法求 解： 取 当 （II）重根情况 设为的重根， 。对于