[工学]matlab 第七-八章.pptVIP

第七章 Matlab数值积分与微分 7.1 Matlab数值积分 7.2 Matlab数值微分 7.1 数值积分 数值积分基本原理 求解定积分的数值方法多种多样： 简单的梯形法、 辛普生(Simpson)法、 牛顿－柯特斯(Newton-Cotes)法等 它们的基本思想都是将整个积分区间[a,b]分成n 个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b. 这样求定积分问题就分解为求和问题。 7.1.2 数值积分的实现方法 变步长辛普生法 Matlab给出了quad、quadl函数来求定积分. 该函数的调用格式为： [I,n]=quad(fname,a,b,tol,trace) [I,n]=quadl(fname,a,b,tol,trace) fname是被积函数名。 a和b分别是定积分的下限和上限。 tol用来控制积分精度，缺省时取tol=10-6。 trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。 返回参数I即定积分值, n为被积函数的调用次数. 例7-1 求定积分。 (1) 建立被积函数文件fesin.m。 function f=fesin(x) f=exp(-0.5*x).*sin(x+pi/6); (2) 调用数值积分函数quad求定积分。 [S,