化工应用数学2第二章数学模型基础介绍.pptVIP

3. Euler法示例 例：对一级反应动力学方程 设定：初始浓度[A]0=0.2000mol/L，反应速率常数k=0.01s-1。由上式可推得其求浓度近似值的递推公式为： 2.4 求微分方程数值解 3. Euler法示例 具体解法如下： 1）A列输入时间，间隔为20s，直到140s。 2）B列用Euler法计算在时间t时的浓度[A]t。在B5单元格输入初始浓度0.2000mol/L。B6输入递推公式:=B5-B5*$D$1*$D$2。 3）选定B6单元格，向下拖拽到B12。由于B5单元格为相对引用，每下一格，浓度值更新一次，得到各个时刻由Euler法计算的浓度。 2.4 求微分方程数值解 3. Euler法示例 4）C列是根据 计算的解析浓度值。在C5单元格输入公式： =$B$5*EXP(-$D$1*A5)。 然后自动填充得其余值。 5）D列为Euler法计算值与解析值的相对误差: =100×(解析值-Euler值)/解析值。 2.4 求微分方程数值解 3. Euler法示例 6) 减小步长h，可改进精度。如步长分别为5 和1时，在140s的误差分别降为3.56％和 0.70％。 2.4 求微分方程数值解 4. Runge-Kutta法 Euler法产生较大误差的原因是f(x,y)为曲线，用Teller公式展开： Euler公式只取了线性项(斜率），忽略了高次项。 2.4 求微分方程数值解 4. Runge-Kutta法 R-K法则包括了Teller展开式的高次项，其中最常用的是四阶R-K公式。在递推公式里x取值为：xi, xi +h/2, xi +h。则微分方程通式的四阶R-K公式为： 2.4 求微分方程数值解 4. Runge-Kutta法 微分方程通式的四阶R-K公式为： 2.4 求微分方程数值解 5. Runge-Kutta法示例 例子同Euler法一样。由于微分方程仅涉及因变量y，R-K法四项表达式简化为： 2.4 求微分方程数值解 5. Runge-Kutta法示例 具体步骤： 1）在F5单元格输入初始浓度0.2000mol/L。 2）根据T1-T4的计算公式， 在B6输入：=-$D$1*$D$2*F5 在C6输入：=-$D$1*$D$2*(F5+B6/2) 在D6输入：=-$D$1*$D$2*(F5+C6/2) 在E6输入：=-$D$1*$D$2*(F5+D6) 2.4 求微分方程数值解 5. Runge-Kutta法示例 具体步骤： 3）根据递推公式， 在F6输入：=F5+(B6+2*C6+2*D6+E6)/6 4）选定区域B6:F6，拖拽填充柄到第12行。 5）通过与解析解的相对误差可以发现，R-K法相当精确，完全可以满足一般反应动力学研究的需要。 2.4 求微分方程数值解 5. Runge-Kutta法示例 2.4 求微分方程数值解 本章小结 2.1 拟合建模方法 2.2 迭代法求解非线性方程 2.3 求解线性方程组 2.3 求微分方程数值解 |Di|是方程组的系数行列式，但其中的第i列元素被常数列阵b所取代。例如|D1|是用常数b列阵取代系数矩阵的第一列元素所得的行列式： 2.3 求解线性方程组 示例演示 例.利用Excel求下面三元一次方程组的解。 2.3 求解线性方程组 Excel操作步骤： 1）输入系数矩阵A，列向量b。利用函数MDETERM( )计算系数行列式|A|=-9，方程组有唯一解。 2）在计算行列式|D1|、|D2|和|D3|的值。 3）由 计算得方程组的解。 2.3 求解线性方程组 2.3 求解线性方程组 矩阵解法 线性方程组的矩阵形式： 经变换可得： 因此，若系数矩阵A可逆，则解线性一次方程组就变成简单的矩阵乘法。 2.3 求解线性方程组 示例演示 同样用矩阵法求解例3. Excel造作步骤： 1）输入系数矩阵A，列向量b。 2）利用函数MINVERSE( )求系数矩阵的逆矩阵。 3）利用函数MMULT( )求得方程组的解。 注意：三键确认。 2.3 求解线性方程组 2.3 求解线性方程组 行列式法与矩阵法 行列式法和矩阵法所能求解的方程组的特点？ 多元一次 非齐次 系数矩阵为方阵 系数矩阵可逆或行列式不为0。 2.3 求解线性方程组 Newton-Raphson迭代法 N-R迭代法原理 N-R迭代法属于间接解非线性方程组的方法。原理类似于一元方程求解。然而因有不止一个自变量，导数变为偏导，迭代公式成为迭代增量的线性方程组。 设有非线性方程组： Excel解非线性方程组 Excel解非线性方程组 对