第5章 MATLAB绘图.doc

1、第 5 章 MATLAB 绘图文库帮手网 免费帮下载 百度文库积分 资料本文由 wzedehao 贡献ppt 文档可能在 WAP 端浏览体验不佳。建议您优先选择 TXT，或下载源文件到本机查看。第 5 章 MATLAB 绘图 章 绘图 5.1 5.2 5.3 5.4 5.5 二维数据曲线图 其他二维图形 隐函数绘图 三维图形 图形修饰处理5.1 二维数据曲线图 5.1.1 绘制单根二维曲线 plot 函数的基本调用格式为： 函数的基本调用格式为： 函数的基本调用格式为 plot(x,y) 其中 x 和 为长度相同的向量 分别用于存储 x 坐标和 为长度相同的向量， 其中 和 y 为长度相同的

2、向量，分别用于存储 坐标和 y 坐标数据。 坐标数据。 坐标数据例 5-1 在 0x2 区间内，绘制曲线 区间内， y=2e-0.5xcos(4x) 程序如下： 程序如下： x=0:pi/100:2*pi; y=2*exp(-0.5*x).*cos(4*pi*x); plot(x,y)绘制曲线。 例 5-2 绘制曲线。 程序如下： 程序如下： t=0:0.1:2*pi; x=t.*sin(3*t); y=t.*sin(t).*sin(t); plot(x,y);plot 函数最简单的调用格式是只包含一个输入参数： 函数最简单的调用格式是只包含一个输入参数： 函数最简单的调用格式是只包含一个输入

3、参数 plot(x) 在这种情况下， 是实向量时， 在这种情况下，当 x 是实向量时，以该向量元素的下 是实向量时 标为横坐标，元素值为纵坐标画出一条连续曲线， 标为横坐标，元素值为纵坐标画出一条连续曲线， 这实际上是绘制折线图。 这实际上是绘制折线图。5.1.2 绘制多根二维曲线 1plot 函数的输入参数是矩阵形式 函数的输入参数是矩阵形式 (1) 当 x 是向量，y 是有一维与 同维的矩阵时，则绘 是向量， 是有一维与 同维的矩阵时， 是有一维与 x 同维的矩阵时 是向量 制出多根不同颜色的曲线。曲线条数等于 y 矩阵的 制出多根不同颜色的曲线。曲线条数等于 矩阵的 另一维数， 被作为这

4、些曲线共同的横坐标 被作为这些曲线共同的横坐标。 另一维数，x 被作为这些曲线共同的横坐标。 (2) 当 x,y是同维矩阵时，则以 对应列元素为横、 是同维矩阵时， 对应列元素为横、 是同维矩阵时 则以 x,y 对应列元素为横 纵坐标分别绘制曲线，曲线条数等于矩阵的列数。 纵坐标分别绘制曲线，曲线条数等于矩阵的列数。(3) 对只包含一个输入参数的 对只包含一个输入参数的 plot函数，当输入参数 函数， 函数 是实矩阵时， 是实矩阵时，则按列绘制每列元素值相对其下标 的曲线，曲线条数等于输入参数矩阵的列数。 的曲线，曲线条数等于输入参数矩阵的列数。 当输入参数是复数矩阵时， 当输入参数是复数矩

5、阵时，则按列分别以元素实部 和虚部为横、纵坐标绘制多条曲线。 和虚部为横、纵坐标绘制多条曲线。2含多个输入参数的 plot 函数 含多个输入参数的 函数 调用格式为： 调用格式为： plot(x1,y1,x2,y2,xn,yn) (1) 当输入参数都为向量时，x1 和 y1，x2 和 y2， 当输入参数都为向量时， 和 ， 和 ， ， xn 和 yn 分别组成一组向量对，每一组向量对的 分别组成一组向量对， 和 分别组成一组向量对 长度可以不同。每一向量对可以绘制出一条曲 长度可以不同。 这样可以在同一坐标内绘制出多条曲线。 线，这样可以在同一坐标内绘制出多条曲线。 (2) 当输入参数有矩阵形

6、式时，配对的 按对应列 当输入参数有矩阵形式时，配对的 x,y 按对应列 元素为横、纵坐标分别绘制曲线， 元素为横、纵坐标分别绘制曲线，曲线条数等 于矩阵的列数。 于矩阵的列数。分析下列程序绘制的曲线。 例 5-3 分析下列程序绘制的曲线。x1=linspace(0,2*pi,100); x2=linspace(0,3*pi,100); x3=linspace(0,4*pi,100); y1=sin(x1); y2=1+sin(x2); y3=2+sin(x3); x=x1;x2;x3; y=y1;y2;y3; plot(x,y,x1,y1-1)3具有两个纵坐标标度的图形 在 MATLAB 中

7、，如果需要绘制出具有不同纵坐标标 中 度的两个图形，可以使用 plotyy 绘图函数。调用 绘图函数。 度的两个图形，可以使用 绘图函数 格式为：格式为： plotyy(x1,y1,x2,y2) 其中 x1,y1 对应一条曲线，x2,y2对应另一条曲线。横 对应一条曲线， 对应另一条曲线。 其中 对应一条曲线 对应另一条曲线 坐标的标度相同，纵坐标有两个， 坐标的标度相同，纵坐标有两个，左纵坐标用于 x1,y1 数据对，右纵坐标用于 x2,y2 数据对。 数据对，右纵坐标用于 数据对。 数据对 数据对用不同标度在同一坐标内绘制曲线 y1=0.2e 例 5-4 用不同标度在同一坐标内绘制曲线 0

8、.5xcos(4x) 和 y2=2e-0.5xcos(x)。 。 程序如下： 程序如下： x=0:pi/100:2*pi; y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); plotyy(x,y1,x,y2);4图形保持 hold on/off 命令控制是保持原有图形还是刷新原有 命令控制是保持原有图形还是刷新原有 图形，不带参数的hold 命令在两种状态之间进行 图形，不带参数的 命令在两种状态之间进行 切换。 切换。采用图形保持， 例 5-5 采用图形保持，在同一坐标内绘制曲线 y1=0.2e-0.5xcos(4x)

9、 和 y2=2e-0.5xcos(x)。 。 程序如下： 程序如下： x=0:pi/100:2*pi; y1=0.2*exp(-0.5*x).*cos(4*pi*x); plot(x,y1) hold on y2=2*exp(-0.5*x).*cos(pi*x); plot(x,y2); hold off5.1.3 设置曲线样式 MATLAB 提供了一些绘图选项，用于确定所绘曲线 提供了一些绘图选项， 提供了一些绘图选项 的线型、颜色和数据点标记符号， 的线型、颜色和数据点标记符号，它们可以组合 使用。例如， 表示蓝色点划线， 使用。例如， “b-.”表示蓝色点划线， “y:d”表示 表示蓝色

10、点划线 表示 黄色虚线并用菱形符标记数据点。当选项省略时， 黄色虚线并用菱形符标记数据点。当选项省略时， MATLAB 规定，线型一律用实线，颜色将根据曲 规定，规定 线型一律用实线， 线的先后顺序依次。 线的先后顺序依次。要设置曲线样式可以在 plot 函数中加绘图选项，其 函数中加绘图选项， 要设置曲线样式可以在 函数中加绘图选项 调用格式为： 调用格式为： plot(x1,y1,选项 选项 1,x2,y2,选项 选项2,xn,yn,选项 选项 n) 选项 选项 选项在同一坐标内， 例 5-6 在同一坐标内，分别用不同线型和颜色绘制曲 ， 线 y1=0.2e-0.5xcos(4x) 和 y

11、2=2e-0.5xcos(x)，标记两 曲线交叉点。 曲线交叉点。 程序如下： 程序如下： x=linspace(0,2*pi,1000); y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); k=find(abs(y1-y2)1e-2); %查找 与 y2 相等点 查找 y1 与 相等点 查找 (近似相等 的下标 近似相等)的下标 近似相等 x1=x(k); %取 y1 与 y2 相等点的 坐标 相等点的 x坐标 取 与 相等点的 y3=0.2*exp(-0.5*x1).*cos(4*pi*x1); %求y1 与 y2

12、值 求 与 值 相等点的 y 坐标 相等点的 坐标 plot(x,y1,x,y2,k:,x1,y3,bp);5.1.4 图形标注与坐标控制 1图形标注 有关图形标注函数的调用格式为： 有关图形标注函数的调用格式为： title(图形名称 图形名称) 图形名称 xlabel(x 轴说明 轴说明) 轴说明 ylabel(y 轴说明 轴说明) 轴说明 text(x,y,图形说明 图形说明) 图形说明 legend(图例 图例 图例 1,图例 图例 图例 2,)函数中的说明文字，除使用标准的 字符外， 函数中的说明文字，除使用标准的 ASCII 字符外， 字符外 还可使用 LaTeX 格式的控制字符，

13、这样就可以在 格式的控制字符， 还可使用 格式的控制字符 图形上添加希腊字母、数学符号及公式等内容。 图形上添加希腊字母、数学符号及公式等内容。 例如， 例如，text(0.3,0.5,sin(omegat+beta)将得到 将得到 标注效果 sin(t+)。 标注效果 。例 5-7 在 0x2 区间内，绘制曲线 区间内，绘制曲线y1=2e-0.5x 和 y2=cos(4x)，并给图形添加图形标注。 ，并给图形添加图形标注。 程序如下： 程序如下： x=0:pi/100:2*pi; y1=2*exp(-0.5*x); y2=cos(4*pi*x); plot(x,y1,x,y2) title(

14、x from 0 to 2pi); %加图形标题 加图形标题 xlabel(Variable X); %加 X 轴说明 加 轴说明 ylabel(Variable Y); %加 Y 轴说明 加 轴说明 text(0.8,1.5,曲线 曲线 y1=2e-0.5x); %在指定位置 曲线 在指定位置 添加图形说明 text(2.5,1.1,曲线 曲线y2=cos(4pix); 曲线 legend(y1, y2) %加图例 加图例2坐标控制 axis 函数的调用格式为： 函数的调用格式为：函数的调用格式为 axis(xmin xmax ymin ymax zmin zmax) axis 函数功能丰富

15、，常用的格式还有： axis 函数功能丰富，常用的格式还有： 函数功能丰富 axis equal：纵、横坐标轴采用等长刻度。 横坐标轴采用等长刻度。 ： axis square：产生正方形坐标系 缺省为矩形 。 缺省为矩形)。 ：产生正方形坐标系(缺省为矩形 axis auto：使用缺省设置。 ：使用缺省设置。 axis off：取消坐标轴。 ：取消坐标轴。 axis on：显示坐标轴。 ：显示坐标轴。给坐标加网格线用 grid 命令来控制。grid on/off 命令 命令来控制。 给坐标加网格线用 命令来控制 命令 控制是画还是不画网格线，不带参数的 grid 命令 控制是画还是不画网格线

16、，不带参数的 命令 在两种状态之间进行切换。 在两种状态之间进行切换。 on/off 命令控 给坐标加边框用 box 命令来控制。box on/off 命令控 给坐标加边框用 box 命令来控制。 命令来控制 制是加还是不加边框线，不带参数的 box 命令在两 制是加还是不加边框线，不带参数的 命令在两 种状态之间进行切换。 种状态之间进行切换。在同一坐标中，可以绘制 3 个同心圆 个同心圆， 例 5-8 在同一坐标中，可以绘制 个同心圆，并加坐 标控制。 标控制。 程序如下： 程序如下： t=0:0.01:2*pi; x=exp(i*t); y=x;2*x;3*x; plot(y) grid

17、 on; %加网格线 加网格线 box on; %加坐标边框 加坐标边框 axis equal %坐标轴采用等刻度 坐标轴采用等刻度5.1.5 图形的可视化编辑 MATLAB 6.1 版本在图形窗口中提供了可视化的图 版本在图形窗口中提供了可视化的图 形编辑工具，利用图形窗口菜单栏或工具栏中的 形编辑工具， 有关命令可以完成对窗口中各种图形对象的编辑 处理。 处理。 在图形窗口上有一个菜单栏和工具栏。 在图形窗口上有一个菜单栏和工具栏。菜单栏包含 File、Edit、View、Insert、Tools、Window 和 、 、 、 、 、 和 Help 共 7 个菜单项，工具栏包含 个命令按钮

18、。 个菜单项， 个命令按钮。 共 个菜单项 工具栏包含 11 个命令按钮5.1.6 对函数自适应采样的绘图函数 fplot 函数的调用格式为：函数的调用格式为： 函数的调用格式为 fplot(fname,lims,tol,选项 选项) 选项 其中 fname 为函数名，以字符串形式出现，lims为 为函数名， 其中 为函数名 以字符串形式出现， 为 x,y 的取值范围，tol 为相对允许误差，其系统默认 的取值范围， 为相对允许误差 为相对允许误差， 的取值范围 值为 2e-3。选项定义与 函数相同。 值为 。选项定义与 plot 函数相同。 函数相同函数绘制 f(x)=cos(tan(x)的

19、曲线。 的曲线。 例 5-9 用fplot 函数绘制 函数绘制 的曲线 命令如下： 命令如下： fplot(cos(tan(pi*x), 0,1,1e-4)5.1.7 图形窗口的分割 subplot 函数的调用格式为： 函数的调用格式为： 函数的调用格式为 subplot(m,n,p) 该函数将当前图形窗口分成 m 个绘图区 个绘图区， 该函数将当前图形窗口分成 n 个绘图区， 即每行 n 个 即每行 个，共 m 行，区号按行优先编号， 行 区号按行优先编号， 且选定第 p 个区为当前活动区 个区为当前活动区。 且选定第 个区为当前活动区。在每一个绘 图区允许以不同的坐标系单独绘制图形。 图区

20、允许以不同的坐标系单独绘制图形。 在图形窗口中， 例 5-10 在图形窗口中，以子图形式同时绘制 多根曲线。 多根曲线。5.2 其他二维图形 5.2.1 其他坐标系下的二维数据曲线图 1对数坐标图形 MATLAB 提供了绘制对数和半对数坐标曲线的函数， 提供了绘制对数和半对数坐标曲线的函数， 提供了绘制对数和半对数坐标曲线的函数 调用格式为： 调用格式为： semilogx(x1,y1,选项 选项 1,x2,y2,选项 选项 2,) 选项 选项 semilogy(x1,y1,选项 选项 1,x2,y2,选项 选项 2,) 选项 选项 loglog(x1,y1,选项 选项 1,x2,y2,选项

21、选项 2,) 选项 选项 例5-11 绘制 y=10x2 的对数坐标图并与直角线性坐标 绘制 的对数坐标图并与直角线性坐标 图进行比较。 图进行比较。2极坐标图 polar 函数用来绘制极坐标图，其调用格式为：函数用来绘制极坐标图， 函数用来绘制极坐标图 其调用格式为： polar(theta,rho,选项 选项) 选项 其中 theta 为极坐标极角，rho为极坐标矢径，选项 为极坐标极角， 为极坐标矢径 为极坐标矢径，其中 为极坐标极角 的内容与 plot 函数相似。 函数相似。 的内容与 函数相似 绘制 r=sin(t)cos(t)的极坐标图 的极坐标图， 例5-12 绘制 r=sin(

22、t)cos(t)的极坐标图，并标记数据 点。 程序如下： 程序如下： t=0:pi/50:2*pi; r=sin(t).*cos(t); polar(t,r,-*);5.2.2 二维统计分析图 在 MATLAB 中，二维统计分析图形很多，常见的有 中 二维统计分析图形很多， 条形图、阶梯图、杆图和填充图等，所采用的函 条形图、阶梯图、杆图和填充图等， 数分别是： 数分别是： bar(x,y,选项 选项) 选项 stairs(x,y,选项 选项) 选项 stem(x,y,选项 选项) 选项 fill(x1,y1,选项 选项1,x2,y2,选项 选项 2,) 选项 选项分别以条形图、阶梯图、 例

23、5-13 分别以条形图、阶梯图、杆图和填充图形式 绘制曲线 y=2sin(x)。 绘制曲线 。 程序如下： 程序如下： x=0:pi/10:2*pi; y=2*sin(x); subplot(2,2,1);bar(x,y,g); title(bar(x,y,g);axis(0,7,-2,2); subplot(2,2,2);stairs(x,y,b); title(stairs(x,y,b);axis(0,7,-2,2); subplot(2,2,3);stem(x,y,k); title(stem(x,y,k);axis(0,7,-2,2); subplot(2,2,4);fill(x,y,

24、y); title(fill(x,y,y);axis(0,7,-2,2);MATLAB 提供的统计分析绘图函数还有很多，例如， 提供的统计分析绘图函数还有很多，例如， 提供的统计分析绘图函数还有很多 用来表示各元素占总和的百分比的饼图、 用来表示各元素占总和的百分比的饼图、复数的 相量图等等。 相量图等等。 绘制图形：例 5-14 绘制图形： (1) 某企业全年各季度的产值 单位：万元 分别为： 某企业全年各季度的产值(单位 万元)分别为 单位： 分别为： 2347,1827,2043,3025，试用饼图作统计分析。 ，试用饼图作统计分析。 (2) 绘制复数的相量图：7+2.9i、2-3i 和

25、-1.5-6i。 绘制复数的相量图： 、 和 。程序如下： 程序如下： subplot(1,2,1); pie(2347,1827,2043,3025); title(饼图 饼图); 饼图 legend(一季度 二季度 三季度 四季度 一季度,二季度 三季度,四季度 一季度 二季度,三季度 四季度); subplot(1,2,2); compass(7+2.9i,2-3i,-1.5-6i); title(相量图 相量图); 相量图5.3 隐函数绘图 MATLAB 提供了一个 提供了一个 ezplot 函数绘制隐函数图形， 函数绘制隐函数图形， 提供了一个 函数绘制隐函数图形 下面介绍其用法。

26、下面介绍其用法。 (1) 对于函数 = f(x)，ezplot 函数的调用格式为： 对于函数 f 函数的调用格式为：， 函数的调用格式为 ezplot(f)：在默认区间 绘制 f 的图形。 ：在默认区间-2x2 绘制 = f(x)的图形。 绘制 的图形 ezplot(f, a,b)：在区间 绘制 f 的图形。 ：在区间 axb 绘制 = f(x)的图形。 绘制 的图形(2) 对于隐函数 = f(x,y)，ezplot 函数的调用格式为： 对于隐函数 f 函数的调用格式为： ， 函数的调用格式为 ezplot(f)：在默认区间 ：在默认区间-2x2 和-2y2 绘制 和 绘制 f(x,y) =

27、0 的图形。 的图形。 的图形 ezplot(f, xmin,xmax,ymin,ymax)：在区间 ： xminxxmax 和 yminyymax绘制 绘制 f(x,y) = 0 的图 和 绘制 的图 形。 ezplot(f, a,b)：在区间 绘制 f(x,y) = 0 ：在区间 axb 和 ay b 绘制 和 绘制 的图形。 的图形。(3) 对于参数方程 = x(t)和 y = y(t)，ezplot 函数的调 对于参数方程 x 和 ， 函数的调 用格式为： 用格式为： ezplot(x,y)：在默认区间 绘制 x=x(t)和 y=y(t) ：在默认区间 0t2 绘制 绘制 和 的图形。

28、 的图形。 ezplot(x,y, tmin,tmax)：在区间 ：在区间 tmin t tmax 绘 绘 的图形。 制 x=x(t)和 y=y(t)的图形。 和 的图形隐函数绘图应用举例。 例 5-15 隐函数绘图应用举例。 程序如下： 程序如下： subplot(2,2,1); ezplot(x2+y2-9);axis equal subplot(2,2,2); ezplot(x3+y3-5*x*y+1/5) subplot(2,2,3); ezplot(cos(tan(pi*x), 0,1) subplot(2,2,4); ezplot(8*cos(t),4*sqrt(2)*sin(t)

29、,0,2*pi)5.4 三维图形 5.4.1 三维曲线 plot3 函数与 函数与 plot 函数用法十分相似，其调用格式为： 函数用法十分相似， 函数与 函数用法十分相似 其调用格式为： plot3(x1,y1,z1,选项 选项1,x2,y2,z2,选项 选项 2,xn,yn,zn,选项 选项 选项 选项 n) 其中每一组 x,y,z 组成一组曲线的坐标参数，选项的 组成一组曲线的坐标参数， 其中每一组 组成一组曲线的坐标参数 定义和 plot函数相同。当 x,y,z 是同维向量时，则 函数相同。 是同维向量时，定义和 函数相同 是同维向量时 x,y,z 对应元素构成一条三维曲线。当 x,y

30、,z 是同维 对应元素构成一条三维曲线。 是同维 矩阵时，则以 x,y,z 对应列元素绘制三维曲线，曲 对应列元素绘制三维曲线，矩阵时，则以 对应列元素绘制三维曲线 线条数等于矩阵列数。 线条数等于矩阵列数。绘制三维曲线。 例 5-16 绘制三维曲线。 程序如下： 程序如下： t=0:pi/100:20*pi; x=sin(t); y=cos(t); z=t.*sin(t).*cos(t); plot3(x,y,z); title(Line in 3-D Space); xlabel(X);ylabel(Y);zlabel(Z); grid on;5.4.2 三维曲面 1产生三维数据 在 MA

31、TLAB 中，利用 中 利用 meshgrid 函数产生平面区域内 函数产生平面区域内 的网格坐标矩阵。其格式为： 的网格坐标矩阵。其格式为： x=a:d1:b; y=c:d2:d; X,Y=meshgrid(x,y); 语句执行后，矩阵 X 的每一行都是向量 的每一行都是向量 x， 语句执行后，矩阵 的每一行都是向量 ，行数等于 向量 y 的元素的个数 矩阵 Y 的每一列都是向量 的元素的个数， 的每一列都是向量 y， 向量 的元素的个数，矩阵 的每一列都是向量 ， 列数等于向量 x 的元素的个数 的元素的个数。列数等于向量 的元素的个数。2绘制三维曲面的函数 surf 函数和 函数和 me

32、sh 函数的调用格式为： 函数的调用格式为： 函数和 函数的调用格式为 mesh(x,y,z,c) surf(x,y,z,c) 一般情况下， 是维数相同的矩阵。 是网格坐 一般情况下，x,y,z 是维数相同的矩阵。x,y 是网格坐 是维数相同的矩阵 标矩阵， 是网格点上的高度矩阵 用于指定在 是网格点上的高度矩阵， 标矩阵，z 是网格点上的高度矩阵，c 用于指定在 不同高度下的颜色范围。 不同高度下的颜色范围。绘制三维曲面图 z=sin(x+sin(y)-x/10。 例 5-17 绘制三维曲面图 。 程序如下： 程序如下： x,y=meshgrid(0:0.25:4*pi); z=sin(x+

33、sin(y)-x/10; mesh(x,y,z); axis(0 4*pi 0 4*pi -2.5 1); 此外，还有带等高线的三维网格曲面函数 meshc 和 此外，还有带等高线的三维网格曲面函数 和 带底座的三维网格曲面函数meshz。其用法与 带底座的三维网格曲面函数 。 mesh 类似，不同的是 类似， 还在 xy 平面上绘制曲 类似 不同的是 meshc 还在 平面上绘制曲 还在 面在 z 轴方向的等高线，meshz 还在 平面上绘制 面在 轴方向的等高线， 还在 xy 平面上绘制 轴方向的等高线 还在 曲面的底座。 曲面的底座。平面内选择区域-8,8-8,8，绘制 种三维曲 例 5

34、-18 在xy 平面内选择区域 平面内选择区域 ，绘制 4 种三维曲 面图。 面图。 程序如下： 程序如下： x,y=meshgrid(-8:0.5:8); z=sin(sqrt(x.2+y.2)./sqrt(x.2+y.2+eps); subplot(2,2,1); mesh(x,y,z); title(mesh(x,y,z) subplot(2,2,2); meshc(x,y,z); title(meshc(x,y,z) subplot(2,2,3); meshz(x,y,z) title(meshz(x,y,z) subplot(2,2,4); surf(x,y,z); title(su

35、rf(x,y,z)3标准三维曲面 sphere 函数的调用格式为： 函数的调用格式为： 函数的调用格式为 x,y,z=sphere(n) cylinder 函数的调用格式为： 函数的调用格式为： 函数的调用格式为 x,y,z= cylinder(R,n) MATLAB 还有一个 还有一个 peaks 函数，称为多峰函数，常 函数，称为多峰函数， 还有一个 用于三维曲面的演示。 用于三维曲面的演示。绘制标准三维曲面图形。 例 5-19 绘制标准三维曲面图形。 程序如下： 程序如下： t=0:pi/20:2*pi; x,y,z= cylinder(2+sin(t),30); subplot(2,2,1); surf(x,y,z); subplot(2,2,2); x,y,z=sphere; surf(x,y,z); subplot(2,1,2); x,y,z=peaks(30); surf(x,y,z);1