MATLAB程序设计教程课后答案.doc

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)

b=n(m);

end

end

max=a

min=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=[::];

for n=1:61

a=b(n);

y(n)=(exp*a)-exp*a))/2*sin(a++log(+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)

j=j+1;

A(j)=s(i);

else continue;

end

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::10;

y=x-x.^3./6;

plot(x,y)

（2）

x=-10::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）