大一上数学实验报告.pdfVIP

以家为家，以乡为乡，以国为国，以天下为天下。——《管子》

西安交通大学实验报告成绩

第1页共12页

课程实验日期2012年12月25日

系别实验报告日期2012年12月28日

专业班级_报告退发(订正、重做)

姓名______学号___教师审批签字

实验I

一、实验问题

单位球面被柱面所截以及单位球面里挖掉柱面

11

x2y2z2x2y22

1()()

22

二、问题分析

该问题应分解为单位球面被柱面所截和单位球面里挖掉柱面两部分考虑，两

者的相同之处在于均需要绘制单位球面，区别在于前者需要将柱面显示出来，而

后者需要将柱面隐去，且需要用到find和nan指令对图像实现精确绘制。

三、程序设计

1.单位球面被柱面所截

x1=-1:0.01:1;

y1=-1:0.01:1;

[X1,Y1]=meshgrid(x1,y1);

Z1=sqrt(1-(X1.^2+Y1.^2));

i1=find(X1.^2+Y1.^21);

Z1(i1)=nan;

x2=0:0.01:1;

z2=-1:0.01:1;

[X2,Z2]=meshgrid(x2,z2);

Y2=sqrt(-X2.^2+X2);

mesh(X1,Y1,Z1)

holdon

mesh(X1,Y1,-Z1)

mesh(X2,Y2,Z2)

mesh(X2,-Y2,Z2)

1

以家为家，以乡为乡，以国为国，以天下为天下。——《管子》

2.单位球面里挖掉柱面

x=-1:0.01:1;

y=-1:0.01:1;

[X,Y]=meshgrid(x,y);

Z=sqrt(1-(X.^2+Y.^2));

i1=find(X.^2+Y.^21);

Z(i1)=nan;

i2=find(X.^2-X+Y.^2=0);

Z(i2)=nan;

mesh(X,Y,Z)