第7章导数与微分的ATLAB求解.pptVIP

由图1可知，二分法每一步执行的操作就是将有根区间一分为二，直至所求得的根达到所要求的精度为止，其执行流程如图2所示。 图2 二分法执行流程 3. 牛顿法及其MATLAB实现 对于方程 ，如果 是线性函数，那么它的求根是容易的。牛顿法实质上就是一种线性化方法，其基本思想是将非线性方程 逐步归结为某种线性方程来求解。 设方程 有近似根 ，将函数 在点 处展开，则有 于是，方程 就可近似地表示为 记该方程的根为 ，则 的计算公式为 上式即称为Newton迭代公式。 由牛顿迭代公式可知， 是点 处 的切线 与 轴的交点的横坐标，如图1所示。 图1 牛顿法几何意义 牛顿法执行流程比较简单，只需按如下流程即可，如图2所示。 图2 牛顿法执行流程 7.12 导数的数值求解 插值型求导公式 若已知函数在一些离散点上的函数值时，则该函数可用插值多项式来近似，然后对多项式进行微分以求得数值导数。该过程可由函数polyfit和polyder实现，相应的MATLAB函数代码如下： 上述函数中，第14行根据实验数据的个数获得插值多项式的最高阶次，第15行对实验数据进行多项式插值，第16~18行依次对多项式进行求导，第19行求解求导后的多项式在求导点处的值。 2.中心差分公式 对于给定的实验数据（或列表函数）的求导问题，最简单的方法就是利用差商来近似代替导数，比如利用向前差商近似代替导数有 类似地，也可以用向后差商近似代替导数 或用中心差商近似代替导数 由上述几个公式可知，步长 越小，计算结果越准确。但从计算角度看， 越小， ， 与 越接近，直接相减会造成有效数字的严重损失。因此需要合理的选取步长，至于步长的合理选取读者可参阅相关参考书。 上面几个公式都是两点微分公式，下面给出一个多点中心差分算法，该算法的一到四阶微分公式为 谢谢大家！ * * * 第7章 导数与微分的MATLAB求解 编者 Outline 7.1 导数概念 7.2 导数的MATLAB符号求解 7.3 函数的微分 7.4 微分中值定理 7.5 洛必达法则 7.6 泰勒公式 7.7 函数的单调性与曲线的凹凸性 7.8 函数的极值与最值 7.9 曲线的渐近线 7.10 曲率 7.11 方程的近似解 7.12 导数的数值求解 7.1 导数概念 1.导数的定义 设函数 在点 的某个邻域内有定义，当自变量 在 处取得增量 （假设点 仍在该邻域内）时，相应的函数取得增量 ；如果 与 之比当 时的极限存在，则称函数 在点 处可导，并称这个极限为函数 在点 处的导数，记为 ，即 也可记作 或 。 将上面导数的定义式中的 换为 即可得到导函数的定义式 根据函数 在点 处的导数 的定义，导数 是一个极限，而极限存在的充分必要条件是左、右极限都存在且相等，因此 存在即 在点 处可导的充分必要条件是左、右极限 及 都存在且相等。这两个极限分别称为函数 在点 处的左导数和右导数，记作 及 ，即 现在可以说，函数 在点 处可导的充分必要条件是左导数 和右导数 都存在且相等。 2.导数的几何意义 函数 在点 处的导数 在几何上表示曲线 在点 处的切线的斜率，即 其中 是切线的倾角。 如果函数 在点 处的导数