实验五MATLAB的科学计算..docVIP

实验五 MATLAB在高等数学中的科学计算 实验目的 了解MATLAB在高等数学相关领域中的应用 实验地点：A404 实验日期： 实验内容 （一）matlab在解析几何中的应用 1、利用Matlab描绘三维空间图形 解析几何教学中，图形是很重要的，空间概念的建立，各种曲面的研究，各种方程的建立都离不开图形，而空间图形既难画又费时。借助Matlab的绘图功能，可以快捷、准确地绘出图形，使教学变得形象、直观、生动，有利于学生观察三维空间图形的形状，掌握图形的性质。 例1：作出球面 x2+y2+z2=a2 和圆柱面x2+y2-ax=0 的交线—— 维维安尼( Viviani )曲线 此题通过联立球面方程x2+y2+z2=a2和圆柱面方程x2+y2-ax=0而得出维维安尼曲线的参数方程： ??? ?x=acos2t ??? ?y=acost * sint? (0≤t2 ∏) ?? z=asint 利用Matlab编制以下程序，可绘出维维安尼曲线的形状，如图1所示： clear; t=0:pi/30:2*pi; a=4; x=a*(cos(t)).^2; y=a*cos(t).*sin(t); z=a*sin(t); plot3(x,y,z) pause(1) grid on 图1 从不同角度观察维维安尼曲线 也可在同一直角坐标系中绘出球面与柱面相交的情形，这更利于我们观察维维安尼曲线的形状。编制的Matlab程序如下： %柱面的Matlab程序 t=0:.1:pi; x=4*(cos(t).^2); y=4*cos(t).*sin(t); n=size(x,2); for i=-5:.01:5 ??? z=i*ones(1,n); plot3(x,y,z)??? hold on end %球面的Matlab程序 u=0:.1:pi; v=-pi:.2:pi/2; [U,V]=meshgrid(u,v); x=4*sin(U).*cos(V); y=4*sin(U).*sin(V); z=4*cos(U); mesh(x,y,z) grid on hold off 运行结果如图2所??????????????????????????????? 图2 从不同角度观察球面与柱面相 交而得的维维安尼曲线 利用Matlab判定空间图形间的位置关系 ??? 掌握空间图形的位置关系是解析几何学习中的难点之一，学生的空间想象能力不强，准确地把握空间图形之间的位置关系更加困难，在同一直角坐标系下形象地显示出各个空间图形是解决这一问题的关键。借助Matlab可准确地显示三维空间图形的位置关系。 例2： ?判断平面 2x-y-3z+11, 球 ??(x-3)2+(y+5)2+(z+2)2=602, 椭球 (x+30)2/402+(y-20)2/1002+(z-100)2/102=1 之间的位置关系。 此题通过计算来判断平面、球、椭球三者之间的位置关系较为繁琐，利用Matlab在三维直角坐标系中同时绘出球、椭球与平面，既可直观地显示它们之间的位置关系，也可帮助我们验证计算结果。 球、椭球与平面的Matlab程序为： ?? u=-100:100; [x,y]=meshgrid(u); z=(2*x-y+11)/3; u=0:pi/20:pi; v=0:pi/20:2*pi; [U,V]=meshgrid(u,v); x1=60*sin(U).*cos(V)+3; y1=60*sin(U).*sin(V)-5; z1=60*cos(U)-2; x2=40*sin(U).*cos(V)-30; y2=100*sin(U).*sin(V)+20; z2=10*cos(U)+100; mesh(x,y,z); hold on surf(x1,y1,z1); surf(x2,y2,z2); view(-40,16) ? ?运行结果如图3 所示。由图H 可知，平面与球相交、与椭球相离，球与椭球相离。 图3 平面、球、椭球的位置关 利用Matlab研究旋转曲面的性质 解析几何学习中，旋转曲面的性质仅仅通过 教师的讲解和板书的静态演示很难让学生掌握， 利用Matlab演示旋转曲面的形成过程则可突破 这一教学难点。 ?? 例3 演示将椭圆 y2/a2+z2/b2=1(ab)x=0 绕短轴（即z 轴)旋转所得的扁形旋转椭球面的形成过程。 因a, b为常数，不妨取a=20,b=2，利用Matlab编制以下程序可实现旋转椭球面形成过程的动态演示： a=20; b=2; t=0:pi/30:2*pi; n=length(t); y=20*cos(t); z=2*sin(t); x=0*ones(1,n); plot3(x,y,z) grid on pause(1