非线性方程(组)的间接求解法.docVIP

第二章 非线性方程（组）的间接求解法 2.0 引言 求n次多项式方程的根：，至今可有对n5的精确求解公式（可查Mathematica），但对一般的代数方程，如：（可视为无穷次），求解变得很复杂，根的位置、重数和数目都是未知的，即使现有的软件也是如此，因此现在问题给出的一般形式为求给定范围内的根。 2.1 求实根的二分法 设f(x)在[a,b]内连续，f(a)f(b)0,则由闭区间上连续函数的性质知，存在，使，不妨设f(a)0, 下面找的近似值：对分[ a,b]得子区间： 若，则得到了一个根； 若，则保留子区间继续讨论； 若，则保留子区间继续讨论，这样的过程进行n次后便得到了方程有根区间的一个序列，满足，且。 于是当n充分大时，可取作为的近似值，此时有。 局限性：只能求一个实根，不能辨识重根或求复根，收敛速度对任何函数均一样且慢，对多元方程的情形须做修改。 2.2 迭代法 对方程f(x)=0,改写成下面的不动点（fixed point）等价形式:，这里要求连续，如取，给定根的初始猜测便可进行迭代：得到，若，则由的连续性便知该极限为根。 几何描述： 收敛定理：设对，，且存在常数0L1,使对任意满足Lipschitz条件，则在[a,b]上存在唯一不动点；且对，迭代公式生成的序列，并有误差估计：。 证明：第一步（的存在唯一性）：构造函数，由定理的条件及f(x)的连续性立得的存在性；再由在[a,b]上满足Lipschitz条件反证唯一性。 第二步（序列的收敛性）： 。 第三步（误差估计）：由于 上式的最后两项分别用于事后估计和事先估计。 问题：若迭代函数不满足Lipschitz条件，则反解出是否应满足？对后面的迭代法解线性方程组的情形也可有此问题。 局部收敛性与收敛阶的概念 定义 若存在的不动点的一个邻域，使得在其中任取初始迭代值均能获得收敛于的迭代序列，则称该迭代公式局部收敛。 定理： 设为的不动点，在的某邻域S内连续，且，则迭代公式局部收敛。 证明：由的连续性知，，继而有，这说明，故由前面的定理知迭代序列收敛于。 注1 直接应用该定理并不方便，因不是预知的。 注2 压缩条件是充分而非必要的，如迭代有不动点=0，但=1，关于其局部收敛性可从并取看出。 例 试用不同迭代法求解。 解：由于方程简单，分别可算出其=-1；0.314；0，因此后两种情况是局部收敛的。 收敛阶（刻画收敛速度的标准之一）：设序列收敛到，记误差，若存在实数及，使，则称序列是p阶收敛的，称为渐近误差常数。P=1时称线性收敛；P=2时称平方收敛；P1时称超线性收敛；=0时称超p阶收敛。 定理：设为的不动点，整数p1，在的邻域连续，且满足，则迭代序列在的邻域是p阶收敛的，且有 。 证明：由条件知迭代序列在的某邻域S是阶收敛的，且再由条件便知结论正确。 Steffensen加速技巧：设是一阶不动点，即有，从而近似有： ，消去C解出。 于是通过构造新的迭代公式（此时的迭代序列记为，的迭代序列仍记为）：。 上面的迭代过程可视为预估校正过程：； 或视为如右新迭代公式：的一步迭代。 Steffensen加速技巧的几何意义由Aitken给出。 定理：若是一阶不动点，，则是二阶不动点。 2.3 Newton迭代法 为求解，把在点附近展成Taylor级数，则得： 取上式的线性部分作为近似方程得： ，于是有: 据此构造不动点形式：，从而有迭代公式： (要求)。 2.3.1 几何意义（切线法）：过以为斜率的直线与x轴的交点横坐标。 注1 若取二次或其他的近似方程，能否构造其不动点形式？ 注2 简化Newton法（少算导数值）： 2.3.2 收敛性与收敛阶： 例 设，从方程利用Newton迭代法求近似。 解：构造迭代公式：, 可以证明半局部收敛性：对，Newton迭代法收敛：由于且，故有极限，再由迭代函数的连续性得证。 2.3.4 Newton下山法 从Newton迭代法的收敛性可见，它仅是局部收敛的，为扩大收敛范围，可改进迭代公式为：，其中---下山因子的选取，要求使，以便于收敛。 例 求的近似解。 解：用Newton迭代法，若由看出迭代不趋收敛，若取（离根较近），迭代3步便得到较好的近似1.32476； 若用Newton下山法，并取（离根较远），，迭代5步仍能得到较好的近似1.32476，由此可见，收敛范围较大。 2.3.5 重根的情形————Newton下山法的应用 设 2.4.1 Newton迭代法的迭代函数，易得：，故为一阶收敛的； 2.4.2 另一种迭代格式：由于f(x)的m重根，相当于的单根，故对用Newton迭代法，则为二阶收敛的，此时的迭代函数为：。 2.4.3 为加速收敛，利用Newton下山法，构造迭代函数：，此时，故为二阶收敛的，但m通常未