第3章 MATLB符号计算.ppt

第3章 MATLAB符号计算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号计算是可以对未赋值的符号对象(可以是常数、变量、表达式)进行运算和处理。MATLAB具有符号数学工具箱(Symbolic Math Toolbox)，将符号运算结合到MATLAB的数值运算环境。符号数学工具箱是建立在Maple软件基础上的。 3.1 符号表达式的建立 Symbolic Math Toolbox2.1版规定在进行符号计算时，首先要定义基本的符号对象然后才能进行符号运算。 sym命令也可以把数值转换成某种格式的符号常量。 语法： sym(常量,参数) %把常量按某种格式转换为符号常量 说明：参数可以选择为’d’、’f’、’e’或’r’ 四种格式，也可省略，其作用如下表。 3.1.2 创建符号变量和表达式 创建符号变量和符号表达式可以使用sym和syms命令。 2.使用syms命令创建符号变量和符号表达式 语法： syms(‘arg1’, ‘arg2’, …,参数) %把字符变量定义为符号变量 syms arg1 arg2 …,参数　 %把字符变量定义为符号变量的简洁形式 3.1.3 符号矩阵 3.2 符号表达式的代数运算 符号运算与数值运算的区别主要有以下几点： 3.2.1符号表达式的代数运算 运算符“.*”，“./”，“.\”，“.^”分别实现符号数组的乘、除、求幂，即数组间元素与元素的运算。 在符号对象的比较中，没有“大于”、“大于等于”、“小于”、“小于等于”的概念，而只有是否“等于”的概念。 运算符“==”、“~=”分别对运算符两边的符号对象进行“相等”、“不等”的比较。当为“真”时，比较结果用1表示；当为“假”时，比较结果则用0表示。 2. 函数运算 复数的共轭conj、求实部real、求虚部imag和求模abs函数与数值计算中的使用方法相同。但注意，在符号计算中，MATLAB没有提供求相角的命令。 3.2.2 符号数值任意精度控制和运算 数值型：MATLAB的浮点运算。 有理数型：Maple的精确符号运算。 VPA型：Maple的任意精度运算。 语法： digits(n) %设定默认的精度 语法： S=vpa(s,n) %将s表示为n位有效位数的符号对象 3. Symbolic Math Toolbox中的三种运算方式的比较 程序分析： 三种运算方式中数值型运算的速度最快。 有理数型符号运算的计算时间和占用内存是最大的，产生的结果是非常准确的。 VPA型的任意精度符号运算比较灵活，可以设置任意有效精度，当保留的有效位数增加时，每次运算的时间和使用的内存也会增加。 数值型变量a1结果显示的有效位数并不是存储的有效位数，显示的有效位数由“format”命令控制。 3.2.3 符号对象与数值对象的转换 sym命令可以把数值型对象转换成有理数型符号对象，vpa命令可以将数值型对象转换为任意精度的VPA型符号对象。 如果不指定，MATLAB将基于以下原则选择一个自由变量： 如果不确定符号表达式中的自由符号变量，可以用findsym函数来自动确定。 语法： findsym(EXPR,n) %确定自由符号变量 3.3.2 符号表达式的化简 同一个数学函数的符号表达式的可以表示成三种形式 【例】三种形式的符号表达式的表示。 1. pretty函数（排版形式输出） 2. collect函数（合并同类项多项式） 3. expand函数（展开为多项式） 4. horner函数（写成嵌套形式） 5. factor函数（因式分解） 6. simplify函数（简化） 7. simple函数（给出多种简化形式） 3.3.3 符号表达式的替换 subexpr函数对子表达式是自动寻找的，只有比较长的子表达式才被置换，比较短的子表达式，即使重复出现多次，也不被置换。 subs函数可用来进行对符号表达式中符号变量的替换。 f(x) g(y) comp