第2章 非线性方程和方程组的数值解法.ppt

弦截法收敛定理 求解方程f(x)=0的快速弦截法 2.5 非线性方程组的牛顿方法 设二阶方程组 其中x,y为自变量。为了方便起见，将方程组写成向量形式： 将 在（x0,y0）附近进行二元泰勒展开，并取其线性部分，得到下面方程组： 2.5 非线性方程组的牛顿方法 令 则有 如果 再将原方程组在u1处进行二元泰勒展开，并取其线性部分 2.5 非线性方程组的牛顿方法 得到下面方程组： 解出 得出 继续做下去，每一次迭代都是一个方程组 为止。 (2-12) 2.5 非线性方程组的牛顿方法 实例 例 求解下面非线性方程组 取初始值 解： 解方程得 2.5 非线性方程组的牛顿方法 实例 继续做下去，直到 时停止。 2.6 应用实例 例 在合成氨生产中，烃类蒸气发生以下转化反应： 已知进料甲烷为1mol，水蒸汽为5mol，反应后总压P=1atm，反应平衡常数为: 试求反应平衡时各组分的浓度。 解：设反应平衡时有x摩尔甲烷转化成CO，同时生成的CO中又有y摩尔转化成CO2，则反应平衡时各组分的摩尔数及分压如下： 将平衡时各组分的分压表达式代入反应平衡常数KP1及KP2的表达式得: 设x的初值为0.1， y的初值为0.05，若采用直接迭代法进行计算，可得： x1=0.999988 y1=1.8707143 若采用直接迭代的方法求解该方程组，结果是发散的，无法得到真实解。但是，若采用松弛迭代法求解，并取松弛因子ω小于0.49，则可得到收敛解。其最后的求解结果为: x*=0.9437 y*=0.6812 * * * 松弛迭代法 实例 例 用（松弛）迭代法求解下面非线性方程组，并分析松弛因子对迭代次数及收敛过程的影响。已知迭代初值x和y均为0，收敛精度ε=0.001 。 解：取以下迭代表达式： 若取松弛因子为1.1，则其迭代过程如表2-2。 若改变松弛因子,迭代过程及迭代所需的次数亦将发生变化，详见表2-3。 表2-2 迭代过程 表2-3 松弛因子及迭代次数的变化 ，则为直接迭代法 2.3 Newton迭代法 设x * 是方程f (x ) = 0的根，又x0 为x * 附近的一个值 ，将f (x ) 在x0附近做泰勒展式 令 ，则 Newton迭代法 去掉 的二次项，有： 即 以x1代替x0重复以上的过程，继续下去得： Newton迭代法 以此产生的序列{Xn}得到f(x)=0的近似 解，称为Newton法，又叫切线法。 Newton迭代法几何解释 几何意义 例题 例2.3.1 用Newton法求 的近似解。 解：由零点定理。 例题 例题 例2.3.2 用Newton法计算 。 解： Newton迭代法算法框图 Newton迭代法算法 Newton迭代法收敛性 定理2.3.1 设函数 ，且满足 若初值 满足 时，由Newton法产生的序列收敛到 在[a,b]上的唯一根。 Newton迭代法收敛性 证明： 根的存在性 根的唯一性 Newton迭代法收敛性 收敛性 Newton迭代法收敛性 Newton迭代法收敛性 Newton迭代法收敛性 推论 在定理2.3.1条件下， Newton迭代法具有平方收敛速度。 代数方程的Newton迭代法 代数方程的Newton迭代法推导 设n次代数方程 用Newton迭代法求有限区间的实根，则要计算 ，一般采用秦九韶算法。 代数方程的Newton迭代法 由Taylor展式 代数方程的Newton迭代法 代数方程的Newton迭代法 同理 代数方程的Newton迭代法 比较x的同次幂系数得： 故代数方程的Newton迭代公式 代数方程的Newton迭代法算法 2.4弦截法 Newton迭代法有一个较强的要求是 且存在。因此，用弦的斜率 近似的替代 。 弦截法 令y=0