第3讲matlab的符号运算..docVIP

第三讲 MATLAB的符号运算 （注：文中红色字体为命令执行的结果，在Command窗口中显示） 3-1 符号对象的创建和使用 1．符号运算入门 符号运算的特点是，运算过程中允许存在非数值的符号变量。先看如下示例： 函数，用MATLAB求它的微积分，命令如下： f’sinx^2’;%定义符号函数fx dfdxdifff%求的指令 intfintf%求的指令 显示的计算结果为： dfdx2*sinx*cosx intf-1/2sinx*cosx+1/2*x 所以，，。此例中，首先定义符号函数f’sinx^2’，然后由符号运算获得的微分和积分。 2．定义符号变量 在使用符号变量之前，应先声明某些要用到的变量是“符号”变量。声明符号变量的语句： syms 变量名列表 或：sym‘变量名’ 其中各个变量名应该用空格分隔，而不能用逗号分隔。如创建符号变量x和a： xsym‘x’ asym‘alpha’ 或用： syms x a %定义符号变量x和a 这里，变量x和a的类型是符号对象，它们被定义后，即可参与符号运算。 3．定义符号表达式和符号方程 符号表达式和符号方程是两种不同的操作对象。区别在于：符号表达式不包含等号（），而符号方程须带等号。它们的创建方式相同。 如：要考虑二次函数fax^2+bx+c，可以创建符号表达式，赋值给符号变量f。 fsym‘a*x^2+b*x+c’ 或： f‘a*x^2+b*x+c’ 此例中，将符号表达式赋给符号变量f，但这不是必需的，引入符号变量是为了以后调用方便。 在这种情况下，没有创建对应于表达式中a、b、c、x项的变量，为了执行符号数学运算（如微分、积分等），必须显式地创建这些变量，可用下列命令创建： syms a b c x 如下例中创建了符号表达式和符号方程，分别赋给相应的符号对象。 syms x a b c f’sinx^2’;%创建符号表达式sinx^2赋给变量f eq‘a*x^2+b*x+c0’ %创建的符号方程赋给变量eq 4．定义抽象函数和符号数学函数 若要创建抽象函数fx，可使用： fsym‘fx’ 则f就象fx一样参与运算。 在Symbolic Math Toolbox中可利用符号表达式创建符号数学函数。如下面的命令将产生符号表达式r、t、f： syms x y z rsqrtx^2+y^2+z^2 %或rsym‘sqrtx^2+y^2+z^2’ tatany/x fsinx*y/x*y 3-2数值与符号的转换 Sym函数有四种选项将数值结果转换为符号表达式。 如对于变量rho，定义为：，在MATLAB命令窗口中定义rho: rho1+sqrt5/2 将这一数值结果转化为符号表达式： （1）symrho,’f’ 返回符号浮点表示形式，结果为： symrho,’f’ ans ‘1.9e3779b97f4a8’*2^0 （2）symrho,’r’ 返回符号有理数表示形式，这是sym的默认设置，当调用sym而没有第二个参数时，相当于sym使用了r选项。结果为： symrho,’r’ ans 座机电话号码*2^-52 例如： sym0.75 ans3/4 （3）symrho,’e’ 返回符号有理数表示形式，同时根据eps给出rho的理论表达式和实际计算的差。 symrho,’e’ ans座机电话号码*2^-52 这里，rho的理论值和实际浮点值相同，对于1/3，可得： sym1/3,’e’ ans1/3-eps/12 （4）symrho,’d’ 返回符号十进制小数表示形式。有效位数由digits定义。Digits的默认值是32位。 symrho,’d’ ans1.座机电话号码0座机电话号码711907 如果想要一个较短的结果，可由digits定义有效位数。 digits5 symrho,’d’ ans 1.6180 上述四种格式中，最后得到的都是符号对象变量。 3-3符号算术运算 1．定义符号矩阵 MATLAB中的数值矩阵不能直接参与符号运算，必须经过转换。Sym函数的一个非常有用的功能就是将数值矩阵转换成符号矩阵。 例如，下面的代码将产生一个3×3的Hilbert矩阵。 Ahilb3 A 1.0000 0.5000 0.3333 0.5000 0.3333 0.2500 0.3333 0.2500 0.2000 对A使用sym命令可获得一个具有无穷精度的3×3的符号形式的矩。 AsymA A [ 1,1/2,1/3] [1/2,1/3,1/4] [1/3,1/4,1/5] 当调用Symbolic Math Toolbox时，默认地调用有理算术运算。当使用sym函数创建符号变量时，有理算术运算将被启动。