MATLAB实用教程第版郑阿奇第三章.ppt

21 ???? 第3章 MATLAB符号计算Using Symbolic Math Toolbox 3.1 符号表达式的建立 3.2符号表达式的代数运算 3.3符号表达式的操作和转换 3.4 符号极限、微积分和级数求和 3.5 符号积分变换 3.6符号方程的求解 MATLAB具有符号数学工具箱(Symbolic Math Toolbox) 。 符号计算是可以对未赋值的符号对象(可以是常数、变量、表达式)进行运算和处理。 与数值运算的区别: ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 符号运算的功能 符号线性代数(linear algebra) 因式分解、展开和简化(simplification and substitution) 符号代数方程求解(solving equations) 符号微积分(Calculus) 符号微分方程 3.1 符号表达式的建立 3.1.1 创建符号常量 符号常量是不含变量的符号表达式。 sym(‘常量’) %创建符号常量 sym(常量,参数) %按某种格式转换为符号常量 说明：参数可以选择为’d’、’f’、’e’或’r’ 四种格式，也可省略。 EX: a=sym(sin(2)) a=sym(sin(2),r) 3.1.2 创建符号变量和表达式 (Creating Symbolic Variables and Expression) 1. 使用sym命令创建符号变量和表达式 sym(‘变量’,参数) %把变量定义为符号对象 2.使用syms命令创建符号变量和符号表达式 syms(‘arg1’, ‘arg2’, …,参数) syms arg1 arg2 …,参数　 例如： f1=sym(‘a*x^2+b*x+c’) %创建表达式 syms a b c x %创建变量 f2=a*x^2+b*x+c %创建表达式 syms(a,b,c,x) f3=a*x^2+b*x+c 符号表达式（）中的参数一定要用 单引号括起来。 3.1.3 符号矩阵 例如，使用sym命令创建的符号矩阵： A=sym([a,b;c,d]) 例如，使用syms命令创建相同的符号矩阵： syms a b c d A=[a b;c d] 比较符号矩阵与字符串矩阵 ： B=‘[a,b;c,d]’ %创建字符串矩阵 C=[a,b;c,d] %创建数值矩阵 ??? Undefined function or variable a. 3.2符号表达式的代数运算 3.2.1符号表达式的代数运算 由于MATLAB采用了重载技术，使得符号表达式的运算符和基本函数都与数值计算中的几乎完全相同 。 例如： A+2 A.’ det(A) A^2 1. 符号运算中的运算符 基本运算符 运算符“＋”，“－”，“*”，“\”，“/”，“^”分别实现符号矩阵的加、减、乘、左除、右除、求幂运算。 运算符“.*”，“./”，“.\”，“.^”分别实现符号数组的乘、除、求幂，即数组间元素与元素的运算。 运算符“′”，“.′”分别实现符号矩阵的共轭转置、非共轭转置。 (2)关系运算符 运算符“= =”、“~=”分别对运算符两边的符号对象进行“相等”、“不等”的比较。 2. 函数运算 三角函数和双曲函数 (2) 指数和对数函数 (3) 复数函数 (4) 矩阵代数命令 3.2.2 符号数值任意精度控制和运算 在Symbolic Math Toolbox中有三种不同的算术运算： 数值型：MATLAB的浮点运算。 有理数型：Maple的精确符号运算。 VPA型：Maple的任意精度运算。 任意精度的VPA型运算可以使用digits和vpa命令来实现。 digits(n) %设定默认的精度 S=vpa(s,n) %将s表示为n位有效位数的符号对象 3.2.3 符号对象与数值对象的转换 将数值矩阵转化为符号矩阵 函数调用格式：sym(A) EX: A=[1/3,2.5;1/0.7,2/5] sym(A) 将符号矩阵转化为数值矩阵 函数调用格式： numeric(A) EX: a=sym(2/3) b=numeric(a) 3.3符号表达式的操作和转换3.3.1符号表达式中自由变量的确定 符号表达式“f=ax2+bx+c” 中只有一个变量是独立变量： 小写字母i和j不能作为自由变量。 符号表达式中如果有多个符号变量，则按照以下顺序选择自由变量：首先选择x作为自由变量；如果没有x，则