讲MATLAB数值计算.pptVIP

* 小结： 一、数据统计与分析 二、多项式计算 三、数值微分与积分 * 4.3 数值微积分 6.2.1 数值微分 1、数值差分与差商 任意函数f(x)在x点的导数定义为： * 以上三式分别称为函数f(x)在x点处以h 为步长的向前差分，向后差分和中心差分。 当步长h充分小时，函数f在点x处的微分接近函数在该点的任意种差分。 定义差分： * 当步长h充分小时，有 分别称为函数在x点处以h为步长的向前差商，向后差商和中心差商。 当步长h充分小时，函数f在点x处的导数接近函数在该点的任意种差商。 * 2、数值微分的实现 MATLAB中，没有直接提供求数值导数的函数，只有计算向前差分的函数diff。可用差商近似代替导数。 DX=diff(X)：计算向量X的向前差分， DX(i)=X(i+1)-X(i)，0in。 DX=diff(X,n)： 计算向量X的n阶向前差分， 例如： diff(X,2)=diff(diff(X))。 * DX=diff(A,n,dim) ： 计算矩阵A的n阶差分， dim=1时(缺省状态)，按列 计算差分， dim=2，按行计算差分。 * 例4.18 设X由[0，2π]间均匀分布的10个点组成，求向量sin(X)的1~3阶差分。 X=linspace(0,2*pi,10) Y=sin(X) DY=diff(Y) %计算Y的一阶差分 D2Y=diff(Y,2) %计算Y的二阶差分， 也可用命令diff(DY)计算 D3Y=diff(Y,3) %计算Y的三阶差分， 也可用diff(D2Y) 或diff(DY,2) * close all h=pi/20;t=0:h:2*pi; y=sin(t); DY=diff(y);DYy=DY/h; %一阶差分和差商 DYY=diff(y,2);DYYy=DYY/h/h; %二阶差分和差商 subplot(3,1,1),stem(t,y,.),title(sin(t)), axis([-1,7,-1.1,1.1]) subplot(3,1,2),stem(t(1:end-1),DYy,.),title(sin(t)), axis([-1,7,-1.1,1.1]) subplot(3,1,3),stem(t(1:end-2),DYYy,.),title(sin(t)), axis([-1,7,-1.1,1.1]) 例如：求sin(t)的一阶、二阶导数。 数值微分的实现：用差商代替导数。 * * 4.3.2数值积分 1、数值积分的基本原理 即， 当ε充分小时，可用I2近似代替I1。 也即，可用p(x)代替f(x)。 选p(x)的原则：选f(x)的插值多项式 * 2、数值积分的实现 (1)被积函数是一个解析式 quad(filename,a,b,tol,trace)： %用Simpson公式计算 quadl(filename,a,b,tol,trace)： %用Lobatto公式计算 filename是被积函数名，a,b为定积分的上、下限，tol是计算精度，缺省值是10-6。trace非0时，展现积分过程,缺省值为0。 * 例4.20 用两种不同的方法求积分。 注意，调用quad 和 quadl 函数时，先要建立一个描述被积函数f(x)的函数文件或语句函数。 * 在MATLAB命令窗口直接输入命令： I=quad(exp(-x.^2),0,1) 或 I=quadl(exp(-x.^2),0,1) 可以完成积分。 * (2)被积函数由一个表格定义 MATLAB中，对由表格形式定义的函数关系的求定积分问题用trapz(X,Y)函数。其中等长向量X、Y定义函数关系Y=f(X)。 I=quadl(sin(x),0,pi) I = 2.0000 例如： * x=[0,pi/4,pi/2,3*pi/4,pi]; y=[0,0.707,1,0.707,0];trapz(x,y) ans =