资源描述：
《matlab课后习题解答第二章.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第2章符号运算习题2及解答1说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？3/7+0.1;sym(3/7+0.1);sym('3/7+0.1');vpa(sym(3/7+0.1))〖目的〗l不能从显示形式判断数据类型，而必须依靠class指令。〖解答〗c1=3/7+0.1c2=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.5286c2=37/70c3=0.528571428571428571
2、42857142857143c4=0.52857142857142857142857142857143Cs1=doubleCs2=symCs3=symCs4=sym2在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量.sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')〖目的〗l理解自由符号变量的确认规则。〖解答〗symvar(sym('sin(w*t)'),1)ans=wsymvar(sym('a*exp(-X)'),1)ans=a8symvar(sym('z*exp(j*th)'),1)ans=z5求符号矩阵的行列式值
3、和逆，所得结果应采用“子表达式置换”简洁化。〖目的〗l理解subexpr指令。〖解答〗A=sym('[a11a12a13;a21a22a23;a31a32a33]')DA=det(A)IA=inv(A);[IAs,d]=subexpr(IA,d)A=[a11,a12,a13][a21,a22,a23][a31,a32,a33]DA=a11*a22*a33-a11*a23*a32-a12*a21*a33+a12*a23*a31+a13*a21*a32-a13*a22*a31IAs=[d*(a22*a33-a23*a32),-d*(a12*a33-a13*a32),d*(a12*a23-a13*a
4、22)][-d*(a21*a33-a23*a31),d*(a11*a33-a13*a31),-d*(a11*a23-a13*a21)][d*(a21*a32-a22*a31),-d*(a11*a32-a12*a31),d*(a11*a22-a12*a21)]d=1/(a11*a22*a33-a11*a23*a32-a12*a21*a33+a12*a23*a31+a13*a21*a32-a13*a22*a31)8（1）通过符号计算求的导数。（2）然后根据此结果，求和。〖目的〗ldiff,limit指令的应用。l如何理解运行结果。〖解答〗symsty=abs(sin(t))8d=diff(y)%求
5、dy/dtd0_=limit(d,t,0,'left')%求dy/dt
6、t=0-dpi_2=limit(d,t,pi/2)%求dy/dt
7、t=pi/2y=abs(sin(t))d=sign(sin(t))*cos(t)d0_=-1dpi_2=09求出的具有64位有效数字的积分值。〖目的〗l符号积分的解析解和符号数值解。l符号计算和数值计算的相互校验。〖解答〗（1）符号积分symsxclearsymsxy=exp(-abs(x))*abs(sin(x))si=vpa(int(y,-10*pi,1.7*pi),64)y=abs(sin(x))/exp(abs(x))si=1.08784949941
8、2904913166671875948174520895458535212845987519414166（2）数值计算复验xx=-10*pi:pi/100:1.7*pi;sn=trapz(exp(-abs(xx)).*abs(sin(xx)))*pi/100sn=1.087710计算二重积分。〖目的〗l变上限二重积分的符号计算法。〖解答〗symsxyf=x^2+y^2;r=int(int(f,y,1,x^2),x,1,2)r=1006/10511在区间，画出曲线，并计算。〖目的〗l在符号计算中，经常遇到计算结果是特殊经典函数的情况。8l如何应用subs获得超过16位有效数字的符号数值结果。l初
9、步尝试ezplot指令的简便。〖解答〗（1）符号计算symstx;f=sin(t)/t;y=int(f,t,0,x)%将得到一个特殊经典函数y5=subs(y,x,sym('4.5'))ezplot(y,[0,2*pi])y=sinint(x)y5=1.6541404143792439835039224868515（2）数值计算复验tt=0:0.001:4.5;tt(1)=eps;yn=trapz
