习题三非线性代数系统.PDFVIP

习题三 非线性代数系统 2005 年 秋10.34，KJ ·比尔 王本，马克·Styczynski 2005年9月28 日 2005 年10 月7 日修改 1：2.A.1 （2 分） 请用如下方程笔算雅可比行列式矩阵： 对雅可比行列式矩阵的计算可用对每一个方程对应的变量的推导，就像每一个雅 可比行列式的元素可以这样定义： 通过以上两个方程我们得到如此下的矩阵： 从一个初步的猜想开始： 我们可以通过用如下方程计算p来推到一个可以表示x[1] 当k=0 时，我们能用插入x[0] 到方程（1）（2 ）和（4 ）计算J和f ： 这是标准的方程Ax=b ，可以用高斯消去法来解。现在，这可以用高斯消去法或 一些笔算来完成。或者，如果我们真得很懒的话，像我一样，在Matlab用“\ ” 操作。我们的结果是： 现在我们来决定修正，我们用把x[1] 带入修正方程中来计算x[1] ： 并且用x[1] 作为新“初始猜想” 用于第二步迭代来估算f(x): 如何确定减少牛顿方法的步骤是否成立？根据两个方程解的2 个准则，确定我们 计算结果是否满足如下关系： 我们的计算值的确满足下面方程： 所以，这样的移动不是必需的执行减少步骤： 此题的评分方法如下： 正确的雅可比行列式0.5 分 得到正确的x[1] 0.5 分 如果x[1]不正确，而del_x 或者p 正确，也得0.5 分 通过准则来讨论减少步骤 0.5 分 2 ：2.A.2 (2 分) 此题是题2.A.2 的数学执行。比尔教授已经同意用FSOLVE 作为工具来完成这道 题。并且此题是一种 FSOLVE 中允许的一种近似。下面的编码是此题的解题方 法： 你的输出会是如下形式： 给出的初步猜想是[2；1]，你应该得到[2；2] 的解。如果你开始有其他不同的猜 想，fsolve 会潜在的寻找其他不同的解。记住牛顿方法可以找到很多的解 （谢谢 大卫·刘对此进行的试验。） 尽管如此，如果我们阅读了题干，他要求我们执行牛顿方法的一个步骤并要检查 牛顿方法中的步骤简化是否能被第一步接受。你可以有效的解释这个问题像你自 己执行牛顿方法，或者至少是牛顿方法的一步。这会使你有这样一种观察力来洞 察一个非线性代数体系的方程并将之转化成线性代数体系的方程，这便是正如此 文件中方程6-9 所述雅可比行列式。如下编码是你可以用来解决此问题的方法： 你的输出应该是这样： 并且如果我们在一步后应[f,x]值，那么我们会得到： 这些数和我们笔算出得数相似。我们把所有连续步骤连起来用迭代作图，得到： 得分点分布： 可执行的编码1 分 如果你得到x[1]或x[solution] 中的一个1 分 3 ：2.B.2 （3 分） 此题关系到建立一个正确的质量守恒方程并且利用 FSOLVE 方法来解这一组线 性无关的代数方程以描述CSTR 反应器。这非常像比尔教授在课堂上所说的。我 们解稳态下的浓度随着物料浓度随时间变化而趋近于零。如需我们应用牛顿方 法。此题最大的不同是我们有5 种物料来计算，3 个反应，不同的进料比，和不 同的温度。 我们假设一个可以忽略不计的体积（VCSTR是一个常数）关系物料从反应中的产 生并假设等温操作。这可以向下面所是一样表示： 我们的方程可以写成： 这些方程式是fsolve 中所称得函数[f]=calc_concentrations( )的一部分。 此题要求一个变化的温度（我们已知两个温度下的速率数据）和一个参量微克（这 是进料A 和出料B 的比值）。最简单的执行方法是套入一些回路，针对不同速率 和微克值，重复执行FSOLVE 。 最后我们要做的是开始的猜想：我们可以从一个简单的猜想初始进料A 和B 浓 度，当C,D,E=0 的时候。所以我们的初始猜想是： 有了这些信息为帮助我们能用如下编码来创作一个有关于 A 随着微克和温度而 转化的图表： 你应该得到如下所示的图： 你们的助教马克得到了满足恒速的阿列纽斯关系的结果，并且现在可以作图操作 CSTR 在许多不同的温度下而非特定的温度下。这是非常必要的如函数微克和温 度的作图转化，但是出图的曲线越多会越好的帮助你理解你的程序和CST