2011秋数学实验基础 作业2 电子版.docVIP

注：图形须对坐标轴、标题、图例加适当的标注，步长、参数和自变量的范围的选取要合适 作数列： 的数据表（保留8位有效数字）,作图观察趋于无穷大时的极限情况。 n 5 2.4883200 2.1613940 -0.3269260 10 2.5937425 2.3871824 -0.2065600 50 2.6915880 2.6403310 -0.0512570 100 2.7048138 2.6784269 -0.0263869 500 2.7155685 2.7101644 -0.0054041 1000 2.7169239 2.7142138 -0.0027101 5000 2.7180101 2.7174667 -0.0005433 10000 2.7181459 2.7178742 -0.0002717 解：clear clc n=[5 10 50 100 500 1000 5000 10000]; a=(1+1./n).^n; b=(1+1./(n+1)).^n; c=b-a; fprintf(a=\n ); fprintf(%8.7f\n ,a); fprintf(b=\n ); fprintf(%8.7f\n ,b); fprintf(c=\n ); fprintf(%8.7f\n ,c); 2．作出平面曲线： （1） 解：x=-pi:pi/1000:pi; y=x.*sin(1./x); plot(x,y) （2）标准正态分布概率密度函数（standard normal distribution density function）: 解：x=-5:0.0005:5; y=1/sqrt(2*pi)*exp(-x.^2/2); plot(x,y) （3）在同一直角坐标系中作出三条曲线： 解：x=-pi:pi/1000:pi; y1=x.*sin(1./x); y2=x; y3=-x; plot(x,y1) hold on plot(x,y2) plot(x,y3) hold off （4）用 ‘subplot’ 命令在6个子窗口中分别画6条曲线: 解：x=-3*pi:pi/1000:3*pi; for i=1:6 y=sin(i*x); subplot(2,3,i) plot(x,y) end （5）用 ‘subplot’ 命令在4个子窗口中分别画4条曲线: 解：t=-pi:pi/1000:pi; for k=1:4 p=cos(k*t); subplot(1,4,k) polar(t,p) end （6）画出参数曲线： 解：t=-pi:pi/1000:pi; a=2; x=a*(cos(t)).^3; y=a*(sin(t)).^3; plot(x,y) 2．作出3-D 曲线: 解：t=-12*pi:pi/1000:12*pi; x=t.*sin(t); y=t.*cos(t); z=t; plot3(x,y,z) 3. 作出下列曲面: (mesh surf contour3 surfc) (1)(2) (3) (4) (5) (6) 解：（1）x=-1:0.01:1; y=-1:0.01:1; [X,Y]=meshgrid(x,y); Z=X.*Y; mesh(X,Y,Z) (2) x=-5:0.1:5; y=-5:0.1:5; [X,Y]=meshgrid(x,y); Z=X.^2/3+Y.^2/2; mesh(X,Y,Z) (3) x=-5:0.3:5; y=-5:0.3:5; [X,Y]=meshgrid(x,y); Z=sin(sqrt(X.^2+Y.^2))/sqrt(X.^2+Y.^2); mesh(X,Y,Z) (4) x=-5:0.3:5; y=-5:0.3:5; [X,Y]=meshgrid(x,y); Z=6+0*X; mesh(X,Y,Z) (5) x=-5:0.1:5; y=-5:0.1:5; [X,Y]=meshgrid(x,y); Z1=sqrt(X.^2+Y.^2); mesh(X,Y,Z1) hold on Z2=-sqrt(X.^2+Y.^2); mesh(X,Y,Z2) hold off (6) u=-5:0.1:5; v=-5:0.1:5; [U,V]=meshgrid(u,v); x=(3+0.3*cos(U)).*cos(V); y=(3+0.3*cos(U)).*sin(V); z=0.3*sin(V); mesh(x,y,z) 年级、专业 级信息与计算科学 姓名 学号 名单序号 实验时间 20年 9月30日