大学数学实验课后习题答案(清华大学出版)

实验名称：MA TLAB 程序设计

(1)作马鞍面：22

,66,8823x y z x y =--≤≤-≤≤

程序: x=-6:0.5:6;y=-8:0.5:8

[X,Y]=meshgrid(x,y);

Z=X.^2./2-Y .^2./3;

mesh(X,Y ,Z)

(2)P441第5题

程序1:n=18;

I(1)=1-exp(-1);%I(1)对应I0

for k=1:n-1

I(k+1)=1-(k+1)*I(k);

end

I

程序2:n=18;

I1=(1/(n+1))*exp(-1);

I2=1/(n+1);

I(18)=(I1+I2)/2;

for k=n:-1:2

I(k-1)=(1-I(k))/n;

end

I

（3）自定义函数：lnsin cos ln tan y x x x =-，并求()?3y π

=

程序：function y=fun(x);

y=log(sin(x))-cos(x)*log(tan(x));

>>fun(pi/3)

（4）P441第10题的（1）、（2）小题。要求建立函数M 文件求解。并求：

201!n T n ==∑

程序1：求!n 自定义函数

function y=fun(n)

A=1;

for k=1:n

A=A*k;

end

A

程序2：求：20

1

!n T n ==∑

s=0;

for n=1:20

A=1;

for k=1:n

A=A*k;

end

s=s+A;

end

s

C

程序3：求n

m

function y=funa(n,m)

A=1;%求

for k=1:n

A=A*k;

end

B=1;

for k=1:m

B=B*k;

end

C=1;

for k=1:n-m

C=C*k;

end

D=A/(B*C) %求组合数

一元函数的图形练习解答：

1．用ezplot画出的图象.

程序：ezplot('asin(x)')

2．用ezplot画出用在(0,)之间的图象.

程序：ezplot('sec(x)',[0 pi])

3．在同一坐标系中画出,,,,的图象.并用gtext加以标记ezplot('sqrt(x)')

hold on

ezplot('x^2')

hold on

ezplot('x^(1/3)')

hold on

ezplot('x^3')

hold on

ezplot('x')

axis([-2 3 -2 2])

gtext('sqrt(x)')

gtext('x^2')

gtext('x^(1/3)')

gtext('x^3')

gtext('x')

4．画出及其反函数的图象.

x=-2:0.01:20;

y=1+log(x+2+eps);

plot(x,y)

hold

plot(y,x,'r')

axis([-4 4 -4 4])

8题：x=100;

y=50;

n=50;

r1=0.2;

r2=0.3;

a1=0.001;

a2=0.002;

for k=1:n

x(k+1)=(1+r1-a1*y(k))*x(k);

y(k+1)=(1-r2+a2*x(k))*y(k);

end

k=0:n;

round([k',x',y'])

plot(k,x,k,y),grid,

2题：function z=exf14(x0,y0,n,r,N,d,a,b);

x=x0;y=y0;

for k=1:n

x(k+1)=x(k)+r*(1-x(k)/N)*x(k)-a*y(k)*x(k)/N;

y(k+1)=(1-d+b*x(k)/N)*y(k);

end

z=[x',y'];

z=exf14(1000,100,100,0.8,3000,0.9,1.6,1.5);

k=0:100;

plot(k,z(:,1),k,z(:,2)),grid