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*| (，即((x) ( [x* –