数值积分与数微分matlab.pptVIP

MATLAB数值积分与数值微分;*;*;为什么要作数值积分;数 值 积 分 的 基 本 思 路;数值积分基本思想;三个求积分公式;矩形法数字积分的演示程序rsums;y=f(x);2.辛普森(Simpson)公式 （抛物线公式）;数值积分的实现方法;用MATLAB 作数值积分;例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;　　　也可不建立关于被积函数的函数文件，而使用匿名求解， 　命令如下： g=@(x)(exp(-0.5*x).*sin(x+pi/6)); %定义语句函数 [S,n]=quad(g,0,3*pi) %不用加引号 　S = 　 0.9008 　n = 　77 ;例2 分别用quad函数和quadl函数求定积分的近似值，并在相同的积分精度下，比较函数的调用次数。 1.调用函数quad求定积分： clc,clear format long; fx=@(x)(exp(-x)); [I,n]=quad(fx,1,2.5,1e-10) I = 0.28579444254766 n = 65;2.调用函数quadl求定积分： clc,clear format long; fx=@(x)(exp(-x)); [I,n]=quad(fx,1,2.5,1e-10) I = 0.28579444254881 n = 18;梯形积分法 在科学实验和工程应用中，函数关系往往是不知道的，只有实验测定的一组样本点和样本值，这时，人们就无法使用quad等函数计算其定积分。在MATLAB中，对由表格形式定义的函数关系的求定积分问题用梯形积分函数trapz。;被积函数由一个表格定义—梯形积分法 　 在MATLAB中，对由表格形式定义的函数关系的求定积分问题用trapz(X,Y)函数。其中向量X,Y定义函数关系Y=f(X)。 例4 用trapz函数计算定积分。 命令如下： clc,clear X=1:0.01:2.5; Y=exp(-X); %生成函数关系数据向量 trapz(X,Y) ans = 0.28579682416393;数值积分;数值积分;数值积分;;数值积分;;矩形域上计算二重积分的命令：;例5 计算二重定积分 (1) 建立一个函数文件fxy.m： function f=fxy(x,y) f=exp(-x.^2/2).*sin(x.^2+y); (2) 调用dblquad函数求解 I=dblquad(fxy,-2,2,-1,1) I =1.57449318974494 ;clear f=@(x,y)(exp(-x.^2/2).*sin(x.^2+y)); I=dblquad(f,-2,2,-1,1) I = 1.5745;例6 计算三重定积分命令如下：fxyz=@(x,y,z)(4*x.*z.*exp(-z.*z.*y-x.*x));triplequad(fxyz,0,pi,0,pi,0,1,1e-7)ans= 1.7328;　　　　　2 数值微分 2.1 数值差分与差商 任意函数f(x)在x点的导数是通过极限定义的：;　　　如果去掉上述等式右端的h→0的极限过程，并引进记号： 　 分别称为函数在x点处以h(h0)为步长的向前差分、向后差分和中心差分。 　;　　当步长h充分小时，有 　 分别称为函数在x点处以h(h0)为步长的向前差商、向后差商和中心差商。当步长h(h0)充分小时，函数f在点x的微分接近于函数在该点的任意一种差分，而f在点x的导数接近于函数在该点的任意一种差商。; 在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，按行计算差分。;;例8