matlb微分与积分.pptVIP

数值微积分以及数值分析 数值微分 数值微分的实现 两种方式计算函数f(x)在给定点的数值导数：1.用多项式或者样条函数 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，按行计算差分。 例子：help diff 拉普拉斯微分算子 MATLAB中的离散拉普拉斯微分算子调用格式dell2(U) Help dell2 数值积分 数值积分基本原理 求解定积分的数值方法多种多样，如简单的梯形法、辛普生(Simpson)法、牛顿－柯特斯(Newton-Cotes)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。这样求定积分问题就分解为求和问题。 数值积分的实现方法 低阶法-自适应递推辛普生法 基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： I=qua