MATLAB基础与应用教程 习题答案 作者 蔡旭晖 刘卫国 蔡立燕 第1-8章答案

第1章

一、思考题

4. （1） B=A(2:5,1:2:5)

（2）A(7)=[]

（3） A=A+30

（4）size(A)

（5） t(find(t==0))=eps

（6）t=reshape(x,3,4)

（7）abs('matlab')

（8）char(93)

5. A=[97,67,34,10;-78,75,65,5;32,5,-23,-59;0,-12,54,7] (1) B=A(1:3,:)

C=A(:,1:2)

D=A(2:4,3:4)

E=B*C

(2)A(find(A>=50 & A<=100))

二、实验题

=-74/180*pi;y=-27/180*pi;

sin(x*x+y*y)/sqrt(tan(abs(x+y)))+pi

=::;

exp*a).*sin(a+

=[2,4;,5];

log(x+sqrt(1+x.*x))/2

4. A=[3,54,2;34,-45,7;87,90,15];B=[1,-2,67;2,8,74;9,3,0];（1）A*B

ans =

129 432 4197

7 -407 -1052

402 591 12489 A.*B

ans =

3 -108 134

68 -360 518

783 270 0

（2） A^3

ans =

-28917 240246 -4368 137883 -259101 27669 171333 252504 38673 A.^3

ans =

27 157464 8 39304 -91125 343 658503 729000 3375（3）A/B

ans =

B/A

ans =

（4）[A,B]

ans =

3 5

4 2 1 -2 67

34 -45 7 2 8 74 87 90 15 9 3 0 [A([1,3],:);B^2]

ans =

3 5

4 2

87 90 15

600 183 -81

684 282 726

15 6 825

=1+2i;b=3+4i;

c=exp(pi*i/6)

c =

+

c+a*b/(a+b)

ans =

+

第2章

一、思考题

=0;

for n=0:63

s=s+2^n;

end

disp(s)

n=0:63;

s=sum(2.^n)

二、实验题

1.

x=input('输入一个四位整数：');

y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)] z=mod((y+7),10)

x=z(3)*1000+z(4)*100+z(1)*10+z(2)

2.

gh=input('输入工号');

h=input('输入工时');

dj=84;

if h>120

gz=dj*120+*dj*(h-120);

elseif h<60

gz=dj*h-700;

else

gz=dj*h;

end

format bank;

display([gh,gz])

3.

循环结构

n=input('input n:');

s=0;

for k=1:n

s=s+1/k^2;

end

display(sqrt(s*6))

向量运算

n=input('input n:');

k=1:n;

display(sqrt(sum(1./k.^2)*6))

4.

y=0;k=0;

while y<3

k=k+1;

y=y+1/(2*k-1);

end

display([k-1,y-1/(2*k-1)])

5.

x0=0;x=1;k=0;

a=input('a=');

b=input('b=');

while abs(x-x0)>=1e-5 && k<500

x0=x;

x=a/(b+x0);

k=k+1;

end

display([k,x]);

display([(-b+sqrt(b^2+4*a))/2,(-b-sqrt(b^2+4*a))/2]); 6.

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

（1）函数文件

function f=fun(n)

f=n+log(n^2+5);

（2）函数文件

function f=fun(n)

a=1:n;

f=sum(a.*(a+1));

第3章

一、思考题

4.

t=0::;

y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3);

5.

x=-10::10;

y=linspace(-6,6,size(x,2))

z=x.^3+3*x.*y.^2;

plot3(x,y,z)

6.

x=100:100:400;

y=100:100:400;

z=[636,697,624,478;698,712,630,478;680,674,598,412;662,626,552,334]; [X,Y]=meshgrid(x,y);

mesh(X,Y,z)

二、实验题

1.

（1）x=-10::10;

plot(x,x-x.^3/3/2)

（2）plot(x,exp(-x.*x/2)/2/pi)

（3）x=-8::8;

plot(x,sqrt((64-x.*x)/2))

（4）t=0:pi/100:6*pi;

plot(t.*sin(t),t.*cos(t))

2. （1）

x1=linspace(0,1,100);

y1=2*;

t=linspace(0,pi,100);

x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);

plot(x1,y1,'r-',x,y,'b:');

text,1,'y=');

text,,'x=sin(3t)cos(t)');

text,,'y=sin(3t)sin(t)');