第一章Matlab基础知识.pptVIP

1.3. 符号工具箱使用简介 早期的MATLAB不具备符号运算功能，在1993年Mathworks公司从加拿大滑铁卢大学购买了Maple的使用权，并在此基础上，利用Maple函数库，开发了符号计算工具箱（Symbolic Toolbox），从此MATLAB便具备了数值计算、符号计算和图形处理三大基本功能于一身，称为在数学计算领域功能最强的语言。 MATLAB实现符号计算功能主要有三种途径： ①调用MATLAB开发的各种功能函数进行常用的符号运算，如符号表达式和符号矩阵的基本操作，符号矩阵的运算，符号微积分，符号线性方程求解，符号微分方程求解，特殊数学符号函数以及符号函数图形等。 ②通过使用函数maple.m和mpa.m与Maple接口，给一些特殊专业的人员提供方便。 ③符号函数计算器（Function Calculator）功能。 1.4.1 建立符号对象 1、建立符号变量和符号常量 Matlab提供了sym和syms两个建立符号对象的函数 1）sym函数，用来建立单个符号变量，调用格式为： 符号变量名=sym(‘符号字符串’) 注：符号字符串可以是常量、变量、函数或者表达式。 2）syms函数，依次可以定义多个符号变量，调用格式为： syms 符号变量1 符号变量2 …… 符号变量n 注意 符号变量名之间不能使用任何标点符号，只能用空格隔开。 例如： f=sym(‘cos(x)’),f=sym(‘sin(x)^2=0’); syms x; f=sin(x)+cos(x), 2、建立符号表达式 利用单引号来生成表达式 利用sym函数建立符号表达式 使用已经定义过的符号变量组成符号表达式 1.4.3 符号表达式运算 1、符号表达式的四则运算 四则运算分别使用symadd, symsub,symmul和symdiv来实现加减乘除运算，使用sympow实现幂运算。（注：6.5以后的版本中已经没有了这些函数，而直接使用+ - * /即可） 2、符号表达式的提取分子和分母 如果符号表达式是一个有理分式或可以展开为有理分式，可利用numden函数来提取符号表达式中的分子或者分母，调用格式为： [n,d]=numden(s) [n,d] = numden(sym(4/5)) returns n = 4 and d = 5. [n,d] = numden(x/y + y/x) returns n = x^2+y^2 , d = y*x 3、符号表达式的各种运算 1）factor(s)分解因式，s可以是正整数或者符号表达式。 2）expand(s)展开S中的每个因式的乘积，多用于多项式函数，三角函数，指数函数和对数函数等 3）collect(s)合并同类项 4）collect(s,v)按变量v进行合并同类项 5）compose(f,g)复合计算函数 6）符号复数的共轭函数conj（X） 7）复数的实部和虚部函数real和imag 8）设置变量的精度函数digits,digits(d)设置当前的可变算数精度的位数为整数d位，9）以控制精度的方式计算函数vpa。 示例：z=1.0e-16; x=1.0e+2; digits(14); y=vpa(1+x*z);digits(15);y=vpa(x*z+1); 10）搜索符号表达式的最简形式函数[r,how]=simple(s)，返回最短的一个 符号表达式的化简函数simplify，化简符号矩阵中的每一个元素。 示例：R1=simplify(sin(x)^4+cos(x)^3);R2=simplify(exp(c*log(sqrt(a+b)))); S=[(x^2+5*x+6/(x+2),sqrt(16))]; R3=simplify(S) 11)矩阵维数的计算函数[m,n]=size(X),m=size(X,dim),dim=1/2分别代表行和列方向。 12）pretty函数 显示数学书写格式的表达式 13）符号表达式与数值表达式之间的转换 （1）利用函数sym可以将数值表达式表示成符号表达式 （2）numeric或者eval函数可以将符号表达式转换成数值表达式 14）多项式的表示 （1）特征多项式函数poly,调用方式：p=poly(A)，若A为数值阵列，返回矩阵A的特征多项式的系数，poly(sym(A))近似等于poly2sym(poly(A)),如果A是符号矩阵，则返回A的变量为x的特征多项式，p=poly(A,v)，返回矩阵A的变量为v的特征多项式。 示例：syms z;A=hilb(4);p=poly(A);q=poly(sym(A));s=poly(sym(A),z) （2）将多项式系数向量转化为带符号变量的多项式系数poly2sym （3）从符号表达式中找出符号变量finds