第3章微积分问题计算机求解.PDFVIP

第3 章 微积分问题的计算机求解 Isaac Newton（1643–1727）和Gottfried Wilhelm Leibniz（1646–1716）创立的微积分学是很 多科学分支的基础。单变量与多元函数微积分、函数极限、级数求和、Taylor 级数展开、Fourier 级 数展开、常微分方程等问题直接求解是微积分学的重要内容。MATLAB 的符号运算工具箱可以 直接求解这样问题的解析解。本章3.1 节中给出基于MATLAB 符号运算工具箱中函数的单边、 多边极限问题及多元函数极限问题的求解方法，3.2 节介绍各种微分问题的计算机求解方法， 3.3 节介绍各种积分问题的解析求解方法。3.4 节将介绍给定单变量函数与多元函数的Taylor 幂 级数展开、给定函数的Fourier 级数逼近方法，并利用MATLAB 的绘图功能研究有限项拟合的 拟合效果和适用范围；还介绍一般级数的求和与求积方法等。3.5 节中将介绍的两类曲线积分和 两类曲面积分及其MATLAB 求解方法补充了微积分学的计算机求解方式，这部分内容大部分 均应该是解析求解和解析推导，属于计算机代数研究的领域，用传统的数值分析方法是不能求 解的。对不熟悉计算机代数系统开发的读者来说，用C 这样的底层语言直接进行解析解推导有 极大难度，必须使用计算机数学语言完成这类问题的分析与求解。通过这几节内容的初步学习， [1] 读者可能会发现借助计算机去求解曾令很多学生望而生畏的吉米多维奇《数学分析习题集》 中的绝大部分计算问题变得轻而易举。 在实际科学与工程研究中，微积分问题解析求解有时也面临困难。例如，若函数本身未知， 只由科学实验测出的一些实验数据，则无法用推导的方式通过数据对其代表的函数求导或求积 分，而需要通过数值的方式进行数值微积分运算。3.6 节中将单变量与多元函数的数值微积分计 算问题。在实际应用中还有很多函数积分的解析解不存在，所以需要通过数值积分的算法进行 近似。3.7 节中将介绍用数值算法求取函数积分及重积分问题的求解方法。 作为本章内容的补充，8.2 节将介绍基于样条插值的数值微积分方法；如果微积分的阶次可 以选择为非整数，还可以引入一个新的学科——分数阶微积分学。本书10.6 节将系统介绍分数 阶微积分学问题及其MATLAB 求解方法。 3.1 极限问题的解析解 应用MATLAB 语言的符号运算工具箱，可以很容易地求解极限问题、微分问题、积分问题 等微积分基本问题。利用本节和后面两节介绍的方法，读者应该能立即具备依赖MATLAB 语言 及其符号运算工具箱中提供的强大函数直接求解一般微积分运算问题的能力。本节主要侧重各 种极限问题的求解方法，包括单变量极限、单边极限和多重极限等问题。 3.1.1 单变量函数的极限 假设已知函数f (x)，则极限问题的一般描述为 L = lim f (x) （3-1-1 ） → ·54 · 高等应用数学问题的MATLAB 求解（第四版） 其物理意义是当自变量x 无限接近x0 时函数f (x) 的取值，其中，x0 可以是一个确定的值，也可 以是无穷大，例如x → ∞。对某些函数来说，还可以如下定义单边极限（或称左右极限）问题。 L1 = lim f (x), 或 L2 = lim f (x) （3-1-2 ） − → → 前者表示x 从左侧趋近于x0 点，所以又称为左极限，后者相应地称为右极限。极限问题在 MATLAB 符号运算工具箱中可以使用limit() 函数直接求出，该函数的调用格式为 L=limit(f ,x ,x0 ) ％求极限 L=