《王小才 绘图实验之马鞍面.doc》由会员分享,可在线阅读,更多相关《王小才 绘图实验之马鞍面.doc(4页珍藏版)》请在金锄头文库上搜索。
1、淮阴工学学院数理学院 数学实验课程 实验报告实验名称 Matlab绘图 姓名 王小才 日期 2012-08-23 【实验目的及意义】1熟悉MATLAB曲面绘图命令;2 掌握MATLAB绘制平面图3 掌握MATLAB绘制空间曲面图;4 掌握MATLAB图形处理命令。 通过该实验的学习,使学生能灵活应用MATLAB软件绘制一些常见的函数图像。【实验内容】一、 画出马鞍面在xy平面投影是单位圆,并且去掉马鞍面的1/4部分。解答: 方法一:用for循环把去掉部分z的值设置为NaN。 方法二:采用极坐标方法直接画出需要画的区域,即把要画的区域用极坐标生成网格,然后极坐标转化为平面坐标。 下面程序采用的是
2、方法二, 方法二程序耗时少。程序命令% 画马鞍面 z=ax2-by2a=3;b=4;rho=0:0.01:1;theta=0:0.01:3/2*pi;Rho,Theta=meshgrid(rho,theta);%用极坐标生成网格x=Rho.*cos(Theta);y=Rho.*sin(Theta);%极坐标转化为平面坐标z=a*x.2-b*y.2;mesh(x,y,z)%马鞍面去掉1/4的部分 hold ont = 0:0.02: 2*pi;u = sin(t);v = cos(t);plot(u, v, m, linewidth, 2.5)s = linspace(0, 1, 50);plo
3、t(s, 0*s, 0*s, -s, b, linewidth, 2)sx = 0:0.05:1;sy = 0*sx;sz=a*sx.2-b*sy.2stem3(sx,sy,sz)sy = -1:0.05:0;sx = 0*sy;sz=a*sx.2-b*sy.2stem3(sx,sy,sz)t = 0:0.02: 2*pi;u1 = sin(t);v1 = cos(t);rz=-4*ones(size(u1);plot3(u1, v1, rz,m, linewidth, 2.5)s = linspace(0, 1, 50);rho=1;theta=0:0.01:3/2*pi;Rho,Theta
4、=meshgrid(rho,theta);x=Rho.*cos(Theta);y=Rho.*sin(Theta);z=a*x.2-b*y.2;n=length(x);for i = 1: n z1=-4:0.01:z(i); x1=x(i)*ones(size(z1); y1=y(i)*ones(size(z1); plot3(x1,y1,z1,y)end view(-30, 30)title(A Part of surface z=ax2-by2, fontsize, 16)xlabel(x,FontName,Times New Roman,FontSize,14)ylabel(y,FontName,Times New Roman,FontSize,14,Rotation,0) zlabel(z,FontName,Times New Roman,FontSize,14)set(gca,FontName,Times New Roman,FontSize,14)运行结果:去掉边界的效果:第- 1 -页 共 4页