微积分线性代数知识在matlab中的实现.docVIP

第一章 微积分、线性代数的基础知识及其在Matlab中的实现 1．1 建模中常用微积分基础知识在Matlab中的实现 在数摸问题的计算中，经常要用到极限、导数、微分和积分等基本运算。利用数学软件Matlab可以使复杂的微积分运算变得很容易。本节介绍微积分的基本运算在计算机上如何实现，包括：微分、积分、极限、级数求和、方程求根以及Taylor展开等功能的实现方法。 下面是MATLAB中的基本符号运算函数（有关符号运算请参考附录A）： 1. x=sym(‘x’) 功能：创建一个符号变量x 2. syms x y z 功能：创建多个符号变量 x , y , z 3. r=collect(S,v) 功能：合并同类项，S是符号表达式,v是变量或表达式，R是合并同类项后的结果。 4. factor(S) 功能：符号计算的因式分解，X是待分解的符号多项式 5. expand(S) 功能：对符号多项式或函数S进行展开 6. r=simple(S) 或 r=simplify(S) 功能：对符号表达式S进行化简 7. subs(S,old,new) 功能：把符号变量中的变量old 用 new 代替，new 可以是一个符号，也可以是具体的数 8. vpa(S) 功能：对符号表达式S计算其任意精度的数值 9. eval(S) 功能：计算符号表达式（或字符串）S 1．1．1 导数、极值和积分、Talyor公式及在Matlab中的实现 一、极限运算 在MATLAB中，计算极限采用下面的命令： 命令 功能 limit(f,x,a) 计算 limit(f,x,inf) 计算 limit(f,x,a,’right’) 计算单侧极限 limit(f,x,a,’left’) 计算单侧极限 注意：在左右极限不相等或左右极限有一个不存在时，Matlab的默认状态是求右极限 例1 求极限与极限 解 Matlab命令为 syms x; y1=(1+4*x)^(1/x); y2=(exp(x)-1)/x; limit(y1,x,0) ans= exp(4) limit(y2,x,0) ans= 1 例2 求极限 解 syms x; y=sqrt(x)-2^(-1/x); limit(y,x,0,’right’) ans= 0 二、求导运算 导数是函数增量与自变量增量之比的极限（），即。在Matlab中求函数的导数及其他一些类似计算均由 diff 命令来完成。 1．一元函数的求导 命令形式1： diff(f) 功能：求函数f的一阶导数，其中f为符号函数。 命令形式2： diff(f,n) 功能：求函数f的n阶导数，其中f为符号函数。 例3 求函数 的二阶导数 解 Matlab的命令为 syms x ; f= 3*x^3+5*x+1 ; diff(f,2) ; ans= 18*x 例4 设，求 解 syms x ; y=3*x^2-2*x+1 ; B=diff(y) , x=1 ; eval(B) 运行结果为 B = 6*x – 2 ans = 4 2．多元函数的偏导数 将函数中的变量看成常量而对变量的导数称为二元函数对变量的偏导数，记为。若把当作常量而对求导的结果称为函数对的偏导数，记为。求偏导的方法和一元函数的求导方法一样，只要把另一个变量看成常量即可。 命令形式1： diff(f,) 功能： 多函数f 对变量的一阶偏导。 命令形式2： diff(f,,n) 功能：多函数f 对变量的n阶偏导。 例5 求关于的偏导数 解 syms x y; z= x^2*sin(2*y); B=diff(z,’x’); 运行结果为 B=2*x*sin(2*y) 3. 全微分、参数方程求导及隐函数求导 (1) 若函数在某点的两个偏导数存在且连续，则函数在该点的全微分为 Matlab中求函数全微分的命令为： diff(z,x)+diff(z,y) (2) 对参数方程所确定的函数，根据公式，连续两次利用diff(f)命令就可以求出结果。 (3) 隐函数求导 方程所确定的隐函数，则其导数；方程确定的隐函数，其导数为，。在Matlab中按照上述公式，分别求出函数的偏导再相除就可得到隐函数的导数。 三、积分运算 1．一元函数的不定积分 命令形式1： int(f) 功能： 求函数f对默认变量的不定积分，用于函数中只有一个变量的情况。 命令形式2：int(f,v) 功能：求符号函数