数学积分与微分-MATLAB学习之三课程方案.pptVIP

MATLAB数值积分与微分; ;1.2 数值积分的实现方法 1．变步长辛普生法 基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： [I,n]=quad(filename,a,b,tol,trace) 其中filename是被积函数名。a和b分别是定积分的下限和上限。tol用来控制积分精度，缺省时取tol=10-6。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数的调用次数。; 例1 求定积分 解：(1) 建立被积函数文件fesin.m： function f=fesin(x) f=exp(-0.5*x).*sin(x+pi/6); (2) 调用数值积分函数quad求定积分： [S,n]=quad(fesin,0,3*pi) S = 0.9008 n = 77;2．被积函数由一个表格定义 在MATLAB中，对由表格形式定义的函数关系的求定积分问题用trapz(X,Y)函数。其中向量X,Y定义函数关系Y=f(X)。 例2 用trapz函数计算定积分 命令如下： X=1:0.01:2.5; Y=exp(-X); %生成函数关系数据向量 trapz(X,Y) ans = 0.28579682416393;1.3 二重定积分的数值求解 使用MATLAB提供的dblquad函数就可以直接求出二重定积分的数值解。该函数的调用格式为： I=dblquad(f(x,y),a,b,c,d,tol,trace) 该函数求f(x,y)在[a,b]×[c,d]区域上的二重定积分。参数tol，trace的用法与函数quad完全相同。;例3 计算二重定积分 解：(1) 建立一个函数文件fxy.m： function f=fxy(x,y) global ki; ki=ki+1; %ki用于统计被积函数的调用次数 f=exp(-x.^2/2).*sin(x.^2+y); (2) 调用dblquad函数求解： global ki;ki=0; I=dblquad(fxy,-2,2,-1,1) ki I = 1.57449318974494 ki = 1038;2.数值微分 ;称 分别为函数在x点处以h(h0)为步长的向前差分、向后差分和中心差分。 称 分别为函数在x点处以h(h0)为步长的向前差商、向后差商和中心差商。 当步长h(h0)充分小时，函数f在点x的微分接近于函数在该点的任意种差分。 f在点x的导数接近于函数在该点的任意种差商。;2.2 数值微分的实现 在MATLAB中，没有直接提供求数值导数的函数，只有计算向前差分的函数diff，其调用格式为： DX=diff(X) 计算向量X的向前差分，DX(i)=X(i+1)-X(i)，i=1,2,…,n-1。 DX=diff(X,n) 计算X的n阶向前差分。例如，diff(X,2)=diff(diff(X))。 DX=diff(A,n,dim) 计算矩阵A的n阶差分，dim=1时(缺省状态)，按列计算差分；dim=2，按行计算差分。 ;例 4 生成以向量V=[1,2,3,4,5,6]为基础的范得蒙矩阵，按列进行差分运算。 命令如下： V=vander(1:6) DV=diff(V) %计算V的一阶差分 DV = 31 15 7 3 1 0 211 65 19 5 1 0 781 175 37 7 1 0 2101 369 61 9 1 0 4651 671 91 11 1 0;例5 设 用不同的方法求函数f(x)的数值导数，并在同一个坐标系中做出f (x)的图像。 用三种方法： 1、用一个5次多项式p(x)拟合函数f(x)，并对p(x)求一般意义下的导数dp(x)，求出dp(x)在假设点的值； 2、直接求f(x)在假设点的数值导数； 3、求出g(x)=f (x)的表达式 然后求f (x)在假设点的数值导数。;程序如下： f=inline(sqrt(x.^3+2*x.^2-x+12)+(x+5).^(1/6)+