MATLAB第3章 高等数学运算（1）.pptVIP

3.1函数与极限 3.1.1函数运算 3.1.2极限运算 3.1.3级数的求和 3.2方程的求解 3.2.1线性方程(组)求解 3.2.2非线性方程求解 3.2.3求函数的区间的零点（根） 3.2.4区间的函数极值 3.3导数与微分 3.3.1导数运算 3.3.2导数的应用 3.4积分及其应用 3.4.1不定积分 3.4.2定积分 3.4.3曲线与曲面积分 第3章 高等数学运算（1） 3.1 函数与极限 本节介绍利用MATLAB对函数和函数极限进行运算与分析的基本方法，包括： 函数的四则运算 函数的极限运算 级数（函数）的求和与收敛分析等。 如前所述，MATLAB对函数有符号型函数和数值型函数两种表示形式。符号型函数可用于对未知变量函数进行运算，数值型函数可用于对赋值变量函数进行运算(章节1.2)。用户可直接利用MATLAB提供的如三角函数等基本函数命令进行函数表达式的编辑和运算(章节2.2)，也可直接调用MATLAB的内定函数和自建函数进行求解(章节2.3)。 3.1.1函数运算 1.函数的加减运算 例3-1 已知 试求1）h(x)=f(x)+g(x)；2）给定变量a4=8;a3=2;a1=9;a0=5;b1=6;b0=10，求h(x)。 解：先定义符号变量及函数求符号多项式函数的和运算；再给系数变量赋值求常系数多项式函数表达式。 2. 函数的乘法运算 事实上是利用了向量卷积(convpolytion)的函数命令格式： hp=conv(fp,gp) 3.1.1函数运算 3. 函数的除法运算 与乘法运算类似，对于符号函数可以利用除法运算符进行函数的运算，对于数值多项式函数可以利用反卷积函数命令格式进行运算： [z,r]=deconv(gp, fp) 3.1.1函数运算 3.1.2极限运算 1. 函数参数的极限 2. 函数的极限 3. 多元函数的极限 注意：在极限不存时, MATLAB可能给出非数NaN或者无结果的结构原式, 但有时也不能就此说明极限是否不存在, 应进一步观察其左右极限情况。如, 设F1=x^n·lnx，在x→0时，F1是0·∞型的未定式，若使用洛必达法则，可等价为G1=-x^n /n，在x→0时有G1=0，即。但用MATLAB直接运算，则不会得到上述结果。 3.1.3 级数的求和 级数可分为常数项级数和函数项级数。 MATLAB所提供的级数求和的函数命令的运算中已将收敛性问题考虑其中，可以由求解结果看到级数的收敛情况，所以级数的收敛性是和极限的概念联系在一起的。MATLAB提供了无穷级数求和的函数命令, 其格式为： symsum(f,v,a,b) 其中，f为级数的通项表达式；v是求和变量；a, b为求和的下限和上限值，缺省时a=0，b=v-1。 对于有限项级数的求和，也可利用命令sum()及数组点运算格式来实现求和。 3.2 方程的求解 求解方程是数学分析的基础内容之一，也是MATLAB符号函数运算的一个重要功能。方程求解分为一元、多元线性方程和非线性方程的求解问题。本节结合高等数学教程中的相关内容介绍方程求解的一般方法, 其它方法将在线性代数中进行介绍。 3.2.1 线性方程(组)求解 命令格式： X=solve(‘eq1’,’eq2’,...,’eqn’,val1,val2,...,valn) %解符号方程、方程组、超越方程 X=linsolve(A,B) %解线性方程组AX=B的解X，中A、B为符号矩阵 X=A\B %解线性方程组AX=B的解X, 求解速度快、精度高、适应性较强 X=inv(A)*B %解线性方程组AX=B的解X。由于计算速度慢，一般不采用。 * * 例3-20,3-21： 说明：命令solve()和linsolve()是一种符号型函数的求解运算，所以给出的结果是符号解的数组结构的形式解。因此例3-21的两种方法所得结果是有区别的，即方法1给出的是符号解的形式解，方法2给出的是解析解。实际上方法2是线性代数的求解方法。另外，还可以看出X=inv(A)*B比X=A\B所需的运算时间要长近8倍或者X=A\B的运行效率提高近8倍。 3.2.1 线性方程(组)求解 3.2.2 非线性方程求解 命令格式: X=fsolve(functoin name,x0,options) 其中，functoin name: 非线性方程组的函数名； x0: 初始值向量或矩阵； options: 设置命令函数fsolve求解过程的参数选项，通常设