天津大学化工数学教程三单元非线性方程组的数值解课件.pptVIP

§3 非线性方程组的数值解 问题1：对串联的油换热器组进行最优化设计时， 得到如下方程组： 求解最佳设计时，油换热器进出口温度。 T1 T2 一、高斯－雅可比迭代法 迭代格式： i=1,n 基本思想：构造一种迭代格式，然后从初始值出发 逐次迭代，直至收敛。 二、高斯－赛德尔迭代 迭代格式： 与非线性方程的收敛条件一致！！ 三、迭代的收敛性分析 即： 1、 收敛条件 2、 收敛准则 绝对收敛准则 相对收敛准则 i=1,n 3、 收敛速度 P=1 例1 用雅可比迭代求解问题1 解：首先将方程改写: （1） （2） 取初值为：T(0)1=180, T (0)2=292 收敛性判断 结论：迭代格式可能不收敛！！ 在初值处 结论：迭代格式收敛 重新改写迭代格式为: 在初值处 （1） （2） 迭代7次达到收敛精度要求 次数 T1 T2 -------------------------------------------------------- 1 180.00000 295.11912 2 181.74554 295.53602 3 181.98080 295.59234 4 182.01262 295.59996 5 182.01693 295.60099 6 182.01751 295.60113 7 182.01759 295.60115 构造雅可比迭代格式: 取ε＝10-4 四、松弛法 第一步作赛德尔迭代 第二步引入松弛因子ω作线性加速 迭代格式： 松弛迭代收敛的必要条件： 0ω2 取1 ω 2 用于加速某收敛的迭代过程 ，超松弛 取0 ω 1 用于非收敛迭代过程使其收敛 ，亚松弛 ω＝1？？？ 高斯赛德尔迭代 展成X（k）附近的台劳级数，取线性近似 五、牛顿拉夫森法 设 上述方程组可改写为： 雅可比矩阵 解方程组，可得： (i=1, n) 则： 牛顿拉夫森法迭代格式 迭代计算步骤如下： （1） （2） 取初值为：（180,292） 解 例2 用牛顿拉夫森法求解问题1 迭代3次收敛 从初值 出发进行迭代 牛顿-拉夫森法求解结果 ------------------------------------------------------ No T1 T2 ------------------------------------------------------ 0 180.00000 292.00000 1 181.95986 295.52774 2 182.01757 295.60112 3 182.01760 295.60115 P=2 应用牛顿拉夫森法注意