电力-授课顾问学院-matlab

程序学习前篇学习好一门程序语言的建议：掌握基本的语法结构，熟悉软件使用环境；养成养好的编程习惯，如文件的命名，函数的命名，程序的注释；培养良好的“逻辑结构”；数学多动手，不怕出问题，重点学会寻找问题的方法Matlab技术与应用2012-2013学年课程安排Matlab的基本知识基于Matlab语言进行实例开发Matlab技术与应用2012-2013学年第一讲Matlab概述Matlab的由来NewMexico大学计算机系，系主任CleverMoler博士，1980年Matlab：MatrixLaboratory，Fortran1984年，MathWorks公司，CMatlab技术与应用2012-2013学年第一讲Matlab概述Matlab的特点“智能化”的语言，一种直译式的高级语言，语法宽松，设计自由大；运算符、函数库丰富；图像功能强大；功能全面，应用广泛，现已成为大学工科类专业的必修课。Matlab技术与应用2012-2013学年第一讲Matlab概述Matlab的功能数据运算：矩阵为基本单位，丰富的数值计算函数，高质量的数值计算功能符号运算：解析解，结合符号计算语言Maple绘图功能：二维、三维图形语言体系：程序结构控制、函数调用、数据结构、输入输出等工具箱：基本和可选工具箱Matlab技术与应用2012-2013学年第一讲Matlab概述Matlab工具箱Matlab技术与应用2012-2013学年MATLAB主工具箱符号数学工具箱SIMULINK仿真工具箱控制系统工具箱信号处理工具箱图象处理工具箱通讯工具箱系统辨识工具箱神经元网络工具箱金融工具箱Matlab在工程中的应用Matlab技术与应用2012-2013学年应用系统的开发流程：算法设计与验证算法在硬件系统的实现Matlab是首选汇编语言，C语言，C++Matlab也可与C,C++混合编程Matlab在工程中的应用Matlab技术与应用2012-2013学年Matlab在工程中的应用Matlab技术与应用2012-2013学年第一讲Matlab概述Matlab的开放性与兼容性Matlab技术与应用2012-2013学年Matlab中的所有函数均对用户开放，用户可以根据自己的需要对函数进行查看、更改；Matlab语言可以与其它语言兼容，如C语言，C++，跨平台兼容。(此处用Amigo平台举例说明)。Matlab安装Matlab技术与应用2012-2013学年MatlabR2008a，Matlab7.6.0第一步：第二步：Matlab使用Matlab技术与应用2012-2013学年Matlab使用Matlab技术与应用2012-2013学年Matlab使用Matlab技术与应用2012-2013学年Matlab使用Matlab技术与应用2012-2013学年Matlab使用Matlab技术与应用2012-2013学年Matlab使用Matlab技术与应用2012-2013学年Matlab使用Matlab技术与应用2012-2013学年Matlab使用Matlab技术与应用2012-2013学年Matlab编程实例Matlab技术与应用2012-2013学年两种编程方式举一个加减法的例子说明。在调试过程中，介绍whos,who,help,clc,clear等命令符。示范matlab的安装与打开。第二讲Matlab数据及其运算-2.1数据的特点Matlab技术与应用2012-2013学年Matlab的基本数据对象：矩阵向量：一行或一列的矩阵单个数据标量：单个元素的矩阵Matlab技术与应用2012-2013学年Matlab的数据类型：双精度，Double，64位，8个字节

单精度，Single，32位，4个字节带符号整数，int无符号整数，uint字符数据，Char逻辑数据：没有专门的逻辑数据，‘1’为‘真’，‘0’为‘假’第二讲Matlab数据及其运算-2.1数据的特点Matlab技术与应用2012-2013学年矩阵中元素的数据类型应该一致结构体（Structure）和单元（Cell）数据类型第二讲Matlab数据及其运算-2.1数据的特点Matlab技术与应用2012-2013学年变量命名变量-内存单元-命名以访问变量以字母开头，后接字母、数字、下划线，最多63字符Speed_12,Var1合法字符12Speed,_Var1非法字符区分大小写：Var与var不同Matlab中提供的函数名和命令都是小写clc,clear,who,whos第二讲Matlab数据及其运算-2.2变量及其操作Matlab技术与应用2012-2013学年预定义变量ans计算结果的默认赋值变量pi圆周率的近似值i,j虚数单位inf,Inf无穷大1/0NaN,nan非数，0/0,inf/inf函数命名时，应尽量避开第二讲Matlab数据及其运算-2.2变量及其操作Matlab技术与应用2012-2013学年数据的保存：方法一save——

将工作空间中所有的变量存到matlab.mat文件中第二讲Matlab数据及其运算-2.2变量及其操作Matlab技术与应用2012-2013学年数据的保存方法二save文件名存储路径：默认/指定saved:\data2aascii省略时以二进制格式存储把所有数据存到指定的文件中第二讲Matlab数据及其运算-2.2变量及其操作Matlab技术与应用2012-2013学年数据的保存方法三save文件名变量名ascii存储路径：默认/指定saved:\data2aascii省略时以二进制格式存储文件名可以不带扩展名，默认‘.mat’格式load文件名变量名ascii第二讲Matlab数据及其运算-2.2变量及其操作Matlab技术与应用2012-2013学年数据的读入loadload文件名load文件名变量名ascii第二讲Matlab数据及其运算-2.2变量及其操作Matlab技术与应用2012-2013学年矩阵的建立A=[123;456],直接建立法利用M文件建立大矩阵C=[A;B]C=[AB]矩阵的修改A(,)=来修改第二讲Matlab数据及其运算-2.3矩阵的表示

Matlab技术与应用2012-2013学年冒号表达式e1:e2:e3,e1为初始值，e2为步长（省略时步长为1），e3为终止值矩阵的拆分A(:,j),A(i,:),A(i,j)A(:,K:K+M),A(K:K+M,:)矩阵中元素的删除A(:,[12])=[]第二讲Matlab数据及其运算-2.3矩阵的表示

Matlab技术与应用2012-2013学年矩阵加减运算A+B,两者的大小要相同；矩阵也可以和标量加减,标量与矩阵的所有元素分别进行加减操作。举例第二讲Matlab数据及其运算-2.3数据的运算

Matlab技术与应用2012-2013学年矩阵乘法C=A*B，A为m*n矩阵，B为n*p矩阵C为m*p矩阵矩阵也可以和标量相乘,标量与矩阵的所有元素分别进行乘法操作。第二讲Matlab数据及其运算-2.3数据的运算

Matlab技术与应用2012-2013学年矩阵除法\和/表示左除和右除标量运算：两者一样1\2,2/1矩阵不同A=[1,2,3;4,2,6;7,4,9];B=[4,3,2;7,5,1;12,7,92]A\B不等于B/A第二讲Matlab数据及其运算-2.3数据的运算

Matlab技术与应用2012-2013学年矩阵乘方A^x每个元素的乘方矩阵点运算点运算表示矩阵的元素的进行相关的运算两者的维数要相同A=[12;34],B=[11;-1-1]A.*B不等于A*B第二讲Matlab数据及其运算-2.3数据的运算

Matlab技术与应用2012-2013学年矩阵的其他运算inv——矩阵求逆det——行列式的值eig——矩阵的特征值diag——对角矩阵’——矩阵转置sqrt——矩阵开方第二讲Matlab数据及其运算-2.3数据的运算

Matlab技术与应用2012-2013学年Name='ShandongUniversity‘Name(1:3)第二讲Matlab数据及其运算-2.4字符串

Matlab技术与应用2012-2013学年结构矩阵的建立a(1).x1='Tom';a(1).x2='[001,99]'a(2).x1='Devi';a(2).x2='[002,69]'结构成员的修改a(1).x2='[001,80]'第二讲Matlab数据及其运算-2.5结构数据

Matlab技术与应用2012-2013学年单元数据的建立b={10,'liu';12,'wang'}单元数据的引用b{1,2}第二讲Matlab数据及其运算-2.6单元数据

Matlab技术与应用2012-2013学年结构数据和单元数据的比较相同点：两者均可把不同属性的数据放在一个变量中不同点：两者对成员的引用方式不同

结构变量名.成员名第二讲Matlab数据及其运算Matlab技术与应用2012-2013学年Matlab语言把多项式表达成一个行向量，该向量中的元素是按多项式降幂排列的。

f(x)=anxn+an-1xn-1+……+loa0可用行向量p=[anan-1

……a1+a0]表示第二讲Matlab数据及其运算-2.7多项式运算Matlab技术与应用2012-2013学年conv例子：a(x)=x2+2x+3;b(x)=4x2+5x+6;c=(x2+2x+3)(4x2+5x+6)a=[123];b=[456];c=conv(a,b)c=4.0013.0028.0027.0018.00p=poly2str(c,'x')p=4x^4+13x^3+28x^2+27x+18第二讲Matlab数据及其运算-2.7多项式乘法Matlab技术与应用2012-2013学年deconv例子：a=[123];c=[4.0013.0028.0027.0018.00][d,r]=deconv(c,a)d:c除a后的整数r:余数第二讲Matlab数据及其运算-2.7多项式除法Matlab技术与应用2012-2013学年ployder例子：a=[12345];poly2str(a,'x')ans=x^4+2x^3+3x^2+4x+5b=polyder(a)b=4664poly2str(b,'x')ans=4x^3+6x^2+6x+4第二讲Matlab数据及其运算-2.7多项式微分Matlab技术与应用2012-2013学年[]：产生空矩阵zeros(m,n):产生全0矩阵zeros(size(A))ones(m,n);产生全1矩阵eye:单位矩阵rand:产生0~1间均匀分布的随机矩阵randn:产生均值为0，方差为1的标准分布随机矩阵第三讲Matlab矩阵分析与处理-3.1特殊矩阵

Matlab技术与应用2012-2013学年例子1：在区间[20,50]内均匀分布的5阶随机矩阵x=20+(50-20)*rand(5)例子2:均值为0.6，方差为0.1的5阶正态分布随机矩阵。y=0.6+sqrt(0.1)*randn(5)第三讲Matlab矩阵分析与处理-3.1特殊矩阵

Matlab技术与应用2012-2013学年矩阵求逆运算：对于一个方阵A，如果存在一个与其同阶的方阵B，使得：A*B=B*A=I(I为单位矩阵)则称B为A的逆矩阵，当然A也为B的逆矩阵Matlab中的调用函数：inv(A)实例：A=[1,-1,1;5,-4,3;2,1,1],求A的逆矩阵B第三讲Matlab矩阵分析与处理

3.2线性方程求解

Matlab技术与应用2012-2013学年广义逆矩阵：若矩阵A不是一个方阵，或者A是一个非满秩的方阵时，矩阵A没有逆矩阵，此时可以找到一个与A的转置矩阵同型的矩阵B，使得：A*B*A=AB*A*B=B此时，矩阵B为矩阵A的伪逆，也称为广义逆矩阵。Matlab中的调用函数：pinv(A)实例：A=[3,1,1,1;1,3,1,1;1,1,3,1]求A的广义逆矩阵第三讲Matlab矩阵分析与处理

3.2线性方程求解

Matlab技术与应用2012-2013学年线性方程组的矩阵表示第三讲Matlab矩阵分析与处理

3.2线性方程求解

Matlab技术与应用2012-2013学年线性方程组的矩阵表示第三讲Matlab矩阵分析与处理

3.2线性方程求解

a为m*n的系数矩阵。m=n,该方程为“恰定”方程m>n,该方程为“超定”方程m<n,该方程为“欠定”方程Matlab技术与应用2012-2013学年第三讲Matlab矩阵分析与处理

3.2线性方程求解

恰定方程求解方程ax+b(a为非奇异)x=a-1

b两种解:x=inv(a)b—采用求逆运算解方程x=a\b—采用左除运算解方程方程ax=ba=[12;23];b=[8;13]x=inv(a)*b

x=a\b

=

ax=b例:x1+2x2=82x1+3x2=13Matlab技术与应用2012-2013学年第三讲Matlab矩阵分析与处理

3.2线性方程求解

方程ax=b,m>n时此时不存在唯一解。方程解(a'a)x=a'bx=(a'

a)-1

a'b——求逆法x=a\bmatlab用最小二乘法找一个准确地基本解。第三讲Matlab矩阵分析与处理

3.2线性方程求解

超定方程求解Matlab技术与应用2012-2013学年例:x1+2x2=12x1+3x2=23x1+4x2=3a=[12;23;34];b=[1;2;3];

解1x=a\b

解2x=inv(a'*a)*a'*b

=

ax=bMatlab技术与应用2012-2013学年第三讲Matlab矩阵分析与处理

3.2线性方程求解

3.欠定方程组的解当方程数少于未知量个数时,即不定情况,有无穷多个解存在。matlab可求出两个解：（1）用除法求的解x是具有最多零元素的解（2）是具有最小长度或范数的解，这个解是基于伪逆pinv求得的。第三讲Matlab矩阵分析与处理

3.2线性方程求解

欠定方程求解Matlab技术与应用2012-2013学年x1+2x2+3x3=12x1+3x2+4x3=2a=[123;234];b=[1;2];x=a\bx=pinv(a)bx=x=1.000.8300.330-0.17=ax=b第三讲Matlab矩阵分析与处理

3.2线性方程求解

Matlab技术与应用2012-2013学年六、微分方程求解微分方程求解的仿真算法有多种，常用的有Euler（欧拉法）、RungeKutta(龙格-库塔法。Euler法称一步法，用于一阶微分方程。龙格-库塔法：实际上取两点斜率的平均斜率来计算的，其精度高于欧拉算法。第三讲Matlab矩阵分析与处理

3.3微分方程求解

Matlab技术与应用2012-2013学年六、微分方程求解实例：第三讲Matlab矩阵分析与处理

3.3微分方程求解

Matlab技术与应用2012-2013学年六、微分方程求解第三讲Matlab矩阵分析与处理

3.4多项式拟合多项式拟合x0=0:0.1:1；y0=[-.4471.9783.115.255.024.664.014.583.455.359.22];p=polyfit(x0,y0,3)p=56.6915-87.117440.0070-0.9043xx=0:0.01:1;yy=polyval(p,xx);figure;plot(xx,yy,'-b',x0,y0,'or')Matlab技术与应用2012-2013学年谢谢大家！Matlab技术与应用2012-2013学年六、微分方程求解第四讲Matlab绘图Matlab语言丰富的图形表现方法，使得数计算结果可以方便的、多样性地实现了可视化，这是其它语言所不能比拟的。Matlab技术与应用2012-2013学年一、二维绘图（一）plot——最基本的二维图形指令plot的功能：用直线连接相邻两数据点来绘制图形根据图形坐标大小自动缩扩坐标轴，将数据标尺及单位标注自动加到两个坐标轴上，可自定坐标轴，可把x,y轴用对数坐标表示Matlab技术与应用2012-2013学年如果已经存在一个图形窗口，plot命令则清除当前图形，绘制新图形可单窗口单曲线绘图；可单窗口多曲线绘图；可单窗口多曲线分图绘图；可多窗口绘图可任意设定曲线颜色和线型可给图形加坐标网线和图形加注功能一、二维绘图Matlab技术与应用2012-2013学年plot的调用格式

plot(x)——缺省自变量绘图格式，x为向量,以x元素值为纵坐标，以相应元素下标为横坐标绘图plot(x,y)——基本格式，以y(x)的函数关系作出直角坐标图，如果y为n×m的矩阵，则以x为自变量，作出m条曲线plot(x1,y1,x2,y2)——多条曲线绘图格式Matlab技术与应用2012-2013学年plot(x,y,’s’)——开关格式，开关量字符串s设定曲线颜色和绘图方式，使用颜色字符串的前1~3个字母，如yellow—yel表示等。或plot(x1,y1,’s1’,x2,y2,’s2’，…)Matlab技术与应用2012-2013学年S的标准设定值如下：

字母颜色标点线型

y黄色·

点线

m粉红○圈线

c亮蓝××线

r大红＋＋字线

g绿色－实线

b蓝色星形线

w白色：虚线

k黑色－·(--)点划线Matlab技术与应用2012-2013学年1.单窗口单曲线绘图例1：x=[0,0.48,0.84,1,0.91,0.6,0.14][x1,x2,x3,x4,x5,x6,x7,]plot(x)Matlab技术与应用2012-2013学年2.单窗口多曲线绘图例2：t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y,t,y1,t,y2)yy1y2Matlab技术与应用2012-2013学年3.单窗口多曲线分图绘图subplot——子图分割命令调用格式：

subplot(m,n,p)——按从左至右，从上至下排列

Matlab技术与应用2012-2013学年行列绘图序号subplot(1,3,1);plot(t,y)

subplot(1,3,2);plot(t,y1)

subplot(1,3,3);plot(t,y2)Matlab技术与应用2012-2013学年subplot(3,1,1);

plot(t,y)

subplot(3,1,2);

plot(t,y1)

subplot(3,1,3);

plot(t,y2)

Matlab技术与应用2012-2013学年4.多窗口绘图

figure(n)——

创建窗口函数，n为窗口顺序号。t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y)——

自动出现第一个窗口figure(2)plot(t,y1)——

在第二窗口绘图figure(3)plot(t,y2)——在第三窗口绘图Matlab技术与应用2012-2013学年5.可任意设置颜色与线型例4(例2)：plot(t,y,'r-',t,y1,'g:',t,y2,'b*')Matlab技术与应用2012-2013学年6.图形加注功能将标题、坐标轴标记、网格线及文字注释加注到图形上，这些函数为：

title——

给图形加标题xlable——

给x轴加标注ylable——

给y轴加标注text——

在图形指定位置加标注gtext——

将标注加到图形任意位置gridon(off)——

打开、关闭坐标网格线legend——

添加图例axis——

控制坐标轴的刻度Matlab技术与应用2012-2013学年例：t=0:0.1:10;y1=sin(t);y2=cos(t);figure;plot(t,y1,'r',t,y2,'b--');x=[1.7*pi;1.6*pi];y=[-0.3;0.8];s=['sin(t)';'cos(t)'];text(x,y,s);title('正弦和余弦曲线');legend('正弦','余弦')xlabel('时间t'),ylabel('正弦、余弦')gridMatlab技术与应用2012-2013学年

axis的用法还有：

axis([xminxmaxyminymax])——

用行向量中给出的值设定坐标轴的最大和最小值。如axis([-2205])axis(equal)——

将两坐标轴设为相等axison(off)——

显示和关闭坐标轴的标记、标志axisauto——

将坐标轴设置返回自动缺省值Matlab技术与应用2012-2013学年（二）fill––––

基本二维绘图函数

fill的功能：绘制二维多边形并填充颜色例：x=[12345];y=[41514];fill(x,y,'r')Matlab技术与应用2012-2013学年（三）特殊二维绘图函数bar––––

绘制直方图polar––––

绘制极坐标图

hist––––

绘制统计直方图

stairs––––

绘制阶梯图

stem––––

绘制火柴杆图

rose––––

绘制统计扇形图

comet––––

绘制彗星曲线Matlab技术与应用2012-2013学年errorbar––––

绘制误差棒图compass––––

复数向量图(罗盘图)feather––––

复数向量投影图(羽毛图)

quiver––––

向量场图

area––––

区域图

pie––––

饼图

convhull––––

凸壳图

scatter––––

离散点图Matlab技术与应用2012-2013学年例，绘制阶梯曲线x=0:pi/20:2*pi;y=sin(x);stairs(x,y)例：阶梯绘图h2=[11;1-1];h4=[h2h2;h2-h2];h8=[h4h4;h4-h4];t=1:8;subplot(8,1,1);stairs(t,h8(1,:));axis('off')subplot(8,1,2);stairs(t,h8(2,:));axis('off')subplot(8,1,3);stairs(t,h8(3,:));axis('off')subplot(8,1,4);stairs(t,h8(4,:));axis('off')subplot(8,1,5);stairs(t,h8(5,:));axis('off')subplot(8,1,6);stairs(t,h8(6,:));axis('off')subplot(8,1,7);stairs(t,h8(7,:));axis('off')subplot(8,1,8);stairs(t,h8(8,:));axis('off')h2=[11;1-1];h4=[h2h2;h2-h2];h8=[h4h4;h4-h4];t=1:8;fori=1:8subplot(8,1,i);stairs(t,h8(i,:))axis('off')end例：绘制极坐标绘图t=0:2*pi/90:2*pi;y=cos(4*t);polar(t,y)例：绘制火柴杆绘图

t=0:0.2:2*pi;y=cos(t);stem(y)例：绘制直方图

t=0:0.2:2*pi;y=cos(t);bar(y)例：绘制彗星曲线图

t=-pi:pi/500:pi;y=tan(sin(t))-sin(tan(t));comet(t,y)x=magic(6);area(x)x=[1234567];y=[0001000];

pie(x,y)pie(x,y,{'North','South','East','West',

'middle','fa','white'})loadseamount

scatter(x,y,50,z)a=rand(200,1);b=rand(200,1);

c=rand(200,1);

scatter(a,b,100,c,'p')二、三维绘图三维绘图的主要功能：绘制三维线图绘制等高线图绘制伪彩色图绘制三维网线图绘制三维曲面图、柱面图和球面图绘制三维多面体并填充颜色（一）三维线图plot3——

基本的三维图形指令调用格式：plot3(x,y,z)——x,y,z是长度相同的向量plot3(X,Y,Z)——X,Y,Z是维数相同的矩阵plot3(x,y,z,s)——

带开关量plot3(x1,y1,z1,'s1',x2,y2,z2,'s2',…)二维图形的所有基本特性对三维图形全都适用。定义三维坐标轴大小

axis([xminxmaxyminymaxzminzmax])gridon(off)绘制三维网格

text(x,y,z,‘string’)三维图形标注子图和多窗口也可以用到三维图形中例：绘制三维线图t=0:pi/50:10*pi;plot3(t,sin(t),cos(t),'r:')(二)三维饼图

pie3([43689])（三）三维多边形fill3=fill——

三维多边形的绘制和填色与二维多边形完全相同调用格式：

fill3(x,y,z,‘s’)——

与二维相同例：用随机顶点坐标画出5个粉色的三角形，并用黄色的○表示顶点y1=rand(3,5);y2=rand(3,5);y3=rand(3,5);fill3(y1,y2,y3,'m');holdon;plot3(y1,y2,y3,'yo')（四）三维网格图mesh——

三维网线绘图函数调用格式：

mesh(z)——z为n×m的矩阵，x与y

坐标为元素的下标

mesh(x,y,z)——x,,y,z分别为三维空间的坐标位置例，矩阵的三维网线图z=rand(6);

0.88080.07290.41680.76940.37750.47760.33810.71010.09640.63520.38260.70860.18950.87910.67470.89650.68760.23800.74310.35940.56260.87840.12170.39100.71890.08990.81300.48650.07680.97590.87920.16100.87820.11310.14330.6288z=round(z)

100100010101011110101100101001101001mesh(z)；mesh(z)例：8阶hadamard矩阵的网线图h2=[11;1-1];h4=[h2h2;h2-h2]h8=[h4h4;h4-h4]111111111-11-11-11-111-1-111-1-11-1-111-1-111111-1-1-1-11-11-1-11-1111-1-1-1-1111-1-11-111-1mesh(h8)三维网线图作图要领生成坐标——[X,Y]=meshgrid(x,y)表达式点运算——Z=X.^2+Y.^2X，Y是n×m的矩阵，维数可任定X—nY—mZ—默认方位角：37。5º，俯角30º。n×m维meshgrid——

网线坐标值计算函数

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)坐标矩阵坐标向量纵坐标矩阵绘图函数colormap([R,G,B])——

色图设定函数matlab的颜色数据集合为红、绿、兰三颜色矩阵[R，G，B]，维数m×3r,g,b在[01]区间连续取值，理论上颜色种类可达无穷多种matlab使用三维向量表示一种颜色，常用颜色数据见下表饱和色

[000]—

黑色

[001]—

兰色

[010]—

绿色

[011]—

浅兰

[100]—

红色

[101]—

粉红

[110]—

黄色

[111]—

白色

调和色[0.50.505]—

灰色[0.500]—

暗红色[10.620.4]—

铜色[0.4910.8]—

浅绿[0.4910.83]—宝石兰（五）三维曲面图surf——

三维曲面绘图函数，与网格图看起来一样与三维网线图的区别：网线图：线条有颜色，空挡是黑色的（无颜色）曲面图：线条是黑色的，空挡有颜色（把线条之间的空挡填充颜色，沿z轴按每一网格变化调用格式：surf(x,y,z)——

绘制三维曲面图，x,y,z

为图形坐标向量例：[X,Y,Z]=peaks(30)※peaks为matlab自动生成的三维测试图形surf(X,Y,Z)

surfc(X,Y,Z)—

带等高线的曲面图

[X,Y,Z]=peaks(30);surfc(X,Y,Z)surfl(X,Y,Z)——

被光照射带阴影的曲面图[X,Y,Z]=peaks(30);surfl(X,Y,Z)

cylinde(r,n)—

三维柱面绘图函数

r为半径；n为柱面圆周等分数例：绘制三维陀螺锥面t1=0:0.1:0.9;