第8章 MATLAB符号计算.pptVIP

MATLAB程序设计教程（第二版） 刘卫国 主编 中国水利水电出版社 所谓符号计算是指在运算时,无须事先对变量赋值,而将所得到结果以标准的符号形式来表示。 MathWorks公司以Maple的内核作为符号计算引擎（Engine），依赖Maple已有的函数库，开发了实现符号计算的两个工具箱：基本符号工具箱和扩展符号工具箱。 第8章 MATLAB符号运算 8.1 符号对象 8.1.1 建立符号对象 定义符号变量 参与符号运算的对象可以是符号变量、符号表达式或符号矩阵。符号变量要先定义，后引用。可以用sym函数、syms函数将运算量定义为符号型数据。引用符号运算函数时，用户可以指定函数执行过程中的变量参数；若用户没有指定变量参数，则使用findsym函数默认的变量作为函数的变量参数。 1．建立符号变量和符号常量 MATLAB提供了两个建立符号对象的函数：sym和syms，两个函数的用法不同。 (1) sym函数 sym函数用来建立单个符号量，一般调用格式为： 符号量名=sym(符号字符串) 该函数可以建立一个符号量，符号字符串可以是常量、变量、函数或表达式。 应用sym函数还可以定义符号常量，使用符号常量进行代数运算时和数值常量进行的运算不同。下面的命令用于比较符号常量与数值常量在代数运算时的差别。 a=sym(a); %定义符号变量a w=a*a*a %符号运算 w=a^3 %符号运算结果 x=3; v=x*x*x v=27 whos %查看变量 作符号计算： a,b,x,y均为符号运算量。在符号运算前，应先将a,b,x,y定义为符号运算量 a=sym(‘a’); %定义‘a’为符号运算量，输出变量名为a y =2/bb=sym(‘b’); x=sym(‘x’); y=sym(‘y’); [x,y]=solve(a*x-b*y-1,a*x+b*y-5,x,y) %以a,b为符号常数，x,y为符号变量 %即可得到方程组的解： x =3/a y =2/b (2) syms函数 函数sym一次只能定义一个符号变量，使用不方便。MATLAB提供了另一个函数syms，一次可以定义多个符号变量。 syms函数的一般调用格式为： syms 符号变量名1 符号变量名2 … 符号变量名n 用这种格式定义符号变量时不要在变量名上加字符串分界符(‘)，变量间用空格而不要用逗号分隔。 syms a b c d 2．建立符号表达式 含有符号对象的表达式称为符号表达式。建立符号表达式有以下3种方法： (1)利用单引号来生成符号表达式。 (2)用sym函数建立符号表达式。 (3) 使用已经定义的符号变量组成符号表达式。 y=1/sqrt(2*x) %生成1/sqrt(2*x) x要再定义 U=sym(a+b*c) %此时不需要定义变量a b c sym(x); v=1/sqrt(2*x) 8.1.2 符号表达式运算 1．符号表达式的四则运算 符号表达式的四则运算与数值运算一样，用+、-、*、/、^ 运算符实现，其运算结果依然是一个符号表达式。 syms x y z; f=2*x+3/y+z; g=4/z+3*y; f+g 2．符号表达式的提取分子和分母运算 如果符号表达式是一个有理分式或可以展开为有理分式，可利用numden函数来提取符号表达式中的分子或分母。其一般调用格式为： [n,d]=numden(s) 该函数提取符号表达式s的分子和分母，分别将它们存放在n与d中。 a=sym(0.333) a=333/1000 [n,d]=numden(a) 如果是符号矩阵numden返回两个新矩阵n、d 3．符号表达式的因式分解与展开 MATLAB提供了符号表达式的因式分解与展开的函数，函数的调用格式为： factor(s)：对符号表达式s分解因式。 expand(s)：对符号表达式s进行展开。 collect(s)：对符号表达式s合并同类项。 collect(s,v)：对符号表达式s按变量v合并同类项。 4．符号表达式系数的提取 如果符号表达式是一个多项式，可利用coeffs函数来提取符号表达式中的系数。其一般调用格式为 c=coeffs(s[,x]) 该函数返回多项式中按指定变量升幂顺序排列的系数，若没有指定变量，则返回所有项的常系数，且按离字符“x”近原则确定主变量。 5．符号表达式的化简 MATLAB提供的对符号表达式化简的函数有： simplif