matlb第6章.pptVIP

第6章数值微积分与常微分方程求解 在许多实际问题中要采用数值方法来求函数的微分或积分。 而在实际问题中遇到的常微分方程往往很复杂，在许多情况下得不出一般解，所以，一般是要求获得解在若干个点上的近似值。 【本章学习目标】 ● 掌握微分与积分的数值计算方法。 ● 掌握常微分方程的数值求解方法。 6.1 数 值 微 分 6.1.1 数值差分与差商 6.1.2 数值微分的实现 数值微分的基本思想是先用逼近或拟合等方法将已知数据在一定范围内的近似函数求出，再用特定的方法对此近似函数进行微分。 1．用diff函数计算差分 用f(x)在点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，按行计算差分。 对于求向量的微分，函数diff计算的是向量元素间的差分，故所得输出比原向量少了一个元素。 【例6.1】设f(x)?=?sinx，用不同的方法求函数f(x)的数值导数，并在同一个坐标