MATLAB课件 第5章 MATLAB符号运算.docVIP

第五章 MATLAB符号运算 5.1符号矩阵的生成 在MATLAB中输入符号向量或者矩阵的方法和输入数值类型的向量或者矩阵的方法在形式上很相像，只不过要用到符号矩阵定义函数sym，或者是用到符号定义函数syms。 1．用命令sym定义矩阵： 这时的函数sym实际是在定义一个符号表达式，这时的符号矩阵中的元素可以是任何的符号或者是表达式，而且长度没有限制，只是将方括号置于用于创建符号表达式的单引号中。如下例：注意：标点符号的区别 【例5.1 】 A = sym([b c d;Tom,HelpMe,NOWAY]) A = [ b, c, d] [ Tom, HelpMe, NOWAY] %符号矩阵的每一行都用[]将其界定，而数值矩阵没有。 B = sym([10 20 30;aa bb cc; sin(x) cos(x) tan(x)]) B = [10 20 30] [aa bb cc] [sin（x）cos（y）tan（z）] 2．用命令syms定义矩阵 先定义矩阵中的每一个元素为一个符号变量，而后像普通矩阵一样输入符号矩阵。 【例5.2】 syms a b c ;s1 = sym(red);s2 = sym( gree);s3 = sym(Blue); sy = [a b c; s1 s2 s3;2 3 8] sy = [ a, b, c] [ red, gree, Blue] [ 2, 3, 8] syms a b c d %符号变量之间用空格隔开，不能用逗号。 A=[a b;c d];D=det(A) D = a*d - b*c 3把数值矩阵转化成相应的符号矩阵。 数值型和符号型在MATLAB中是不相同的，它们之间不能直接进行转化。将数值型转化成符号型的命令：sym()。 【例5.3】 a1=1/3,a2=sym(1/3) a1 = 0.3333 a2 = 1/3 A = [1/3,sqrt(2),3.42;exp(0),log10(100), 23^(-1)] B= sym(A) A = 0.3333 1.4142 3.4200 1.0000 2.0000 0.0435 B = [ 1/3, sqrt(2), 171/50] [ 1, 2, 1/23] 注意：不论矩阵是用分数形式还是浮点形式表示的，将矩阵转化成符号矩阵后，都将以最接近原值的有理数形式表示或者是函数形式表示。 5.2符号矩阵的运算 5.2.1 算术符号操作 命令 + - * .* \ .\ / ./ ^ .^ 功能 符号矩阵的算术操作 用法如下：（和数值矩阵的用法基本是一致的） A+B、A-B 符号阵列（矩阵）的加法与减法。即符号矩阵和符号阵列的加减法一致。 若A与B为同型阵列时，A+B、A-B分别对对应分量进行加减；若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行加减。 【例5.4】 syms a b c d;A=[a b;c d],B=[d c;b a],c1=A+B,c2=A+1 A = [ a, b] [ c, d] B = [ d, c] [ b, a] c1 = [ a + d, b + c] [ b + c, a + d] c2 = [ a + 1, b + 1] [ c + 1, d + 1] A*B 符号矩阵乘法。 A*B为线性代数中定义的矩阵乘法。按乘法定义要求必须有：矩阵A的列数等于矩阵B的行数。即：若An*k*Bk*m=(aij)n*k.*(bij)k*m=Cn*m=(cij)n*m，则，i=1,2,…,n；j=1,2,…,m。或者至少有一个为标量时，方可进行乘法操作，否则将返回一出错信息。 【例5.5】 syms a b c d;A=[a b;c d] ;B=[d c;b a] ;c1=A*B,c2=B*A,c3=2*A c1 = [ b^2 + a*d, a*b + a*c] [ b*d + c*d, c^2 + a*d] c2 = [ c^2 + a*d, b*d + c*d] [ a*b + a*c, b^2 + a*d] c3 = [ 2*a, 2*b] [ 2*c, 2*d] A.*B 符号数组的乘法。 A.*B为按参量A与B对应的分量进行相乘。A与B必须为同型阵列，或至少有一个为标量。即：An*m.*Bn*m=(aij)n*m.*(bij)n*m=Cn*m=(c