化工系统第2章精编.pptVIP

* * 4. 求雅可比矩阵的逆矩阵时工作量大 在求解过程中，雅可比矩阵有可能是奇异矩阵。两种原因： 建模型时出错，此时需对模型进行检查修正； 迭代过程中出现的（相当于单变量情况下 时的情形），此时可重设初始值。 * * 5. 解析求导困难时，可用下式进行数值求导： Δxi 为差分步长，xi 一般可取为0.001 * * 例：对串联的油换热器组进行最优设计时，得到如下方程组 求油换热器进出口的温度T1 和 T2，已知初始值是（180，292）。要求精度为0.01 解：（1）写出 即： * * （2）计算 （3）求出 J * * （4）计算 （5）求逆 （6）迭代 （7）重复 2 至 6 步，直到满足计算精度要求 计算 若：J 的行列式为非奇异（满秩）方阵，即 则 J 的逆矩阵为： * * J 的伴随矩阵 矩阵求逆 求 J*：先求 中元素的代数余子式： 则 * * 三、布洛伊顿（Broyden）拟牛顿法 牛顿-拉夫森法优点：有一定理论基础，收敛速度较快，效果较好。 缺点：迭代公式中含有一个待解方程组中函数向量 的雅可比矩阵的逆矩阵 雅可比矩阵的求取，是同方程组数学形式直接有关的运算，当方程组的规模较大时，相当麻烦。 用解析求导，工作量很大，有时还有困难；用数值求导也并不简单。 * * 比照着牛顿-拉夫森法的迭代公式把迭代公式写成如下形式： （2-31） 如果 ，（2-31）式就是牛顿-