prob_ans_a【DOC精选】.doc

符号运算 习题2及解答 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？ 3/7+0.1; sym(3/7+0.1); sym(3/7+0.1); vpa(sym(3/7+0.1)) 〖目的〗 不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+0.1 c2=sym(3/7+0.1) c3=sym(3/7+0.1) c4=vpa(sym(3/7+0.1)) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = 0.5286 c2 = 37/70 c3 = 0.52857142857142857142857142857143 c4 = 0.52857142857142857142857142857143 Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量. sym(sin(w*t)),sym(a*exp(-X)),sym(z*exp(j*th)) 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym(sin(w*t)),1) ans = w symvar(sym(a*exp(-X)),1) ans = a symvar(sym(z*exp(j*th)),1) ans = z 求以下两个方程的解 （1）试写出求三阶方程正实根的程序。注意：只要正实根，不要出现其他根。 （2）试求二阶方程在时的根。 〖目的〗 体验变量限定假设的影响 〖解答〗 （1）求三阶方程正实根 reset(symengine) %确保下面操作不受前面指令运作的影响 syms x positive solve(x^3-44.5) ans = (2^(2/3)*89^(1/3))/2 （2）求五阶方程的实根 syms a positive %注意：关于x的假设没有去除 solve(x^2-a*x+a^2) Warning: Explicit solution could not be found. In solve at 81 ans = [ empty sym ] syms x clear syms a positive solve(x^2-a*x+a^2) ans = a/2 - (3^(1/2)*a*i)/2 a/2 + (3^(1/2)*a*i)/2 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。 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))等来观察。 〖目的〗 理解准确符号数值的创建法。 高精度误差的观察。 〖解答〗 （1）x=7/3 x=7/3;a=x,b=sym(x),c=sym(x,d),d=sym(7/3), a = 2.3333 b = 7/3 c = 2.3333333333333334813630699500209 d = 7/3 v1=vpa(abs(a-d)),v2=vpa(abs(b-d)),v3=vpa(abs(c-d)) v1 = 0. v2 = 0. v3 = .1480297366166876e-15 （2）x=pi/3 x=pi/3;a=x,b=sym(x),c=sym(x,d),d=sym(pi/3), a = 1.0472 b = pi/3 c = 1.0471975511965976313177861811710 d = pi/3 v1=vpa(abs(a-d)),v2=vpa(abs(b-d)),v3=vpa(abs(c-d)) v1 = 0. v2 = 0. v3 = .1148364282799222e-15 （3）x=pi*3^(1/3) x=pi*3^(1/3);a=x,b=sym(x),c=sym(x,d),d=sym(pi*3^(1/3)) a = 4.5310 b = 5101408179057732*2^(-50) c = 4.5309606547207899041040946030989 d = pi*3^(1/3) v1=vpa(abs(a-d)),v2=vpa(abs(b-