第一章代数方程及代数方程组的求解在化学中的应用.pptVIP

例：解方程组 迭代公式 取初值 第一次迭代 第二次 迭代 Jacobi Gauss —Seidel 直至 收敛 1-2-2-2 Gauss—Seidel迭代法——方法原理 GSD(A,B,N,X,EPS,MAXIT) 初值X(I)=0 I=1,2,…,N IT=0 IT=IT+1 ITMAXIT RETURN no yes yes yes no no 1-2-2-3 Gauss—Seidel迭代法——程序框图 开始 输入：物质名称 给定不同波长下各纯物质及混合物的摩尔吸光度Eij, Ai,即为线性方程组的系数A(N,M)；M＝N＋1 选主元精度EPS 调用Gauss—Seidel迭代法子程序求解线性方程组EC＝A 计算混合物中各组分的浓度Cj 输出：混合物中各组分的浓度Cj 结束 显示程序 显示输出 1-2-2-4 Gauss—Seidel迭代法——应用示例 例：根据分光光度法测得的下列数据，确定五组分混合物中各组分的物质的量浓度。 0.1350 120 16 8 4 2 5 0.2992 24 140 18 4 1 4 0.2700 2 10 90 30 30 3 0.2215 0.9 5 10 120 12 2 0.1135 0.5 1 2 10 100 1 5 4 3 2 1 实验测得的总吸光度 各组分的摩尔吸光系数 波长编号 解：由题意，得到方程组为： 1-2-2-4 Gauss—Seidel迭代法——应用示例 * * 第一章 代数方程及代数方程组的求解在化学中的应用 一元N次非线性 方程的求解 线性方程组 的求解 直接法 间接法 直接法 间接法 N3 二分法 弦截法 插值法 牛顿迭代法 Gauss消去法 Gauss-Seidel迭代法 ★ ★ ★ ★ 例：压力对化学反应中反应物平衡转化率的影响 任何气相化学反应 达平衡时有 起始nB,0/mol a b 0 0 平衡nB/mol a(1-x) b(1-x) cx dx 1-1-1–1 二分法——问题的提出 温度一定，K0为常数， 平衡时，K0与总压p之间关系为： 即 令 上式整理为： 1-1-1–1 二分法——问题的提出 令 此方程为含x的高次方程，x:0.0~1.0之间 用二分法 1-1-1–1 二分法——问题的提出 有函数f(x)=0，求在区间[a,b]内的实根。 设：①函数f(x)在区间[a,b]内连续。 ②函数在a、b两点之值f(a)与f(b)异号，说明在 区间[a,b]中至少有一实根。 a b x0 x1 可得一系列区间套： 它们必须满足： x y f(x) x2 1-1-1–2 二分法——方法原理 1.取x的初值 计算函数f(x0)值，若| f(x0) |EPS，（EPS为10-5），则 x0为所求之根。 （对分） 2.若| f(x0) |EPS，计算并比较函数值，取与f(x0)异号的区间 为新的区间，（可能为[x0,b]或[a,x0]），如图） a [a,x0] [x0,b] b x0 （EPS取10-5） 判据： 1-1-1–2 二分法——方法原理（计算步骤 ） A1=A B1=B F0=F(A) X=(A1+B1)/2 F=F(X) A1=X RETURN HALF(A,B,EPS,F,X) F=0 F0*F0 no B1=X |(B1-A1)/X|EPS yes no yes yes no 改进的二分法：将区间 [a,b]划分成若干小区间 1-1-1–3 二分法 程序框图 开始 输入：M,N(I),K0,P 计算SR,SP,KN A=0.0 B=1.0 A1=A B1=B F1=F(A1) F2=F(B1) 调用二分法子程序求解方程F＝0，解出转化率x 输出：转化率x 结束 定义函数SN，F=F(x) 显示程序 1-1-1–4 二分法——应用示例 计算反应 在773K时，不同压力下的平衡转化率，结果如下：（已知773K时反应平衡常数K0=0. 3007） p/MPa 转化率x 0.1 0.986547 5.0 0.998097 10.0 0.998655 30.0 0.999223 1-1-1–4 二分法——应用示例 例：已知醋酸HAc浓度为c的溶液，