绘图命令
（一）二维绘图命令
二维绘图的基本命令有plot，loglog，semilogx，semilogy和polar。它们的使用方法基本相同，其不同点是在不同的坐标中绘制图形。
plot命令使用线性坐标空间绘制图形；
loglog命令在两个对数坐标空间中绘制图形；
而semilogx（或semilogy）命令使用x轴（或y轴）为对数刻度，另外一个轴为线性刻度的坐标空间绘制图形；
polar使用极坐标空间绘制图形。
二维绘图命令plot为了适应各种绘图需要，提供了用于控制线色、数据点和线型的3个基本参数。它的使用格式如下：
plot(x, y, 'color_point_linestyle')该命令是绘制y对应x的轨迹的命令。y与x均为向量，且具有相同的元素个数。用字符串'color_point_linestyle'完成对上面3个参数的设置。
线色（r-red，g-green, b-blue，w-white，k-balck，i-invisible，y-yellow）
数据点（.，o，x，+，*，S，H，D，V，^，>，<，p）
线型（-，-.，--，:）
都可以根据需要适当选择。
当plot(x, y)中的x和y均为mxn矩阵时，plot命令将绘制n条曲线。
如果多重曲线对应不同的向量绘制，可使用命令
plot(t1, x1, t2, x2, t3, x3)其中：x1对应t1；x2对应t2；等等。在这种情况下，t1，t2和t3可以具有不同的元素个数，但要求x1、x2和x3必须分别与t1、t2和t3具有相同的元素个数。
也可以使用另外一种形式，将多张图在一张图上显示。
x = [0:0.5:7]
y1 = sin(x)
y2 = cos(x)
y3 = (1/2)*x-0.5
hold on
plot(x, y1)
plot(x, y2)
plot(x, y3)
hold off
plot(x, y1, 'r', x, y2, 'g', x, y3, 'b')
subplot命令使得在一个屏幕上可以分开显示n个不同坐标系，且可分别在每一个坐标系中绘制曲线。其命令格式如下：
subplot(r,c,p)该命令将屏幕分成r x c个子窗口，而p表示激活第p个子窗口。窗口的排号是从左到右，自上而下。
subplot(2, 3, 2)
subplot(2, 3, 3)
subplot(2, 3, 4)
在图形绘制完毕后，执行如下命令可以再在图中加入标题、标号、说明和分格线等。
这些命令有title，xlabel，ylabel，text，gtext等。它们的命令格式如下：
title('My Title')
xlabel('My X-axis Label')
ylabel('My Y-axis Label')
text(x, y,'Text for annotation')
gtext('Text for annotation')
gridgtext命令是使用鼠标定位的文字注释命令。当输入命令后，可以在屏幕上得到一个光标，然后使用鼠标控制 它的位置。单击鼠标左键，即可确定文字设定的位置。
hold on是图形保持命令，可以把当前图形保持在屏幕上不变，同时在这个坐标系内绘制另一个图形。hold命令是一个交替转换命令，即执行一次，转变一个状态（相当于hold on，hold off）
y = x
subplot(2, 3, 2)
title('my data1')
xlabel('my data xlabel')
ylabel('my data ylabel')
text(x, y,'data')
gtext('data symbol')
grid
（二）显函数、符号函数或隐函数的绘图
fplot（fun，lims）绘制由字符串fun指定函数名的函数在x轴区间为lims=[xmin, xmax]的函数图。若lims=[xmin, xmax, ymin, ymax]，则y轴也被限制。
例1：
function y = Afun1(x);
if x<1y = x + 1;
elsey = 1+ 1/x;
end然后在命令窗口输入
>> fplot('Afun1', [-3, 3])就可以画出函数f(x)的图形。
也可以使用匿名函数：
>> Afun1 = @(x) (x + 1)*(x < 1) + (1 + 1/x)*(x >= 1);
>> fplot(Afun1, [-3, 3])>> ezplot(Afun1)%函数Afun1的图形，横轴的近似范围为[-2π，2π]，也可以使用参数来代替默认的参数例2
>> ezplot('x^2 + (y^2)/4 = 1')
（三）三维图形
在实际工程计算中，最常用的三维绘图是三维曲线图、三维网格图和三维曲面图3中基本类型。与此对应，Matlab也提供了一些三维基本绘图命令，如三维曲线命令plot3，三维网格图命令mesh和三维表面图命令surf。
1.三维曲线
plot3（x, y, z）通过描点连线画出曲线，这里x，y，z都是n维向量，分别表示该曲线上点集的横坐标、纵坐标、竖坐标。
例3
在区间[0, 10π]画出参数曲线x=sint，y=cost，z=t，并分别标注。
t = 0:pi/50:10*pi;
plot3(sin(t), cos(t), t)
xlabel('sin(t)'), ylabel('cos(t)'), zlabel('t')
2.网格图
命令mesh（x, y, z）画网格曲面，这里x，y，z是三个同维数的数据矩阵，分别表示数据点的横坐标、纵坐标、竖坐标，命令mesh(x, y, z)将该数据点在空间中描出，并连成网格。
例4
x = -3:0.1:3; y = -5:0.1:5;
[x, y] = meshgrid(x, y); %生成网格数据
z = (sin(x.*y) + eps)./(x.*y + eps);
mesh(x, y, z)
3.表面图
命令surf(x, y, z)画三维表面图，这里x，y，z是三个同维数的数据矩阵，分别表示数据点的横坐标、纵坐标、竖坐标。
例5
[x, y] = meshgrid([-3:0.2:3]); %生成网格数据
z = (sin(x.*y) + eps)./(x.*y + eps);
surf(x, y, z)
4.旋转曲面
例6
x = 0:10:600;
[X, Y, Z] = cylinder(30*exp(-x/400).*sin((x + 25*pi)/100) + 130);
surf(X, Y, Z)
例7
画出x^2 + (y-5)^2= 16绕x轴旋转一周所形成的旋转曲面。
解：因为这里的函数是隐函数，化成显函数后有两支，必须使用参数方程，旋转面的参数方程为
x = 4cosα，
y = (5 + 4sinα)cosβ，
z = (5 + 4sinα)sinβ. 式中：α，β∈[0, 2π]。
alpha = [0:0.1:2*pi]';
beta = 0:0.1:2*pi;
x = 4*cos(alpha)*ones(size(beta));
y = (5 + 4*sin(alpha))*cos(beta);
z = (5 + 4*sin(alpha))*sin(beta);
surf(x, y, z)也可以使用匿名函数
x = @(alpha, beta) 4*cos(alpha);
y = @(alpha, beta) (5 + 4*sin(alpha)).*cos(beta);
z = @(alpha, beta) (5 + 4 *sin(alpha)).*sin(beta);
ezsurf(x, y, z)
5.其他二次曲面
Matlab中使用绘图命令ezmesh或ezsurf也很方便，只需要把曲面方程写成两个变量的显函数方程或参数方程即可。
例8
（1）旋转单叶双曲面 (x^2 + y^2)/9 - (z^2)/4 = 1;
x = @(s, t) 3 *sec(s).*cos(t);
y = @(s, t) 3*sec(s).*sin(t);
z = @( s, t) 2*tan(s);
ezmesh(x, y, z)
（2）旋转双叶双曲面 (x^2)/9 - (y^2 + z^2)/4 = 1;
x = @(s, t) 3 *sec(s);
y = @(s, t) 2*tan(s).*cos(t);
z = @( s, t) 2*tan(s).*sin(t);
ezmesh(x, y, z)
（3）抛物柱面 y^2 = x;
ezsurf(@(y, z) y.^2, 50)
（4）椭圆锥面 (x^2)/9 + (y^2)/4 = z^2;
z = @(s, t) tan(s);
x = @(s, t) 3*tan(s).*cos(t);
y = @( s, t) 2*tan(s).*sin(t);
ezsurf(x, y, z)
（5）椭球面 (x^2)/9 + (y^2)/4 + (z^2)/6 = 1;
ellipsoid(0, 0, 0, 3, 2, sqrt( 6))
（6）马鞍面 z = xy;
ezsurf(@(x, y) x.*y)
（7）椭圆柱面 (x^2)/9 + (y^2)/4 = 1;
x = @(s, t) 3*cos(s);
y = @( s, t) 2*sin(s);
z = @(s, t) t;
ezmesh(x, y, z)
（四）3 - D可视化图形
例9
[x, y, z, v] = flow;
isosurface(x, y, z, v);
例10
画出v = (x^2)*y*(z + 1)的示意图。
x = 1:20;
y = 1:10;
z = -10:10;
[x, y, z] = meshgrid(x, y, z);
v = x.^2.*y.*(z + 1);
isosurface(x, y, z, v)
