化工过程模仿与分析(第四章非线性方程组的求解).pptVIP

三、弧长参数微分法(Path parameterization by arc length) 曲线x (t)的弧长 x t dt dx ds 对于n维向量函数xi=xi(t)，下式成立： 此即广义Pythagorean定理 根据广义勾股定理，可利用弧长参数微分法求取同伦曲线 悟查坛划骤孝象烧椿乱你茫绊鸳哺艺丸慢刚软担砌令涪找危尼甩骏昔欢冲化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 记 则有： n个常微分方程，n+1个自变量 与广义勾股定理联立即可求解 该问题如何求解最方便？ 数值方法每求解出dx及dt后，可应用牛顿－拉佛逊法修正一次以使(x (t), t)更接近同伦曲线。 初始条件为： 该方程组方程数比自变量数少1 联立求解则可获得修正值 痈护枯街瘦客底潜酋狄蔑宝睬缺襄捧震俺疟燃凡爵淮针仑荒持功傅铸卑哮化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 弧长参数微分法步骤： （1）给定初始值，循环（2）~（4）直至收敛； （2）求解弧长常微分方程获得dx, dt, 更正x及t; （3）使用牛顿－拉佛逊法获得δx, δt，再次更正x和t; 聪偶婆扇火浇炉悟妙坛怜忻鳃捉杠孩郎坞受恭煮掣根铅陋玫仑您吮粟谰皋化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 万锡人 计算非线性方程组所有解的同伦连续法 南京化工学院学报，vol 16, No. 3，88～94，1994. 同伦算法在化工中的应用： 多重定态问题、反应蒸馏、全局优化 T.L. Wayburn J.D. Seader. Homotopy continuation method for computer-aided process design. Comput. Chem. Engng, 11 (1), 7~25, 1987. J.D. Seader etal. Mapped continuation methods for computing all solutions to general systems of nonlinear equations. Comput. Chem. Engng, 14 (1), 71~85, 1990. 剂恕嗓荒屋肾陨上骚箭曰掀诱纶喧插律舆剔苟彭吮季痈辑柄言剩萍揣七澎化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 作业1 用牛顿－拉佛逊法求解下列非线性方程组，分别以初始点x1 = 0, x2 = 0及x1 = 2, x2 = 2，检测是否收敛到同一个解。 作业2 用Wegstein法求解下列显示非线性方程组 注：求解区域为0≤ x1, x2≤ 1，自行确定初始解 端锭菱个续语幸榨舱墅其诲奉樱符旺办五敛寺偏颊考锦茨涡页妒渡屹胁甜化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 作业3 用数值延拓法求解下列非线性方程组，其中将同伦参数t离散化为8个区间。 初始点：x1 = 1, x2 = 0 注：该方程组以牛顿－拉佛逊法求解时不能收敛。 汀谗仅涟遥绷妊谤邯绽宙聪阉忍粥喊伙徐森蛋批辗壹娇吧对烁选又丑键覆化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 邹斥滓痴肠珠讨因旷浮拷劣粕嘛烙恼践箔窜诣冰伟茵难铡排帜泥拐抵掣对化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 第四章 非线性方程组的求解 4.1 概述 一、非线性方程组的一般形式 解法： 迭代逼近的方法：牛顿法等 最优化方法 或 毋名圈托汐咎役翠言毅略柔斥斡大掂颈烧苫蚂吴绚械简棘蔑诡却踌喻怀况化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 二、解的存在性—压缩映射原理 不动点 压缩映射 压缩映射原理 若G在闭集D0上是压缩映射，且G D0 ? D0 ，则G (x)在D0上存在唯一的不动点 排妈逆捂碧质呕抒赠优赚蛆滥诱慕孰洗任挎按红狙括蝶尤僧捐迫涎江巡瑰化工过程模拟与分析(第四章非线性方程组的求解)化工过程模拟与分析(第四章非线性方程组的求解) 收敛性 对于某个非线性方程组，若其解为 ，某一算法产生 迭代序列x (k)，且有下式成立，则称该算法收敛。 收敛速度 对于一个收敛于解的序列x (k)，若存在一个正实数α及一个与迭代步数k无关的正常数q，由某k0开始下式成立，则称该迭代序列具有α阶收敛速度。 溅犬竞蕴赛壮奸亢汽幽浊魄扑砰绸膏尸街蜗惧成荚否整喝佑湃歹牡酋彤捐化工过程模拟与分析(第四章非线性方程组的求解)化工过程模