(完整word版)MATLAB程序設(shè)計(jì)教程(第二版)課后答案

1、MATLAB 第二版課后答案 unit3-8unit3實(shí)驗(yàn)指導(dǎo)1、 n=input( 請(qǐng)輸入一個(gè)三位數(shù)： ); a=fix(n/100);b=fix(n-a*100)/10); c=n-a*100-b*10; d=c*100+b*10+a2(1)n=input( 請(qǐng)輸入成績(jī) );switch ncase num2cell(90:100)p=A;case num2cell(80:89)p=B;case num2cell(70:79)p=C;case num2cell(60:69)p=D; otherwisep=E;endprice=p( 2) n=input( 請(qǐng)輸入成績(jī) );if n=90&n

2、=80&n=70&n=60&naa=n(m);elseif n(m)bb=n(m);endendmax=amin=b法2 n=1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6;min=min(n)max=max(n)4 b=-3.0:0.1:3.0;for n=1:61a=b(n);y(n)=(exp(0.3*a)-exp(-0.3*a)/2*sin(a+0.3)+log(0.3+a)/2);endy5y1=0;y2=1;n=input( 請(qǐng)輸入 n 的值： );for i=1:ny仁 y1+1/2;y2=y2*(4*i*i)/(2*i-1)*

3、(2*i+1);endy1y26A=1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6; n=input(請(qǐng)輸入n的值:);if n=0disp(A(n,:);elseif n0disp(lasterr);else disp(A(6,:);disp(lasterr);end7( 1 )f=;for n=1:40f(n )=n+10*log( n2+5);endy=f(40)/(f(30)+f(20)(2)f=;a=0;for n=1:40 f(n)=a+n*(n+1); a=f(n);endy=f(40

4、)/(f(30)+f(20)8y=0;m=input(輸入 m 的值:); n=input( 輸入 n 值： );for i=1:n y=y+iAm;endy*function s=shi8_1(n,m)s=0;for i=1:n s=s+iAm;end*shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2) 思考練習(xí)2N=1,2,3,4,5;2.*NN./21./N1./N.A23s=fix(100*rand(1,20)*9/10+10) y=sum(s)/20 j=0;for i=1:20if s(i)y&rem(s(i),2)=0j=j+1;A(j)=s(i)

5、;else continue;endendA4y1=0;y2=0;n=input( 請(qǐng)輸入 n 的值： ); for i=1:ny仁 y1+-(-1)9/(2*i-1); y2=y2+1/49;endy1y2unit4 實(shí)驗(yàn)指導(dǎo)1(1)x=-10:0.05:10;y=x-x.A3./6;plot(x,y)(2)x=-10:0.5:10;ezplot(xA2+2*yA2-64,-8,8);grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t);subplot(2,2,1);bar(t,y);title( 條形圖 (t,y);axis(-pi,pi,0,1);subplot

6、(2,2,2);stairs(t,y,b);title( 階梯圖 (t,y);axis(-pi,pi,0,1);subplot(2,2,3);stem(t,y,k);title(桿圖(t,y);axis(-pi,pi,0,1);subplot(2,2,4);loglog(t,y,y);title( 對(duì)數(shù)坐標(biāo)圖 (t,y);3( 1 )t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title(rho=5*costheta+4);(2)t=-pi/3:pi/50:pi/3;r=5.*(s in (t)A2)./cos(t);polar(t,r);4(1)t=0:p

7、i/50:2*pi;x=exp(-t./20).*cos(t); y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);grid on;(2) x,y=meshgrid(-5:5);z=zeros(11)+5;mesh(x,y,z);shading interp;5x,y,z=sphere(20); surf(x,y,z);axis off;shading interp; m=moviein(20);for i=1:20axis(-i,i,-i,i,-i,i) m(:,i)=getframe;endmovie(m,4); 思考練習(xí) 2(1) x=-5:0.1:5;y=(

8、1./(2*pi).*exp(-(x.A2)/2); plot(x,y);(2)t=-2*pi:0.1:2*pi; x=t.*sin(t); y=t.*cos(t);plot(x,y);grid on;3t=0:pi/1000:pi;x=sin(3.*t).*cos(t);y1=sin(3.*t).*sin(t);y2=2.*x-0.5;plot(x,y1,k,x,y2);hold on;k=find(abs(y1-y2)0.5)n=n+1;endendny=n/300005p=45,74,54,55,14;78,98,45,74,12;87,98,85,52,65M,S=max(p)N,H=

9、min(p)junzhi=mean(p,1) fangcha=std(p,1,1) zong=sum(p,2);Max,wei=max(zong)Min,wei=min(zong)zcj,xsxh=sort(zong,descend)6x=1:10:101;y=0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043; p,s=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,:o,a,y1,-*)unit6實(shí)驗(yàn)指導(dǎo)1A=1/2,1/3,1/4;1/3,1/4,

10、1/5;1/4,1/5,1/6;p=0.95,0.67,0.52;x=ApA=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;p=0.95,0.67,0.53;x=Apcon d(A)2 (1)x1=fzero(fu nx1,-1)fun ctio n fx=fu nx1(x)fx=xA41+xA3+1;(2)x2=fzero(fu nx2,0.5)fun ctio n fx=fun x2(x)fx=x-s in( x)/x;(3)optio ns=optimset(Display,off);x=fsolve(fu n3,1,1,1,optio ns)q=fu n3(x

11、)function q=fun 3(p)x=p(1);y=p(2);z=p(3);q(1)=si n(x)+yA2+log (z)-7;q(2)=3*x+2Ay-zA3+1;q(3)=x+y+z-5;3 (1)t0=0;tf=5;y0=1;t,y=ode23(fu n4,t0,tf,y0);tyfunction yp=fu n4(t,y)yp=-(1.2+si n(10*t)*y;(2)t0=0;tf=5;y0=1;t,y=ode23(fu n5,t0,tf,y0);tyfunction yp=fu n5(t,y)yp=cos(t)-y/(1+L2);4x=fm inbn d(m ymin,

12、0,2);-mymi n(x)fun ctio n fx=mymi n(x)fx=-(1+x.A2)/(1+x.A4);5optio ns=optimset(Display,off);x,fval=fmi neon (fu n6,0,0,0,a,b,lb,ub)-fvalfun ctio n f=fun 6(x)f=-(sqrt(x(1)+(400-x(1)*1.1+(sqrt(x(2)+(400-x(1)*1.1-x(2)*1.1+sqrt(3)+(400-x(1)*1.1-x(2)*1.1-x (3)*1.1+sqrt(x(x4);思考練習(xí)1 (1)A=2,3,5;3,7,4;1,-7,1

13、;B=10,3,5;C1=i nv(A)*BC2=ABL,U=lu(A);x=U(LB)(2)A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;B=-4,13,1,11;C1=i nv(A)*BC2=ABL,U=lu(A);x=U(LB)2 (1)x仁 fzero(fu nx1,1.5)fun ctio n fx=fu nx1(x)fx=3*x+si n( x)-exp(x);(2)x仁 fzero(fu nx2,1)fun ctio n fx=fun x2(x)fx=x-1/x+5;(3)optio ns=optimset(Display,off);x=fsolv

14、e(fu n3,3,0,optio ns)q=fu n3(x)function q=fun 3(p)x=p(1);y=p(2);q(1)=xA2+yA2-9;q(2)=x+y-1;3 (1)to=o；tf=5;yo=o,i;t,y=ode45(vdpol,t0,tf,y0);t,yfunction ydot=vdpol(t,y);ydot(1)=(2-3*y(2)-2*t*y(1)./(1+tA2);ydot(2)=y(1);ydot=ydot;(2)t0=0;tf=5;y0=i;0;2;t,y=ode45(vdpoll,t0,tf,y0);t,yfunction ydot=vdpoll(t,

15、y);ydot(1)=cos(t)-y(3)./(3+si n( t)+5*y(1).*cos(2*t)/(t+1).A2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot;4x=fmi nbn d(mymi n,0,pi);-mymi n(x)fun ctio n fx=mymi n(x)fx=-s in( x)-cos(x.A2);5x,y1=fmi nbn d(mymax,0,1.5);-yifun ctio n fx=mymax(x);fx=-(9*x+4*x.A3-12*x.A2);unit7實(shí)驗(yàn)指導(dǎo)1 (1)format lo ngfx=i nl

16、i ne(s in (x)./x);I, n=quadl(fx,0,2,1e-10)(2)format lo ngfx=i nlin e(1./(x-0.3).A2+0.01)-1./(x-0.9).A2+0.04)-6);I,n=quad(fx,0,1,1e-10)2(1)global ki;ki=0;I=dblquad(fxy,0,1,0,1)ki(2)f=inline(abs(cos(x+y),x,y);I=dblquad(f,0,pi,0,pi)3X=0.3:0.2:1.5;F=0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325;trapz

17、(X,F)4p=0:pi/5:2*pi;for n=1:3nDX=diff(sin(p),n)end5f=inline(sin(x)./(x+cos(2.*x);g=i nlin e(cos(x).*(x+cos(2*x)-si n(x).*(1-2.*si n(2*x)/(x+cos(2.*x).A2); x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x);%求 dp 在假設(shè)點(diǎn)的函數(shù)值dx=diff(f(x,3.01)/0.01;%直接對(duì) f(x)求數(shù)值導(dǎo)數(shù)gx=g(x);%求函數(shù)f的導(dǎo)函數(shù)g在假設(shè)點(diǎn)的導(dǎo)數(shù)plo

18、t(x,dpx,x,dx,.,x,gx,-);%作圖思考練習(xí)2format longfx=i nlin e(1./(1+x.A2);I,n=quad(fx,-Inf,Inf,1e-10)I,n=quadl(fx,-Inf,Inf,1e-10) x=-100000:100000;y=1./(1+x.A2);trapz(x,y)format short3(1)format longfx=inline(log(1+x)./(1+x.A2);I,n=quad(fx,0,1,1e-10)(2)format longfx=i nli ne(sqrt(cos(t.A2)+4*si n( (2*t)92)+1

19、);I,n=quad(fx,0,2*pi,1e-10)4f=inlin e(4.*x.*z.*exp(-z.A2.*y-x.A2);I=triplequad(f,0,pi,0,pi,0,1)5f=inline(sin(x); g=inline(cos(x);x=0:0.01:2*pi;p=polyfit(x,f(x),5); dp=polyder(p);dpx=polyval(dp,x);dx=diff(f(x,2*pi+0.01)/0.01; gx=g(x);plot(x,dpx,x,dx,.,x,gx,-)unit8實(shí)驗(yàn)指導(dǎo)1syms x y;s=xA4-yA4;factor(s)fact

20、or(5135)2syms x; f=(x-2)/(xA2-4);limit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)/sqrt(x+1); limit(f,x,-1,right)3sym x;f=sin(1/x); diff(f,x) diff(f,x,2) sym x;f=(1-cos(2*x)/x;diff(f,x)diff(f,x,2)4sym x;f=sqrt(exp(x)+1);int(f,x)syms x y;f=x/(x+y);int(f,y)sym x;f=exp(x)*(1+exp(x)A2;int(f,x,0,log(2)sym x;f=x

21、*log(x); int(f,x,1,exp(1)5sym x;s=symsum(-1)A(x+1)/x,1,l nf)sym y;z=symsum(yA(2*y-1)/(2*y-1),1,l nf)6sym x;f1=(exp(x)+exp(-x)/2;f2=sqrt(xA3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0)7syms x y a;x=solve(xA3+a*x+1=0,x)x=solve(sin(x)+2*cos(x)-sqrt(x)=0,x)x y=solve(log(x/y)=9,exp(x+y)=3,x,y)8syms n;x,y=dsolve(x*(D2y)+(1-n)*(Dy)+y=0,y(0)=0,Dy(0)=0,x) 思考練習(xí)2syms x B1 B2 a bs1=2*(cos(x)A2)-(si n(x)2; s2=sin(B1)*cos(B2)-cos(B1)*sin(B2); s3=sqrt(a+sqrt(aA2-b)/2)+sqrt(a-sqrt(aA2-b)/2); s4=(4*xA2+8*x+3)/(2*x+1);h1=simplify(s1)h2=simplify(s2)