大上数学实验报告.docVIP

西安交通大学实验报告 第 1 页 共 12 页 课程 实 验 日 期 2012 年 12 月 25 日 系 别 实 验 报 告 日 期 2012 年 12 月 28 日 专业班级_ 报 告 退 发 ( 订正 、 重做 ) 姓 名___ ___学号_ __ 教 师 审 批 签 字 实验I 一、实验问题 单位球面被柱面所截以及单位球面里挖掉柱面 二、问题分析 该问题应分解为单位球面被柱面所截和单位球面里挖掉柱面两部分考虑，两者的相同之处在于均需要绘制单位球面，区别在于前者需要将柱面显示出来，而后者需要将柱面隐去，且需要用到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.^2 1 ); 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) hold on mesh(X1,Y1,-Z1) mesh(X2,Y2,Z2) mesh(X2,-Y2,Z2) 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.^2 1 ); Z(i1)= nan; i2=find(X.^2-X+Y.^2=0); Z(i2)=nan; mesh(X,Y,Z) hold on mesh(X,Y,-Z) 四、问题求解结果与结论 1.单位球面被柱面所截 2.单位球面里挖掉柱面 五、问题的进一步拓展与实验 本实验中最显著的问题在于单位球面的绘制：步长大，运算快，但图像缺陷多；步长小，图像细腻，但运算太慢。于是需要寻找一种既运算快又图像细腻的单位球面程序设计方法。百度并作部分修改后得到利用参数方程的方法，程序和结果如下： 1.单位球面被柱面所截 r=1; theta=0:0.01:pi; phi=0:0.01:2*pi; [u,v]=meshgrid(theta,phi); x1=r.*sin(u).*cos(v); y1=r.*sin(u).*sin(v); z1=r.*cos(u); x2=0:0.01:1; z2=-1:0.01:1; [X2,Z2]=meshgrid(x2,z2); Y2=sqrt(-X2.^2+X2); mesh(x1,y1,z1) hold on mesh(X2,Y2,Z2) mesh(X2,-Y2,Z2) 2.单位球面里挖掉柱面 r=1; theta=0:0.01:pi; phi=0:0.01:2*pi; [u,v]=meshgrid(theta,phi); x=r.*sin(u).*cos(v); y=r.*sin(u).*sin(v); z=r.*cos(u); i=find(x.^2-y+y.^2=0); z(i)=nan; mesh(x,y,z) 六、实验总结与体会 第三、四部分中呈现的程序和结果较为容易实现，但在改进过程中，发现各种各样、各不相同的方法均能绘制出球体，其中使用的函数指令更是五花八门，十分令人费解。实验方法改进过程中最大的收获就是学会了利用参数方程绘制曲面图形。 并且，由于这是第一个实验，在此过程中对Matlab的操作方法、文件命名规则和存储方法的熟悉程度也在摸索中有了较大的提高，操作方面障碍的解除，为接下来的几个实验提供了极大的便利，也节省了许多时间。 实验II 一、实验问题 水仙花数 若一个三位自然数的各位数字的立方和等于该数本身，则称该自然数为水仙花数，例如153=13+53+33，所以153就是一个水仙花数，编程计算出所有的水仙花数 二、问题分析 首先需要利用三重循环语句表示出所有的三位数，再验证是否满足各位数字的立方和等于该数本身的条件，满足者输出，不满足者跳过。 三、程序设计 for a=1:1:9 for b=0:1:9 for c=0:1:9 D=a^3+b^3+c^3; E=100*a+10*b+c;