Matlab 第5章 符号运算要点.pptVIP

【例5-4】 使用syms函数定义符号变量。 syms a b c x f=sym(a*x^2 + b*x + c) f = a*x^2 + b*x + c g=f^2+4*f-2 g = (a*x^2+b*x+c)^2+4*a*x^2+4*b*x+4*c-2 注意：使用这种简写形式定义符号变量时，各个变 量要用空格来进行分隔，而不能使用逗号或分号等。 比如syms a,b,c,x 则为错误的格式。 ② 用子矩阵创建符号矩阵 这种方法是指首先创建出符号矩阵中的子矩 阵，再利用这些子矩阵组合成所需矩阵。需要注 意的是：创建子矩阵时，要保证各子矩阵相对应 的各列长度一致。 【例5-24】 用子矩阵创建符号矩阵。 B=[[a ,b ,c ];[100,200 ,300];[3*x,cos(y),1/z]] B = [a ,b ,c ] [100,200 ,300] [3*x,cos(y),1/z] * 在上面的例子中可以看到，3个子矩阵的第1列长度都为3个符号长度，长度不足的部分用空格符补足。比如第1个子矩阵中元素a后还须添加两个空格符以保证各子矩阵相对应的列长度一致。第2列、第3列也是如此。 * ③ 由数值矩阵转换为符号矩阵 方法很简单，直接调用sym指令即可。 【例5-25】 把数值矩阵转换为符号矩阵。 C=[1 2 3 4;11 13 15 17;24 26 28 22;3.3 5.5 0 0] C = 1.0000 2.0000 3.0000 4.0000 11.0000 13.0000 15.0000 17.0000 24.0000 26.0000 28.0000 22.0000 3.3000 5.5000 0 0 * D=sym(C) D = [ 1, 2, 3, 4] [ 11, 13, 15, 17] [ 24, 26, 28, 22] [ 33/10, 11/2, 0, 0] * 2. 符号矩阵的运算 符号矩阵的常用运算形式有四则运算、 逆运算、幂运算、求秩和求行列式值等，在 MATLAB中，这些运算的指令都很简单。下 面通过几个例子对上述运算分别进行介绍。 【例5-26】 符号矩阵的四则运算。 A=sym([1 2 3 4;1/5 1/6 1/7 1/8;60 70 80 90;a b c d]); B=sym([1 3 5 7;2 4 6 8;x x^2 y y^2;a b c d]); C=A+B * C = [ 2, 5, 8, 11] [ 11/5, 25/6, 43/7, 65/8] [ 60+x, 70+x^2, 80+y, 90+y^2] [ 2*a, 2*b, 2*c, 2*d] D=A-B D = [ 0, -1, -2, -3] [ -9/5, -23/6, -41/7, -63/8] [ 60-x, 70-x^2, 80-y, 90-y^2] [ 0, 0, 0, 0] * E=A*B E = [ 5+3*x+4*a, 11+3*x^2+4*b, 17+3*y+4*c, 23+3*y^2+4*d] [8/15+1/7*x+1/8*a,19/15+1/7*x^2+1/8*b,2+1/7*y+1/8*c,41/15+1/7*y^2+1/8*d] [ 200+80*x+90*a, 460+80*x^2+90*b, 720+80*y+90*c,980+80*y^2+90*d] [ a+2*b+c*x+d*a,3*a+4*b+c*x^2+d*b,5*a+6*b+c*y+d*c,7*a+8*b+c*y^2+d^2] * 【例5-27】求符号矩阵的行列式值、特征根、逆运算和秩。 syms x y; A=[1 2;x y]; D=det(A) %求矩阵A的行列式值 D =y-2*x I=inv(A) %求矩阵A的逆矩阵 I = [ -y/(-y+2*x), 2/(-y+2*x)] [ x/(-y+2*x), -1/(-y+2*x)] E=eig(A) %求矩