MATLAB程序设计-关于常用绘图函数的运用方法及说明.docVIP

MATLAB程序设计作业 陈杰杰201309030207 2014-11-3 MATLAB具有强大的图形处理功能。下面给出了3个m脚本文件,请在MA TLAB环境下运行,观察其输出。 要求根据每个m文件输出的图形(共18个,用中文翻译并解释产生每个输出图形的函数具体是什么?其功能是什么? 文件1: clear all %清除工作区间所有的变量 clf %清除图形窗口的内容 mfilename(fullpath %返回当前正在运行的函数所在文件的文件名(全部路径 echo on %显示M文件执行的每一条命令 subplot(2,3,1 %使(2*3幅子图中第一个子图成为当前图 t = 0:0.1:10; %将以0为起点、以10为终点、以0.1为步长的一维矩阵赋值给t z = impulse(1, [1 1 1], t; %动力系统的脉冲响应数据,以时间t步长0.1为单位stairs(t(1:5:end,z(1:5:end %绘制阶梯状图,从第1行开始,间隔5行取1行,到最后1行为止 hold on %保持当前坐标轴和图形,并接受下一次绘制 plot(t,z,r %用红线绘制横轴为t、纵轴为z的二维函数图 plot([0 t(end], [0 0], k: %用黑色虚线绘制函数图像,要求经过原点平行于横轴、取值范围为0到t的最后一个值 title(Impulse Response - (STAIRS %将此图命名为Impulse Response - (STAIRS (脉冲响应-(阶梯图 subplot(2,3,2 %使(2*3幅子图中第二个子图成为当前图 theta = 2*pi*(0:74/75; %将以0为起点、以2*pi*74/75为终点、2*pi/75为步长的一维矩阵赋值给theta x = cos(theta; %计算cos(theta的值,并赋给x y = sin(theta; %计算sin(theta的值,并赋给y z = abs(fft(ones(10,1, 75; %ones(10,1生成十行一列的全一矩阵;fft(ones(10,1, 75进行快速傅里叶变换;z = abs(fft(ones(10,1, 75取幅值并转置 stem3(x, y, z %绘制三维杆状图 title(Polar FFT - (STEM3 %将此图命名为Polar FFT - (STEM3(极坐标下快速傅里叶变换-(三维针状图 subplot(2,3,3 %使(2*3幅子图中第二个子图成为当前图 [X,Y,Z] = peaks(-2:0.25:2; %产生-2为起点、2为终点、0.25为步长的guassian分布矩阵,返回峰函数的三个坐标轴空间上的数值,X表示在x轴,Y表示在y轴,Z表示在z轴,这样每个点就对应一个(X,Y, [U,V] = gradient(Z, 0.25; %返回二维数值梯度的U、V部分,这里的0.25指定了沿着梯度的方向取点的间隔为0.25 contour(X,Y,Z,10; %绘制矩阵Z的等高线,绘制的等高线被限定在由X、Y指定的区域内,等高线条数为10 hold on %保持当前坐标轴和图形,并接受下一次绘制 quiver(X,Y,U,V; %使用箭头来直观的显示矢量场,该调用格式表示通过在(X, Y 指定的位置绘制小箭头来表示以该点为起点的向量(U,V title(Surface Gradient - (CONTOUR QUIVER %将此图命名为Surface Gradient - (CONTOUR QUIVER(表面梯度-(等高线图和矢量场图 theta = 0:0.1:4*pi; %将以0为起点、以4*pi为终点、0.1为步长的一维矩阵赋值给theta [x,y] = pol2cart(theta(1:5:end, theta(1:5:end; %把极坐标(theta(1:5:end, theta(1:5:end转换为对应的二维笛卡尔坐标(x,y,theta(1:5:end表示从第一行到最后一行,每五行取一行 subplot(2,3,4 %使(2*3幅子图中第四个子图成为当前图 polar(theta,theta %绘制极坐标图像,第一个theta是用弧度制表示的角度,第二个theta是对应的半径 axis([-13 13 -12.5 14.5] %横坐标范围为-13到13,纵坐标范围为-12.5到14.5 title(Spiral Plot - (POLAR %将此图命名为Spiral Plot - (POLAR(螺旋图-极坐标图 subplot(2,3,5 %使(2*3幅子图中第五个子图成为当前图 compass(x,y %绘制罗盘图 axis([-13 13 -12.5 14.5]