MATLAB教程R2007习题答案张志涌2007.docVIP

MATLAB教程（张志涌 杨祖樱编2006.8第一版）习题答案 第二章 %p92_1 a=7/3;class(a) b=sym(7/3);class(b) c=sym(7/3,d),class(c) d=sym(7/3),class(d) vpa(abs(a-d)), vpa(abs(b-d)), vpa(abs(c-d)) %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_4 syms a11 a12 a13 a21 a22 a23 a31 a32 a33 W; A=[a11 a12 a13;a21 a22 a23;a31 a32 a33] adet=det(A), [adet1 W]=subexpr(adet,W) ainv=inv(A), [ainv1 W]=subexpr(ainv,W) %p92_6 clear, syms x k, x=abs(x); 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=subs(df,t,0), 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 %p92_11 clear, syms x n; f=sin(x).^n; yn=int(f,x,0,pi/2), class(yn) %y(5)=? yn1=subs(yn,n,sym(5)) %p93_19 clear, syms x y S, S=solve(x*x+y*y=1,x*y=2,x,y), double(S.x), double(S.y) [x,y]=solve(x*x+y*y=1,x*y=2,x,y), vpa(x,6) %p93_20 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_21 clear, f=dsolve(Dx=a*t*t+b*t,x(0)=2,t) %p93_22 [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,g(0)=1,x) 第三章 %122_1 a=linspace(0,2*pi,10),b=0:2*pi/9:2*pi %122_2 rand(state,0),A=rand(3,5),B=A0.5,C=find(B),[ci,cj]=ind2sub(size(A),C)(此法太繁) [Ci,Cj]=find(B) %122_3 rand(state,123),A=rand(1,1000) B=(A=0.5)*(+1),C=(~(A=0.5))*(-1),D=C+B,[ri,rj]=size(find(C)) %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)=-exp(-0.5*t(k))*cos(2*t(k)), end, plot(t,yy) %122_6 clear,format long, rand(state,1),A=rand(3,3),B=diag(diag(A)),C=A.*(~B) %p122_7 x=-3*pi:pi/15:3*pi; % delete 0 element % x=x+(x==0)*eps; y=x; [X,Y]=meshgrid(x+eps,y+eps);warning off; Z=sin(X)./X.*sin(Y)./Y; %1 how many NaN a=size(find(X==0)); a1=a(1)*a(2); b=size(find(Y=