数值微分与数值积分.pptxVIP

2025/6/61第5章

数值微分与数值积分

微积分是高等数学中的重要内容，在化学工程上有许多非常重要的应用微积分的数值方法，不同于高等数学中的解析方法，尤其适合求解没有或很难求出微分或积分表达式的实际化工问题的计算，例如：列表函数求微分或积分引言1---数值微积分方法不同于解析方法

数值微分和数值积分与插值和拟合往往是密不可分的如在进行数值微分时，常针对离散的数据点，利用插值和拟合可以减少误差；而数值积分的基本思路也来自于插值法。例如如果所积函数的形式比较复杂或以表格形式给出，则可通过构造一个插值多项式来代替原函数，从而使问题简化1引言2---数值微分和数值积分与插值和拟合关系密切2

§5.1数值微分42025/6/6化工领域的实际问题中时常需要求列表函数在节点和非节点处的导数值，这是数值微分所要解决的问题。数值微分方法可近似求出某点的导数值例如在反应动力学的研究中，根据实验数据确定反应的动力学方程：这里实验测得一批离散点，要计算只能借助数值微分求导解决表5-1反应动力学实验数据0导入

0.1建立数值微分公式的三种思路52025/6/6常用三种思路建立数值微分公式：从微分定义出发，通过近似处理，得到数值微分的近似公式从插值近似公式出发，对插值公式的近似求导可得到数值微分的近似公式先用最小二乘拟合方法根据已知数据得近似函数，再对此近似函数求微分可得到数值微分的近似公式。然后对各方法数值微分后得到的多项式求值，即可求出任意点处的任意阶微分

方法概述在微积分中，一阶微分的计算可以在二相邻点x+h和x间函数取下列极限求得：注：高阶微分项可以利用低阶微分项来计算，如二阶微分式可以表示为：5.1.1差分近似微分取其达到极限前的形式，就得到以下微分的差分近似式：对应的差分式有：上式中三种不同表示形式依次是一阶前向差分、一阶后向差分和一阶中心差分来近似表示微分。其中一阶中心差分的精度较高。

2差分的Matlab实现72025/6/6在Matlab中，可用diff函数进行离散数据的近似求导01调用形式：Y=diff(X,n)02其中：X表示求导变量，可以是向量或矩阵。如是矩阵形式则按各列作差分；n表示n阶差分，即差分n次Y是X的差分结果03注：用diff函数进行离散数据的近似求导与前向差分近似，但误差较大。最好将数据利用插值或拟合得到多项式，然后对近似多项式进行微分04

例5.1：丁二烯的气相二聚反应如下：实验在一定容器的反应器中进行，3260C时，测得物系中丁二烯的分压(mmHg)与时间的关系如表5-2所示。。表5-2丁二烯二聚反应实验数据用数值微分法计算所列时刻每一瞬间的反应速率

92025/6/6解：程序如下：t=[0:5:90];pA=[632.0590.0552.0515.0485.0458.0435.0414.0396.0378.0362.0348.0336.0325.0314.0304.0294.0284.0274.0];dt=diff(t);％求时间t的差分dpA=diff(pA);％求压力的差分q=dpA./dt％q为数值微分结果执行结果：q=Columns1through8-8.4000-7.6000-7.4000-6.0000-5.4000-4.6000-4.2000-3.6000Columns9through16-3.6000-3.2000-2.8000-2.4000-2.2000-2.2000-2.0000-2.0000Columns17through18-2.0000-2.0000

§5.1.2三次样条插值函数求微分若三次样条插值函数收敛于，那么导数收敛于，因此用样条插值函数作为函数，不但彼此的函数值非常接近，而且导数值也很接近。用三次样条插值函数求数值导数是可靠的，这是化工计算中求数值微分的有效方法。的近似

方法概述用三次样条插值函数建立的数值微分公式为：01求导得：02（5-2）03式（5-1）和（5-2）不但适用于求节点处的导数，而且可求非节点处的导数。04（5-1）05其中，06

122025/6/62三次样条插值函数求微分的Matlab函数Matlab求离散数据的三次样条插值函数微分方法分三个步骤：Step1:对离散数据用csapi函数（或spline函数）得到其三次样条插值函数调用形式pp=csapi(x,y)其中：x,y分别为离散数据对的自变量和因变量；