MATLAB习题答案_张志涌2008.docVIP

MATLAB教程（张志涌 杨祖樱编2006.8第一版）习题答案 第一章 5,6,7 第二章 2,3,6,7,8,9,10,11,22,24 第三章 2,5,6 第四章 2,5,11,13,17,18 第五章 1,2,5,6,7 第六章 1,2,3 第七章 1,2 第一章5,6,7 %p41_5 直接计算时得到处于第一象限的方根 a=-8; a^(2/3) %-2.0000 + 3.4641i (a^2)^(1/3) %等价于64^(1/3) 得4.0000 (a^(1/3))^2 %等价于a^(1/3)也就是1.0000 + 1.7321i的平方 得-2.0000 + 3.4641i 第二章2,3,6,7,8,9,10,11,22,24 %p92_3 A1=sym(sin(w*t)); findsym(A1) A2=sym(a*exp(-X)); findsym(A2) A3=sym(x*exp(j*theta)); findsym(A3) %p92_6 clear, syms x k, x=abs(x);%或syms x positive; f1=sym(2./(2*k+1).*((x-1)./(x+1)).^(2*k+1)) s1=symsum(f1,k,0,inf), simple(s1), simple(ans) %p92_7 clear t, syms t df=diff(abs(sin(t))),class(df) df1=limit(df,t,0,left) df2=subs(df,t,sym(pi/2)) %p92_8 clear, syms x; f=exp(-abs(x)).*abs(sin(x)); fint=int(f,x,-5*pi,1.7*pi), digits(64), vpa(fint) %p92_9 clear,syms x y,f=x.^2+y.^2, fint=(int(int(f,y,1,x.^2),x,1,2)), double(fint) %p92_10 clear, syms t x; f=sin(t)./t, yx=int(f,t,0,x), ezplot(yx,[0 2*pi]) fint=subs(yx,x,4.5),%或yxd=int(f,t,0,4.5),fint=double(yxd) hold on, plot(4.5,fint,*r) %p92_11 clear, syms x n; f=sin(x).^n; yn=int(f,x,0,pi/2), class(yn) %y(1/3)=? yn1=subs(yn,n,sym(1/3)),vpa(yn1) %或yn=limit(yn,n,1/3),vpa(yn) %或yy=int(sin(x).^(1/3),x,0,pi/2) ,vpa(yy) %p93_22 clear, syms x y S S = dsolve(Dy*y/5+x/4=0,x) ezplot(subs(S(1),C1,1),[-2,2 -2,2],1), hold on ezplot(subs(S(2),C1,1),[-2,2 -2,2],1) %p93_24 [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,g(0)=1,x) 第三章 2,5,6 %122_1 a=linspace(0,2*pi,10),b=0:2*pi/9:2*pi %122_2 rand(twister,0),A=rand(3,5),B=A0.5,C=find(B) [ci,cj]=ind2sub(size(A),C) %此法太繁 [Ci,Cj]=find(B) %122_5 t=linspace(1,10,100) (1) y=1-exp(-0.5.*t).*cos(2.*t),plot(t,y) (2) L=length(t) for k=1:L, yy(k)=1-exp(-0.5*t(k))*cos(2*t(k)), end, plot(t,yy) %122_6 clear,format long, rand(twister,1),A=rand(3,3), B=diag(diag(A)),C=A.*(~B)%或C=A-B%或C=triu(A,1)+tril(A,-1) 第四章2,5,11,13,17,18 %p171_2 d=0.5; tt=0:d:10; t=tt+(tt==0)*eps; y=sin(t)./t; s=d*trapz(y) % 计算出积分值 ss=d*(cumtrapz(y)) %计算梯形法累计积分并绘积 分曲线 pl