符号计算答案.docVIP

符号运算作业 观察一个数（在此用@记述）在以下四条不同指令作用下的异同： a = @ b = sym( @ ) c = sym( @ ,d ) d = sym( @ ) %这给出完全准确值 在此，@ 分别代表具体数值 7/3 , pi/3 , pi*3^(1/3) ；而异同通过vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d))等来观察。 程序： a = 7/3 , b = sym(7/3), c = sym( 7/3 ,d), d = sym( 7/3 ) vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d)) a =pi/3 , b = sym( pi/3 ), c = sym( pi/3,d), d = sym( pi/3 ) vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d)) a = pi*3^(1/3) , b = sym( pi*3^(1/3) ), c = sym( pi*3^(1/3),d), d = sym( pi*3^(1/3) ) vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d)) 〖答案〗 （1）x=7/3 a = 2.3333 b = 7/3 c = 2.3333333333333334813630699500209 d = 7/3 v1 = 0. v2 = 0. v3 = .1480297366166876e-15 （2）x=pi/3 a = 1.0472 b = pi/3 c = 1.0471975511965976313177861811710 d = pi/3 v1 = 0. v2 = 0. v3 = .1148364282799222e-15 （3）x=pi*3^(1/3) a = 4.5310 b = 5101408179057732*2^(-50) c = 4.5309606547207899041040946030989 d = pi*3^(1/3) v1 = .2660111416629094e-15 v2 = .2660111416629094e-15 v3 = .2660111416629094e-15 说出以下三条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”对象？ 3/7+0.1, sym(3/7+0.1), vpa(sym(3/7+0.1)) 程序： class(3/7+0.1), class(sym(3/7+0.1)), class(vpa(sym(3/7+0.1),4)) 〖答案〗 c1 = 0.5286 c2 = 37/70 c3 = .52857142857142857142857142857143 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是独立自由变量。 sym(sin(w*t)) , sym(a*exp(-X) ) , sym(z*exp(j*theta)) 程序： findsym(sym(sin(w*t)),1) findsym(sym(a*exp(-X)),1) findsym(sym(z*exp(j*th)),1) 〖答案〗 ans = w ans = a ans = z 求符号矩阵的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 程序： syms a11 a12 a13 a21 a22 a23 a31 a32 a33 A = [a11 a12 a13;a21 a22 a23;a31 a32 a33] a=det(A) B=inv(A) C=subexpr(B) [RS,w]=subexpr(B,w) 〖答案〗 A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22 IAs = [(a22*a33-a23*a32)/d, -(a12*a33-a13*a32)/d, -(-a12*a23+a13*a22)/d] [-(a21*a33-a23*a31)/d