第08章节MATLAB符号计算.pptVIP

第8章 MATLAB符号计算 8.1 符号对象 8.2 符号微积分 8.3 级 数 8.4 符号方程求解 8.1 符号对象的创建和使用方法 8.1.1 符号变量和表达式的创建 1．建立符号变量 MATLAB提供两种建立符号对象的方法：sym和syms。 (1) sym函数 sym函数用来建立单个符号量，格式为： 符号量名=sym(符号字符串) 该函数可以建立一个符号量，符号字符串可以是常量、变量、函数或表达式。 x=sym(‘x’) y=sym(‘y’) a=sym(s,real) a = s a*a ans = s^2 a*a*a*a/2 ans = 1/2*s^4 (2) syms命令 函数sym一次只能定义一个符号变量，使用不方便。MATLAB提供了另一个函数syms，一次可以定义多个符号变量。syms函数的一般调用格式为： syms 符号变量名1 符号变量名2 … 符号变量名n 不要在变量名上加字符串分界符(‘) 变量分隔用空格而不要用逗号分隔 注意： syms一次可以定义多个符号变量； syms命令sym函数的缩写形式； syms x y real 等价于 x=sym(‘x’,‘real’); y=sym(‘y’,‘real’) MatLab提倡使用syms命令，因为书写更方便、清楚。 x=sym(a*b*c) x = a*b*c f=x^2-x+1 f = a^2*b^2*c^2-a*b*c+1 syms x y kk y=sym(x^2) y = x^2 kk=sym(x^2+x^4+y^2) kk = x^2+2*x^4 2．建立符号表达式 含有符号对象的表达式称为符号表达式。建立符号表达式有以下2种方法： (1)用sym函数建立符号表达式。 f=sym(cos(x)+sin(x)) f = cos(x)+sin(x) (2)使用已经定义的符号变量组成符号表达式。 syms a b c y=a^2+2*b+c y = a^2+2*b+c S=sym(x^2+y)中x和y未定义。可使用： syms x y S=sym(x^2+y) S = x^2+y 或 syms x y S=x^2+y S = x^2+y 8.1.2 符号矩阵 符号矩阵也是一种符号表达式。符号表达式运算都可以在矩阵意义下进行。 符号矩阵的创建方法： syms a b c d A=[a b;c d] B=[a c;b d]; A*B ans = [ a^2+b^2, a*c+b*d] [ a*c+b*d, c^2+d^2] MATLAB还有一些专用于符号矩阵的函数，这些函数作用于单个的数据无意义。例如 transpose(s)：返回s矩阵的转置矩阵。 det (s)：返回s矩阵的行列式值。 其实，曾介绍过的许多应用于数值矩阵的函数，如diag、triu、tril、inv、rank、eig等，也可直接应用于符号矩阵。 transpose(A) ans = [ a, c] [ b, d] det(A) ans = a*d-b*c diag(A) ans = a d triu(A) ans = [ a, b] [ 0, d] tril(A) ans = [ a, 0] [ c, d] inv(A) ans = [ d/(a*d-b*c), -b/(a*d-b*c)] [ -c/(a*d-b*c), a/(a*d-b*c)] eig(A) ans = 1/2*a+1/2*d+1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2) 1/2*a+1/2*d-1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2) 8.1.3 符号表达式运算 1．符号表达式的四则运算 符号表达式的加、减、乘、除运算在低版本的MatLab中可分别由函数symadd、symsub、symmul和symdiv来实现，在6.5以后的版本中已经废弃，可直接用运算符运算。 A*B ans = [ a^2+b^2, a*c+b*d] [ a*c+b*d, c^2+d^2] A/B ans = [ (-b^2+a*d)/(a*d-b*c), a*(b-c)/(a*d-b*c)] [ -d*(b-c)/(a*d-b*c), (a*d-c^2)/(a*d-b*c)] 2．符号表达式的因式分解与展开 MATLAB提供了符号表达式的因式分解与展开的函数，函数的调用格式为： factor(s)：对符号表达式s分解因式。 expand(s)：对符号表达式s进行展开。 horner：将符号表达式转化为嵌套形式 collect(s)：对符号表达式