数学应用软件作业2 Matlab作图

1、注意：上机作业文件夹以自己的姓名学号命名，文件夹中包括如下上机报告和Matlab程序。佛山科学技术学院上 机 报 告课程名称 数学应用软件 上机项目 Matlab作图 专业班级 姓 名 学 号 一. 上机目的本章学习了有关Matlab作图的有关知识，主要有：1. 二维绘图函数：（1）曲线图：plot(X,Y,S) （2）ezplot，fplot等函数作图2. 三维图形：（1）三维曲线：plot3(x,y,z,S) （2）空间曲面surf(X,Y,Z) ,mesh(X,Y,Z);3. 处理图形：（1）subplot,grid on,grid off ,xlabel,ylabel,tiltle,g

2、text ,hold on, figure等4. 特殊二、三维图形polar(theta,rho,s) ,contour3(x,y,z,n), scatter3等二. 上机内容1、在同一坐标系中，取自变量区间15，200，用不同的颜色或线型分别画出下面四个函数的图形：观测当增大时图形的走向。注：可选用r,g,b颜色区分各曲线。2、用ezplot绘制函数在-3 ,3上的图形。3、用ezplot绘制摆线的图形。4、将一幅图分成四块，在各子图中分别绘制函数 和的图像，区间取0,5。注：可用上subplot,grid on,xlabel,ylabel,tiltle,gtext等命令。5、用surf，m

3、esh绘制曲面，。6、用polar绘制阿基米德螺线和三叶玫瑰线。（a=100）7、在同一平面的两个窗口中分别画出心形线和马鞍面。三. 上机方法与步骤给出相应的问题分析及求解方法，并写出Matlab程序。并有上机程序显示。第1题：要在同一坐标系中分别画出四个函数的图形，可以利用plot来将这四条曲线画在一起。Matlab程序：x=linspace(15,200,50);y1=(1+1./x).x;y2=(1+1./x).(x+1);y3=1+1./(x+1).x;y4=2.7183;plot(x,y1,'r',x,y2,'g',x,y3,'b',x

4、,y4,'co')第2题：直接用ezplot命令绘制函数的图形。Matlab程序：ezplot('（exp(x*y)-sin(x+y)）',-3,3)第3题：取a=1，直接用ezplot命令绘出这两条曲线。Matlab程序：ezplot('t-sin(t)','1-cos(t)',0,2*pi)第4题：利用subplot命令将这四个函数的图像画在两行两列的同一个图形上。Matlab程序：x=0:0.1:5;subplot(2,2,1);plot(x,sin(5*x);grid on;title('plot-y=sin(5*

5、x)'); subplot(2,2,2);plot(x,cos(3*x);grid on;title('plot-y=cos(3*x)'); subplot(2,2,3);plot(x,sin(5*x).*cos(3*x);grid on;title('plot-y=sin(5*x).*cos(3*x)'); subplot(2,2,4);plot(x,sin(5*x)+cos(3*x);grid on;title('plot-y=sin(5*x)+cos(3*x)');第5题：分别用命令surf（X,Y,Z）、命令mesh（X,Y,Z）

6、将该数据在空间上描出，并连成网络。Matlab程序：（1）x=-3:0.1:3;y=-3:0.1:3;X,Y=meshgrid(x,y);Z=2.*(X.2)+Y.2;surf(X,Y,Z)（2）x=-3:0.1:3;y=-3:0.1:3;X,Y=meshgrid(x,y);Z=2.*(X.2)+Y.2;mesh(X,Y,Z)第6题：先考虑用polar绘制阿基米德螺和三叶玫瑰线，再利用命令subplot将这两个图像画在同一个图形上。Matlab程序：theta=linspace(0,2*pi);rho1=100*theta;rho2=100*cos(theta);subplot(1,2,1);

7、polar(theta,rho1,'g');title('polar plot of 100*theta');subplot(1,2,2);polar(theta,rho2,'g');title('polar plot of 100*cos(theta)');第7题：要在同一平面的两个窗口中分别画出心形线和马鞍面，可以用命令subplot绘出一行两列的图形。Matlab程序：theta=linspace(0,2*pi);rho=100*1-cos(theta);subplot(1,2,1);polar(theta,rho,'g');title('心形线');x=-3:0.1:3;y=-5:0.1:5;X,Y=meshgrid(x,y);Z=X.2-Y.2;subplot(1,2,2);mesh(X