第6章_MATLAB数值计算_part2.pptVIP

第6章 MATLAB数值计算 6.1 数据处理与多项式计算 6.2 数值微积分 6.3 线性方程组求解 6.4 最优化问题求解 6.5 常微分方程的数值求解 6.2 数值微积分 6.2.1 数值微分（导数） 不关心微分的形式和性质，只关心该微分在一串离散点的近似值以及所计算的近似值有多大的误差。 MATLAB下求数值导数的两种方法： 1）用多项式对任意函数进行拟合，再利用多项式求导函数进行求导，从而得出某点的导数值； 2）用向前差商的方法求任意函数在某点的导数值。 用向前差分方法求数值微分 (1) 在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，按行计算差分。 用向前差商方法求数值导数 向前差分除以步长：diff(x)/h 称为函数在所求的点处以h(h0)为步长的向前差商。 只要步长足够小，向前差分即无限接近于函数在该点的微分，向前差商的结果无限接近于函数在该点的数值导数。 还可以用向后差分、中