西安邮电大学matlab考试,MATLAB实验报告西安邮电大学.doc

MATLAB实验报告西安邮电大学

西安邮电学院

《Matlab》

实验报告

(三)

2011- 2012 学年第 1 学期

专业：自动化班级：自动0903学号：姓名：

2011 年 10 月 20 日

实验三 MATLAB图形系统

一、实验目的

1．．．．、实验仪器和设备实验原理

1．二维数据曲线图

绘制单根二维曲线 plot(x,y);

绘制多根二维曲线 plot(x,y) 当x是向量，y是有一维与x同维的矩阵时，则绘制多根不同颜色的曲线。当x，y是同维矩阵时，则以x，y对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。

含有多个输入参数的plot函数 plot(x1,y1,x2,y2,…,xn,yn)

具有两个纵坐标标度的图形 plotyy(x1,y1,x2,y2)

2．图形标注与坐标控制

1)title (图形名称)；

2)xlabel(x轴说明)

3)ylabel(y轴说明)

4)text(x，y图形说明)

5)legend(图例1，图例2，…)

6)axis ([xmin xmax ymin ymax zmin zmax])

3．图形窗口的分割

subplot(m,n,p)

4．三维曲线

plot3(x1,y1,z1,选项1，x2,y2，选项2，…,xn,yn,zn,选项n)

5．三维曲面

mesh(x,y,z,c) 与surf(x,y,z,c)。一般情况下，x，y，z是维数相同的矩阵。X，y是网格坐标矩阵，z是网格点上的高度矩阵，c用于指定在不同高度下的颜色范围。

6．图处理

imread和imwrite函数 这两个函数分别用于将图象文件读入工作空间，以及将图象数据和色图数据一起写入一定格式的图象文件。

image和imagesc函数 这两个函数用于图象显示。为了保证图象的显示效果，一般还应使用colormap函数设置图象色图。四、预习要求

1．复习二维与三维图形的绘图函数。

2．复习图形辅助操作。五、实验内容及步骤

1．设，在x＝0～2π区间取101点，绘制函数曲线。x=0:2*pi/100:2*pi;

y=(0.5+3.*sin(x)./(1+x.*x))./cos(x);

plot(x,y);

grid

title('x-y')

xlabel('x')

ylabel('y')

2．已知y1=x2,y2=cos(2x),y3=y1*y2,完成下列操作：

(1)在同一坐标系下用不同的颜色和线型绘制三条曲线；x=0:0.01:2*pi;

y1=x.*x;

y2=cos(2.*x);

y3=y1.*y2;

plot(x,y1,x,y2,x,y3)

grid

(2)分别用条形图、阶梯图、杆图和填充图绘制三条曲线。y3=y1.*y2;

subplot(2,2,1);

bar(x,y1,'r');

title('bar(x,y1,"r")');

grid;

subplot(2,2,2);

fill(x,y1,'r');

title('fill(x,y1,"r")');

grid;

subplot(2,2,3);

stairs(x,y1,'r');

title('stairs(x,y1,"r")');

grid;

subplot(2,2,4);

stem(x,y1,'r');

title('stem(x,y1,"r")');

grid

同理 y2，y3如下图 所示

3．已知

在－5<=x<=5区间绘制函数曲线。x=-5:0.1:5;

y=((x+sqrt(pi))./exp(2)).*(x<=0)+1/2.*log(x+sqrt(1+x.^2)).*(x>0);

plot(x,y);

grid;

4．绘制函数的曲面图和等高线

其中x的21个值均匀分布在[-5,5]范围，y的31个值均匀分布在[0,10]，要求使用subplot(2，1，1)和subplot(2，1，2)将产生的曲面图和登高图画在同一个窗口上。x=-5:(10/20):5;

y=0:(10/30):10;

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

z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)./4)

subplot(2,1,1);

mesh(x,y,z);

subplot(2,1,2);

meshc(x,y,z);

5．，即一个周期的曲线。

2)正割曲线为红色点划线输出，余割曲线为蓝色实线输出。

3)图形上面表明正割和余割公式，横轴标x，纵轴标y。

4)将图形窗口分成两格，正割在上，余割在下。

代码 ：

x=0:0.01:2*pi;

y1=sec(x);

y2=csc(x);

subplot(2,1,1);

plot(x,y1,'r');

title('y=sec(x