方程解数值法精要.pptVIP

第二章 非线性方程的数值解法  简介（Introduction） 我们知道在实际应用中有许多非线性方程的例子，例如 （1）在光的衍射理论（the theory of diffraction of light)中,我们需要求x-tanx=0的根 （2）在行星轨道（ planetary orbits）的计算中,对任意的a和b，我们需要求x-asinx=b的根 （3） 在数学中，需要求n次多项式xn + a1 xn-1+...+an-1 x + an ＝0的根 方程根的数值计算步骤 判断根的存在 确定根的分布范围 根的精确化 Iteration: The act or an instance of iterating; repetition. Mathematics: A computational procedure in which the desired result is approached through a repeated cycle of operations, each of which more closely approximates the desired result. Computer Science: The process of repeating a set of instructions a specified number of times or until a specific result is achieved. 迭代法流图 几何意义: 收敛性分析P32 定义 若存在常数?(0≤ ?1),使得对一切x1,x2∈［a,b］, 成立不等式 |g(x1)-g(x2)|≤ ? |x1-x2|， 则称g(x)是［a,b］上的一个压缩映射， ?称为压缩系数 考虑方程 x = g(x), g(x)?C[a, b], 若 ( I ) 当 x?[a, b] 时， g(x)?[a, b]； ( II )在［a,b］上成立不等式：|g(x1)-g(x2)|≤ ? |x1-x2| 。 则（1）g在［a,b］上存在惟一不动点x* （2）任取 x0?[a, b]，由 xk+1 = g(xk) 得到的序列 {xk}（?［a,b]） 收敛于x* 。 （3）k次迭代所得到的近似不动点xk与精确不动点x*有有误差估计式： Remark： （迭代法的局部收敛定理） 定理条件非必要条件，而且定理中的压缩条件不好验证，一般来讲， 若知道迭代函数g(x)∈C1『a,b］(1),并且满足|g′(x)|≤ ? ≤1(2),对任意的x∈［a,b］(1), 则g（x)是［a,b］上的压缩映射 例题 已知方程2x-7-lgx＝0，求方程的含根区间，考查用迭代法解此方程的收敛性。 在这里我们考查在区间[3.5,4]的迭代法的收敛性 很容易验证：f(3.5)0,f(4)0 将方程变形成等价形式：x＝(lgx+7)/2 举例 用一般迭代法求x3-x-1=0的正实根x* 例题 用一般迭代法求方程x-lnx＝2在区间（2，?）内的根，要求|xk-xk-1|/|xk|=10-8 将方程化为等价方程：x＝2＋lnx k xi 0 3.000000000 1 3.098612289 2 3.130954362 3 3.141337866 4 3.144648781 5 3.145702209 6 3.146037143 另一种迭代格式： 定理 设x*为g的不动点，p≥2为正整数，g在x*的某邻域Ｕ(x*)内p阶连续可微，且 g′(x*)=g″(x*)=…=g (p-1)(x*)=0，而g(p)(x*)≠0， 则存在?0，当x0∈ ［x* - ? ，x*+ ?］ (x0≠x*)时，由迭代法产生的序列｛xk｝以p阶收敛速度收敛于x*. Prove: (1)由g′(x*)=0?必存在?0，当x0∈ ［x* - ? ，x*+ ?］ ? U(x)时，由Picard迭代产生的序列{xk}收敛于x*,并有xk∈ ［x* - ? ，x*+ ?］ (2)由泰勒公式有xk+1=g(xk)=g(x* ）＋g′(x*)(xk- x*)+…+g (p-1) (x*)(xk-x*) p-1/(p-1)! + g (p)(x*+ ?(xk-x*))(xk-x*) p /p! ，0?1.  利用g在x*的各阶导数条件及g(x*)=x*,上式可改写成 (3)由于g在x*处p阶连续可微且g(p)(x*)≠0，知必存在x*的某邻域Ｕ(x*)，当x∈U(x*)时，有g (p) (x)≠0. 由于x*+ ?(xk-x*