4计算函数零点和极值点的迭代法.docVIP

第4章 计算函数零点和极值点的迭代法 本章讨论非线性方程（组）的求解问题 4.1 不动点迭代法及其收敛性 1．不动点 设非线性方程组f(x) = 0 (4.1-1) 等价： x = ?(x) (4.1-2) 则有迭代格式：x(k+1) = ?(x(k))，k = 0，1，2，… 若?连续，且迭代序列{x(k)}收敛到x*，则两边取极限得x* = ?(x*)，即x*满足(4.1-2)，从而满足(4.1-1)，即x*为f零点。 称x*为?(x)的不动点。 注：(1) 求零点求不动点 (2) ?(.)称为迭代函数，{x(k)}称为迭代序列 (3) 不同方法构造迭代函数，得不同的迭代序列 2．迭代法的基本问题 (1) 如何构造适当的迭代函数?(.)使迭代序列{x(k)}收敛 (2) 收敛的速度和误差 (3) 如何加速 4.1.1 解一元方程的迭代法 1. 根的隔离 设一元方程f(x) = 0，f连续，其实根可能有很多，需将各根隔离，即f在[a，b]内有且仅有一根。 方法：设f ?C[a，b]，f(a)f(b) 0，且f在[a，b]上单调，则f在[a，b]内有且仅有一根。 2. 迭代序列的收敛性 因为可以有多种迭代函数，所产生的迭代序列{x(k)}有可能： (1) 收敛快 (2) 收敛慢 (3) 不收敛 例1 f(x) = x3 – x – 1 = 0，求f在x = 1.5附近的根，初值取x(0) = 1.5。(p328) 取——收敛 取?(x) = x3 – 1——不收敛 定理4.1-1 (1) 设?(x) ? C[a，b]，且对于任意x ? [a，b]有? (x) ? [a，b]，则?在[a，b]上必有不动点。 (2) 进一步，若?(x) ? C1[a，b]，且存在L 1，使对于任意的x ? [a，b]有|? (x)| ? L 1 (4.1-4) 则对于任意的x(0) ? [a，b]，x(k+1) = ?(x(k))收敛于唯一不动点x* ? [a，b]。且有 (4.1-5) 证明：(1) 若?(a) = a或?(b) = b，则结论显然成立 现设a ?(a)，?(b) b，令?(x) = ?(x) – x，则?(x) ? C[a，b]，且?(a) = ?(a) – a 0，?(b) = ?(b) – b 0，故存在x* ? [a，b]，使?(x*) = 0，即?(x*) – x* = 0 ? ?(x*) = x* (2) 设?(x) ? C1[a，b]，且满足(4.1-4)，若有x1*，x2* ? [a，b]，x1* ? x2*，?(xi*) = xi* i = 1，2 由微分中值定理，|x1* – x2*| = |?(x1*) – ? (x2*)| = |?(ξ)| |x1* – x2*| ? L|x1* – x2*| |x1* – x2*| 矛盾，所以不动点唯一。 由|x(k) – x*| = |?( x(k-1)) – ? (x*)| ? L|x(k-1) – x*| ? … ? L k|x(0) – x*| 及0 L 1知 即x(k)收敛于x*。 并且由|x(k) – x*| ? L|x(k-1) – x*|得 |x(k) – x*| ? L|x(k-1) – x(k) + x(k) – x*| ? L|x(k-1) – x(k)| + L|x(k) – x*| 从而有 又因 |x(k) – x(k-1)| = |?(x(k-1)) – ? (x(k-2))| ? L| x(k-1) – x(k-2)| ? …? L k-1| x(1) – x(0)|，代入上式的右边，即得 注：(1) 若L ? 1，则收敛很慢，须考虑加速问题 (2) (4.1-5)中第一式为后验公式—迭代停止准则 第二式为先验公式—预测迭代次数 (3) 定理是以[a，b]中任一点作初值，迭代都收敛，称为全局收敛。 （此要求较难满足，故考虑在） x*的某一邻域内—局部收敛性 定理4.1-2 设x*为?的不动点，? 在x*的某邻域内连续，且|? (x*)| 1，则存在? 0，只要x(0) ? [x* –?，x* + ?]，迭代法x(k+1) = ?(x(k))收敛。 证明：∵|?(x*)| 1，及? (x)在U(x*)内连续，?存在? 0，使[x* –?，x* + ?] ? U(x*)，且 ? x ? [x* –?，x* + ?]有|? (x)| ? q 1 ? ? x ? [x* –?，x* + ?]，|?(x) –x*| = |?(x) – ?(x*)| = |?(ξ)| |x – x*| ? q|x – x*| ?，即?(