数值微分与数值积分.pptxVIP

会计学 1 数值微分与数值积分 2 2022-10-22 微积分是高等数学中的重要内容，在化学工程上有许多非常重要的应用 微积分的数值方法，不同于高等数学中的解析方法，尤其适合求解没有或很难求出微分或积分表达式的实际化工问题的计算，例如：列表函数求微分或积分 引言1---数值微积分方法不同于解析方法 第1页/共47页 3 2022-10-22 数值微分和数值积分与插值和拟合往往是密不可分的 如在进行数值微分时，常针对离散的数据点，利用插值和拟合可以减少误差；而数值积分的基本思路也来自于插值法。例如如果所积函数的形式比较复杂或以表格形式给出，则可通过构造一个插值多项式来代替原函数，从而使问题简化 引言2---数值微分和数值积分与插值和拟合关系密切 第2页/共47页 4 2022-10-22 §5.1 数值微分 化工领域的实际问题中时常需要求列表函数在节点和非节点处的导数值，这是数值微分所要解决的问题。数值微分方法可近似求出某点的导数值 例如在反应动力学的研究中，根据实验数据确定反应的动力学方程 ： 这里实验测得一批离散点，要计算 只能借助数值微分求导解决 表5-1 反应动力学实验数据 t1 t2 … tn pA1 pA2 … PAn 0 导入 第3页/共47页 5 2022-10-22 0.1 建立数值微分公式的三种思路 常用三种思路建立数值微分公式： 从微分定义出发，通过近似处理，得到数值微分的近似公式 从插值近似公式出发，对插值公式的近似求导可得到数值微分的近似公式 先用最小二乘拟合方法根据已知数据得近似函数，再对此近似函数求微分可得到数值微分的近似公式。然后对各方法数值微分后得到的多项式求值，即可求出任意点处的任意阶微分 第4页/共47页 6 2022-10-22 1 方法概述 在微积分中，一阶微分的计算可以在二相邻点x+h和x间函数取下列极限求得： 取其达到极限前的形式，就得到以下微分的差分近似式： 注：高阶微分项可以利用低阶微分项来计算，如二阶微分式可以表示为： 对应的差分式有： §5.1.1 差分近似微分 上式中三种不同表示形式依次是一阶前向差分、一阶后向差分和一阶中心差分来近似表示微分。其中一阶中心差分的精度较高。 第5页/共47页 7 2022-10-22 2 差分的Matlab实现 在Matlab中，可用diff函数进行离散数据的近似求导 调用形式：Y = diff(X,n) 其中：X表示求导变量，可以是向量或矩阵。如是矩阵形式则按各列作差分； n表示n阶差分，即差分n次 Y是X的差分结果 注：用diff函数进行离散数据的近似求导与前向差分近似， 但误差较大。最好将数据利用插值或拟合得到多项式，然后对近似多项式进行微分 第6页/共47页 8 2022-10-22 例5.1：丁二烯的气相二聚反应如下： 实验在一定容器的反应器中进行，3260C时，测得物系中丁二烯的分压 (mmHg)与时间的关系如表5-2所示。 。 表 5-2 丁二烯二聚反应实验数据 t(min) (mmHg) t(min) (mmHg) 0 632.0 50 362.0 5 590.0 55 348.0 10 552.0 60 336.0 15 515.0 65 325.0 20 485.0 70 314.0 25 458.0 75 304.0 30 435.0 80 294.0 35 414.0 85 284.0 40 396.0 90 274.0 45 378.0 用数值微分法计算所列时刻每一瞬间的反应速率 第7页/共47页 9 2022-10-22 解：程序如下： t=[0:5:90]; pA=[632.0 590.0 552.0 515.0 485.0 458.0 435.0 414.0 396.0 378.0 362.0 348.0 336.0 325.0 314.0 304.0 294.0 284.0 274.0]; dt=diff(t); ％ 求时间t的差分 dpA=diff(pA); ％ 求压力的差分 q=dpA./dt ％ q为数值微分结果 执行结果： q = Columns 1 through 8 -8.4000 -7.6000 -7.4000 -6.0000 -5.4000 -4.6000 -4.2000 -3.6000 Columns 9 through 16 -3.6000 -3.2000 -2.8000 -2.4000 -2.2000 -2.2000 -2.0000 -2.0000