北京科技大学数学实验第二次实验二维绘图报告.docVIP

姓名+学号第一次实验报告.doc 《数学实验》报告 实验名称 二维绘图 学 院 计通学院 专业班级????通信15xx 姓??? 名??? 刘xx 学??? 号?4152xxxx? ?? 2017年5月 【实验目的】 掌握二维绘图基础知识P79 1,3,5；及补充习题：Matlab 在重积分中的应用 【实验程序】 1.绘制y=e^pi/3sin(3x)(x属于[0,4*pi])的图像，要求用蓝色的星号画图；并且画出其包络线y=+-e^pi/3的图像，用红色的点划线画图。 x=0:pi/50:4*pi; y=exp(x./3).*sin(3*x); plot(x,y,b*); hold on; y1=exp(x./3); y2=-exp(x./3); plot(x,y1,r.-,x,y2,r.-); 实验结果 2.在同一图形窗口画三个子图，要求使用指令gtext、axis、legend、title、xlabel和ylabel： （1）y=xcosx，-pi=x=pi; （2）y=xtan1/xsinx^3,pi=x=4pi; （3）y=e^1/xsinx,1=x=8; x1=-pi:pi/50:pi; y1=x1.*cos(x1); x2=pi:pi/50:4*pi; y2=x2.*tan(1./x2).*sin(x2.^3); x3=1:0.1:8; y3=exp(1./x3).*sin(x3); xlabel(x?á); ylabel(y?á); subplot(1,3,1),plot(x1,y1,r*),grid on,title(y1=x1.*cos(x1)),gtext(y1=x1.*cos(x1)), legend(y1=x1.*cos(x1)),grid on; subplot(1,3,2),plot(x2,y2,k.-),grid on,title(y2=x2.*tan(1./x2).*sin(x2.^3)), gtext(y2=x2.*tan(1./x2).*sin(x2.^3)),legend(y2=x2.*tan(1./x2).*sin(x2.^3)), grid on; subplot(1,3,3),plot(x3,y3,bo),grid on,title(y3=exp(1./x).*sin(x)), gtext(y3=exp(1./x).*sin(x)),legend(y3=exp(1./x).*sin(x)),grid on; 结果： 3.绘制圆锥曲线的图像并加各种标注，圆锥螺线的参数方程为： X=t cos ( pi/6*t) Y=ts sin (pi/6*t) Z=2t 其中0=t=20*pi; t=0:pi/50:20*pi; x=t.*cos(pi/6.*t); y=t.*sin(pi/6.*t); z=2.*t; plot3(x,y,z); xlabel(X?á); ylabel(Y?á); zlabel(Z?á); title(?2×??ú??í???); 结果： 第4题补充习题：Matlab 在重积分中的应用 4. 交换 的 D2： 程序： y1=0:0.1:2; y2=-2:0.1:0; y3=-2:0.01:2; y4=-2:0.1:2; x=-1:0.01:3; x1=(y1.*y1-1)/2; x2=(y2.*y2-1)/2; x3=3/2*(y3./y3); x4=(16-y4.*y4)/8; y=0*x; plot(x1,y1,g,x2,y2,g,x3,y3,r,x4,y4,r); hold on; plot(x,y,--); axis([-0.6,2.1,-2.1,2.1]); gtext(y= sqrt(2x+1)); gtext(y= -sqrt(2x+1)); gtext(x=3/2); gtext(y= sqrt(16-8x)); gtext(y= -sqrt(16-8x)); 结果： 7