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

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&a=0&t=1&t=2&tb disp(A.3=); disp(帕萨卡矩阵P性能更好); disp(A.3); elseif a=90&s=80&s=70&s=60&s70); B=0.95,0.67,0.52 disp(D); X1=AB else B(3)=0.53 disp(E); X2=AB end disp(系数矩阵A的条件数:); T2-switch: disp(cond(A); s=input(please enter the score:); T6: while (s100) a=1:25; 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)&(x(i)120 y(i)=x(i)2-5*x(i)+6; w=120*84+(t-120)*1.15*84; else elseif t60 y(i)=x(i)2-x(i)-1; w=t*84-700; end else end w=84*t; y end T2-if: disp(应发工资为:); s=input(please enter the score:); disp(w); while (s100) 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=+&s=-&s=*&s=/) y=y+1/(2*i-1); 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&i500) i=i+1; c x(i+1)=a/(x(i)+b); T5: A=rand(5,6) end n=input(请输入n:); r1=(-b+sqrt(b2+4*a)/2; r2=(-b-sqrt(b2+4*a)/2; while (n5 ; 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) T1-2: f0=length(find(f=0) n=input(请输入n:); T5: m=1./(1:n).2; j=0;s=0; s=sum(m); for i=3:50 pii=sqrt(6*s) m=i*(i-1)-1; T2: if isprime(m)=1 n=1000;y=0; m 3 s=s+m; s(n)=10*rem(m(t),10)+fix(m(t)/10); %j=j+1; 挑选出的素数进行位置交换 if isprime(s(n)=1 %再挑选出交换 end end 绝对素数 disp(亲密数的对数:,num2str(j); r(k)=m(t); disp(所有亲密素数之和:,num2str(s); k=k+1; 实验五: end fn1: n=n+1; function f=fn1(n) end f=n+10*log(n.2+5); disp(所有两位绝对素数:,num2str(r) fn2: T5: function f=fn2(n) y1=fn1(40)/(fn1(30)+fn1(20) %函数调用 f=0; y2=fn2(40)/(fn2(30)+fn2(20) for i=1:n 实验六: f=f+i*(i+1); T1: end x=linspace(0,2*pi,101); fx: y=(0.5+3*sin(x)./(1+x.2).*cos(x); function f=fx(x) plot(x,y); f=1./(x-2).2+0.1)+1./(x-3).4+0.01); T2: T1: x=linspace(0,2*pi,100); x=input(请输入一个复数:); y1=x.2; fe=exp(x) y2=cos(2*x); flg=log(x) y3=y1.*y2; fsin=sin(x) subplot(2,2,1); fcos=cos(x) plot(x,y1,r:,x,y2,y-,x,y3,b-.); T2: subplot(2,2,2); m1=input(请输入m1:); plot(x,y1,r:); m2=input(请输入m2:); subplot(2,2,3); t=input(请输入theta:); plot(x,y2,y-); A=m1*cos(t),-m1,-sin(t),0;m1*sin(t),0,. subplot(2,2,4); cos(t),0;0,m2,-sin(t),0;0,0,-cos(t),1 plot(x,y3,b-.); B=0,m1*9.8,0,m2*9.8 T2-3: disp(a1,a2,N1,N2=); x=linspace(0,2*pi,10); disp(AB); y1=x.2; T3: y2=cos(2*x); j=1;n=1;k=1; y3=y1.*y2; for i=10:99 % 挑选出1099中所subplot(3,4,1); 有的素数 bar(x,y1,r);axis(0,7,0,40); if isprime(i)=1 subplot(3,4,2); m(j)=i; stairs(x,y1,y);axis(0,7,0,40); j=j+1; subplot(3,4,3); end stem(x,y1,b);axis(0,7,0,40); end subplot(3,4,4); for t=1:length(m); fill(x,y1,g);axis(0,7,0,40); 4 subplot(3,4,5); theta=0:0.001:2*pi; bar(x,y2,r);axis(-1,7,-1,1); y=a.*sin(b+n.*theta); polar(theta,y); subplot(3,4,6); stairs(x,y2,y);axis(-1,7,-1,1); T5: subplot(3,4,7); x=linspace(-5,5,21); stem(x,y2,b);axis(-1,7,-1,1); y=linspace(0,10,31); subplot(3,4,8); x,y=meshgrid(x,y); fill(x,y2,g);axis(-1,7,-1,1); z=cos(x).*cos(y).*exp(-sqrt(x.2+y.2)/4); subplot(3,4,9); subplot(1,2,1); bar(x,y3,r);axis(0,7,-30,50); surf(x,y,z); subplot(3,4,10); subplot(1,2,2); stairs(x,y3,y);axis(0,7,-30,50); contour3(x,y,z); subplot(3,4,11); T6: stem(x,y3,b);axis(0,7,-30,50); s=0:pi/2; subplot(3,4,12); t=0:3*pi/2; fill(x,y3,g);axis(0,7,-30,50); s,t=meshgrid(s,t); T3: x=cos(s).*cos(t); x=linspace(-5,5,100); y=cos(s).*sin(t); y=; z=sin(s); for x0=x; surf(x,y,z); if x0=0; shading interp; y=y,(x0+sqrt(pi)/exp(2); 实验八: else T1: y=y,1/2*log(x0+sqrt(1+x02); A=rand(1,30000); end disp(均值为:); end mean(A) subplot(1,2,1); disp(标准方差为:); plot(x,y); std(A) disp(最大元素为:); max(A) disp(最小元素为:); for i=1:length(x) min(A) if x(i)0.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(x)+5*sin(x)-78.5=0) det(B) x4,y