2022年Matlab程序设计与应用第二版刘卫国课后实验答案 .pdf

《2022年Matlab程序设计与应用第二版刘卫国课后实验答案 .pdf》由会员分享，可在线阅读，更多相关《2022年Matlab程序设计与应用第二版刘卫国课后实验答案 .pdf（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1 实验一：T1：%第一小题z1=2*sin(85*pi/180)/(1+exp(2) %第二小题x=2,1+2i;-0.45,5; z2=1/2*log(x+sqrt(1+x.2); z2 %第三小题a=-3.0:0.1:3.0; z3=1/2*(exp(0.3*a)-exp(-0.3*a).*sin(a+0.3)+log(0.3+a)/2) %第四题t=0:0.5:2.5 z4=(t=0&t=1&t=2&t=A&ab disp(帕萨卡矩阵P性能更好 ); elseif ab disp(希尔伯特矩阵H 性能更好 ); else disp(两个矩阵性能相同); end T3: a=1:25;

2、A=reshape(a,5,5) disp(行列式的值： ); disp(det(A); disp(矩阵的秩： ); disp(rank(A); disp(矩阵的迹： ); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 7 页 - - - - - - - - - 2 disp(trace(A); disp(矩阵的范数：); disp(norm(A); T4: A=-29,6,18;20,5,12;-8,8,5 V,D=eig(A) T5: A=1/2,1/3,1/4;1

3、/3,1/4,1/5;1/4,1/5,1/6 B=0.95,0.67,0.52 X1=AB B(3)=0.53 X2=AB disp(系数矩阵A 的条件数： ); disp(cond(A); T6: a=1:25; A=reshape(a,5,5) disp(B1=sqrtm(A); B1=sqrtm(A) disp(B2=sqrt(A); B2=sqrt(A) disp(B1*B1); B1*B1 disp(B2.*B2); B2.*B2 实验三：T1: x=-5.0:2:5.0 for i=1:length(x); if (x(i)=0)&(x(i)5)&(x(i)=2)&(x(i)=3)

4、; y(i)=x(i)2-5*x(i)+6; else y(i)=x(i)2-x(i)-1; end end y T2-if: s=input(please enter the score:); while (s100) disp(the score is not reasonable); s=input(please enter the score:); end if (s=90&s=80&s=70&s=60&s70); disp(D); else disp(E); end T2-switch: s=input(please enter the score:); while (s100) d

5、isp(the score is not reasonable); s=input(please enter the score:); end switch fix(s/10) case9,10 disp(A); case8 disp(B); case7 disp(C); case6 disp(D); case0,1,2,3,4,5 disp(E); end T3: t=input( 请输入工时time=); if t120 w=120*84+(t-120)*1.15*84; elseif t60 w=t*84-700; else w=84*t; end disp(应发工资为： ); disp

6、(w); T4: a=10+floor(rand(1)*89) b=10+floor(rand(1)*89) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 7 页 - - - - - - - - - 3 s=input( 请输入 +或-或* 或/,s); while(s=+&s=-&s=*&s=/) disp(输入的符号错误，请重新输入); s=input( 请输入 +或-或*或/,s); end switch s case+ c=a+b; case- c=a-b; c

7、ase* c=a*b; case/ c=a/b; end c T5: A=rand(5,6) n=input( 请输入 n:); while (n5 B=A(5,:) else B=A(n,:) End 实验四：T1: n=input( 请输入 n:); m=0; for i=1:n; m=m+1/i2; end pii=sqrt(6*m) T1-2: n=input( 请输入 n:); m=1./(1:n).2; s=sum(m); pii=sqrt(6*s) T2: n=1000;y=0; for i=1:n; y=y+1/(2*i-1); if y=3 disp(最大 n 值为： );

8、i-1 disp(对应的 y 值为： ); y-1/(2*i-1) break; end end T3: format long a=input( 请输入 a:); b=input( 请输入 b:); x(1)=a/(1.0+b);i=1;x(2)=a/(x(i)+b); while (abs(x(i+1)-x(i)0.00001&i500) i=i+1; x(i+1)=a/(x(i)+b); end r1=(-b+sqrt(b2+4*a)/ 2; r2=(-b-sqrt(b2+4*a)/ 2; disp(x,num2str(i+1),=,num2str(x(i+1),8); disp(x,n

9、um2str(i),=,num2str(x(i),8); disp(r1=,num2str(r1,8),r2=,num2str(r2,8); format short T4: f(1)=1;f(2)=0;f(3)=1;i=4; while(i0) fn=length(find(f0) f0=length(find(f=0) T5: j=0;s=0; for i=3:50 m=i*(i-1)-1; if isprime(m)=1 m 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3

10、页，共 7 页 - - - - - - - - - 4 s=s+m; j=j+1; end end disp(亲密数的对数：,num2str(j); disp(所有亲密素数之和：,num2str(s); 实验五：fn1: function f=fn1(n) f=n+10*log(n.2+5); fn2: function f=fn2(n) f=0; for i=1:n f=f+i*(i+1); end fx: function f=fx(x) f=1./(x-2).2+0.1)+1./(x-3).4+0.01); T1: x=input( 请输入一个复数：); fe=exp(x) flg=lo

11、g(x) fsin=sin(x) fcos=cos(x) T2: m1=input( 请输入 m1:); m2=input( 请输入 m2:); t=input( 请输入 theta:); A=m1*cos(t),-m1,-sin(t),0;m1*sin(t),0,. cos(t),0;0,m2,-sin(t),0;0,0,-cos(t),1 B=0,m1*9.8,0,m2*9.8 disp(a1,a2,N1,N2=); disp(AB); T3: j=1;n=1;k=1; for i=10:99 % 挑选出1099 中所有的素数if isprime(i)=1 m(j)=i; j=j+1; e

12、nd end for t=1:length(m); s(n)=10*rem(m(t),10)+fix(m(t)/10); %挑选出的素数进行位置交换if isprime(s(n)=1 %再挑选出交换绝对素数r(k)=m(t); k=k+1; end n=n+1; end disp(所有两位绝对素数：,num2str(r) T5: y1=fn1(40)/(fn1(30)+fn1(20) %函数调用y2=fn2(40)/(fn2(30)+fn2(20) 实验六：T1: x=linspace(0,2*pi,101); y=(0.5+3*sin(x)./(1+x.2).*cos(x); plot(x,

13、y); T2: x=linspace(0,2*pi,100); y1=x.2; y2=cos(2*x); y3=y1.*y2; subplot(2,2,1); plot(x,y1,r:,x,y2,y-,x,y3,b-.); subplot(2,2,2); plot(x,y1,r:); subplot(2,2,3); plot(x,y2,y-); subplot(2,2,4); plot(x,y3,b-.); T2-3: x=linspace(0,2*pi,10); y1=x.2; y2=cos(2*x); y3=y1.*y2; subplot(3,4,1); bar(x,y1,r);axis(

14、0,7,0,40); subplot(3,4,2); stairs(x,y1,y);axis(0,7,0,40); subplot(3,4,3); stem(x,y1,b);axis(0,7,0,40); subplot(3,4,4); fill(x,y1,g);axis(0,7,0,40); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 7 页 - - - - - - - - - 5 subplot(3,4,5); bar(x,y2,r);axis(-1,7,-1,1

15、); subplot(3,4,6); stairs(x,y2,y);axis(-1,7,-1,1); subplot(3,4,7); stem(x,y2,b);axis(-1,7,-1,1); subplot(3,4,8); fill(x,y2,g);axis(-1,7,-1,1); subplot(3,4,9); bar(x,y3,r);axis(0,7,-30,50); subplot(3,4,10); stairs(x,y3,y);axis(0,7,-30,50); subplot(3,4,11); stem(x,y3,b);axis(0,7,-30,50); subplot(3,4,1

16、2); fill(x,y3,g);axis(0,7,-30,50); T3: x=linspace(-5,5,100); y=; for x0=x; if x0=0; y=y,(x0+sqrt(pi)/exp(2); else y=y,1 /2*log(x0+sqrt(1+x02); end end subplot(1,2,1); plot(x,y); for i=1:length(x) if x(i)0.5); b=length(k); disp(大于0.5 的随机数个数占总数的百分比： ); sprintf(%f%, b/ 300) T2: P=45+50*rand(100,5) disp

17、(最高分及序号：); Y,U=max(P,1) disp(最低分及序号：); Y1,U1=min(P,1) disp(每门课平均分：); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 7 页 - - - - - - - - - 6 mean(P,1) disp(每门课的标准方差：); std(P,0,1) s=sum(P ,2); disp(总分最高分：); Y2,U2=max(s,1) disp(总分最低分：); Y3,U3=min(s,1) zcj,xsxh=sor

18、t(s,1,descend) %descend 降序T3: h=6:2:18; t1=18.0,20.0,22.0,25.0,30.0,28.0,24.0; t2=15.0,19.0,24.0,28.0,34.0,32.0,30.0; t=6.5:2:17.5; T1=interp1(h,t1,t,spline) T2=interp1(h,t2,t,spline) T4: x=1:10:101; y=log10(x); p=polyfit(x,y,5) x1=1:0.01:101; y1=log10(x1); y2=polyval(p,x1); plot(x1,y1,r-,x1,y2,b-.)

19、; T5: P1=1,2,4,0,5; P2=0,0,0,1,2; P3=0,0,1,2,3; P4=conv(P2,P3); sp1=length(P1); sp2=length(P4); P=zeros(1,sp2-sp1),P1+P4 x=roots(P) A=-1,1.2,-1.4;0.75,2,3.5;0,5,2.5; Ax=polyval(P,A) AA=polyvalm(P,A) 实验九：T1: f=inline(det(x x2 x3;1 2*x 3*x2;0 2 6*x); i=1;x=1; while x=3.01 g(i)=f(x); i=i+1; x=x+0.01; e

20、nd dx=diff(g)/0.01; dx(1) dx(101) dx(length(g)-1) T2: f1=inline(sqrt(cos(t.2)+4*sin(2*t).2+1) f2=inline(log(1+x)./(1+x.2) I1=quad(f1,0,2*pi) I2=quad(f2,0,2*pi) T3: A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2; B=-4,13,1,11; x1=AB L,U=lu(A); x2=U(LB) Q,R=qr(A); x3=R(QB) fun: function F=fun(X); x=X(1); y=X(

21、2); z=X(3); F(1)=sin(x)+y2+log(z)-7; F(2)=3*x+2y-z3+1; F(3)=x+y+z-5; T5: f=inline(3*x+sin(x)-exp(x); fzero(f,1.5) X=fsolve(fun,1 1 1,optimset(Display,off) 实验十：T1: x=sym(6); y=sym(5); z=(x+1)/(sqrt(3+x)-sqrt(y) T2: syms x y; t=sym(5135); factor(x4-y4) factor(t) T3: syms beta1 beta2 x; 名师资料总结 - - -精品资

22、料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 7 页 - - - - - - - - - 7 simplify(sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2) simplify(4*x2+8*x+3)/(2*x+1) T4: syms a b c d e f g h k; p1=0,1,0;1,0,0;0,0,1; p2=1,0,0;0,1,0;1,0,1; A=a,b,c;d,e,f;g,h,k; B=p1*p2*A B1=inv(B) B1*B; tril

23、(B) det(B) T5: syms x t a y f1=(x*(exp(sin(x)+1)-2*(exp(tan(x)-1)/sin(x)3; limit(f1,x,0) %(1) f2=(sqrt(pi)-sqrt(acos(x)/sqrt(x+1); limit(f2,x,-1,right) %(2) f3=(1-cos(2*x)/x; diff(f3,x,1) %(3) diff(f3,x,2) A=ax,t3;t*cos(x),log(x); diff(A,x,1) &(4) diff(A,t,2) diff(diff(A,x,1),t,1) f5=(x2-2*x)*exp(-x

24、2-y2-x*y); yx=-diff(f5,x,1)/diff(f5,y,1) %(5) fxy=diff(diff(f5,x,1),y,1) x=sym(0);y=sym(1); eval(fxy) T6: x=sym(x); f1=1/(1+x4+x8); int(f1,x) f2=1/(asin(x)2/sqrt(1-x2); int(f2,x) f3=(x2+1)/(x4+1); int(f3,x,0,inf) f4=exp(x)*(1+exp(x)2; int(f4,x,0,log(2) 实验十一：T1：syms n x ; S1=symsum(1 /(2*n-1),n,1,10)

25、 S2=symsum(n2*x(n-1),n,1,inf) S3=symsum(n2/5n,n,1,inf) T2：x=sym(x); f=log(x); taylor(f,x,6,1) T3：x1=solve(log(1+x)-5/(1+sin(x)=2) x2=solve(x2+9*sqrt(x+1)-1=0) x3=solve(3*x*exp(x)+5*sin(x)-78.5=0) x4,y4=solve(sqrt(x2+y2)-100=0,3*x+5*y-8=0) T4：dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x) T5：x,y,z=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z,t) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 7 页 - - - - - - - - -