马鞍面等曲线曲面画图第四章

《数学实验》报告实验名称:姓名马鞍面等曲线曲面画图第四章学号2010年11月22一．【实验目的】掌握常微分方程求解和曲线拟合的方法，能够应用Matlab解决概率中的一些基本运算．二．【实验任务】79页359三．【实验程序】79页课后习题：3.在同一图形窗口画三个子图，要求使用指令gtext、legend、title、xlabel和ylabel：(1)y=xcosx,x(-,)(2)y=xtan(1/x)sinx3,x(,4)(3)y=e1/xsinx,x[1,8]5.绘制圆锥螺线的图像并加各种标注，圆锥螺线的参数方程为：x=tcos(/6*t)y=tsin(/6*t),(0=t=20)z=2t9.画三维曲面z=5-x2-y2(-1=x=1)与平面z=3的交线.一、【实验程序】79页的3题：x1=-pi:pi/50:pi;y1=x1.*cos(x1);x2=pi:pi/50:4*pi;y2=x2.*tan(1./x2).*sin(x2.^3);3x3=1:1/50:8;y3=exp(1./x3).*sin(x3);subplot(1,3,1),plot(x1,y1,'m.'),gridon,title('y1=xcosx')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:1/50:8;y3=exp(1./x3).*sin(x3);subplot(1,3,1),plot(x1,y1,'m.');legend('y1=xcosx'),gridon,title('y1=xcosx'),xlabel('x轴'),ylabel('y轴')subplot(1,3,2),plot(x2,y2,'rp');legend('y2=xtan(1/x)sin(x^3)'),gridon,title('y2=xtan(1/x)sin(x^3)'),xlabel('x轴'),ylabel('y轴')subplot(1,3,3),plot(x3,y3,'b-');legend('y3=exp(1/x)sinx'),gridon,title('y3=exp(1/x)sinx'),xlabel('x轴'),ylabel('y轴')gtext('y1=xcosx'),gtext('y2=xtan(1/x)sin(x^3)'),gtext('y3=exp(1/x)sinx')79页的5题：t=1:pi/50:20*pi;x=t.*cos(pi/6.*t);y=t.*sin(pi/6.*t);z=2.*t;plot3(x,y,z);gridon;title('圆锥螺线')xlabel('x轴');ylabel('y轴');zlabel('z轴');axissquare79页的9题：t=-1:0.1:1;[x,y]=meshgrid(t);z1=5-x.^2-y.^2;subplot(1,3,1),surf(x,y,z1),title('三维曲面')a=input('a=(-50a50)'),z2=a*ones(size(x));subplot(1,3,2),surf(x,y,z2),title('平面')r0=abs(z1-z2)=1;zz=r0.*z2;yy=r0.*y;xx=r0.*x;subplot(1,3,3),plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'x'),title('交线')二、【实验结果】79页的3题：479页的5题：579页的9题：