matlab实验报告一.doc

20112012学年第一学期合肥学院数理系实验报告课程名称： MATLAB实验 实验项目： 特殊函数与图形 实验类别：综合性 设计性 验证性 专业班级： 09信息与计算科学 姓 名： 韦志强 学 号： 0907011039 实验地点： 7#604 实验时间： 2011-10-19 指导教师： 李旭 成 绩： 一.实验目的：1学会绘制matlab中的一些符号函数和各种简单图形；2初步了解关于简单图形的简单程序的编写；3懂得脚本文件和运用函数文件的使用方法。二.实验内容：1.ezmesh(6+2*cos(u)*cos(v),(6+2*cos(u)*sin(v),2*sin(u),0,2*pi,0,2*pi) axis equal hold on ezmesh(2*cos(u)*cos(v),(2*cos(u)*sin(v),2*sin(u),0,2*pi,0,2*pi) axis equal title(圆环面和球面)2.ezmesh(2*cos(u)*cos(v),(2*cos(u)*sin(v),2*sin(u),0,2*pi,0,2*pi) axis equal title(球面) grid off3.ezmesh(2*cos(u)*cos(v),(2*cos(u)*sin(v),6*sin(u),0,2*pi,0,2*pi) axis equal title(椭球面) grid off4.ezsurf(3*tan(u)*cos(v),3*tan(u)*sin(v),5*sec(u),-pi/2,3*pi/2,0,2*pi); axis auto title(双叶双曲面)5.ezsurf(3*sec(u)*cos(v),3*sec(u)*sin(v),5*tan(u),-pi/2,pi/2,0, 2*pi); axis auto title(单叶双曲面)6.x,y=meshgrid(-25:1:25,-25:1:25); z=x.2/9-y.2/4; surf(x,y,z) title(马鞍面) grid on meshz(x,y,z)3. 实验要求1.学会并了解Matlab中一些符号函数的绘图功能和简单图形的绘制，学会熟练运用各种绘图语句；2.在实验中要求领会各条语句的意义，要求将每个所需要画出的图形的代码写成函数文件或者脚本文件。 四 实验步骤及结果1.环面和球面的组合图： u,v=meshgrid(0:1/100:2*pi); x=(6+2.*cos(u).*cos(v); y=(6+2.*cos(u).*sin(v); z=2.*sin(u); mesh(x,y,z); axis equal; hold on a,b,c=sphere(20); surf(2*a,2*b,2*c) shading flat2.1，球面： x,y,z=sphere(50); surf(4.*x+1,4.*y+1,4.*z+1)2.2,椭球面： ezsurf(2*cos(u)*cos(v),(2*cos(u)*sin(v),6*sin(u), 0,2*pi,0,2*pi) axis equal title(椭球面)2.3,双叶双曲面： ezsurf(3*tan(u)*cos(v),3*tan(u)*sin(v),5*sec(u),-pi/2, 3*pi/ 2, 0,2*pi); axis auto title(双叶双曲面)2.4,单叶双曲面：ezsurf(3*sec(u)*cos(v),3*sec(u)*sin(v),5*tan(u),-pi/2,pi/2, 0, 2*pi); axis auto title(单叶双曲面)5,马鞍面： x,y=meshgrid (-20:1:20,-20:1:20); z=x.2/9-y.2/4; meshz(x,y,z) colormap(gray) title(马鞍面) grid off6,黎曼函数： function liman(n) % n=input (Enter n:); x=; y=; k=1; for i=2:n for j=1:i-1 if gcd(i,j)=1 x(k)=j/i; y(k)=1/i; k=k+1; end end end plot(x,y,.b) axis(0,1,0,1)五 实验总结：通过这次实验我学会了如何熟练运用matlab中各种简