实验六MATLAB绘图并且附有答案
上机习题6 MATLAB7.0三维绘图
实验六MATLAB7.0三维绘图实验目的：①掌握绘制三维图能形的方法；②掌握图形修饰处理方法；③知道图像处理方法，了解动画制作方法。
实验要求:给出程序和实验结果。
实验内容：一、绘制三维曲线sin()cos()sin()cos()x ty tz t t t=⎧⎪=⎨⎪=⎩二、绘制三维曲面图z=sin(x+sin(y))-x/10。
三、绘制z=x2+y2的三维网线图形；四、绘制三维陀螺锥面；（仅供参考：t1=0:0.1:0.9;t2=1:0.1:2;r=[t1 -t2+2];[x,y,z]=cylinder(r,30);surf(x,y,z);grid）五、在xy平面内选择区域[-8,8]×[-8,8]，利用mesh、meshc、meshz和surf绘制z=六、绘制光照处理后的球面，取三个不同的光照位置进行比较。
（提示：可以利用函数sphere和 light）七、利用peaks产生数据，绘制多峰曲面图。
八．22yxxez--=，当x和y的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图。
九绘制peaks函数的表面图，用colormap函数改变预置的色图，观察色彩的分布情况。
十、用sphere函数产生球表面坐标，绘制不透明网线图、透明网线图、表面图和带剪孔的表面图。
十一、将5.9题中的带剪孔的球形表面图的坐标改变为正方形，以使球面看起来是圆的而不是椭圆的，然后关闭坐标轴的显示。
实验结果：（1）（2）（3）（4）（5）（6）（7）（8）（9）（10）（11）。
MATLAB全部实验及答案解析
实验项目序号1 MATLAB基本操作一、实验目的熟悉MATLAB软件环境，掌握命令窗口的使用。
二、实验基本知识1、熟悉MATLAB的桌面工具和环境2、Start菜单三、主要仪器及耗材计算机四、实验内容和步骤1、MATLAB的启动及操作界面认识说明有哪几种启动的方法；操作界面有哪些窗口；对常用的窗口作些简单的介绍。
➢使用Windows的“开始”菜单。
➢运行MATLAB系统的启动程序matlab.exe。
➢利用桌面上的快捷方式。
➢MATLAB主窗口➢命令窗口(Command Window)➢当前目录窗口(Current Directory)➢工作空间窗口(Workspace)➢命令历史窗口(Command History)2、下列符号中哪些是MATLAB的合法变量名？用给变量赋值的方法在机器上验证你的答案（3vars、global、help、My_exp、sin、X+Y、_input、E-4、AbCd、AB_C_。
）help，My_exp，sin，AbCd，AB_C_；3、命令窗口的简单使用1）简单矩阵的输入1 2 3A = 4 5 67 8 92）求[12+2×(7-4)]÷32的算术运算结果4、有关向量、矩阵或数组的一些运算1）设A=15；B=20；求C=A+B与c=a+b？2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B与A.*B？A*B遵循矩阵运算，A.*B遵循每个对应元素之间进行相乘3）设A=10，B=20；求C=A/B与D=A\B？4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全下标的形式），并将其单下标转换成全下标>> a=[1 -2 3;4 5 -4;5 -6 7];/%定义矩阵>> b = a(find(a<0));%找到a中小于零的元素b =-2-6-4>> c = find(a<0) %找到小于零元素的单下标c =468>> [x,y] = find(a<0) %找到小于零元素的双下标>> [x,y] = ind2sub(size(a),find(a<0)) %单下标转换为双下标5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行A=[1,2;3,4]+i[5,6;7,8]，结果又如何？6）请设计出程序完成下列计算：a=[1 2 3;3 4 2;5 2 3]，求a^2=？,a.^2=？解：>> a=[1 2 3;3 4 2;5 2 3];>> a^2 % 矩阵运算>> a.^2 %数组与矩阵的乘法7）有一段程序设计如下，请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : )解：转化X为列向量8）使用两种方法，创建一稀疏矩阵解：1）>>a = eye(5); %创建5阶单位矩阵>>sparse(a) %将a构造为稀疏矩阵ans =(1,1) 1(2,2) 1(3,3) 1(4,4) 1(5,5) 12）>> data = [1 2 2 2];>> ir = [1 3 2 4];>> jc = [1 2 4 1];>> s=sparse(ir,jc,data,4,4)s =(1,1) 1(4,1) 2(3,2) 2(2,4) 29）写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^AC =3 9 27 >> D = A.^BD =1 32 7295、已知⎪⎭⎫⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。
matlab绘图试题及答案
matlab绘图试题及答案MATLAB绘图试题及答案1. 绘制正弦曲线题目：使用MATLAB绘制函数 \( y = \sin(x) \) 在区间 \([0, 2\pi]\) 上的图像。
答案：```matlabx = 0:0.01:2*pi;y = sin(x);plot(x, y);title('Sin(x) Function');xlabel('x');ylabel('y');```2. 绘制二次函数图像题目：绘制函数 \( y = ax^2 + bx + c \) 在区间 \([-10, 10]\) 上的图像，其中 \( a = 1, b = 2, c = 3 \)。
答案：```matlabx = -10:0.1:10;y = x.^2 + 2*x + 3;plot(x, y);title('Quadratic Function y = x^2 + 2x + 3');xlabel('x');ylabel('y');```3. 绘制柱状图题目：给定一组数据 [10, 20, 30, 40, 50]，使用MATLAB绘制柱状图。
答案：```matlabdata = [10, 20, 30, 40, 50];bar(data);title('Bar Chart');xlabel('Index');ylabel('Value');```4. 绘制散点图题目：给定两组数据 \( x = [1, 2, 3, 4, 5] \) 和 \( y = [2, 4, 6, 8, 10] \)，使用MATLAB绘制对应的散点图。
答案：```matlabx = [1, 2, 3, 4, 5];y = [2, 4, 6, 8, 10];scatter(x, y);title('Scatter Plot');xlabel('x');ylabel('y');```5. 绘制饼图题目：给定一组数据 [10, 20, 30, 40]，使用MATLAB绘制饼图。
Matlab实验及答案-推荐下载
ans =
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置各试时类卷，管调需路控要习试在题验最到；大位对限。设度在备内管进来路行确敷调保设整机过使组程其高1在中正资，常料要工试加况卷强下安看与全22过，22度并22工且22作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
MATLAB实验指导书(附答案)
MATLAB基础实验指导书漳州师范学院物电系2010年10月目录实验一MATLAB环境的熟悉与基本运算 (2)实验二MATLAB数值运算 (8)实验三MATLAB语言的程序设计 (12)实验四MATLAB的图形绘制 (16)实验五采用SIMULINK的系统仿真 (20)实验六MATLAB在电路中的应用 (25)实验七MATLAB在信号与系统中的应用 (27)实验八MATLAB在控制理论中的应用 (29)实验一 MATLAB环境的熟悉与基本运算一、实验目的：1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识：1.熟悉MATLAB环境:MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器文件和搜索路径浏览器。
2.掌握MATLAB常用命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英文字母开头（3）长度不大于31个（4）区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。
MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算4.MATLAB的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式5.MATLAB的基本运算表7 两种运算指令形式和实质内涵的异同表6.MATLAB的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）2、学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
3、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、exerc2、exerc3……），学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。
MATLAB)课后实验答案-精简版.do
实验一1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4)2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.52. 已知：34787,2033657327A B ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A(5) [A,B]和[A([1,3],:);B^2] 3.设有矩阵A 和 B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
ch='ABC123d4e56Fg9';则要求结果是：实验二1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；E=eye(3); R=rand(3,2); O=zeros(2,3);S=diag([1 2]);A=[E R a=[E,R+R*SO S] O ,S^2] A^2=a 2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
matlab课后习题答案（附图）
matlab课后习题答案（附图）习题2.1画出下列常见曲线的图形y （1）⽴⽅抛物线3x命令：syms x y;ezplot('x.^(1/3)')(2)⾼斯曲线y=e^(-X^2);命令：clearsyms x y;ezplot('exp(-x*x)')(3)笛卡尔曲线命令：>> clear>> syms x y;>> a=1;>> ezplot(x^3+y^3-3*a*x*y)(4)蔓叶线命令：>> clear>> syms x y;>> a=1ezplot(y^2-(x^3)/(a-x))(5)摆线:()()tsin-=,=-by1命令：>> clear>> t=0:0.1:2*pi;>> x=t-sin(t);>>y=2*(1-cos(t)); >> plot(x,y)7螺旋线命令：>> clear >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t;>>plot3(x,y,z)(8)阿基⽶德螺线命令：clear>> theta=0:0.1:2*pi;>> rho1=(theta);>> subplot(1,2,1),polar(theta,rho1)(9) 对数螺线命令：cleartheta=0:0.1:2*pi;rho1=exp(theta);subplot(1,2,1),polar(theta,rho1)(12)⼼形线命令：>> clear >> theta=0:0.1:2*pi; >> rho1=1+cos(theta); >> subplot(1,2,1),polar(theta,rho1)练习2.21. 求出下列极限值（1）nnn n3→命令：>>syms n>>limit((n^3+3^n)^(1/n)) ans = 3（2）)121(lim n n n n ++-+∞→命令：>>syms n>>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf) ans = 0（3）x x x 2cot lim 0→命令：syms x ；>> limit(x*cot(2*x),x,0) ans = 1/2 （4）)(coslimcm xx ∞→命令：syms x m ； limit((cos(m/x))^x,x,inf) ans = 1（5）)111(lim 1--→exx x命令：syms x>> limit(1/x-1/(exp(x)-1),x,1) ans =(exp(1)-2)/(exp(1)-1) （6）)(2lim x x xx -+∞>> limit((x^2+x)^(1/2)-x,x,inf)ans = 1/2练习2.41. 求下列不定积分，并⽤diff 验证：（1）+x dxcos 1>>Clear >> syms x y >> y=1/(1+cos(x)); >> f=int(y,x) f =tan(1/2*x) >> y=tan(1/2*x); >> yx=diff(y ,x); >> y1=simple(yx) y1 =1/2+1/2*tan(1/2*x)^2 (2)+exdx1clear syms x yy=1/(1+exp(x));f=int(y,x) f =-log(1+exp(x))+log(exp(x)) syms x yy=-log(1+exp(x))+log(exp(x)); yx=diff(y,x); y1=simple(yx) y1 = 1/(1+exp(x)) (3)dx x x ?sin 2syms x yy=x*sin(x)^2; >> f=int(y,x) f =x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2 clearsyms x y y=x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2; yx=diff(y,x); >> y1=simple(yx) y1 = x*sin(x)^2 (4)xdx ?sec3syms x y y=sec(x)^3;f=int(y,x) f =1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)) clear syms x yy=1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)); yx=diff(y,x); y1=simple(yx) y1 =1/cos(x)^32. 求下列积分的数值解 1）dx x-10clearsyms xy=int(x^(-x),x,0,1) y =int(x^(-x),x = 0 .. 1) vpa(y,10) ans =1.291285997 2）xdx e x cos3202?πclearsyms xy=int(exp(2*x)*cos(x)^3,x, clear syms xy=int((1/(2*pi)^(1/2))*exp(-x^2/2),x,0,1) y =7186705221432913/36028797018963968*erf(1/2*2^(1/2))*2^(1/2)*pi^(1/0,2*pi) y =22/65*exp(pi)^4-22/65vpa(ans,10)(3)dx xe21221-π>> clear >> syms x>> y=int(1/(2*pi)^(1/2)*exp(-x^2/2),0,1); >> vpa(y,14) ans =.341344746068552（4）>> clear >> syms x>> y=int(x*log(x^4)*asin(1/x^2),1,3); Warning: Explicit integral could not be found. > In sym.int at 58 >> vpa(y,14) ans = 2.45977212823752（5） >> clear >> syms x1判断下列级数的收敛性，若收敛，求出其收敛值。
MatLab练习题-绘图参考解答
MatLab 练习题-绘图参考解答学号：班级：姓名：注意: 1,请按要求作题,直接输入答案,并保存文件.2,仔细填写“学生信息”.3,把此word文档的文件名更名为你的“学号第三次”.在下课前15分钟内上传到http://10.1.9.91.例：“MatLab练习册.doc”——-> “200504101第三次.doc” ) 1,在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注“y=sin(x)”, “y=cos(x)” ,x轴,y轴,标题为“正弦余弦函数图象”.答案：函数文件如下：t=0:0.02:4*pi;plot(t,sin(t),'m*',t,cos(t),'rO')text(pi/3,sin(pi/3),'y=sin(x)')text(pi/3,cos(pi/3),'y=cos(x)')title('正线函数和余弦函数')xlabel('x轴')ylabel(' y轴')结果显示：2, 画出椭圆221 169x y+=。
答案：函数文件如下：%画椭圆tt=0:0.02:2*pi;xx=4*cos(tt);yy=3*sin(tt); axis equalplot(xx,yy,'m*')title('椭圆x^2/16+y^2/9=1')答案2：函数文件：xx=-4:0.02:4;yy=3*sqrt(1-xx.^2/16);plot(xx,yy,'m*')axis equalhold onplot(xx,-yy,'r+')title('椭圆x^2/16+y^2/9=1')3,绘制螺旋线.(12cos)(12sin)3x uy uz u=+⎧⎪=+⎨⎪=⎩[0,30]u∈答案：u=0:0.01:30;x=(1+2*cos(u));y=1+2*sin(u);z=3*u; plot3(x,y,z,'r*')4,绘制四叶玫瑰线sin2ρθ=。
matlab绘图参考答案
matlab绘图参考答案matlab绘图参考答案Matlab是一种强大的数学软件，被广泛应用于科学计算、数据分析和绘图等领域。
在进行数据可视化时，Matlab提供了丰富的绘图函数和工具，可以轻松地创建各种类型的图表。
本文将介绍一些常见的Matlab绘图函数和技巧，帮助读者更好地掌握绘图的基本方法和应用。
一、折线图折线图是最常见的一种图表类型，用于显示数据随时间或其他变量的变化趋势。
在Matlab中，可以使用plot函数来绘制折线图。
例如，下面的代码演示了如何绘制一条简单的折线图：```matlabx = 1:10;y = [1 2 3 4 5 4 3 2 1 2];plot(x, y);```上述代码中，x是自变量，表示时间或其他变量的取值范围，y是因变量，表示对应时间或变量的取值。
plot函数将x和y作为输入参数，绘制出对应的折线图。
二、散点图散点图用于显示两个变量之间的关系，通常用于探索变量之间的相关性。
在Matlab中，可以使用scatter函数来绘制散点图。
例如，下面的代码演示了如何绘制一组随机生成的散点图：x = rand(1, 100);y = rand(1, 100);scatter(x, y);```上述代码中，x和y分别是两个变量的取值，scatter函数将x和y作为输入参数，绘制出对应的散点图。
三、柱状图柱状图用于比较不同类别或组之间的数据差异。
在Matlab中，可以使用bar函数来绘制柱状图。
例如，下面的代码演示了如何绘制一组随机生成的柱状图：```matlabx = 1:5;y = rand(1, 5);bar(x, y);```上述代码中，x是表示不同类别或组的标签，y是对应类别或组的数据。
bar函数将x和y作为输入参数，绘制出对应的柱状图。
四、饼图饼图用于显示不同类别或组占总体的比例关系。
在Matlab中，可以使用pie函数来绘制饼图。
例如，下面的代码演示了如何绘制一组随机生成的饼图：```matlabx = rand(1, 5);```上述代码中，x是表示不同类别或组的比例，pie函数将x作为输入参数，绘制出对应的饼图。
Matlab实验指导书(含答案)汇总
实验Matlab 操作环境熟悉、实验目的1．初步了解Matlab 操作环境。
2．学习使用图形函数计算器命令funtool 及其环境。
二、实验内容熟悉Matlab 操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format 命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who 和whos 命令查看内存变量信息；学会使用图形函数计算器funtool ，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1) y=sin(x);(2) y=(1+x)A3*(2-x); 求下列函数的符号积分(1) y=cos(x);(2) y=1/(1+x^2);(3) y=1/sqrt(1-x^2);(4) y=(x-1)/(x+1)/(x+2);求反函数(1) y=(x-1)/(2*x+3);(2) y=exp(x);(3) y=log(x+sqrt(1+x^2));代数式的化简(1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2) sin (x)A2+cos(x)A2;(3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。
从y=xA2通过参数的选择去观察下列函数的图形变化(1) y1=(x+1)A2⑵ y2=(x+2)A2(3) y3=2*x^2⑷ y4=x^2+2(5) y5=x^4⑹ y6=x^2/23．两个函数之间的操作求和(1) sin(x)+cos(x)(2) 1+x+x^2+x^3+x^4+x^5乘积(1) exp(-x)*sin(x)(2) sin(x)*x商(1) sin(x)/cos(x);⑵ x/(1+x^2);(3) 1/(x-1)/(x-2);求复合函数(1) y=exp(u) u=sin(x)(2) y=sqrt(u) u=1+exp(xA2)(3) y=sin(u) u=asin(x)(4) y=sinh(u) u=-x三、设计提示1．初次接触Matlab 应该注意函数表达式的文本式描述。
matlab实验指导答案详解（非常详细正确）
matlab实验指导答案详解（⾮常详细正确）实验⼀ MATLAB ⼯作环境熟悉及简单命令的执⾏⼀、实验⽬的：熟悉MATLAB 的⼯作环境，学会使⽤MATLAB 进⾏⼀些简单的运算。
⼆、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌⾯（Desktop ），包括菜单（Menu ）、⼯具条（Toolbar ）、命令窗⼝(Command Window)、历史命令窗⼝、⼯作空间(Workspace)等；完成⼀些基本的矩阵操作；学习使⽤在线帮助系统。
三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌⾯。
2、在命令窗⼝执⾏命令完成以下运算，观察workspace 的变化，记录运算结果。
（1）（365-52?2-70）÷3 >>(365-52*2-70)/3 ans = 63.6667（2）>>area=pi*2.5^2 area = 19.6350（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -= >>x=3 >>y=4>>z = x ^2 * y ^3 / (x - y) ^2 z = 576（4）将下⾯的矩阵赋值给变量m1,在workspace 中察看m1在内存中占⽤的字节数。
m1=11514412679810115133216 执⾏以下命令>>m1 =[16 2 3 13 ; 5 11 10 8 ; 9 7 6 12 ; 4 14 15 1 ] >>m1( 2 , 3 ) ans = 10 >>m1( 11 ) ans = 6>>m1( : , 3 ) ans =3 10 6 15>>m1( 2 : 3 , 1 : 3 ) ans =5 11 10 9 7 6>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) ans = 34（5）执⾏命令>>help abs查看函数abs 的⽤法及⽤途，计算abs( 3 + 4i ) （6）执⾏命令>>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y)（6）运⾏MATLAB 的演⽰程序，>>demo ，以便对MATLAB 有⼀个总体了解。
MATLAB数学实验练习题附答案PPT课件
%2、求100~999之间的水仙花数
clear all；
clc;
for n=100:999
n1=floor(n/100); %取出百位数字n1
ቤተ መጻሕፍቲ ባይዱ
n2=mod(floor(n/10),10); %取出十位数字n2
2
n3=mod(n,10) ; %取出个位数字n3
if n1^3+n2^3+n3^3==n
fprintf('%d是“水仙花数”\n', n)
让我们共同进步
2021/3/9
5
===============================================================
%5.作xoy面上的曲线(x-a).^2+y.^2=r.^2绕z轴旋转所得到的图形(0<r<a),取r=1,a=2.
clear all,clc;
x=1:0.1:3;
3
[x,y,z]=cylinder(sqrt(1-(x-2).^2));
surf(x,y,z)
2021/3/9
%6. 作x2/2-y2/3=2z图形. clear,clf,clc; [X,Y,Z]=meshgrid(-1:0.1:1); V=X.^2./2+Y.^2./2-2*Z; fv=isosurface(X,Y,Z,V,0);%显示V=0等值面 p=patch(fv); set(p,'FaceColor','blue','EdgeColor','none');%设置图形的相关属性 view(3) hold on camlight
end
202e1n/3d/9
MATLAB实验报告（绘图）
MATLAB实验报告（绘图）MATLAB实验报告——MATLAB绘图学号：学院：班级：姓名：1、绘制y = e x 3sin 3x (x ∈ 0,4π ) 的图像，要求用蓝色的星号画图；并且画出其包络线y = ±e x 3 的图像，用红色的点划线画图. MATLAB Code: x = 0:pi/50:4*pi;y = exp(x/3).*sin(3*x);z(1,:) = exp(x/3);z(2,:) = -exp(x/3);plot(x,y,'b*',x,z(1,:),'r-.',x,z(2,:),'r-.');2、用fplot 和ezplot 命令绘出函数y = e2t 3sin ?(1+2t ) 在区间[1,10]上的图像.MATLAB Code:subplot(1,2,1),fplot(@(t)exp(-2*t/3).*sin(1+2*t),[1,10]); title('fplot');xlabel('t');subplot(1,2,2),ezplot('exp(-2*t/3).*sin(1+2*t)',[1,10]); title('ezplot');3、在同一图形窗口画三个字图，要求使用指令gtext,axis,legend,title,xlabel,ylabel:(1)y=xcosx,x∈(?π,π)(2) y=xtan1xsin x3,x∈(π,4π)(3) y=e 1sin x,x∈[1,8]MATLAB Code:subplot(2,2,1),x1 = -pi:pi/50:pi;y1 = x1.*cos(x1);plot(x1,y1,'r'),axisequal,legend('x*cosx'),title('f1'),xlabel('x'),ylabel('y'); subplot(2,2,2),x2 = pi:pi/50:4*pi;y2 = x2.*tan(1./x2).*sin(x2.^3);plot(x2,y2,'b'),axisequal,legend('x*tan(1/x)*sin(x^3)'),title('f2'),xlabel('x') ,ylabel ('y');subplot(2,2,3:4),x3 = 1:0.01:8;y3 = exp(1./x3).*sin(x3);plot(x3,y3,'g'),legend('exp(1/x)*sinx'),title('f3'),xlabel( '时间'),ylabel('位移');gtext('x*cosx');4、使用合适的单轴对数坐标函数绘制函数y=e x2的图像（其中1≤x≤10）MATLAB Code:x = 1:0.01:10;y = exp(x.^2);semilogy(x,y,'b-.');xlabel('x'),ylabel('y');5、绘制圆锥螺线的图像并添加各种标注，圆锥螺线的参数方程为：x =tcos πt y =tsin π6t z =2t(0≤t ≤20π) MATLAB Code:t = 0:pi/50:20*pi;x = t.*cos(pi/6*t);y = t.*sin(pi/6*t);z = 2*t;plot3(x,y,z,'b'),grid on;title('圆锥螺线');xlabel('x = tcost');ylabel('y = tsint');zlabel('z = 2t');6、在同一个图形窗口画半径为1的球面、柱面x2+y2=1以及极sin4t,t∈[0,2π].坐标图形ρ=12MATLAB Code:subplot(1,2,1),sphere(100);hold on;cylinder;hold on;t = 0:pi/50:2*pi;r = 0.5*sin(4*t);subplot(1,2,2);polar(t,r);title('r = 0.5*sin4t');7、用mesh与surf命令绘制三维曲面z=x2+3y2的图像，并使用不同的着色效果及光照效果.MATLAB Code:t = -3:0.1:3;[x,y] = meshgrid(t);z = x.^2+3*y.^2;subplot(1,2,1),mesh(x,y,z),title('网格z = x^2+3y^2'),shading flat;light('position',[10,10,2]);subplot(1,2,2),surf(x,y,z),title('表面z = x^2+3y^2');shading interp; light('position',[5,-4,-2]);8、绘制由函数x 29+y216+z24=1形成的立体图，并通过改变观测点获得该图形在各个坐标平面的投影.MATLAB Code:t = 0:pi/50:2*pi;[x,y] = meshgrid(t,t);X = 3*sin(y).*cos(x);Y = 4*sin(y).*sin(x);Z = 2*cos(y); subplot(2,2,1);mesh(X,Y,Z);title('x^2/9+y^2/16+z^2/4=1');subplot(2,2,2);mesh(X,Y,Z);view(90,0);title('x^2/9+y^2/16+z^2/4=1在yoz面投影'); subplot(2,2,3);mesh(X,Y,Z);view(0,0);title('x^2/9+y^2/16+z^2/4=1在xoz面投影'); subplot(2,2,4);mesh(X,Y,Z);view(0,90);title('x^2/9+y^2/16+z^2/4=1在xoy面投影');9、画三维曲面z=5?x2?y2?2≤x,y≤2与平面z=3的交线. MATLAB Code:t = -2:0.01:2;[x,y] = meshgrid(t);z1=(5-x.^2-y.^2);z2 = 3*ones(length(t));r0 = (abs(z1-z2)<=0.03);zz = r0.*z2;yy = r0.*y;xx = r0.*x;plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'b.');xlabel('x'),ylabel('y'),zlabel('z');title('z = 5-x^2-y^2 与 z = 3的交线');10(附加)、利用迭代复函数z n+1=z n2+c,c∈C构造Mandelbrot集. Mandelbrot集定义为：M: =c∈C| |c|≤2,|c2+c|≤2,|(c2+c)2+c|≤2…….MATLAB Code:%最大迭代数网格细度%本次作图取maxloopNum = 100,Eps = 0.05function drawMandelbrot(maxloopNum,Eps)t = -2:Eps:2;[x,y] = meshgrid(t);z = x+i*y;isMandel = zeros(length(t));for j = 1:length(z)for k = 1:length(z)iter = 0;current = z(j,k);temp = current;while abs(temp)<=2 && iter<maxloopnum< p="">temp = temp^2+current;iter = iter + 1;endif iter == maxloopNumisMandel(j,k) = 1;endendendplot(z.*isMandel,'.');axis equal;title('MandelBrot set');legend('loop = 100,eps = 0.05','Location','northwest'); xlabel('x'),ylabel('iy');end</maxloopnum<>。
MATLAB实验报告实验六
MATLAB实验报告实验六实验六1.x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x.^2));plot(x,y)2.x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'r:',x,y2,'g-',x,y3,'b--')x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(1,3,1);plot(x,y1,'r:');title('y1=x.^2'); subplot(1,3,2);plot(x,y2,'g-');title('y2=cos(2*x)'); subplot(1,3,3);plot(x,y1,'b--');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2);bar(x,y1,'r');title('y1=x^2');subplot(2,2,3);bar(x,y2,'g');title('y2=cos(2*x)');subplot(2,2,4);bar(x,y3,'b');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2); stairs(x,y1,'r:');title('y1=x^2');subplot(2,2,3);stairs(x,y2,'g-');title('y2=cos(2*x)');subplot(2,2,4);stairs(x,y3,'b--');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2); stem(x,y1,'r:');title('y1=x^2');subplot(2,2,3);stem(x,y2,'g-');title('y2=cos(2*x)');subplot(2,2,4);stem(x,y3,'b--');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2);fill(x,y1,'r:');title('y1=x^2');subplot(2,2,3);fill(x,y2,'g-');title('y2=cos(2*x)');subplot(2,2,4);fill(x,y3,'b--');title('y3=y1.*y2');3.x=-5:0.01:5;y=(x+sqrt(pi))./(exp(2)).*(x<=0)+0.5*log(x+sqrt(1+x.*x)).*(x>0); plot(x,y4.theta=0:pi/100:2*pi; a=input('请输⼊a:');b=input('请输⼊b:');n=input('请输⼊n:'); rho=a*sin(b+n*theta); polar(theta,rho,'g')a=b=n=2a=b=2,n=4a=n=2,b=4a=4,b=n=2采⽤控制变量法的办法，固定两个参数，变动第三个参数观察输出图象的变化。
实验06_高层绘图操作(第5章)
x=linspace(-5,5,21);
y=linspace(0,10,31);
[X,Y]=meshgrid(x,y);
Z=cos(X).*cos(Y).*exp(-sqrt(X.^2+Y.^2)/4);
subplot(2,1,1);surf(X,Y,Z);
subplot(2,1,2);contour3(X,Y,Z);
subplot(2,2,4);fill(x,y3,'g');%填充图
图形:
3. 绘制分段函数的曲线
已知
在-5≤x≤5区间绘制函数曲线。（注意：本曲线不连续！）
程序：
fplot('(x+sqrt(pi))/exp(2)',[-5,0]);
holdon;
fplot('log(x+sqrt(1+x^2))/2',[0,5]);
x
\0
∅
\rceil
ù
\surd
√
\midபைடு நூலகம்
附参考答案：
实验06高层绘图操作
（第5章MATLAB绘图）
一、实验目的
1. 掌握绘制二维图形的常用函数。
2. 掌握绘制三维图形的常用函数。
3. 掌握绘制图形的辅助操作。
二、实验内容
1. 绘制函数的曲线
设 ，在x=0~2π区间取101点，绘制函数的曲线。
程序：
x=linspace(0,2*pi,101);
例
%一个单位圆，一个复数参数
t=0:0.01:2*pi;
x=exp(i*t);%cos(t)+i*sin(t)
plot(x); axisequal;
MATLAB实验报告绘图
68 54 35;
45 25 12;
48 68 45;
68 54 69];
x=sum(t);
h=pie(x);
textobjs=findobj(h,'type','text');
str1=get(textobjs,{'string'});
val1=get(textobjs,{'extent'});
输出图像:
3.画出所表示的三维曲面（图1.3-4）。的取值范围是。
输入程序: x=-8:0.5:8;
y=x';
X=ones(size(y))*x;
Y=y*ones(size(x));
R=sqrt(X.^2+Y.^2)+eps;
Z=sin(R)./R;
surf(X,Y,Z);
colormap(cool)
实验名称
MATLAB绘图
1、实验目的和要求:
2、了解MATLAB的绘图指令。
3、掌握绘图函数的用法、简单的图形标注,简单颜色设定。
4、掌握MATLAB二维曲线绘图、三维曲线绘图,以及一些特殊图形的绘制。
5、掌握MATLAB二维曲线绘图、三维曲线绘图，以及一些特殊图形的绘制。
实验内容和步骤:
1.作多条曲线
plot3(y1,y2,t);
title('helix'),text(0,0,0,'origin');
xlable('sin(t)'),ylable('cos(t)'),zlable('t');
gridon;
matlab实验 三维画图
实验（六）项目名称：三维绘图一、实验目的：熟悉MATLAB中几个常用的绘图命令，掌握集中常见三维图形的画法。
二、实验原理三维绘图命令：Plot（X,Y,Z）//曲线；mesh(X,Y,Z)//网状；surf(X,Y,Z)//表面；contour(X,Y,Z)//等高线。
三、实验环境1.硬件：PC机2. 软件：Windows操作系统、matlab2015四、实验内容、步骤以及结果4.1.1实验要求：用plot函数画出的三维曲线。
4.1.2实验步骤（1）启动matlab，新建一个M文件；（2）输入程序，如图1；（3）保存文件；（4）编译源程序，观察屏幕上显示的编译信息，修改出现的错误，直到编译成功；图1：plot函数画三维曲线4.1.3运行结果如下：图2：三维曲线4.2.1实验要求：用mesh函数画出的三维网状图。
4.2.2实验步骤（5）启动matlab，新建一个M文件；（6）输入程序，如图3；（7）保存文件；（8）编译源程序，查看运行结果，如图4。
图3：mesh函数画三维网状图图4：运行结果4.3.1实验要求：用surf函数画出的三维表面图。
4.3.2实验步骤（9）启动matlab，新建一个M文件；（10）输入程序，如图5；（11）保存文件；（12）编译源程序，查看运行结果，如图6。
图5：surf函数画三维表面图图6：运行结果4.3.1实验要求：用contour函数画出的等高线图。
4.3.2实验步骤（13）启动matlab，新建一个M文件；（14）输入程序，如图7；（15）保存文件；（16）编译源程序，查看运行结果，如图8。
图7：contour函数画等高线图8：运行结果五、实验总结MATLAB具有强大的图形功能，能够将它们直观的表现出来，解决很多的问题。
MATLAB实验六（精品）
第二节 逐点比较法逐点比较法的基本原理是，在刀具按要求轨迹运动加工零件轮廓的过程中，不断比较刀具与被加工零件轮廓之间的相对位置，并根据比较结果决定下一步的进给方向，使刀具沿着坐标轴向减小偏差的方向进给，且只有一个方向的进给。
也就是说，逐点比较法每一步均要比较加工点瞬时坐标与规定零件轮廓之间的距离，依此决定下一步的走向，如果加工点走到轮廓外面去了，则下一步要朝着轮廓内部走；如果加工点处在轮廓的内部, 则下一步要向轮廓外面走，以缩小偏差，周而复始，直至全部结束，从而获得一个非常接近于数控加工程序规定轮廓的刀具中心轨迹。
逐点比较法既可实现直线插补，也可实现圆弧插补。
其特点是运算简单直观，插补过程的最大误差不超过一个脉冲当量，输出脉冲均匀，而且输出脉冲速度变化小，调节方便，但不易实现两坐标以上的联动插补。
因此，在两坐标数控机床中应用较为普遍。
一般来讲，逐点比较法插补过程每一步都要经过如图3-1所示的四个工作节拍：（1）偏差判别 判别刀具当前位置相对于给定轮廓的偏差情况，即通过偏差值符号确定加工点处在理想轮廓的哪一侧，并以此决定刀具进给方向。
（2）坐标进给 根据偏差判别结果，控制相应坐标轴进给一步，使加工点向理想轮廓靠拢，从而减小其间的偏差。
（3）偏差计算 刀具进给一步后，针对新的加工点计算出能反映其偏离理想轮廓的新偏差，为下一步偏差判别提供依据。
（4）终点判别 每进给一步后都要判别刀具是否达到被加工零件轮廓的终点，若到达了则结束插补，否则继续重复上述四个节拍的工作，直至终点为止。
一、逐点比较法I 象限直线插补（一）基本原理设第一象限直线OE ，起点为坐标原点O(0,0)，终点为E （X e ，Y e ），另有一个动点为N （X i ，Y i ），如图3-2所示。
其中，各个坐标值均是以脉冲当量为单位的整数，以便于后面的推导与讲解，并且在脉冲增量式插补算法中都是这样约定的。
ee i i X Y X Y = （3-1a ） 即 X e Y i —X i Y e =0 （3-1b ） 当动点N 处于直线OE 的下方N ′处时，直线N O '的斜率小于直线OE 的斜率，从而有ii X Y ＜e e X Y （3-2a ）即 X e Y i —X i Y e ＜0 （3-2b ） 当动点N 处于直线OE 的上方N ″处时，直线N O ''的斜率大于直线OE 的斜率，从而有ee i i X Y X Y ＞ （3-3a ） 即 X e Y i —X i Y e ＞０ （3-3b ） 由上述关系可以看出，表达式（X e Y i —X i Y e ）的符号就能反映出动点N 相对直线OE 的偏离情况，为此取偏差函数F 为F =X e Y i —X i Y e （3-4）根据上述过程可以概括出如下关系：当F =0时，动点N （X i ，Y i ）正好处在直线OE 上；当F ＞0时，动点N （X i ，Y i ）落在直线OE 上方区域；当F ＜0时，动点N （X i ，Y i ）落在直线OE 下方区域。
第5章--MATLAB绘图-习题答案
第5章 MATLAB绘图习题5一、选择题1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。
DA．12 B．7 C．4 D．32．下列程序的运行结果是（）。
Ax=0:pi/100:2*pi;for n=1:2:10plot(n*sin(x),n*cos(x))hold onendaxis squareA．5个同心圆B．5根平行线C．一根正弦曲线和一根余弦曲线D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。
CA．{\alpha}+{\beta} B．{\α}+{\β} C．α+βD．\α+\β4．subplot(2,2,3)是指（）的子图。
AA．两行两列的左下图B．两行两列的右下图C．两行两列的左上图D．两行两列的右上图5．要使函数y=2e x的曲线绘制成直线，应采用的绘图函数是（）。
CA．polar B．semilogx C．semilogy D．loglog6．下列程序的运行结果是（）。
B[x,y]=meshgrid(1:5);surf(x,y,5*ones(size(x)));A．z=x+y平面B．与xy平面平行的平面C．与xy平面垂直的平面D．z=5x平面7．下列函数中不能用于隐函数绘图的是（）。
DA．ezmesh B．ezsurf C．ezplot D．plot38．下列程序运行后，看到的图形（）。
Ct=0:pi/20:2*pi;[x,y]=meshgrid(-8:0.5:8);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);授课：XXXsurf(x,y,z)view(0,90);axis equalA．像墨西哥帽子B．是空心的圆C．边界是正方形D．是实心的圆9．下列程序运行后得到的图形是（）。
A[x,y]=meshgrid(-2:2);z=x+y;i=find(abs(x)<1 & abs(y)<1);z(i)=NaN;surf(x,y,z);shading interpA．在一个正方形的正中心挖掉了一个小的正方形B．在一个正方形的正中心挖掉了一个小的长方形C．在一个正方形的上端挖掉了一个小的正方形D．在一个正方形的下端挖掉了一个小的正方形10．在使用MA TLAB“绘图”选项卡中的命令按钮绘图之前，需要（）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验六、MA TLAB绘图
一、实验目的
1、掌握绘制MA TLAB二维、三维和特殊图形的常用函数
2、掌握绘制图形的辅助操作
二、实验内容
1、二维图形的绘制
1)设计程序，在一个图形窗体中画出函数y=e x和y=logx的
特性曲线。
其中y=e x中x取值为（-2:0.1:2），y=logx中x
取值为（0.1:0.1:5）。
要求：要有图例、标题（'二维图'）、坐标轴标签（'x轴数
据'和'y轴数据'）；每条曲线的clm自己设置。
x1=-2:0.1:2;
y1=exp(x1);
x2=0.1:0.1:5;
y2=log(x2);
plot(x1,y1,'r',x2,y2,'g');
title('二维图');
legend('y=exp(x)','y=logx');
xlabel('X轴数据');ylabel('Y轴数据');
grid on;
2)绘制tan图形，其中x的取值为（0:pi/100:pi/2）。
要求：要有图例、标题、坐标轴标签（'x轴数据'和'y轴数
据'）；曲线的clm自己设置。
绘制出的图形应如下右图所
示。
x=0:pi/100:pi/2;
y=tan(x);
plot(x,y,'b+-');
axis([0,pi/2,0,5]);
legend('y=tan x');
title('正切函数');
xlabel('X轴数据');ylabel('Y轴数据');
3）绘制曲线
x=tcos(3t)
y=tsin2t 其中-π≤t≤π，步长取π/100。
要求：要有图例、标题、坐标轴标签；每条曲线的clm 自己设置。
t=-pi:pi/100:pi;
x=t.*cos(3*t);
y=t.*sin(t).*sin(t);
plot(x, y);
title(date);
legend(strvcat('x=tcos(3t)','y=tsin2t'));
xlabel('T轴数据');ylabel('X,Y轴数据');
三维图形的绘制
分别用线框图和曲面图表现函数z=cos(x)sin(y)/y，其中x
的取值为[-1.5pi,1.5pi]，y=x，要求：要有标题、坐标轴标
签
x=-1.5*pi:1.5*pi;
y=x;
[X,Y]=meshgrid(x,y);
Z=cos(X).*sin(Y)./Y; %使用大写了
subplot(1,2,1);
surf(Z) %曲面图
subplot(1,2,2)
mesh(Z) %线框图
2、特殊图形的绘制
自己举例使用子图的方式绘制出条状图、面积图、饼图、火柴杆图四种图形式。
a=[10 2 3 5;5 8 10 3;9 7 6 1 ;3 5 7 2 ;4 7 5 3 ];
subplot(2,2,1);bar(a);
subplot(2,2,2);area(a);
subplot(2,2,3);pie(a);
subplot(2,2,4);stem(a);
3、思考题：要绘制出如图所示的图形，请正确填写下列空格subplot( 6 , 4 , [ 2，3，6，7，] );plot(1:10);grid on; subplot(6,4 ,[10,11,12,14,15,16] , , );plot(peaks) ;grid on;
subplot(6, 4 , [ 5，9] , );plot(membrane);grid
on;
subplot( 6 , 4 , [17 18 19 21 22
23 ] );surf(membrane);grid on;
4、
其中：plot(membrane)plot(peaks)
surf(membrane)。