大学数学实验课后习题答案.docVIP

实验名称：MATLAB程序设计 (1)作马鞍面： 程序: 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）自定义函数：，并求 程序：function y=fun(x); y=log(sin(x))-cos(x)*log(tan(x)); fun(pi/3) （4）P441第10题的（1）、（2）小题。要求建立函数M文件求解。并求： 程序1：求自定义函数 function y=fun(n) A=1; for k=1:n A=A*k; end A 程序2：求： s=0; for n=1:20 A=1; for k=1:n A=A*k; end s=s+A; end s 程序3：求 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