Matlab二维图形和三维图形的创建
数学2-用MATLAB绘制二维-三维图形(lq)
x=0:0.1*pi:2*pi; y=sin(x); z=cos(x); plot(x,y,'b-',x,z,'k.-','linewidth',3,'markersize',15)
axis([-0.2*pi 2.2*pi -1.2 1.2]) %重置坐标范围 grid xlabel('variable\it{x}') %标记横坐标,\it{x}表示x斜体 ylabel('variable\it{y}') %标记纵坐标轴
例:在图形窗口极坐标方程 r=2sin2θ*cos2θ的图形
h2=figure; %打开第二个图形窗口 theta=linspace(0,2*pi); %linspace()函数等分角,默认100等分 rho=sin(2*theta).*cos(2*theta); %生成相应极坐标方程的极径rho向量 polar(theta,rho,‘r’) %绘制相应的极坐标方程图形 title(‘polar plot’) %添加标题 如果想对第二个图形加粗的话,可以用如下命令 set(h2,'linewidth',3)
其中,点线的颜色代码与线型代码 表2-1 颜色代码表
表2-2 点型和线型代码
meshz在matlab中的用法
MATLAB中的Meshz函数用法简介MATLAB是一种强大的科学计算软件,其中的meshz函数是在三维图形中创建网格形状的常用函数之一。
它可以将一个二维或三维的矩阵数据在三维空间中绘制出来,并以网格形式进行可视化展示。
本文将详细介绍meshz函数的语法、用法和一些应用示例。
语法meshz函数的基本语法格式如下:meshz(X, Y, Z)其中:•X代表自变量x的数值矩阵,可以为一维或二维矩阵。
•Y代表自变量y的数值矩阵,可以为一维或二维矩阵。
•Z代表因变量z的数值矩阵,维度应与X和Y相匹配。
使用示例下面通过一些具体的示例来说明meshz函数的使用方法。
示例一:二维数据的网格绘制首先,我们生成一个二维数据,使用[X, Y] = meshgrid函数生成x轴和y轴的网格点坐标,然后计算每个网格点上的z轴数据。
最后,调用meshz函数进行可视化展示。
[X, Y] = meshgrid(-2:0.2:2, -2:0.2:2);Z = X.*exp(-X.^2 - Y.^2);meshz(X, Y, Z);运行上述代码后,将得到一个以x和y为坐标轴,z为高度的网格形状图。
示例二:三维数据的网格绘制如果要绘制三维数据,则需要在生成网格点坐标后,计算每个网格点上的z轴数据。
下面的示例中,我们使用[X, Y, Z] = meshgrid函数生成x、y和z轴的网格点坐标,然后根据一定的规则计算每个网格点上的z轴数值。
[X, Y, Z] = meshgrid(-2:0.2:2, -2:0.2:2, -2:0.2:2);V = X.*exp(-X.^2 - Y.^2 - Z.^2);meshz(X, Y, Z, V);运行上述代码后,将得到一个以x、y和z为坐标轴,V为高度的网格形状图。
注意事项1.meshz函数仅适用于二维或三维的数值数据,对于其他数据类型会导致错误。
2.数据矩阵的维度需要与网格点的坐标大小相匹配,否则会引发尺寸不匹配错误。
实验二MATLAB绘制图形
grid on %在所画出的图形坐标中加入栅格
绘制图形如下
50
10
1
0.8
40
10
0.6
0.4
30
10
0.2
0
1020
-0.2
-0.4
1010
-0.6
-0.8
0
10
-1
-2
0
2
-2
0
2
10
10
10
10
10
10
如果在图中不加栅格
程序如下:
clear x=logspace(-1,2);%在10^(-1)到10^2之间产生50个 对数等分的行向量 subplot(121); loglog(x,10*exp(x),'-p') subplot(122); semilogx(x,cos(10.^x))
(2)plot(x,y): 基本格式,x和y可为向量或矩阵. 1. 如果x,y是同维向量,以x元素为横坐标,以y元素 为纵坐标绘图. 2. 如果x是向量,y是有一维与x元素数量相等的矩阵, 则以x为共同横坐标, y元素为纵坐标绘图,曲线数目 为y的另一维数. 3. 如果x,y是同维矩阵,则按列以x,y对应列元素为 横、纵坐标绘图,曲线数目等于矩阵列数.
y=2*exp(-0.5*x).*cos(4*pi*x);
2
plot(x,y)
1.5
1
0.5
0
-0.5
-1
-1.5
-2
0
1
2
3
4
5
6
7
例4 绘制曲线
t=(0:0.1:2*pi);
x=t.*sin(3*t);
y=t.*sin(t).*sin(t);
实验五 MATLAB二维、三维图形的绘制
实验五 MATLAB二维、三维图形的绘制一、实验目的1.掌握二维、三维图形的绘制;2.掌握特殊二维图形的绘制;3.掌握绘图参数的设置;4.了解并学习简单动画的制作。
二、实验内容1.运行下列程序,学会并掌握标题、坐标轴标签和网格线的设置方法x=0:1:10;y=x.^2-10*x+6;plot(x,y);title ('Plot of y=x.^2-10*x+6');xlabel ('x');ylabel ('y');grid on;2.运行下列程序,学会并掌握线型、点型、颜色的设置方法x = -pi:pi/20:pi;y1 = sin(x);y2 = cos(x);plot(x,y1,'bo',x,y2,'r:');title('线型、点型和颜色');xlabel('时间'),ylabel('Y');grid on;3.同一坐标系内多条曲线的绘制1)使用 plot(x,[y1;y2;…])x = -pi:pi/20:pi;y1 = sin(x);y2 = cos(x);plot(x,[y1;y2]);legend('sin x','cos x');2)使用hold命令x = -pi:pi/20:pi;y1 = sin(x);y2 = cos(x);plot(x,y1);hold on;plot(x,y2,‘r’);3)在plot后使用多输入变量x = -2*pi:pi/20:2*pi;y1 = 2*sin(x);y2 = 2*cos(x);plot(x,y1,'ro',x,y2,'b:');title('线型、点型和颜色');xlabel('时间'),ylabel('Y');4) 使用plotyy命令x = -pi:pi/20:pi;y1 = sin(x);y2 = 5*cos(x);plotyy(x,y1,x,y2);grid on;gtext(‘sinx’) ; gtext(‘5cosx’) ;4.子图形窗口的绘制subplot(2,1,1);x= -pi:pi/20:pi;y=sin(x);plot(x,y) ; grid on;title('正弦曲线');subplot(2,1,2);x= -pi:pi/20:pi;y=cos(x);plot(x,y); grid on;title('余弦曲线');5.对数坐标图形x=0:0.1:10;y=x.^2 -10.*x +25;subplot(2,2,1);plot(x,y); grid on;xlabel('a) x、y轴线性刻度');subplot(2,2,2);semilogx(x,y); grid on;xlabel('b) x轴对数刻度、y轴线性刻度');subplot(2,2,3);semilogy(x,y); grid on;xlabel('c) x轴线性刻度、y轴对数刻度');subplot(2,2,4);loglog(x,y); grid on;xlabel(‘d) x、y轴对数刻度');6.极坐标下的绘图theta = 0:pi/20:2*pi;r = 0.5+cos(theta);polar(theta,r);7.复数的绘图,并比较下面几种情况的不同1)t = 0:pi/20:6*pi;y = exp(0.1*t) .* (cos(t) + i * sin(t));plot(y);grid on ;title('Plot of Complex Function vs Time');xlabel('Real Part');ylabel('Imaginary Part');2)t = 0:pi/20:6*pi;y = exp(0.1*t) .* (cos(t) + i * sin(t));plot(t, y);grid on ;title('Plot of Complex Function vs Time');xlabel('t');ylabel('y(t)');3)t = 0:pi/20:6*pi;y = exp(0.1*t) .* (cos(t) + i * sin(t));plot(t, real(y),'b-');grid on;hold on;plot(t, imag(y),'r-');title('Plot of Complex Function vs Time');xlabel('t');ylabel('y(t)');legend('real','imaginary');hold off;4)t = 0:pi/20:6*pi;y = exp(0.1*t) .* (cos(t) + i * sin(t));polar(angle(y),abs(y));title('Plot of Complex Function');8.特殊二维图形的绘制1)x = [1 2 3 4 5 6];y = [2 6 8 7 8 5];stem(x,y);title('Example of a Stem Plot');xlabel('x');ylabel('y');axis([0 7 0 10]);将上述程序中的stem语句换为stairs、bar、barh和compass,即可实现阶梯图、条形图、罗盘图的绘制。
在MATLAB中用数据文件制作二维、三维图形
由于该 组数据 中 。 n 际上分别是 11。 2 向 巾, 实 8 1 0的
量 根 据作图的维数要求 。 z应为 2 " 8的矩阵 作 图如图 01
3 示 所 对应. m文件 为:
利 用 p毗 函数 即可 完成 作图 . l 如图 1 示 。编 写. 件如 所 m文
下:
fr _ ::。 0 o _i 2 l 0
de r a
c ) Ⅱk ) 【J= 0 a= (k: k )c 1 b
ed n
l d m. t 产 装载文 本数据 - o  ̄ a t x / 一 将 l 18 x 00的矩阵 a 转换 为 3 0 3的矩阵 b . - a 6x l n q三 个向量 备 占一列 。将 3 0 l 6 x 的矩 阵 z 即 h ( 的第 3 ) 列 转换为 2x 8 0 1 的矩
据 转换 。 ( 作图。 4)
图。 用于 局部放电模式 的识别 。 例如 , 二维 Q N 巾一 — 、 N谱 图。
三维 N O — 一巾 谱图 。
如要 将格 式 为【 n 巾 , 一 。 。 巾 n】接有 10 0 0个 数 据。 文件 名 为 x y的文 本文 件数据 读 出 , 并做 出相应 的二 维 巾一 N谱 图 我们 先将该 1 1 0 " 0 0的矩 阵转 换为 5 0 2的矩 0*
n h w t rc s aa f o te tx e a d po D a d 3 f ue i te D na i - i s i h e Ab la t T i a e o u ̄ o o o p o es d t rm h e tfl ,n lt2 n D g rs wt h a ild sh sr c h s p p rfc s
第5章matlab绘制二维图形及三维图形的方法
实验四
专业:电子信息工程2班姓名:李书杰学号:3121003210
一、实验目的
1.掌握绘制二维图形及三维图形的方法。
2.掌握图形控制与修饰处理的方法。
3.了解图像处理及动画制作的基本方法。
二、实验内容
1.绘制下列图形曲线。
(1)y=x-x^3/3! (2)x^2+2Y^2=64
解:程序如下
2.设y=1/(1+e^-t),-pi<=t<=pi,在同一个图形窗口中采用子图的形式绘制条形图、阶梯图、杆图和对数坐标等不同图形,并对不同图形加标注说明。
解:程序如下
3.绘制下列极坐标图。
(1)ρ=5cosθ+4 (2)γ=5sin^2φ/cosφ,-π/3<φ<π/3 解:程序如下
思考练习:
2.绘制下列曲线
(1)y=1/2πe^(-x^2/2) (2)x=tsint y=tcost
解:程序如下
(1)
结果如下:
(2)
结果如下:
3.在同一坐标中绘制下列两条曲线并标注两曲线交叉点。
(1)y=2x-0.5
(2)x=sin(3t)cost
Y=sin(3t)sint
解:程序如下
4.分别用plot和fplot函数绘制y=sin(1/x)的曲线,分析两曲线的差别。
解:程序如下
结果如下:
5.绘制下列极坐标图:
(1)p=12/sqrt(θ) (2)γ=3asinφcosφ/(sin^3φ+cos^3φ)解:程序如下
结果如下:。
MATLAB图形及三维动画仿真设计
[X,Y]=meshgrid(x,y);
p=sqrt(4-X.^2/9-Y.^2/4);
subplot(3,2,1);mesh(p); %有网格 格子图,色调方向:有上
subplot(3,2,2);surf(p); %默认的方向:色调方向:有上到
subplot(3,2,3);surfc(p); %带等高线;色调方向:有上到下
axis([0 6 8 18])
z=f(x,y) — 根据x,y坐标找出z的高度 例:绘制z=x2+y2的三维网线图形
>>x=-5:5; y=x; >>[X,Y]=meshgrid(x,y); >>Z=X.^2+Y.^2 ; >>mesh(X,Y,Z)
(五)三维曲面图
■ surf —— 三维曲面绘图函数,与网格图看起来一样 与三维网线图的区别: 网线图:线条有颜色,内部是黑色的(无颜色) 曲面图:线条是黑色的,内部有颜色(把线条之间的空
■
■
■ fill3(x,y,z,"w"); ■ hold on; ■ plot3(x,y,z,"ro") ■ grid on ■ xlabel("x轴"),ylabel("y轴"),zlabel("z轴") ■ axis([0 10 0 10 0 10])
(四)三维网格图
所谓网格图,是指把相邻的数据点连接起来形成的网状曲 面。利用在X-Y平面的矩形网格点上的Z轴坐标值,Matlab 定义了一个网格曲面。三维网格图的形成原理为:在X-Y平 面上指定一个长方形区域,采用与坐标轴平行的直线将其分 格;计算矩形网格点上的函数值,即Z轴的值,得到三维空 间的数据点;将这些数据点分别用处于X-Z或者平行面内的 曲线和处于Y-Z或者平行面的曲面连接起来,即形成网格图。 网格图对于显示大型的数值矩阵很有用处。
使用matlab绘制三维图形的方法
使用matlab绘制三维图形的方法要使用MATLAB绘制三维图形,首先需要了解MATLAB中的三维绘图函数和绘图选项。
下面将介绍一些常用的绘制三维图形的方法。
1.绘制基本的三维图形要绘制基本的三维图形,可以使用以下函数:- plot3(函数:用于在三维坐标系中绘制线条。
- scatter3(函数:用于在三维坐标系中绘制散点图。
- surf(函数:用于绘制三维曲面图。
- mesh(函数:用于绘制三维网格图。
- bar3(函数:用于绘制三维条形图。
- contour3(函数:用于绘制三维等高线图。
例如,下面的代码演示了如何使用plot3(函数绘制一个三维线条图:```x = linspace(0, 2*pi, 100);y = sin(x);z = cos(x);plot3(x, y, z, 'LineWidth', 2);xlabel('X');ylabel('Y');zlabel('Z');title('3D Line Plot');```2.添加颜色和纹理在绘制三维图形时,可以使用颜色和纹理来增加图形的信息。
MATLAB 提供了一系列函数来处理颜色和纹理,如:- colormap(函数:用于设置颜色映射。
- caxis(函数:用于设置坐标轴范围。
- shading(函数:用于设置颜色插值方法。
- texturemap(函数:用于设置纹理映射方法。
例如,下面的代码展示了如何使用纹理映射来绘制一个球体:```[X, Y, Z] = sphere(50);C = colormap('jet');surface(X, Y, Z, 'FaceColor', 'texturemap', 'CData', C);axis equal;```3.绘制多个数据集要在同一张图中绘制多个数据集,可以使用hold on和hold off命令。
国家开放大学《Matlab语言及其应用》实验报告(第三章--绘制二维和三维图形)
——绘制二维和三维图形
姓名:学号:
实验名称
绘制二维和三维图形
实验目标
利用Matlab常见函数完成二维图形的绘制和图形的标注;实现三维曲线和曲面图形的绘制。
实验要求
熟悉Matlab基本绘图函数、图形处理函数,了解三维曲线和曲面图形的绘制方法。
实验步骤
1、用Matlab基本绘图函数绘制二维图形:根据已知数据,用plot函数画出正弦函数曲线,并进行相应标注。
enon
实验内容
1.二维曲线绘图
例:精细指令实例
2.三维曲线绘图
【例】三维曲线绘图基本指令演示一:plot3
t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);
plot3(x,y,z,'b-',x,y,z, 'rd')三维曲线绘图(蓝实线和红菱形)
box on
legend('链','宝石')在右上角建立图例
subplot(121);
surf(x1,y1,z1);
subplot(122);
[x2,y2,z2]=sphere (30);
surf(x2,y2,z2);
clear;clf;
z=peaks;
subplot(1,2,1);mesh(z);% 透视
hidden off
subplot(1,2,2);mesh(z);%不透视
2、用三维曲线绘图基本指令plot 3绘制三维曲线图:t=0~2pi;x=sin(t);y=cos(t);z=cos(2*t);用plot3函数画出关于x,y,z的三维曲线图,并适当加标注。
如何在Matlab中进行二维和三维绘图
如何在Matlab中进行二维和三维绘图在科学研究和工程领域,数据可视化是一项十分重要的任务,而Matlab作为一种功能强大的数值计算和数据分析软件,自然也提供了丰富的绘图功能。
本文将介绍如何在Matlab中进行二维和三维绘图,并探讨一些常见的绘图技巧和应用。
一、二维绘图Matlab中的二维绘图是最常见和基础的绘图任务之一。
在绘制二维图形时,我们通常会用到plot函数。
这个函数可以接受单个向量作为输入,将这个向量的值作为y轴上的数据点,自动生成与该向量长度相同的x轴坐标。
例如,我们可以用以下代码绘制一个简单的二维折线图:```x = 0:0.1:2*pi;y = sin(x);plot(x, y);```上述代码中,x参量取从0到2π的均匀间隔的值,而y则是根据x计算得到的sin函数值。
plot函数会自动根据输入绘制折线图,并添加相应的轴标签和图例。
在实际应用中,我们经常需要绘制多条曲线在同一个坐标系中进行对比分析。
可以通过在plot函数中传入多个x和y向量实现这一功能。
例如,我们可以通过以下代码绘制一个简单的双曲线图:```x = 0:0.1:2*pi;y1 = sin(x);y2 = cos(x);plot(x, y1, x, y2);```这样,就会在同一个坐标系中同时绘制sin曲线和cos曲线。
除了折线图,Matlab还支持其他常见的二维绘图类型,如散点图、柱状图和面积图等。
这些绘图类型可以通过不同的函数实现,例如scatter、bar和area等。
这里不再一一赘述,读者可以通过Matlab的帮助文档或官方网站了解更多的用法和示例。
二、三维绘图除了二维绘图,Matlab也提供了丰富的三维绘图功能,用于可视化更为复杂的数据和模型。
在绘制三维图形时,我们通常会用到surf函数。
这个函数可以接受两个二维矩阵作为输入,将这两个矩阵的值分别作为x、y轴上的坐标,而将第三个二维矩阵的值作为z轴上的数据点。
MATLAB画图(包括二维三维画图)
Matlab教程2_绘图_ 二维0Matlab绘图Matlab二维图形的绘制本章主要介绍2维图形对象的生成函数及图形控制函数的使用方法,还将简单地介绍一些图形的修饰与标注函数及操作和控制MAT LAB各种图形对象的方法.第一节图形窗口与坐标系一.图形窗口1.MATLAB在图形窗口中绘制或输出图形,因此图形窗口就像一张绘图纸.2. 在MATLA B下,每一个图形窗口有唯一的一个序号h,称为该图形窗口的句柄.MATLAB通过管理图形窗口的句柄来管理图形窗口;3.当前窗口句柄可以由MA TLAB函数gcf获得;4.在任何时刻,只有唯一的一个窗口是当前的图形窗口(活跃窗口);figure(h)----将句柄为h的窗口设置为当前窗口;5.打开图形窗口的方法有三种:1)调用绘图函数时自动打开;2)用File---New---Figure新建;3)figure命令打开,close命令关闭.在运行绘图程序前若已打开图形窗口,则绘图函数不再打开,而直接利用已打开的图形窗口;若运行程序前已存在多个图形窗口,并且没有指定哪个窗口为当前窗口时,则以最后使用过的窗口为当前窗口输出图形.6.窗口中的图形打印:用图形窗口的File菜单中的Print项.7.可以在图形窗口中设置图形对象的参数.具体方法是在图形窗口的Edit菜单中选择Pr opert ies项,打开图形对象的参数设置窗口,可以设置对象的属性.二.坐标系1.一个图形必须有其定位系统,即坐标系;2.在一个图形窗口中可以有多个坐标系,但只有一个当前的坐标系;3.每个坐标系都有唯一的标识符,即句柄值;4.当前坐标系句柄可以由MATLAB函数gca获得;5.使某个句柄标识的坐标系成为当前坐标系,可用如下函数:axes(h) h为指定坐标系句柄值.6.一些有关坐标轴的函数:1)定义坐标范围:一般MATL AB自动定义坐标范围,如用户认为设定的不合适,可用:axis([Xmin, Xmax, Ymin, Ymax]) 来重新设定;2) 坐标轴控制:MATLAB的缺省方式是在绘图时,将所在的坐标系也画出来,为隐去坐标系,可用axis off;axis on则显示坐标轴 (缺省值).3)通常MATL AB的坐标系是长方形,长宽比例大约是 4:3,为了得到一个正方形的坐标系可用:axis square4)坐标系横纵轴的比例是自动设置的,比例可能不一样,要得到相同比例的坐标系,可用:axis equal第二节二维图形的绘制一. plot函数plot 函数是最基本的绘图函数,其基本的调用格式为:1.plot(y)------绘制向量y对应于其元素序数的二维曲线图, 如果y为复数向量, 则绘制虚部对于实部的二维曲线图.例:绘制单矢量曲线图.y=[0 0.6 2.3 5 8.3 11.7 15 17.7 19.4 20];plot(y)由于y矢量有10个元素,x坐标自动定义为[1 2 3 4 5 6 7 8 9 10].图形为:2.plot(x,y)------ 绘制由x,y所确定的曲线.1)x,y是两组向量,且它们的长度相等,则plot(x,y)可以直观地绘出以x为横坐标,y为纵坐标的图形.如画正弦曲线:t=0:0.1:2*pi; y=sin(t);plot(t,y)2)当 plot(x,y)中,x是向量,y是矩阵时,则绘制y矩阵中各行或列对应于向量x的曲线.如果y阵中行的长度与x向量的长度相同,则以y的行数据作为一组绘图数据;如果y阵中列的长度与x向量的长度相同,则以y的列数据作为一组绘图数据;如果y阵中行, 列均与x向量的长度相同,则以y的每列数据作为一组绘图数据.例:下面的程序可同时绘出三条曲线.MATLAB在绘制多条曲线时,会按照一定的规律自动变化每条曲线的的颜色.x=0:pi/50:2*pi;y(1,:)=sin(x);y(2,:)=0.6*sin(x);y(3,:)=0.3*sin(x);plot(x,y)或者还可以这样用:x=0:pi/50:2*pi;y=[ sin(x); 0.6*sin(x); 0.3*sin(x)]; plot(x,y)3) 如果x,y是同样大小的矩阵,则plot(x,y)绘出y中各列相应于x中各列的图形.例:x(1,:)=0:pi/50:2*pi;x(2,:)=pi/4:pi/50:2*pi+pi/4;x(3,:)=pi/2:pi/50:2*pi+pi/2;y(1,:)=sin(x(1,:));y(2,:)=0.6*sin(x(2,:));y(3,:)=0.3*sin(x(3,:));plot(x,y)x=x'; y=y';plot(x,y)在这个例子中,x------ 3x101,y------3x101,所以第一个p lot按列画出 101条曲线,每条3个点;而x'------101x3,y'------ 101x3,所以第二个p lot 按列画出3条曲线,每条101个点.3. 多组变量绘图:plot(x1, y1, 选项1, x2, y2, 选项2, ……)上面的plo t格式中,选项是指为了区分多条画出曲线的颜色,线型及标记点而设定的曲线的属性.MATLAB在多组变量绘图时,可将曲线以不同的颜色,不同的线型及标记点表示出来.这些选项如下表所示:各种颜色属性选项'r' 红色'm' 粉红'g' 绿色'c' 青色'b' 兰色'w' 白色'y' 黄色'k' 黑色各种线型属性选项'-' 实线'--' 虚线':' 点线'-.' 点划线'.' 用点号绘制各数据点'^' 用上三角绘制各数据点'+' 用'+'号绘制各数据点'v' 用下三角绘制各数据点'*' 用'*'号绘制各数据点'>' 用右三角绘制各数据点' .' 用'.'号绘制各数据点'<' 用左三角绘制各数据点's'或squar用正方形绘制各数据点'p' 用五角星绘制各数据点'd'或diamo nd用菱形绘制各数据点'h' 用六角星绘制各数据点这些选项可以连在一起用,如:'-.g'表示绘制绿色的点划线,'g+'表示用绿色的 '+'号绘制曲线.注意:1)表示属性的符号必须放在同一个字符串中;2)可同时指定2~3个属性;3)与先后顺序无关;4)指定的属性中,同一种属性不能有两个以上.例:t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1, '-r',t,y2, ':g',t,y3, '*b')该程序还可以按下面的方式写:t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1, '-r')hold onplot(t,y2, ':g')plot(t,y3, '*b')hold off注:在MATLA B中,如画图前已有打开的图形窗口,则再画图系统将自动擦掉坐标系中已有的图形对象,但设置了ho ld on后,可以保持坐标系中已绘出的图形.还可以进一步设置包括线的宽度(LineWi dth),标记点的边缘颜色(Marker EdgeC olor),填充颜色 (Marker FaceC olor)及标记点的大小(Marker Size)等其它绘图属性.例:设置绘图线的线型,颜色,宽度,标记点的颜色及大小.t=0:pi/20:pi;y=sin(4*t).*sin(t)/2;plot(t,y,'-bs','LineWi dth',2,... %设置线的宽度为2'Marker EdgeC olor','k',... %设置标记点边缘颜色为黑色'Marker FaceC olor','y',... %设置标记点填充颜色为黄色'Marker Size',10) %设置标记点的尺寸为104. 双Y轴绘图:plotyy()函数.其调用格式为:plotyy(x1,y1,x2,y2)绘制由x1,y1和x2,y2确定的两组曲线,其中x1, y1的坐标轴在图形窗口的左侧,x2,y2的坐标轴在图形窗口的右侧.Plotyy(x1,y1,x2,y2, 'functi on1','functi on2')功能同上,functi on是指那些绘图函数,如plot,semilo gx, loglog等.例如:在一个图形窗口中绘制双Y轴曲线.x=0:0.3:12;y=exp(-0.3*x).*sin(x)+0.5;plotyy(x,y,x,y,'plot','stem')stem: 绘制stem形式的曲线(上端带圈的竖线).绘图结果:两条图线自动用不同的颜色区分,两个坐标的颜色与图线的颜色相对应,左边的Y轴坐标对应的是plot形式的曲线,右边的Y坐标对应的是stem形式的曲线.二.对数坐标图绘制函数:在对数坐标图的绘制中,有三种绘图函数:semilo gx,semilo gy和lo glog函数.1) semilo gx( )绘制以X轴为对数坐标轴的图. 调用格式为:semilo gx(x,y,'属性选项') 其中属性选项同plot函数.该函数只对横坐标进行对数变换,纵坐标仍为线性坐标.2) semilo gy( )绘制以Y轴为对数坐标轴的图. 调用格式为:semilo gy(x,y,'属性选项')该函数只对纵坐标进行对数变换,横坐标仍为线性坐标. 35 3) loglog( )绘制X,Y轴均为对数坐标轴的图形. 其调用格式为:loglog(x,y,'属性选项') 该函数分别对横,纵坐标都进行对数变换.例:x=0:0.1:6*pi;y=cos(x/3)+1/9;subplo t(221), semilo gx(x,y);subplo t(222), semilo gy(x,y);subplo t(223), loglog(x,y);4) logspa ce( )函数可按对数等间距地分布来产生一个向量,其调用格式为: x=logspa ce(x1,x2,n)这里,x1 表示向量的起点;x2表示向量的终点;n表示需要产生向量点的个数(一般可以不给出,采用默认值50).在控制系统分析中一般采用这种方法来构成频率向量w.关于它的应用后面还要讲到.三.极坐标图的绘制函数:绘极坐标图可用pola r( )函数.其调用格式如下:polar(theta, rho,'属性选项')theta:角度向量,rho:幅值向量,属性内容与plo t函数基本一致.例如:极坐标模型为:3145/)/)cos((+ =θρ,],[πθ80∈则绘出极坐标图的程序为:theta=0:0.1:8*pi;p=cos((5*theta)/4)+1/3;polar(theta,p)四. 绘制多个子图:subplo t( )函数MATLAB允许在一个图形窗口上绘制多个子图(如对于多变量系统的输出),允许将窗口分成nxm个部分.分割图形窗口用subp lot函数来实现,其调用格式为:subplo t(n,m,k) 或subpl ot(nmk)n,m分别表示将窗口分割的行数和列数,k表示要画图部分的代号,表示第几个图形,nmk三个数可以连写,中间不用符号分开.例如:将窗口划分成2x2=4个部分,可以这样写:subplo t(2,2,1),plot(……)subplo t(2,2,2),……subplo t(2,2,3),……subplo t(2,2,4),……注:subplo t函数没有画图功能,只是将窗口分割.第三节图形的修饰与标注MATLAB提供了一些特殊的函数修饰画出的图形,这些函数如下: 1)坐标轴的标题:title函数其调用格式为:title('字符串')字符串可以写中文如:title('My own plot')2)坐标轴的说明:xlabel和ylab el函数格式:xlabel('字符串') ylabel(' 字符串')如:xlabel('This is my X axis') ylabel('My Y axis')3)图形说明文字:text和g text函数A.text函数:按指定位置在坐标系中写出说明文字.格式为:text(x1, y1, '字符串', '选项')x1,y1为指定点的坐标;'字符串'为要标注的文字;'选项'决定x1,y1 的坐标单位,如没有选项,则x1,y1的坐标单位和图中一致;如选项为'sc', 则x1,y1表示规范化窗口的相对坐标,其范围为0到1. (1,1)规范化窗口(0,0)如:text(1,2, '正弦曲线')B.gtext函数:按照鼠标点按位置写出说明文字.格式为:gtext('字符串')当调用这个函数时,在图形窗口中出现一个随鼠标移动的大十字交叉线, 移动鼠标将十字线的交叉点移动到适当的位置,点击鼠标左键,gtext参数中的字符串就标注在该位置上.5) 给图形加网格:grid函数在调用时直接写grid即可.上面的函数的应用实例:例: 在图形中加注坐标轴标识和标题及在图形中的任意位置加入文本.t=0:pi/100:2*pi;y=sin(t);plot(t,y),grid, axis([0 2*pi -1 1]);xlabel('0 leq itt rm leq pi','FontSi ze',12);ylabel('sin(t)','FontSi ze',12);title('正弦函数图形','FontNa me',' 隶书' ,'FontSi ze',16); text(pi,sin(pi),'←sin(t)=0','FontSi ze',12);text(3*pi/4,sin(3*pi/4),'←sin(t)=0.707', 'FontSi ze',12); text(5*pi/4,sin(5*pi/4),' sin(t)=-0.707→','FontSi ze',12,'Horizo ntalA lignm ent','right')5) 在图形中添加图例框:legend函数其调用格式为:A.legend('字符串1', '字符串2', ……)以字符串1,字符串2……作为图形标注的图例.B.legend('字符串1', '字符串2', ……, pos)pos指定图例框显示的位置.图例框被预定了6个显示位置:0------ 取最佳位置;1------右上角(缺省值);2------左上角;3------左下角;4------右下角;-1------图的右侧.例:在图形中添加图例.x=0:pi/10:2*pi;y1=sin(x);y2=0.6*sin(x);y3=0.3*sin(x);plot(x,y1,x,y2,'-o',x,y3,'-*') legend('曲线1','曲线2','曲线3')6) 用鼠标点选屏幕上的点:ginput函数格式为: [x, y, button]=ginput(n)n为所选择点的个数;x,y均为向量,x为所选n个点的横坐标; y为所选n个点的纵坐标.button为n维向量,是所选n个点所对应的鼠标键的标号: 1----左键; 2----中键; 3----右键. 可用不同的鼠标键来选点,以区别所选的点.此语句可以放在绘图语句之后,它可在绘出的图形上操作,选择你所感兴趣的点,如峰值点,达到稳态值的点等,给出点的坐标,可求出系统的性能指标.第四节 MATLAB下图形对象的修改MATLAB图形对象是指图形系统中最基本,最底层的单元,这些对象包括:屏幕(Root),图形窗口(Figure s),坐标轴(Axes),控件(Uicont rol),菜单(Uimenu),线(Lines),块(Patche s),面(Surfac e),图像 (Images),文本(Text)等等.根据各对象的相互关系,可以构成如下所示的树状层次:RootFigure sAxes Uicont rol Uimenu Uicont extme nu (对象菜单) Images Line PatchSurfac e Text对各种图形对象进行修改和控制,要使用MAT LAB的图形对象句柄 (Handle).在MATLA B中, 每个图形对象创立时,就被赋予了唯一的标识,这个标识就是该对象的句柄. 句柄的值可以是一个数,也可以是一个矢量.如每个计算机的根对象只有一个,它的句柄总是0,图形窗口的句柄总是正整数,它标识了图形窗口的序号等. 利用句柄可以操纵一个已经存在的图形对象的属性,特别是对指定图形对象句柄的操作不会影响同时存在的其它图形对象,这一点是非常重要的.一.对图形对象的修改可以用下面函数:1)set函数:用于设置句柄所指的图形对象的属性.Set函数的格式为:set(句柄, 属性名1, 属性值1, 属性名2, 属性值2, ……)例:h=plot(x,y)set(h, 'Color', [1,0,0])------将句柄所指曲线的颜色设为红色.2)get函数:获取指定句柄的图形对象指定属性的当前值.格式为:get(句柄, '属性名')如: get(gca, 'Xcolor')------获得X轴的当前颜色属性值.执行后可返回X轴的当前颜色属性值[0,0,0](黑色).3)如果没有设置句柄,则可以使用下列函数获得:gcf:获得当前图形窗口的句柄;gca:获得当前坐标轴对象的句柄;gco:获得当前对象的句柄.如:A. 要对图形窗口的底色进行修改,可用set(gcf, 'Color', [1,1,1])------将图形窗口底色设为白色B. 要把当前X轴的颜色改为绿色,可用set(gca, 'Xcolor', [0,1,0])C. 还可对坐标轴的显示刻度进行定义:t=-pi:pi/20:pi;y=sin(t);plot(t,y)set(gca,'xtick',[-pi:pi/2:pi],'xtickl abel',['-pi','-pi/2','0','pi /2','pi'])本例中用'xtick'属性设置x轴刻度的位置(从-pi~pi,间隔pi/2,共设置 5个点),用'xtickl abel'来指定刻度的值,由于通常习惯于用角度度量三角函数, 因此重新设置['-pi','-pi/2','0','pi/2','pi']5个刻度值.二.一些常用的属性如下:1)Box属性:决定图形坐标轴是否为方框形式,选项为'on'(有方框), 'off'(无方框);2)'ColorO rder'属性:设置多条曲线的颜色顺序,默认值为:[1 1 0;1 0 1;0 1 1;1 0 0;0 1 0;0 0 1]黄色粉色天蓝红色绿色兰色颜色向量还有:[1 1 1]------白色;[0 0 0]------黑色.3)坐标轴方向属性:'Xdir','Ydir','Zdir',其选项为:'normal'------正常'revers e'------反向4) 坐标轴颜色和线型属性:'Xcolor','Ycolor','Zcolor'------ 轴颜色, 值为颜色向量'LineWi dth'------ 轴的线宽,值为数字'Xgrid','Ygrid','Zgrid'------坐标轴上是否加网格,值为'on'和'off'.5)坐标轴的标尺属性:'Xtick','Ytick','Ztick'------ 标度的位置,值为向量'Xtickl abel','Ytickl abel','Ztickl abel'------ 轴上标度的符号,它的值为与标度位置向量同样大小(向量个数相同)的向量.6)字体设置属性:'FontAn gle'------ 设置字体角度,选项为: 'normal'------ 正常;'italic'------ 斜体;'obliqu e'------ 倾斜;'FontNa me'------ 字体名称;'FontSi ze'------ 字号大小'FontWe ight'------ 字体的轻重,选项为:'light', 'normal','bold'Matlab绘制三维图形三维曲线plot3函数与plo t函数用法十分相似,其调用格式为:plot3(x1,y1,z1,选项1,x2,y2,z2,选项2,…,xn,yn,zn,选项n)其中每一组x,y,z组成一组曲线的坐标参数,选项的定义和plot函数相同。
MATLAB实用教程第三讲
补充:医学图像处理
• MATLAB图像处理工具箱支持四种图像类型,分别为真 彩色图像(RGB)、索引色图像、灰度图像(I)和二值 图像(BW)。许多图像处理工作对图像类型有特定的要 求,比如要对一幅索引图滤波,首先要把它转换成真彩色 图像或灰度图像,直接滤波的结果是毫无意义的。 • 在图像处理工具箱中有类型转换函数可以实现各种图像类 型之间的相互转换,如rgb2gray、ind2gray、im2bw 等。 • 在MATLAB中,改变图像文件格式的方法是:先用阅读 函数将文件读入工作空间,然后用写入函数以适当的格式 保存文件。由此可以方便的实现医学图像格式与其他格式 之间的相互转换,以利于医学图像的传输和应用。
绘制多边形
任意多边形及多边形区域是用Patch对 象表示的,该对象由patch函数创建。 使用patch函数生成面片,有两种方 法,一种是顶点/颜色法,通过指定每个小 面的顶点坐标和颜色来创建;另一种是顶 点/小面法,通过指定所有顶点的坐标和各 小面的顶点和颜色来创建。
绘制文本
文本用Text对象表示,该对象由text函数创 建,text函数是创建文本图形对象的低级函数。使 用它将文本字符串放在指定的位置上。 下面创建多种形式的文本:
例:把一个TIF格式的图像li2.tif转换 成DICOM格式的图像li2.dcm。
x=imread('li2.tif'); %将转换文件读入工作空间 imshow(x) %如图(a)显示的TIF格式图像 st=dicomwrite(x,‘li2.dcm’); %转换为DICOM格式 I=dicomread('li2.dcm'); %读入DICOM文件 imshow(I) %如图(b)显示的DICOM格式图像
使用matlab绘制三维图形地方法
使用matlab绘制三维图形地方法使用matlab绘制三维图形的方法三维曲线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,z是同维矩阵时,则以x,y,z对应列元素绘制三维曲线,曲线条数等于矩阵列数。
例绘制三维曲线。
程序如下:t=0:pi/100:20*pi;x=sin(t);y=cos(t);z=t.*sin(t).*cos(t);plot3(x,y,z);gridtitle('Line in 3-D Space');xlabel('X');ylabel('Y');zlabel('Z');如下图:XLine in 3-D SpaceYZ三维曲面1.产生三维数据在MATLAB 中,利用meshgrid 函数产生平面区域内的网格坐标矩阵。
其格式为:x=a:d1:b; y=c:d2:d; [X,Y]=meshgrid(x,y);语句执行后,矩阵X 的每一行都是向量x ,行数等于向量y 的元素的个数,矩阵Y 的每一列都是向量y ,列数等于向量x 的元素的个数。
2.绘制三维曲面的函数surf 函数和mesh 函数的调用格式为:mesh(x,y,z,c):画网格曲面,将数据点在空间中描出,并连成网格。
surf(x,y,z,c):画完整曲面,将数据点所表示曲面画出。
一般情况下,x,y,z 是维数相同的矩阵。
x,y 是网格坐标矩阵,z 是网格点上的高度矩阵,c 用于指定在不同高度下的颜色范围。
例绘制三维曲面图z=sin(x+sin(y))-x/10。
程序如下:[x,y]=meshgrid(0:0.25:4*pi); %在[0,4pi]×[0,4pi]区域生成网格坐标 z=sin(x+sin(y))-x/10; mesh(x,y,z);axis([0 4*pi 0 4*pi -2.5 1]); 如下图:-2.5-2-1.5-1-0.500.51此外,还有带等高线的三维网格曲面函数meshc 和带底座的三维网格曲面函数meshz。
MATLAB实验二:二维图形与三维图形的绘制
.实验报告(201 /201 学年第学期)课程名称实验名称二维图形与三维图形的绘制实验时间年月日实验室指导教师学生姓名学号班级专业实验报告三、实验内容及原理(包括硬件原理图、算法、逻辑框图,关键代码等,可续页)(一)二维图形的绘制1、绘制二维曲线的基本函数:○1plot函数plot函数的基本调用格式为:plot(x,y);其中x和y为长度相同的向量,分别用于存储x坐标和y坐标数据。
例,绘制参数方程曲线。
程序如下:含多个输入参数的plot函数调用格式为:p lot(x1,y1,x2,y2,…,xn,yn);含选项的plot函数调用格式为:plot(x1,y1,选项1,x2,y2,选项2,…,xn,yn,选项n);例,用不同线型和颜色在同一坐标内绘制曲线y=2e-0.5x sin(2πx)及其包络线。
程序如下:○2双纵坐标函数plotyyplotyy函数是MATLAB 5.X新增的函数。
它能把函数值具有不同量纲、不同数量级的两个函数绘制在同一坐标中。
调用格式为:plotyy(x1,y1,x2,y2);其中x1-y1对应一条直线,x2-y2对应另一条曲线。
横坐标的标度相同,纵坐标有两个,左纵坐标用于x1-y1数据对,右纵坐标用于x2-y2数据对。
2、绘制二维图形的其他函数在线性直角坐标系中,其他形式的图形有条形图、阶梯图、杆图和填充图等,所采用的函数分别是:bar(x,y,选项);stairs(x,y,选项);stem(x,y,选项);fill(x1,y1,选项1,x2,y2,选项2,…);例,分别以条形图、填充图、阶梯图和杆图形式绘制曲线y=2e-0.5x。
程序如下:(二)三维图形的绘制1、绘制三维曲线的基本函数plot3函数与plot函数用法十分相似,其调用格式为:plot3(x1,y1,z1,选项1,x2,y2,z2,选项2,…,xn,yn,zn,选项n); 例,绘制空间曲线。
程序如下:2、绘制三维曲面的函数surf函数和mesh函数的调用格式为:surf(x,y,z,e);mesh(x,y,z,e);例,绘制两个直径相等的圆管的相交图形。
matlab的meshgrid函数用法
matlab的meshgrid函数用法一、概述Matlab的meshgrid函数是用于生成网格坐标的函数,它可以将二维或三维数组作为输入,生成对应的网格坐标矩阵。
这个函数在绘制三维图形、进行数值分析等方面有着广泛的应用。
二、基本语法语法格式:meshgrid(x, y) 或 meshgrid(X, Y, Z)其中,x和y可以是任何一维或二维数组,分别表示x轴和y轴的坐标;X和Y可以是任何一维或二维数组,表示x轴和y轴坐标的范围;Z是可选的,如果提供,则生成三维图形的网格坐标矩阵。
三、使用方法1. 生成二维网格坐标矩阵:使用meshgrid函数可以方便地生成二维网格坐标矩阵,用于绘制二维图形。
例如,假设有一个二维数组A,可以使用以下代码生成对应的网格坐标矩阵:A = [1 2 3; 4 5 6; 7 8 9];[X, Y] = meshgrid(A);plot(X, Y);输出结果为:图形:三条线段组成的网格图2. 生成三维网格坐标矩阵:使用meshgrid函数可以方便地生成三维网格坐标矩阵,用于绘制三维图形。
例如,假设有三个二维数组A、B和C,可以使用以下代码生成对应的网格坐标矩阵:A = [1 2; 3 4];B = [5 6; 7 8];C = [9 10; 11 12];[X, Y, Z] = meshgrid(A, B, C);surf(X, Y, Z);输出结果为:图形:三维曲面图四、注意事项1. meshgrid函数生成的网格坐标矩阵是一个矩阵,而不是一个数组。
因此,在使用plot函数绘制图形时,需要将生成的网格坐标矩阵转换为向量或矩阵的形式。
2. 在使用meshgrid函数时,需要注意输入数据的维度和范围。
如果输入的数据不符合要求,可能会导致错误或无效的结果。
3. 在生成三维图形时,如果只使用两个维度生成网格坐标矩阵,可能会出现重复或缺失的情况。
因此,在使用meshgrid函数生成三维图形时,建议同时提供三个维度数据。
MATLAB二维与三维绘图
4. 图形窗口的分割 subplot函数的调用格式为: subplot(m,n,p)
n
p
m
4. 图形窗口的分割 subplot函数的调用格式为: subplot(m,n,p) 例5.7 在一个图形窗口中以子图形式同时绘制 正弦、余弦、正切、余切曲线。 程序如下:
例5.4 用不同标度在同一坐标内绘制曲 线y1=e-0.5xsin(2πx)及曲线y2=1.5e0.1xsin(x)。
5.1.2 绘制图形的辅助操作 1. 图形标注 有关图形标注函数的调用格式为: title(图形名称) xlabel(x轴说明) ylabel(y轴说明) text(x,y,图形说明) legend(图例1,图例2,…)
t=0:pi/50:2*pi; x=8*cos(t); y=4*sqrt(2)*sin(t); z=-4*sqrt(2)*sin(t); plot3(x,y,z,'p'); title('Line in 3-D Space');text(0,0,0,'origin'); xlabel('X'),ylabel('Y'),zlabel('Z');grid;
5.1.3 绘制二维图形的其他函数 1. 其他形式的线性直角坐标图 在线性直角坐标系中,其他形式的图形有条形图 、阶梯图、杆图和填充图等,所采用的函数分 别是: bar(x,y,选项) 与plot参数相同 stairs(x,y,选项) stem(x,y,选项) fill(x1,y1,选项1,x2,y2,选项2,…)
教你用MATLAB快速作一维、二维、三维图
subplot(2,2,3);
• ct=cos(x)./(sin(x)+eps);
plot(x,t);
title('tangent(x)');
• subplot(2,2,1); 分成2×2区域且指定1号为活动区 axis ([0 2*pi -40 40]);
• plot(x,y);
subplot(2,2,4);
• 下述程序段绘制一正方形并以黄色填充:
精选可编辑ppt
20
• x=[0 1 1 0 0]; 正方形顶点坐标向量
• y=[0 0 1 1 0]; • fill(x,y,'y');绘制并以黄色填充正方形图
• 再如:
• x=[0:0.025:2*pi];
• y=sin(3*x);
• fill(x,y,[0.5 0.3 0.4]); 颜色向量
24
2、多条曲线
plot3(x,y,z)
其中x,y,z都是m*n矩阵,其对应的每一列表示一条曲线.
例 画多条曲线观察函数Z=(X+Y).^2.
解 x=-3:0.1:3;y=1:0.1:5; [X,Y]=meshgrid(x,y); Z=(X+Y).^2; plot3(X,Y,Z)
列程序段将绘制条形图形
x=[-2.5:0.25:2.5];
y=exp(-x.*x); bar(x,y); 绘制条形图命令
精选可编辑ppt
19
• 6.3 填充图形
• fill(x,y,’c’)函数用来绘制并填充二维多
边图形,x和y为二维多边形顶点坐标向 量。字符 ’c’ 规定填充颜色,其取值前 已叙述。
-0.6
精选可编辑ppt
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
二维图形和三维图形的创建
1.生成1×10维的随机数向量a,在同一幅图片上分别用红、黄、蓝、绿色绘出其连线图、脉冲图、阶梯图和条形图,并分别标出标题“连
线图”、“脉冲图”、“阶梯图”、“条形图”。
2.绘制向量x=[1 3 0.5 2.5 2]的饼形图,并把3对应的部分分离出来。
3.用hold on命令在同一个窗口绘制曲线y=sin(t),y1=sin(t+0.25) y2=sin(t+0.5),其中t=[0 10]。
t=0:1/100:10;
y=sin(t);
plot(t,y,'-r')
hold on
y1=sin(t+0.25);
plot(t,y1,'--b')
y2=sin(t+0.5);
plot(t,y2,'-y')
4.绘制曲线 x=tcos(3t)
y=tsin2t 其中-π≤t≤π,步长取π/100。
要求:要图形注解、标题、坐标轴标签, 并在曲线上截取一点,将相对应的坐标值文本标注出来(ginput())。
;
5.在三个子图像中,分别绘制三维曲线,三维曲面,三维网格的半径为6,坐标为(6,7,6)的由900个面构成的球面(sphere()),对每个图形标注标题
6.(1)绘一个圆柱螺旋线(形似弹簧)图。
圆柱截面直径为10,高度为5,每圈上升高度为1。
如左图所示。
(2)利用(1)的结果,对程序做少许修改,得到如右图所示图形。
思考题:如果要绘制出如图所示的图形,请先指出这四个图形分别对应哪副图,以及请正确填写下列空格
subplot(6 ,4 ,[2 3 4 6 7 8] );plot(1:10);grid on;
subplot( 6 ,4 ,[10 11 14 15] );plot(peaks);grid on;
subplot(6 ,4 ,[5 9]);plot(membrane);grid on;
subplot(6, 4,[17 18 21 22]);surf(membrane);grid on;。