数值微分_47259.pptVIP

数学模型与数学建模方法 Slide * 第八章 数值微分与积分 1 数值微分 数值微分是用离散方法近似地计算函数在某点的 导数值。 在实验和工程应用中，常常得到的离散形式 的数据， 这时，数值微分显得尤为重要。 若 在 x=a 可导， 根据导数的定义，可以用差 商近似代替微商（导数），有 (1) (2) （h0且足够小）分别称为向前差商和向后差商。 两式平均得 (3) 称为中心差商。 中心差商精度较高。 和差商一样，称 分别为向前差分、向后差分、中心差分。 高阶导数也可用差商法求得， 例如二阶导数公式为 (4) 2 数值微分的MATLAB命令 diff(x) 向量x的向前差分；% diff(x)/h 即为计 算数值导数 diff(x, n) 向量x的n阶向前差分； 例1： ，求函数 的数值导数，并画出 的图像。 解：f=inline(sqrt(x.^3+2*x.^2-x+12)+(x+5).^(1/6)+5*x+2); x=-3:0.01:3; df=diff(f([x,3.01]))/0.01; plot(x,df) 除差商公式外，还有三点公式等其它方法计算 数值导数，具体可参考《数值分析》教材。 3 数值积分 在对已知求积分时，理论上可以利用牛顿－莱布 尼兹公式求解， 但在工程中却不实用， 到的大多数函数都不能找到其积分函数， 因为实际中遇 有些函数的 表达式非常复杂，用牛顿－莱布尼兹公式求解会变得 很繁琐，给计算精确带来不变。 数值积分便显得方便而精确。 在对函数求积分时用 3.1 定积分 用数值方法近似地求一个函数 在区间[a, b]上 的定积分，其基本思路可归结为定积分的定义 (5) 其中 显然， 取值不同，数值积分 结果就不同。 的 而各种不同的数值方法就在于，研究选 择什么样的简单函数代替 ，使得既能保证一定 精度，计算量又小。 在每个小区间上用常数函数代替 ，分别取小 区间的左端点和右端点处的函数值为常数函数的值， 得到如下的矩形公式： (6) (7) 从几何意义上说，是在小区间上和矩形面积近似 代替相应曲边梯形面积。 若用分段线性函数作为 的近似， 上用梯形面积近似代替曲边梯形面积， 即每个小区间 得到（复化） 梯形公式： (8) 为进一步提高精度，可用二次函数代替 要用到相邻两个小区间端点的三个函数值， ，每段 数目必须是偶数， 小区间的 记 ，由此可得到（复化） 辛普森公式 3.2 重积分 重积分的数值计算可通过若干次单积分的组合实现， 原则上计算定积分的梯形公式和辛普森公式都可以 推广到重积分， 但当积分域复杂时难以实现， 矩形积分域为例说明。 下面就 考察矩形区域 上的二重积分 将区间 分别 m 等分和 n 等分， 则有 其中 是下面 矩阵 Q 的元素： 上式称为二重积分的复化梯形公式。 4 数值积分的MATALB命令 sum(x) 输入数组x , 输出为x的和，可用于按矩形公式计算积分; cumsum(x) 输入数组x , 输出数组为x的依次累加和，也可用于按矩形公式计算积分; trapz(y) 输入数组y , 输出为按梯形公式计算的y的积分（单位步长）; trapz(x,y) 输入数组x ，表示积分区间；y是与x同维数的数组，表示被积函数； 输出为按梯形公式计算的y对x的积分; quad(fun,a,b) 用自适应辛普森算法计算区间 (a,b)上的积分，缺省误差为10-6. quad(fun,a,b,tol) 同上，指定误差为tol. quadl (fun,a,b,tol) 用自适应勒贝陶算法计算函数在区间 (a,b)上的积分，精度更高. dblquad(fun, xmin, xmax, ymin, ymax) 计算函数在矩形区域上的二重积分. triplequad(fun, xmin, xmax, ymin, ymax, zmin, zmax) 计算函数在长方形闭体上的三 重积分 例2：用不同方法计算 解 用以下几种方法计算： h=0.01; x=1:0.01:2; y=exp(1./x); z1=sum(y(1:100))