MATLAB程序设计教程（第二版）课后答案

MATLAB第二版课后答案unit3-8

unit3

实验指导

1、 n=input('请输入一个三位数：'); a=fix(n/100); b=fix((n-a*100)/10); c=n-a*100-b*10; d=c*100+b*10+a 2（1）

n=input('请输入成绩'); switch n

case num2cell(90:100) p='A';

case num2cell(80:89) p='B';

case num2cell(70:79) p='C';

case num2cell(60:69) p='D'; otherwise

p='E'; end

price=p

（2）n=input('请输入成绩'); if n>=90&n<=100 p='A';

elseif n>=80&n<=89 p='B';

elseif n>=70&n<=79 p='C';

elseif n>=60&n<=69 p='D'; else p='E'; end price=p （3）try n;

catch

price='erroe' end

3

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=2:20 if n(m)>a a=n(m); elseif n(m)

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:61 a=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:n

y1=y1+1/i^2;

y2=y2*((4*i*i)/((2*i-1)*(2*i+1))); end y1 y2

6

A=[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>=0 disp(A([n],:)); elseif n<0

disp(lasterr); else disp(A([6],:)); disp(lasterr); end

7（1） f=[];

for n=1:40

f(n)=n+10*log(n^2+5); end

y=f(40)/(f(30)+f(20))

（2）

f=[];a=0; for n=1:40

f(n)=a+n*(n+1); a=f(n); end

y=f(40)/(f(30)+f(20)) 8

y=0;

m=input('输入m的值:'); n=input('输入n值：'); for i=1:n y=y+i^m; end y

************************************************************ function s=shi8_1(n,m) s=0; for i=1:n

s=s+i^m; end

************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2) 思考练习 2

N=[1,2,3,4,5]; 2.*N N./2 1./N 1./N.^2

3

s=fix(100*rand(1,20)*9/10+10) y=sum(s)/20 j=0;

for i=1:20

if s(i)

end

A 4

y1=0;y2=0;

n=input('请输入n的值：'); for i=1:n

y1=y1+-(-1)^i/(2*i-1); y2=y2+1/4^i; end y1 y2 unit4 实验指导 1（1） x=-10:0.05:10; y=x-x.^3./6; plot(x,y) （2） x=-10:0.5:10;

ezplot('x^2+2*y^2-64',[-8,8]); grid on;

2

t=-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(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('对数坐标图(t,y)'); 3（1） t=0:pi/50:2*pi; r=5.*cos(t)+4;

polar(t,r);

title('\\rho=5*cos\\theta+4'); （2）

t=-pi/3:pi/50:pi/3; r=5.*((sin(t)).^2)./cos(t); polar(t,r); 4（1） t=0:pi/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; 5

[x,y,z]=sphere(20); surf(x,y,z); axis off; shading interp; m=moviein(20); for i=1:20

axis([-i,i,-i,i,-i,i]) m(:,i)=getframe; end

movie(m,4); 思考练习 2（1）

x=-5:0.1:5;

y=(1./(2*pi)).*exp((-(x.^2))/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; 3

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;