第6讲 MATLAB数据分析.ppt

例: f = 1/(5+4*cos(x))’ g = int(int(diff(f,2))) e1= f-g e2=f/g [n,d]=numden(e2) simple(e1) 9.1.4 符号矩阵 符号矩阵也是一种符号表达式，所以前面介绍的符号表达式运算都可以在矩阵意义下进行。但应注意这些函数作用于符号矩阵时，是分别作用于矩阵的每一个元素。 由于符号矩阵是一个矩阵，所以符号矩阵还能进行有关矩阵的运算。MATLAB还有一些专用于符号矩阵的函数，这些函数作用于单个的数据无意义。例如 transpose(s)：返回s矩阵的转置矩阵。 determ(s)：返回s矩阵的行列式值。 其实，曾介绍过的许多应用于数值矩阵的函数，如diag、triu、tril、inv、det、rank、eig等，也可直接应用于符号矩阵。 符号微积分 1 符号极限 limit函数的调用格式为： (1) limit(f,x,a)：求符号函数f(x)的极限值。即计算当变量x趋近于常数a时，f(x)函数的极限值。 (2) limit(f,a)：求符号函数f(x)的极限值。由于没有指定符号函数f(x)的自变量，则使用该格式时，符号函数f(x)的变量为函数findsym(f)确定的默认自变量，即变量x趋近于a。 (3) limit(f)：求符号函数f(x)的极限值。符号函数f(x)的变量为函数findsym(f)确定的默认变量；没有指定变量的目标值时，系统默认变量趋近于0，即a=0的情况。 (4) limit(f,x,a,right)：求符号函数f的极限值。right表示变量x从右边趋近于a。 (5) limit(f,x,a,‘left’)：求符号函数f的极限值。‘left’表示变量x从左边趋近于a。 例9-1 求下列极限。 极限1： syms a m x; f=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/(x+a); limit(f,x,a) ans = (1/2*a*exp(sin(a))+1/2*a-exp(tan(a))+1)/a 极限2： syms x t; limit((1+2*t/x)^(3*x),x,inf) ans = exp(6*t) 极限3： syms x; f=x*(sqrt(x^2+1)-x); limit(f,x,inf,left) ans = 1/2 极限4： syms x; f=(sqrt(x)-sqrt(2)-sqrt(x-2))/sqrt(x*x-4); limit(f,x,2,right) ans = -1/2 符号导数 diff函数用于对符号表达式求导数。该函数的一般调用格式为： diff(s)：没有指定变量和导数阶数，则系统按findsym函数指示的默认变量对符号表达式s求一阶导数。 diff(s,v)：以v为自变量，对符号表达式s求一阶导数。 diff(s,n)：按findsym函数指示的默认变量对符号表达式s求n阶导数，n为正整数。 diff(s,v,n)：以v为自变量，对符号表达式s求n阶导数。 例: syms x t diff(sin(x^2)) diff(t*sin(x^2),t) 符号积分 符号积分由函数int来实现。该函数的一般调用格式为： int(s)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分。 int(s,v)：以v为自变量，对被积函数或符号表达式s求不定积分。 int(s,v,a,b)：求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf)。当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。当a,b中有一个是inf时，函数返回一个广义积分。当a,b中有一个符号表达式时，函数返回一个符号函数。 例9-3 求下列积分。 例: int(-2*x/(1+x^2)^2) int(x/(1+z^2),z) int(x*log(1+x),0,1) int(2*x, sin(t), 1) 级 数 １ 级数符号求和 求无穷级数的和需要符号表达式求和函数symsum，其调用格式为： symsum(s,v,n,m) 其中s表示一个级数的通项，是一个符号表达式。v是求和变量，v省略时使用系统的默认变量。n和m是求和的开始项和末项。 例 symsum(k^2,0,10) symsum(x^k/sym(k!), k, 0,inf) ２函数的泰勒级数 MATLAB提供了taylor函数将函数展开为幂级数，其调用格式为： taylor(f,v,n,a) 该函数将函数f按变量v展开为泰勒级数