王小才 绘图实验之圆柱相交.doc

淮阴工学学院 数理学院 数学实验课程 实验报告 实验名称 Matlab绘图 姓名 王小才 日期 2012-08-23 【实验目的及意义】 [1]熟悉MATLAB；] 掌握MATLAB [3] 掌握MATLAB绘制空间曲面图； [4] 掌握MATLAB。 通过该实验的学习，使学生能灵活应用MATLAB软件，相交部分的图形。 解答：用for循环把去掉部分z的值设置为不定值（NaN），耗时较大，建议用find语句实现。 程序命令 clc; clear; t=[0:0.01:pi+0.01]; s=[0:0.01:2]; x=cos(t)*(0*s+1); y=sin(t)*(0*s+1); z=(0*t+1)*s; xx=cos(t)*(0*s+1); zz=sin(t)*(0*s+1); yy=(0*t+1)*s; z1=x.^2+z.^2; [k1,k2]=find(z11); for i=1:length(k1) z(k1(i),k2(i))=NaN; end z2=xx.^2+yy.^2; [k1,k2]=find(z21); for i=1:length(k1) zz(k1(i),k2(i))=NaN; end figure(color,[1 1 1]) h=surf(x,y,z) %利用对称性作图 hold on h1=surf(x,-y,-z); h2=surf(x,-y,z); h3=surf(x,y,-z); h4=surf(xx,yy,zz) %利用对称性作图 h5=surf(xx,-yy,-zz); h6=surf(xx,-yy,zz); h7=surf(xx,yy,-zz); axis equal light(position,[ 1 1 1]); lighting phong shading interp; title(两个圆柱x^2+y^2=1,x^2+z^2=1 相交部分的图形, fontsize, 16) axis off box on camlight(headlight) 运行结果： 第- 3 -页 共 3页