第5章MATLAB符号运算汇编.ppt

第5章 MATLAB符号运算;本章主要包括： 5.1 符号对象及其表达式 5.2 符号算数运算 5.3 独立变量与表达式化简 5.4 符号微积分运算 5.5 符号积分变换 5.6 方程的解析解;5.1 符号对象及其表达方式 ;5.1.1定义符号常量 ;其中flag为可选参数，有4种形式，分别是r、d、e或f。 r：用有理数格式表达符号量（默认格式） d：十进制数 e：带有机器浮点误差的有理数 f：浮点数 例5.1 将一组数值常量定义成符号常量。 log(2) %数值常量 f1=sym(log(2))　　%符号常量，注意f1在工作空间中的类型 ;2、定义符号变量 一般有两种方法 （1）使用sym( )函数： sym(x) sym(x,real) sym(x,unreal) （2）? 使用syms命令 syms arg1 arg2 ... syms arg1 arg2 ... real syms arg1 arg2 ... Unreal ;例子5.3,5.4 sym(x) %定义符号变量x sym(y,real) %定义符号变量y，且同时为实型符号量 syms a b c syms m n real syms x y z unreal ;5.1.2 符号表达式 ;5.1.3 符号矩阵 ; 再通过sym命令可直接将数值矩阵转换为符号矩阵，例如： S=sym(M) 例5.8 用size()函数求符号矩阵的大小 s=size(M) [s_r,s_c]=size(M) s_r=size(M,1) s_c=size(M,2) 例5.9 抽取符号矩阵中的元素 B=sym([a,b,c;d,e,f;g,h,k]) B(2,3),B(6) ;5.2 符号算术运算 ;5.2.1符号对象的加减;A*B、A/B可分别用来求A和B两个符号矩阵的乘法与除法。A.*B则用来实现两个符号数组的乘法。其中矩阵除法也可用来求解符号线性方程组的解。 例子5.12 syms a b c d e f g h; A = [a b; c d]; B = [e f; g h]; C1 = A.*B D1=A*B;例5-13已知多项式，，求两个多项式的积和商。 syms x fx gx fx = 3*x^5-x^4+2*x^3+x^2+3 gx= 1/3*x^3+x^2-3*x-1 fx*gx expand(fx*gx) %展开积的符号表达式 ;5.3 独立变量与表达式化简 ;findsym( )查找MATLAB的符号表达式或矩阵中哪一个变量是独立变量 （1）r = findsym(S)：以字母表的顺序返回表达式S中的所有符号变量（注：符号变量由除了i与j的字母与数字构成的、字母打头的字符串）。若S中没有任何的符号变量，则findsym返回一空字符串。 （2）r = findsym(S,n)：返回表达式S中依接近x的顺序而排列的n个符号变量 syms a b c d x y z w; findsym(a+b+c+x+y+z+w) ans =a, b, c, w, x, y, z findsym(a+b+c+x+y+z+w,1) ans =x findsym(a+b+c+x+y+z+w,2) ans =x,y findsym(a+b+c+x+y+z+w,3) ans =x,y,w;例子5.14 syms a x y z t alpha beta S1 = findsym(sin(pi*t*alpha+beta),1) S2 = findsym(x+i*y-j*z+eps-nan) S2 = NaN, x, y, z S3 = findsym(a+y,2) S3 = y,a ;5.3.2 表达式化简 ;1. 合并同类项（collect） ;2. 多项式展开（expand） ;3. 因式分解（factor） ;例5-17 因式分解示例 syms a b x y F1 = factor(x^4-y^4) F3 = factor(sym(12345678901234567890)) F3 = (2)*(3)^2*(5)*(101)*(3803)*(3607)*(27961)*(3541) ;4. 一般化简（simplify） ;5. 不定化简（simple） ;7. 书写格式美化（pretty） ;5.4 符号微积分运算 ;例5-