实验名称：MA TLAB 程序设计
(1)作马鞍面：z =x 2y 2
2-3, -6≤x ≤6, -8≤y ≤8
程序: 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）自定义函数：y =lnsin x -cos x ln tan x ，并求y (π
3) =?
程序：function y=fun(x);
y=log(sin(x))-cos(x)*log(tan(x));
>>fun(pi/3)
（4）P441第10题的（1）、（2）小题。要求建立函数M 文件求解。并求：
程序1：求n ! 自定义函数
function y=fun(n)
A=1;
for k=1:n
A=A*k;
end
A
20
程序2：求：T =∑n !
n =1
s=0;
for n=1:20
A=1; 20T =∑n ! n =1
for k=1:n
A=A*k;
end
s=s+A;
end
s
程序3：求C 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,) 之间的图象. 的图象. n
程序：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
实验名称：MA TLAB 程序设计
(1)作马鞍面：z =x 2y 2
2-3, -6≤x ≤6, -8≤y ≤8
程序: 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）自定义函数：y =lnsin x -cos x ln tan x ，并求y (π
3) =?
程序：function y=fun(x);
y=log(sin(x))-cos(x)*log(tan(x));
>>fun(pi/3)
（4）P441第10题的（1）、（2）小题。要求建立函数M 文件求解。并求：
程序1：求n ! 自定义函数
function y=fun(n)
A=1;
for k=1:n
A=A*k;
end
A
20
程序2：求：T =∑n !
n =1
s=0;
for n=1:20
A=1; 20T =∑n ! n =1
for k=1:n
A=A*k;
end
s=s+A;
end
s
程序3：求C 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,) 之间的图象. 的图象. n
程序：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