Matlab程序设计与应用第二版刘卫国课后实验解答

Matlab程序设计与应用第二版刘卫国课后实验答案实验一: C=A*B T1: D=C(3:5,2:3) T4-1: %第一小题 z1=2*sin(85*pi/180)/(1+exp(2) a=100:999; %第二小题 b=find(rem(a,21)=0); x=2,1+2i;-0.45,5; c=length(b) z2=1/2*log(x+sqrt(1+x.2); T4-2: z2 a=input(请输入一个字符串:,s); %第三小题 b=find(a=A disp(帕萨卡矩阵P性能更好); disp(A.3); elseif a=90 disp(the score is not reasonable); A=reshape(a,5,5) s=input(please enter the score:); disp(B1=sqrtm(A); end B1=sqrtm(A) switch fix(s/10) disp(B2=sqrt(A); case9,10 B2=sqrt(A) disp(A); disp(B1*B1); case8 B1*B1 disp(B); disp(B2.*B2); case7 B2.*B2 disp(C); 实验三: case6 T1: disp(D); x=-5.0:2:5.0 case0,1,2,3,4,5 for i=1:length(x); disp(E); if (x(i)=0) w=120*84+(t-120)*1.15*84; else elseif t100) T4: disp(the score is not reasonable); a=10+floor(rand(1)*89) s=input(please enter the score:); b=10+floor(rand(1)*89) 2 s=input(请输入+或-或*或/,s); for i=1:n; while(s=+ if y=3 disp(输入的符号错误，请重新输入); s=input(请输入+或-或*或/,s); disp(最大n值为:); end i-1 switch s disp(对应的y值为:); case+ y-1/(2*i-1) c=a+b; break; case- end end c=a-b; T3: case* c=a*b; format long a=input(请输入a:); case/ b=input(请输入b:); c=a/b; x(1)=a/(1.0+b);i=1;x(2)=a/(x(i)+b); end while (abs(x(i+1)-x(i)0.00001 B=A(5,:) format short else T4: B=A(n,:) f(1)=1;f(2)=0;f(3)=1;i=4; End while(i0) pii=sqrt(6*m) fn=length(find(f0.5); y(i)=(x(i)+sqrt(pi)/exp(2); b=length(k); else disp(大于0.5的随机数个数占总数的百分y(i)=1/2*log(x(i)+sqrt(1+x(i)2); 比:); end sprintf(%f%, b/300) end T2: subplot(1,2,2); P=45+50*rand(100,5) plot(x,y); disp(最高分及序号:); T4: Y,U=max(P,1) a=input(请输入a:); disp(最低分及序号:); b=input(请输入b:); Y1,U1=min(P,1) n=input(请输入n:); disp(每门课平均分:); 5 mean(P,1) i=i+1; disp(每门课的标准方差:); x=x+0.01; end std(P,0,1) s=sum(P,2); dx=diff(g)/0.01; disp(总分最高分:); dx(1) Y2,U2=max(s,1) dx(101) disp(总分最低分:); dx(length(g)-1) Y3,U3=min(s,1) T2: zcj,xsxh=sort(s,1,descend) %descend降f1=inline(sqrt(cos(t.2)+4*sin(2*t).2+1) 序 f2=inline(log(1+x)./(1+x.2) T3: I1=quad(f1,0,2*pi) h=6:2:18; I2=quad(f2,0,2*pi) t1=18.0,20.0,22.0,25.0,30.0,28.0,24.0; T3: t2=15.0,19.0,24.0,28.0,34.0,32.0,30.0; A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2; t=6.5:2:17.5; B=-4,13,1,11; T1=interp1(h,t1,t,spline) x1=AB T2=interp1(h,t2,t,spline) L,U=lu(A); T4: x2=U(LB) x=1:10:101; Q,R=qr(A); y=log10(x); x3=R(QB) p=polyfit(x,y,5) fun: x1=1:0.01:101; function F=fun(X); y1=log10(x1); x=X(1); y2=polyval(p,x1); y=X(2); plot(x1,y1,r-,x1,y2,b-.); z=X(3); T5: F(1)=sin(x)+y2+log(z)-7; P1=1,2,4,0,5; F(2)=3*x+2y-z3+1; P2=0,0,0,1,2; F(3)=x+y+z-5; P3=0,0,1,2,3; T5: P4=conv(P2,P3); f=inline(3*x+sin(x)-exp(x); sp1=length(P1); fzero(f,1.5) sp2=length(P4); X=fsolve(fun,1 1 1,optimset(Display,off) P=zeros(1,sp2-sp1),P1+P4 实验十: x=roots(P) T1: A=-1,1.2,-1.4;0.75,2,3.5;0,5,2.5; x=sym(6); Ax=polyval(P,A) y=sym(5); AA=polyvalm(P,A) z=(x+1)/(sqrt(3+x)-sqrt(y) 实验九: T2: T1: syms x y; f=inline(det(x x2 x3;1 2*x 3*x2;0 2 t=sym(5135); 6*x); factor(x4-y4) i=1;x=1; factor(t) while x=3.01 T3: g(i)=f(x); syms beta1 beta2 x; 6 simplify(sin(beta1)*cos(beta2)-cos(beta1)*sinsyms n x ; (beta2) S1=symsum(1/(2*n-1),n,1,10) S2=symsum(n2*x(n-1),n,1,inf) simplify(4*x2+8*x+3)/(2*x+1) T4: S3=symsum(n2/5n,n,1,inf) syms a b c d e f g h k; T2: p1=0,1,0;1,0,0;0,0,1; x=sym(x); p2=1,0,0;0,1,0;1,0,1; f=log(x); A=a,b,c;d,e,f;g,h,k; taylor(f,x,6,1) B=p1*p2*A T3: B1=inv(B) x1=solve(log(1+x)-5/(1+sin(x)=2) B1*B; x2=solve(x2+9*sqrt(x+1)-1=0) tril(B) x3=solve(3*x*exp(