第5章节 符号数学基础.pptVIP

物理与电气工程学院 物理与电气工程学院 * 第5章 符号数学基础 　　MATLAB以MAPLE为内核,作为符号计算的“引擎”，开发了相关的符号数据工具箱，采用了一个全新的数据结构：符号对象。 　　符号计算过程：首先要定义基本的符号对象（常数、变量或表达式等），然后利用这些基本符号对象去构成新的表达式，最后进行所需的符号运算。其中,符号表达式也是符号对象。 第1节　符号对象的创建 一、创建符号变量和表达式 创建和定义符号对象的函数：sym和syms 1、sym函数定义符号变量、符号表达式 x=sym(‘x’) 定义符号变量 x=sym(‘expression’)　定义符号表达式 f=sym(‘expression=0’) 定义符号方程 【例5-1】使用sym函数创建符号变量。 a=sym(‘a’) 　　% sym object 1×1 b=sym(‘hello’)　 % sym object 1×1 c=sym(‘(1+sqrt(5))/2’)　 % sym object 1×1 cstr= (1+sqrt(5))/2 % char array 1×13 y=sym(x^3+5*x^2+12*x+20) % sym object 1×1 2、syms函数定义符号变量 在符号变量和变量值相同时，可以用sym的简捷方式建立符号变量，即表达在符号变量和变量值相同时，可以用sym的简捷方式建立符号变量，即：syms函数来定义符号变量 syms arg1 arg2 … 　　等价于arg1=sym(‘arg1’); 　　　　　arg2=sym(‘arg2’) 区别：sym 中的参数要用单引号括起来，且只能定义一个符号变量；而syms的变量参数是用空格分隔，参量不用单引号，同时可以定义多个符号变量 【例5-2】使用syms函数创建符号变量。 syms a b c d 物理与电气工程学院 二、创建符号矩阵 1、用syms函数创建 选用syms创建符号矩阵所用的变量；再创建符号矩阵 【例5-3】创建一个循环矩阵。 syms a b c d n=[a b c d;b c d a;c d a b;d a b c] 2、用sym函数创建 矩阵元素是不带等号的符号表达式，各矩阵元素的长度可以不同，行与行之间用分号隔开，元素之间用逗号或空格分开 【例5-4】将3阶Hilbert矩阵转换为符号矩阵。 h=hilb(3) h1=sym(h) h2=sym([h1 h1+h1/2;3*h1 h1]) 三、默认符号变量 　　MATLAB根据数学约定，自变量一般都是小写字母，并且在字母表的后面(如x,y,z)。可以用findsym函数找到自变量的符号。符号表达式中的缺省变量为靠近小写字母x的优先，在x后面的字母优先。 　　i和j为虚数单位，在符号运算中不能作为自变量 【例5-5】求符号函数在不同自变量情况下的结果 syms x n f=x^n; diff(f); diff(f,n) findsym(expression,n) 　 　　　查找符号表达式中前n个符号变量 【例5-6】查询符号函数中的默认自变量。 syms a b n t x f=a*x^n+b*t findsym(f,1) findsym(f,2) findsym(f,5)%列出前5个变量 findsym(f)%列出所有变量 第2节　符号表达式的化简和替换 一、符号表达式的化简 1、因式分解 格式：factor(s) 若s是个多项式、系数为有理数，则会将s表示成系数有理数的低阶多项式相乘； 若s不能被分解成有理式，则返回s自身； 若s是矩阵，则对矩阵s中的每个元素进行分解； factor命令也可对符号整数 进行因式分解 【例5-7】对表达式 进行因式分解。 syms x f=factor(x^9-1) pretty(f) %按手写格式显示 【例5-8】对大整数12345678901234567890进行因式分解。 factor(sym(12345678901234567890)) 2、符号表达式的展开 格式： expand(s) 　将多项式s展开 【例5-9】展开表达式f=(x+1)^5和f=sin(x+y)。　syms x y　f=(x+1)^5;　expand(f)　f=sin(x+y);　expand(f) 3、同类项合并函数collect MATLAB的同类项合并函数为collect(s)。 collect(s) s为关于符号变量的多项式形式，它将s中的x(缺省)的同幂项系数进行合并。 collect(s,v) 将s中的以指定变量v的同幂项系