MATLAB程序設(shè)計教程(第二版)課后答案

1、MATLAB第二版課后答案unit3-8unit3請輸入一個三位數(shù):n=inpu t('a=fix(n/100);b=fix(n-a*100)/10);c=n-a*100-b*10;d=c*100+bT0+a)2( 1');請輸入成績 n=input(' 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'end p rice=p請

2、輸入成績）（2 n=input（' elseif n>=7 0&nv=79');if n>=90&nv=100p='A'elseif n>=80&nv=89p 二B;P=C;elseif n>=6 0&nv=69p='D'elseP='E'endp rice=pn;trycatchp rice='erroe'end n=1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908, 6;a=n(1);b=n(1); for m

3、=2:20if n(m)>aa=n( m);elseif n(m)vbb=n(m); end endmax=a min=b 法 2n=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) 4b=-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); end y 5 y1=0;y2=1;n=input(' 請輸入 n 的值： ');for i=1:nendy1=y1+

4、1/i2;y2=y2*(4*i*i)/(2*i-1)*(2*i+1);y1 y2 6A=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(' 請輸入 n 的值 :');if n<=5&n>=0dis p(A( n,:);elseif n<0 dis p(lasterr); else dis p(A(6,:);disp(lasterr);end 7 ( 1)f=;for n=1:40f(n)=n+10*log(n2+5);endy=f(40)/(

5、f(30)+f(20)f=;a=0; for n=1:40f(n)二a+n*(n+1);a=f( n);endy=f(40)/(f(30)+f(20)y=0;:');輸入 m 的值 m=input('');值：輸入 n=input('n for i=1:ny=y+im; end y*s=0;* function s=shi8_1(n,m)for i=1:ns=s+im;end*shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)考練習2N=1,2,3,4,5;2.*NN./21./N1./N.23s=fix(100*rand(1,2

6、0)*9/10+10) y=sum(s)/20j=0; for i=1:20if s(i)<y&rem(s(i),2)=0j=j+1;A(j)=s(i); else continue; endend2t=-pi:p i/10: pi;y=1./(1+ex p(-t);y1=0;y2=0;n=input（' 請輸入 n的值：'）;for i=1:ny仁y1+-(-1)Ai/(2*i-1);y2=y2+1/4i;end yiy2 unit4 實驗指導(dǎo) （1 1x=-10:0.05:10;y=x-x.A3./6;plot(x,y)x=-10:0.5:10;ezpl ot

7、('xA2+2*yA2-64',-8,8);grid on;sub plot(2,2,1);bar(t,y);(t,y)'); title(' 條形圖axis(卜 pi,p i,0,1);sub plot(2,2,2);stairs(t,y,'b');(t,y)');階梯圖 title('axis(-pi,p i,0,1);sub plot(2,2,3); stem(t,y,'k');(t,y)'); title(' 桿圖axis(-pi,p i,0,1);sub plot(2,2,4);loglo

8、g(t,y,'y');(t,y)'); title('對數(shù)坐標圖t=0: pi/50:2* pi;r=5.*cos(t)+4;grid on; ( 2)x,y=meshgrid(-5:5);polar(t,r);title('rho=5*costheta+4');t=-p i/3: pi/50: pi/3;r=5.*(sin(t)A2)./cos(t);polar(t,r); 4 ( 1 )t=0: pi/50:2* pi;x=ex p(-t./20).*cos(t);y=ex p(-t./20).*sin(t);z=t;plot3(x,y,z)

9、;axis(-i,i,-i,i,-i,i)m(:,i)=getframe;end movie(m,4);思考練習)(2 1 x=-5:0.1:5;z=zeros(11)+5; mesh(x,y,z); shading inter p;x,y,z=s phere(20);surf(x,y,z); axis off;shading inter p;m=moviein(20); for i=1:20y=(1./(2* pi).*ex 卩(心八2)/2);Plot(x,y);(2)t=-2* pi:0.1:2* pi;x=t.*sin(t);y=t.*cos(t); p lot(x,y);grid o

10、n;t=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)v1e-2);x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');x=-2:0.01:2;y=sin(1./x);sub plot(2,1,1);plot(x,y);sub plot(2,1,2);fplot('sin(1./x) ',-2,2,1e-4);i=-4* pi:0.1:10;j=

11、12./sqrt(i);polar(i,j);title('rho=12/sqrt(theta)')a=-p i/6:0.01: pi/6;b=3.*asin(a).*cos (a)./(sin(a).3+(cos (a).3);polar(a,b); 6 ( 1 )u,v=meshgrid(-4:0.1:4);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4.*u.2;sub pl ot(2,1,1);mesh(x,y,z);sub plot(2,1,2); surf(x,y,z);x,y=meshgrid(-3:0.2:3);z=-5./(1+x.2+y.

12、2);sub plot(1,2,1);mesh(x,y,z);sub pl ot(1,2,2);surf(x,y,z);unit5實驗指導(dǎo)1A=randn(10,5) x=mean(A) y=std(A)Max=max(max(A)Min=min(min(A)Sumhang=sum(A,2)SumA=sum(Sumhang)B=sort(A);C=sort(B,2,'descend');2 ( 1 )a=0:15:90;b=a./180.* pi;s=sin(b)c=0:15:75;d=c./180.* pi;t=tan(d)');e=in put('請輸入想計

13、算的值:S=sin(e/180* pi)T=tan(e/180* pi)S1=interp1(a,s,e,'s pline')T1=inter p1(c,te's pline')P1= polyfit(a,s,5);P 2=polyfit(c,t,5);S2=polyval( P1,e)T2=polyval (P 2,e) n=1,9,16,25,36,49,64,81,100;N=sqrt (n);x=inpu t('ji suan zhi :');interp1(n,N,x,'cubic')3N=64;T=5;t=lins p

14、ace(O,T,N);h=exp(-t);dt=t (2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;p lot(f,abs(F),'-*')4P=2,-3,0,5,13;Q=1,5,8;p=p olyder( P)q=p olyder( P,Q)a,b=polyder(P ,Q)P 1=1,2,4,0,5; P 2=0,1,2;P 3=1,2,3;P=P 1+conv( P2, P3)X=roots (P)p=p olyval( P ,A)A=-1,12-1.4;0.75,2,3.5;0,5,2.5;思考練習A=rand(

15、1,30000);a=mean(A)b=std(A)Max=max(A)Min=min(A) n=0;for i=1:30000if(A(i)>0.5)n=n+1; endendy=n/30000 5p=45,74,54,55,14;78,98,45,74,12;87,98,85,52,65M,S=max(p) N,H=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')x=1:10:101

16、;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= po lyval (p ,a); p Iot(x,y,':o'ay1,'-*')unit6實驗指導(dǎo)A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;p=0.95,0.67,0.52'x=ApIl tfl «IV WIIA=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;p=0.95,

17、0.67,0.53'x=Apcond(A)x1=fzero(funx1,-1)function fx=funx1(x)fx=x41+x3+1;x2=fzero(funx2,0.5)function fx=funx2(x)fx=x-sin(x)/x;op tions=op timset('Dis play','off);x=fsolve(fun3,1,1,1',o ptions)q=fun3(x)function q=fun3( p) x=p (1)；y=P( 2); z=p (3);q(1)=sin(x)+yA2+log (z)-7;q(2)=3*x+2

18、y-z3+1;q(3)=x+y+z-5;3 (1 )t0=0;tf=5;y0=1;t,y=ode23(fun4,t0,tf,y0);t'y'function yp=fun4(t,y) yp=-(1.2+sin(10*t)*y;t0=0;tf=5;y0=1;t't,y=ode23(fun5,t0,tf,y0);y'function yp=fun5(t,y)Il *iIlliABiiyp=cos(t)-y/(1+t2);x=fminbnd(mymin,0,2);4-mymin(x)function fx=mymin(x)fx=-(1+xA2"(1+xA4)

19、;options=optimset('Display','off');x,fval=fmincon(fun6,0,0,0,a,b,lb,ub)-fvalfunction f=fun6(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);思考練習 1 （ 1 ）A=2,3,5;3,7,4;1,-7,1;B=10,3,5'C1=inv(A)*B C2=ABL,U=lu(A);x=U

20、(LB)2）C2=AB L,U=lu(A);x=U(LB)2( 1 )A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;B=-4,13,1,11'C1=inv(A)*Bq(2)=x+y-1;3 ( 1 )x1=fzero(funx1,1.5)function fx=funx1(x)fx=3*x+s in( x)-ex p( x);x1=fzero(funx2,1)function fx=funx2(x)fx=x-1/x+5;op tions=op timset('Dis play','off);x=fsolve(fun3,3,0

21、9;,o ptions)q=fun3(x)function q=fun3( p)x=p (1);y=p(2);q(1)=x2+y2-9;to=o；tf=5;y0=0,1;t,y=ode45(vd pol,t0,tf,y0);t,yfunction ydot=vd pol(t,y);ydot (2)=y(1);ydot(1)=(2-3* y( 2)-2*t*y(1)./(1+t2);ydot=ydot't0=0;tf=5;y0=1;0;2;t,y=ode45(vd poll,t0,tf,y0);t,yfunction ydot=vd poll(t,y);ydot(1)=cos(t)-y

22、(3)./(3+sin(t)+5*y(1).*cos(2*t)/(t+1)A2)-y(2);ydot(2)=y(1);ydot(3)=y(2); ydot=ydot' 4x=fminbnd(mymin,O, pi);-mymin(x)function fx=mymin(x) fx=-sin(x)-cos(x.A2);x,y1=fminbnd(mymax,0,1.5);-y1function fx=mymax(x);fx=-(9*x+4*x.3-12*x.2);unit7實驗指導(dǎo)1 ( 1 )format longfx=inline('sin(x)./x');l,n=qu

23、adl(fx,0,2,1e-10)(2)format longfx=inline('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;325;l=dblquad(fxy,0,1,0,1)ki ( 2)I=dblquad(f,0, pi,0, pi)f=inline('abs(cos(x+y)','x','y');X=0.3:0.2:1.5;F=0.3895,0.6598,0.9147,1.1611,1.397

24、1,1.6212,1.8trap z(X,F)4p=0: pi/5:2* pi;for n=1:3endDX=diff(sin( p),n)f=inline('sin(x)./(x+cos(2.*x)');g=inline('(cos(x).*(x+cos(2*x)-sin(x).*(1-2.*sin(2*x)/(p=p olyfit(x,f(x),5);x+cos(2.*x).2'); x=-pi:0.01: pi;dp=po lyder( p);在假設(shè)點的函數(shù)% dp求dpx=po lyval(d p, x);值直接 %dx=diff(f(x,3.01)/對

25、f(x)求數(shù)值導(dǎo)數(shù)在假設(shè)點的導(dǎo)數(shù)g 的導(dǎo)函 f %求函數(shù) 0.01; gx=g(x);數(shù)作圖 plot(x,dpx,x,dx,'.',x,gx,'-')思考練習2format longfx=inline('1./(1+x.2)');l,n=quad(fx,-lnf,lnf,1e-10)l,n=quadl(fx,-lnf,lnf,1e-10)x=-100000:100000;y=1./(1+x.2);trap z(x,y)format short 3(1)format longfx=inline('log(1+x)丿(1+xA2)'

26、);l,n=quad(fx,0,1,1e-10)format longfx=inline('sqrt(cos(tA2)+4*sin(2*t)A2)+1)');I,n=quad(fx,0,2* pi ,1e-10)f=inline(4*x.*z.*ex p(-z.A2.*y-x.A2)');l=tri plequad(f,0, pi,0, pi ,0,1)x=0:0.01:2* pi;f=inline('sin(x)'); g=inline('cos(x)');p=p olyfit(x,f(x),5);dp=po lyder( p);dpx

27、=po lyval(d p, x);gx=g(x);dx=diff(f(x,2* pi+0.01)/0.01;plot(x,d px,x,dx,'.',x,gx,'-')unit8實驗指導(dǎo)1 syms x y;s=x4-y4; factor(s)factor(5135) 2 syms x; f=(x-2)/(x2-4);Iimit(f,x,2)sym x;f=(sqrt( pi)-sqrt(acos(x)/sqrt(x+1);limit(f,x,-1,'right')sym x;f=sin(1/x); diff(f,'x')dif

28、f(f,'x',2)sym x;f=(1-cos(2*x)/x;diff(f,'x') diff(f,'x',2)4 sym x;f=sqrt(ex p( x)+1);int(f,'x')syms x y;f=x/(x+y);int(f,y)sym x;f=ex p(x)*(1+ex p(x)A2;int(fX0,log(2)sym x; f=x*log(x);int(f,'x',1,ex p(1)5 sym x;s=symsum(-1)A(x+1)/x,1,lnf)sym y;z=symsum(yA(2*y-1)

29、/(2*y-1),1,lnf)6 sym x;f1=(ex p(x)+ex p(-x)/2;f2=sqrt(x3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0)syms x y a;x=solve('x3+a*x+1=0','x')x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')x y=solve('log(x/y)=9','ex p(x+y)=3','x',y)8 syms n;x,y=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0','x') 思考練習syms x B1 B2 a bs1=2*(cos(x)A2)-(sin(x)A2;s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);s3=sqrt(a+sq