数学实验——第三章 多元微积分.pptVIP

* * * 程序Ch2_2_2.m * 关于零点问题的更进一步讨论，有兴趣的读者可参阅同济大学计算数学教研室编写的《现代数值数学和计算》。 * * 程序yy4.m * 程序Ch2_8.m * 程序Ch2_8.m * 程序Ch1_9.m * 程序:Ch2_10 我们知道, 在定积分中有公式 是奇数. 例 计算广义积分 输入命令 计算结果 例 计算积分 输入命令 计算结果 执行命令 返回 再执行命令 返回 作数值积分: 返回值 例 用 计算积分 输入语句: 结果为 例 用上面四个函数计算积分 程序如下: 五、级数 基本格式 其中 表示级数求和中的初始项和终止项. 1.级数求和 例 求下列级数的和. ⑴ ⑵ 输入命令: 结果为 注 上面的结果都是符号解. ⑵函数展开成Taylor级数 基本格式 功能: 例 将 在点 处展开到9阶. 输入语句: 对函数 在点 处展开到 阶. 例 对函数 在点 展开到12阶. 命令如下: 的多项式 一个有意思的例子 欧拉常数的来源. 程序如下: Taylor级数逼近计算器 在MatLab中, 有一个用来图示Taylor级数逼近情况的 在命令窗口中输入: 是一个交互式的Taylor级数逼近计算器. 计算器. 输入函数 展开的阶数 例 求函数 在区间 的 阶Taylor级数. * 程序Ch2_1.m * 程序ezplot1.m * 程序ch2_3.m * 程序Ch2_4.m * 程序plot31.m * 程序:fdhs.m * 程序Ch2_5.m * 程序text1.m * 程序Ch2_6.m * 程序text2.m * 程序Ch1_7.m * 程序Ch1_7.m 程序L1.m * 程序Ch2_2_1.m * * 讨论: 二分法能求出函数在指定区间中的零点, 但其 迭代过程比较长. 计算量比较大. 一个比较好的方法叫 做Newton切线法. 切线法的意义是: 设函数 初始点为 在该点处作函数的切线: 得到该切线与轴的交点 重复该过程, 直到 并判定 是否满足精度要求. … 满足精度要求. 切线法的几何描述: 用切线法求解上面问题. 程序如下: 结果如下: 由此可见收敛速度远快于二分法. 在MatLab中函数求零点的命令是 基本格式 例 求函数 的零点. 定义函数文件: 执行命令: 2.导数 基本格式 功能 对表达式及指定变量求 阶导数. 例 求函数 的导数. 输入命令 输出结果 即: 例 求 的三阶导数并化简. 输入语句: 输出结果 化简 结果 即: 3.极值和最大最小值 MatLab中求函数极值的命令为: 和 例 求解函数 在区间 中 的极值和最小值. 先作出函数图形: 再输入命令 返回值为 执行命令 结果与上相同. 程序如下: 三、求解方程（组） 1.方程（组）的符号解 格式 例 给出一元二次方程 的公式解. 编写程序: 例 求解方程 输入语句: 结果为: 再执行命令 得到第二个根的近似值 例 求解方程组 输入语句: 结果为: 2.方程（组）的数值解 许多方程（尤其是超越方程）的解很难有精确的表示 式, 因而需要用数值的方法求得相应的近似解. 例如, 对于方程 由零点定理知该方程在区间 中可解, 若执行命令: 返回值为 进一步地 从而得到该问题的近似解. 画出该函数的图形, 可以看到根的存在范围. 为此输入语句: 得到曲线为: 局部放大: MatLab中方程求根的数值解函数为 格式1 格式2 对于格式2, 要求函数在区间端点异号. 在上例中, 继续执行命令 得到近似解 再执行: 结果相同: 例 求解方程 首先作出函数的大致图形: 由此得到函数的零点范围. 再输入下面语句: 由此得到零