Matlab程序設(shè)計與應(yīng)用第二版劉衛(wèi)國課后實驗答案

1、Matlab 程序設(shè)計與應(yīng)用第二版劉衛(wèi)國課后實驗答案 實驗一: 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+x42); T4-2: z2 a=input( 請輸入一個字符串 :,s); % 第三小題 b=find(a=A&a=0&t=1 &t=2&tb disp(A43=); disp( 帕薩卡矩陣 P 性能更好);disp(A43)

2、; 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(系數(shù)矩陣 A 的條件數(shù):);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=sq

3、rtm(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)A2-5*x(i)+6; w=120*84+(t-120)*1

4、.15*84; else elseif t60 y(i)=x(i)A2-x(i)-1; w=t*84-700; end else end w=84*t; y end T2-if: disp( 應(yīng)發(fā)工資為 :); 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=i nput(請輸入 + 或-

5、或 * 或 /,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īng)的 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;

6、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(bA2+4*a)/2; r2=(-b-sqrt(bA2+4*a)/2; while (n 5 ; 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=l

7、ength(find(f0) T1-2: f0=length(find(f=0) n=input( 請輸入 n:); T5: m=1./(1:n).A2; 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; 挑選出的素數(shù)進行位置交換 if isprime(s(n)=1 % 再挑選出交換 end end 絕對素數(shù) disp( 親密數(shù)的對數(shù) :,num2str(j)

8、; r(k)=m(t); disp( 所有親密素數(shù)之和 :,num2str(s); k=k+1; 實驗五 : end fn1: n=n+1; 實驗四: function f=fn1(n) end 所有兩位絕對素數(shù) :,num2str(r) fn2: T5: function f=fn2(n) y1=fn1(40)/(fn1(30)+fn1(20) % 函數(shù)調(diào)用 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

9、.A2).*cos(x); function f=fx(x) plot(x,y); f=1./(x-2).A2+0.1)+1./(x-3).A4+0.01); T2: T1: x=linspace(0,2*pi,100); x=input( 請輸入一個復(fù)數(shù) :); y1=x.A2; 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:); p

10、lot(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.A2; T3: y2=cos(2*x)

11、; j=1;n=1;k=1; y3=y1.*y2; for i=10:99 % 挑選出 1 099 中所 subplot(3,4,1); 有的素數(shù) 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); f=n+10*log( n. A2+5); disp( end stem(x,y1,b);axis(0,7,0,40); end subplot(3,4,4); for t=1:length(m); fill

12、(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(

13、x,y); fill(x,y2,g);axis(-1,7,-1,1); z=cos(x).*cos(y).*exp(- sqrt(x.A2+y.A2)/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;

14、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(xO+sqrt(1+xOA2); A=ra nd(1,30000); end disp(

15、均值為 :); end mean(A) subplot(1,2,1); disp( 標(biāo)準(zhǔn)方差為 :); 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 的隨機數(shù)個數(shù)占總數(shù)的百分 y(i)=1/2*log(x(i)+sqrt(1+x(i)A2); 比:); end sprintf(%f%, b/300) end T2: subplot(1,

16、2,2); P=45+50*rand(100,5) plot(x,y); disp( b=input( 請輸入 b:); Y1,U1=min(P,1) n=input( 門課平均分 :); mean(P,1) i=i+1; disp( 每門課的標(biāo)準(zhǔn)方差 :); 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,d

17、escend) %descend f1=i nli ne(sqrt(cos(t42)+4*si n(2*t)42+1) 最高分及序號 :); T4: Y,U=max(P,1) a=input( 請輸入 a:); disp( 最低分及序號 :); 請輸入 n:); disp( 每 序 f2=i nli ne(log(1+x)丿(1+X.A2) 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,

18、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

19、(2); plot(x1,y1,r-,x1,y2,b-.); z=X(3); T5: F(1)=si n(x)+yA2+log(z)-7; P 仁1,2,4,0,5; F(2)=3*x+2Ay-zA3+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),P

20、1+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=inlin e(det(x xA2 xA3;1 2*x 3*xA2;0 2 t=sym(5135); 6*x); factor(xA4-yA4) i=1;x=1; factor(t) while x=3.01 T3: g(i)=f(x); syms beta1 beta2

21、 x; 6 simplify(sin(beta1)*cos(beta2)-cos(beta1)*sinsyms n x ; (beta2) S1=symsum(1/(2*n-1),n,1,10) S2=symsum(nA2*xA(n-1),n,1,inf) simplify(4*xA2+8*x+3)/(2*x+1) T4: S3=symsum(nA2/5An,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(xA2+9*sqrt(x+1)-1=0) tril(B) x3=solve(3*x*exp(x)+5*sin(x)- 78.5=0) det(B) x4,y4=sol