数学建模案例分析MATLAB在电气工程中的应用.pptx

《数学建模案例分析MATLAB在电气工程中的应用.pptx》由会员分享，可在线阅读，更多相关《数学建模案例分析MATLAB在电气工程中的应用.pptx（82页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2023/2/141MATLAB简介对一般的计算语言来说，必须采用两层循环才能得到结果:vari,j:integer;/用pascal语言x,y:array0.10ofarray0.20ofreal;constm=10;n=20;beginfori:=1tomdoforj:=1tondoyij:=exp(xij)*sin(xij);end;MATLAB只用一条指令:y=exp(-2*x).*sin(5*x)。例:已知x是(nm)维数组，计算：第1页/共82页2023/2/142MATLAB简介例：对于求解代数方程问题：当A的列数大于行数时，有无数解。MATLAB只用一条指令:x=Ab。当A是标

2、量时：当A是非奇异矩阵时：当A是行数大于列数的满秩阵时，称超定方程组或矛盾方程组：第2页/共82页2023/2/143MATLAB简介在MATLAB中仅用以下指令：ezsurf(y/(1+x2+y2)例：绘制二元函数曲面：第3页/共82页2023/2/144MATLAB简介MATLAB是专门为科学计算而设计的，具有强大数据图示功能的工具软件。语法结构简明、数值计算高效、图形功能完备、易学易用。在矩阵代数、数值计算、数字信号处理、震动理论、神经网络控制、动态仿真、建模、系统识别、控制系统、非线性系统、模糊控制、优化技术、通讯系统、财政金融等领域有着广泛应用。第4页/共82页2023/2/145教

3、学内容n第一章MATLAB基础n第二章MATLAB基本语法n第三章数组及矩阵运算n第四章数值计算n第五章符号运算n第六章MATLAB计算结果可视化n第七章MATLAB在系统仿真中的应用n第八章MATLAB在信号处理中的应用第5页/共82页2023/2/146通过本课程学习，使学生掌握利用MATLAB进行数值计算的基本方法，熟悉MATLAB编程环境、语 言 语 法、程 序 结 构、编 程 及 调 试 技 术，掌 握MATLAB中M文件、M函数编写方法及调试技术、MATLAB的绘图和图形控制函数等内容，上机练习MATLAB数值解算方法，具备上机操作的技能，学习MATLAB在电气工程学科中的建模与分

4、析方法，为后续专业课程学习奠定基础。课程任务第6页/共82页2023/2/147考试形式考试方式：1、平时成绩(出勤率+课堂习题报告)2、考试第7页/共82页2023/2/148MATLAB的发展历程MATLAB 是 Matrix Laboratory的 缩 写。20世 纪 70年 代 由 CleveMoler博士用FORTRAN语言开发了最初版本的矩阵分析软件。80年代初期，由Moler、JohnLittle和SteveBangert采用C语言改写了MATLAB的内核。并于1984年成立了Mathworks软件开发公司，将MATLAB正式推向市场。1988年推出了MATLAB3.1（DOS）

5、版本。1993年，推出了应用于Windows操作系统的MATLAB4.0版本。1996年12月，推出MATLAB5.0版本(Release8)。1997年5月，推出MATLAB5.1版本(Release9)。1998年3月，推出MATLAB5.2版本(Release10)。1999年11月，推出MATLAB5.3版本(Release11)。2000年11月，推出MATLAB6.0版本(Release12)。2001年6月，推出MATLAB6.1版本(Release12.1)。2002年7月，推出MATLAB6.5版本(Release13)。第8页/共82页2023/2/149MATLAB的发展

6、历程2004年6月，推出MATLAB7.0版本(Release14)。2004年11月，推出版本(Release14SP1)。2005年3月，推出版本(Release14SP2)。2005年9月，推出MATLAB7.1版本(Release14SP3)。2006年3月，推出MATLAB7.2版本(Release2006a)。2006年9月，推出MATLAB7.3版本(Release2006b)。2007年3月，推出MATLAB7.4版本(Release2007a)。2007年9月，推出MATLAB7.5版本(Release2007b)。2008年3月，推出MATLAB7.6版本(Release2

7、008a)。2008年10月，推出MATLAB7.7版本(Release2008b)。2009年3月，推出MATLAB7.8版本(Release2009a)。2009年9月，推出MATLAB7.9版本(Release2009b)。2010年3月，推出MATLAB7.10版本(Release2010a)。第9页/共82页2023/2/1410MATLAB的主要特点 MATLAB 的基本单位为矩阵。MATLAB 语言以解释方式工作。具有非常友好的人机界面。具有强大的作图和数据可视化功能.智能化程度高。具有极强的可扩展性。第10页/共82页2023/2/1411第一章 MATLAB 基础第11页/共

8、82页2023/2/14121.1MATLAB的基本组成MATLAB软件主要由主包、Simulink 动态系统仿真和工具箱（Toolbox）三大部分组成。*主包包括了 MATLAB 语言、数学函数库、工作环境、图形处理系统和应用程序接口五部分。*Simulink是用于系统仿真的交互式图形系统，是MATLAB非常重要的共生产品。第12页/共82页2023/2/1413 *Toolbox实际是用MATLAB的基本语句编成的开放式各种子程序集，用于解决各专业领域的特定问题，例如：控制系统、信号处理、图像处理、模糊集合、神经网络、小波分析等。*MATLAB编译器，可以将M文件转化为C或C+语言代码，增

9、强了MATLAB灵活性。第13页/共82页2023/2/14141.2MATLAB系统界面MATLAB的启动桌面快捷方式启动开始菜单方式启动MATLAB的退出Matlab菜单命令/关闭按钮在Matlab命令窗口输入exit或quit第14页/共82页2023/2/1415MATLAB的用户界面窗口CommandwindowWorkspaceHistorycommandwindowCurrentdirectoryLaunchpadM-FileFigureModelGUI第15页/共82页2023/2/1416历史命令历史命令发射台发射台命令命令窗口窗口标题栏标题栏菜单栏菜单栏常用工具栏常用工具栏

10、工作区间工作区间当前目录当前目录第16页/共82页2023/2/1417第17页/共82页2023/2/1418 MATLAB命令窗口lMATLAB语句形式变量表达式；当键入回车键时，该语句被执行。语句执行之后，窗口自动显示出语句执行的结果。如果希望结果不被显示，则需在语句之后加上分号。此时尽管结果没有显示，但它依然被赋值并在MATLAB工作空间中分配了内存。l 方向键可以编辑修改已输入的命令 ：调出上一行语句；：调出下一行语句；第18页/共82页2023/2/1419l 多行命令（）如果命令语句超过一行或者太长希望分行输入，则可以使用多行命令继续输入。例如：S=1-12+13+4+9-4-1

11、8;等价于：S=1-12+13+4+9-4-18;命令窗口的分页输出 more off：不允许分页 more on：允许分页 more(n)：指定每页输出的行数 回车键显示下一行，空格键显示下一页，q结束当前显示。第19页/共82页2023/2/1420l页面显示的疏密控制formatloose(默认):稀疏显示格式；formatcompact:密集显示格式；l清命令窗口 clc第20页/共82页2023/2/1421显示所有目前内存中的MATLAB变量、数学结构、类型以及所占的空间并可进行变量的编辑等。MATLAB的workspace中的变量管理who：显示当前工作空间中所有变量的一个简单列

12、表；whos：显示出变量的大小、数据格式等详细信息；clear：清除工作空间中所有的变量；clear变量名：清除指定的变量；MATLAB工作空间 第21页/共82页2023/2/1422lMATLAB的workspace中的变量保存和载入save：将工作空间的所有变量存到matlab.mat的二进制文件中。savefilename：将所有变量保存到磁盘文件filename中。savefilenamevariables：将variables所表示的变量列表保存到磁盘文件filename中。各个不同的变量之间只能用空格来分隔。load（filename）（variables）：将数据文件载入MAT

13、LAB工作空间。第22页/共82页2023/2/1423 MATLAB历史命令窗 MATLAB当前目录l各种历史命令输入的记载，便于随时查看和调用，通过M文件运行的命令不显示。l显示和改变当前目录，显示目录下的文件并提供搜索功能。第23页/共82页2023/2/1424helpwin指令在命令窗口中键入helpwin，进入帮助窗口；helpdesk指令在命令窗口中键入helpdesk（或doc，或点击工具条中的？按钮），进入帮助窗口，显示HTML格式的帮助内容。1.3 MATLAB帮助系统第24页/共82页2023/2/1425help命令help：列出所有的帮助主题，每个帮助主题对应于MAT

14、LAB搜索路径中的一个目录；help库名：得到库中全部函数名；help函数名或命令；得到函数或命令的帮助信息；lookfor命令如在命令窗口中键入：lookforsu，则只要有与su匹配的单词就显示出来。其收索条件比help命令宽。demo命令demo：可得到MATLAB及其所以工具箱中的全部演示程序。是学习MATLAB的一种很好的途径。第25页/共82页2023/2/14261.5MATLAB两种工作方式指令行操作方式键入命令行。执行M文件方式键入程序名。1.4 MATLAB常用命令第26页/共82页2023/2/1427第二章 MATLAB 基本语法第二章 数组及矩阵运算第27页/共82页

15、2023/2/1428l任何MATLAB的语句的执行结果都可以在屏幕上显示，同时赋值给指定的变量，没有指定变量时，赋值给一个特殊的变量ans，数据的显示格式由format命令控制。lformat只是影响结果的显示，不影响其计算与存储；MATLAB总是以双字长浮点数（双精度）来执行所有的运算。数的显示格式 2.1 数值、变量及其赋值第28页/共82页2023/2/1429format命令控制数的显示格式：format(short)(默认)：如pi=3.1416;formatlongformatrat：如pi=355/113;formatshorte：如pi=3.1416e+000;formatl

16、onge：如pi=;formatbank：2位十进制，如pi=3.14;formathex：如pi=400921fb54442d18;format+：正数、零、负数分别用“+”、空格、“-”表示第29页/共82页2023/2/1430常量MATLAB的保留常量和变量保留常量不能重复定义！第30页/共82页2023/2/1431MATLAB语言的标识符命名规则标识符：标志常量名、变量名、函数名和文件名的字符串的总称。标识符的合法字符是52个英文字母（区分大小写）、10个数字和下划线；第一个字符必须是英文字母（不能超过63个字符）；标识符必须是不含空格的单个词；第31页/共82页2023/2/14

17、32矩阵和变量的赋值直接赋值语句赋值变量表达式（或常量）表达式由运算符、函数、变量名和数字组成，运 算 符 有+、-、*、/、等。（ex201）函数调用语句返回变量列表函数名（输入变量列表）例：a,b,c=my_fun(d,e,f,c)冒号表达式v=s1:s2:s3其中s1是起始值，s2是步长（若省略步长为1），s3是最大值。（ex202）第32页/共82页2023/2/1433复数i，j为MATLAB中的虚数符号，可以不用乘号，直接连写在数字后面。如果在程序中另外给i，j赋值，它们的虚数意义就失效。conj(z)命令表示对复数矩阵z进行共轭变换。Z命令表示对复数矩阵z进行共轭转置变换。（ex

18、203）第33页/共82页2023/2/1434 特殊矩阵的赋值函数linespace(a,b,N)均分向量，初值为均分向量，初值为a,终值为终值为b，元素个数为，元素个数为Nlogspace(a,b,N)向量初值为向量初值为10a,终值为终值为10b，元素个数为，元素个数为Nones(n,m)产生产生n*m阶全阶全1矩阵矩阵rand(n,m)产生产生n*m阶随机数矩阵阶随机数矩阵randn(n,m)产生产生n*m阶正态随机数矩阵阶正态随机数矩阵zeros(n,m)产生产生n*m阶全零矩阵阶全零矩阵eye(n)产生产生n*n阶单位方阵阶单位方阵（ex204）第34页/共82页2023/2/14

19、35 2.2 矩阵的初等运算*矩阵的加减法矩阵A和B进行加减运算的条件是：size(B)=size(B）*矩阵的乘法1）矩阵A*B的条件是：A的列数等于B的行数。2）不符合交换律，即A*BB*A（ex205）第35页/共82页2023/2/1436*矩阵的除法1）左除若线性方程组表示为A*X=B，如果A为非奇异方阵，即其逆阵存在，则方程组的解 用 MATLAB表 示 为：X=inv(A)*B 或 X=AB符号“”称为左除。2）右除若线性方程组表示为X*C=D，如果C为非奇异方阵，即其逆阵存在，则方程组的解 用 MATLAB表 示 为：X=D*inv(C)或 X=D/C符号“/”称为右除。第36

20、页/共82页2023/2/1437例：求解线性方程组的解 将该线性方程组表示为矩阵相乘为：（ex206）第37页/共82页2023/2/1438*数组运算矩阵运算遵循线性代数中的矩阵运算规则，数组运算将矩阵按逐个元素分别进行运算。MATLAB中几乎所有的初等函数都适合数组运算。运算运算矩阵矩阵数组数组运算符运算符规则规则运算符运算符规则规则乘法乘法*A*B:A的列数的列数=B的行数的行数.*同维数组同维数组幂幂方阵方阵.数数组组的的标标量量乘乘方方A.p：A对对应应元元素素的的p次次方方，结结果果和和A同维同维标量的数组乘方标量的数组乘方p.A：p的数组元素次方，结果和的数组元素次方，结果和A

21、同维同维左除左除BA=inv(B)*A./或或.A./B=B.A右除右除/B/A=B*inv(A)（ex207）第38页/共82页2023/2/14392.3MATLAB语言流程控制关系运算六种关系运算符：=，=，=运算结果是逻辑1（真）或0（假）逻辑运算逻辑运算符：&（与），|（或），（非）运算结果是逻辑1（真）或0（假）逻辑关系函数：（见讲义第三章P26）第39页/共82页2023/2/1440循环结构for结构while结构第40页/共82页2023/2/1441(ex208)、(ex209)(ex210)、(ex211)第41页/共82页2023/2/1442 顺序结构第42页/共82

22、页2023/2/1443表达式语句段if真假表达式语句段1if真假语句段2else第43页/共82页2023/2/1444(ex212)表达式语句段1if真假语句段2elseif表达式语句段3真假else第44页/共82页2023/2/1445 开关结构开关结构第45页/共82页2023/2/1446第46页/共82页2023/2/1447当开关表达式的值等于某表达式，执行该语句后结束该结构。程序的执行结果和各个case顺序无关。当需要在开关表达式满足若干个表达式之一时执行某一程序段，则用单元形式（用大括号把这些表达式括起来，用逗号分隔）case语句中条件不能重复，否则列在后面的条件将不能执行

23、（ex213）第47页/共82页2023/2/1448 其它控制语句其它控制语句1）break语句终止循环语句的运行，对于多层嵌套循环结构，break只能终止最里层的循环。2）return语句中断当前函数运行返回调用处。3）pause语句pause：暂停程序运行，如果敲击键盘，程序继续运行；pause(n):暂停n秒；pauseon：该语句后的pause语句有效；pauseoff：该语句后的pause语句无效；第48页/共82页2023/2/1449例：用循环求解使不等式成立的最小的m。(ex214)第49页/共82页2023/2/1450M文件M文件分类命令式文件Scriptfile函数式文

24、件Functionfile2.4 M文件及程序调试第50页/共82页2023/2/1451函数式文件的构成1）function定义行：function输出参量=函数名(输入参量)2）函数体：实现函数的功能。3）函数说明。函数式文件的特点1）输入、输出参量是用来进行变量传递的，可以有多个。2）程序中的变量均为局部变量，不保存在工作空间中，除非用global声明为全局变量。3）文件名必须与函数名相同。（ex215）第51页/共82页2023/2/1452分分析析：本本题题中中，需需要要在在R1=1k，R1=10k，R1=20k时时分分别别计计算算Uc(t)值值，可可编编写写一一个个函函数数实实现现

25、的计算，传递参数是的计算，传递参数是R1、C1、t，t是初始化后的数组，返回值是初始化后的数组，返回值Uc。第52页/共82页2023/2/1453解：解：(1)编写函数编写函数CalUcfunctionUc=CalUc(R1,C1,t)tao=R1*C1;%计算时间常数计算时间常数Uc=20*(1-exp(-t./tao);%返回曲线值返回曲线值将上述语句存为文件将上述语句存为文件CalUc.m。(2)编写主程序编写主程序t=0:0.001:0.04;%t初始化初始化Uc1=CalUc(1000,1e-6,t);%调用调用CalUc函数计函数计算算R=1k时的函数值时的函数值Uc10=Cal

26、Uc(10000,1e-6,t);Uc20=CalUc(20000,1e-6,t);plot(t,Uc1,t,Uc10,t,Uc20);第53页/共82页2023/2/1454运行结果图第54页/共82页2023/2/1455两类文件的区别 1）命令式文件没有输入输出参量，函数式文件有；2）命令式文件运行中可以调用工作空间的数据，运行中产生的所以变量为全局变量。函数式文件不能调用工作空间的数据，运行中产生的所有变量为局部变量;3）命令式文件可以直接运行，函数式文件不能，只能调用。第55页/共82页2023/2/1456M文件的调试MATLAB程序两种错误：语法和逻辑错误。文件编辑器的使用1）用

27、颜色区分内容；2）自动缩进功能，用于判断循环层次。调试方法1）根据命令窗中显示的出错信息调试。2）通过程序编辑器提供的调试工具调试。(ex216)3）取消语句后的“;”，使程序运行时在命令窗口显示变量的值。第56页/共82页2023/2/14572.5 MATLAB开发环境 与其它软件的接口关系 MATLAB的文件管理系统第57页/共82页2023/2/1458与其它软件的接口关系人机交互命令input和keyboard(以return退出）与磁盘操作系统关系与word的接口与记事本的接口(diary命令)低层输入输出命令库(iofun)与C或FORTRAN子程序的动态链接第58页/共82页2

28、023/2/1459 MATLAB的文件管理系统安装后的MATLAB文件系统MATLAB中特有的文件格式1）M文件2）数据文件3）可执行文件4）图形文件5）模型文件、仿真文件等文件管理和搜索路径1）MATLAB管理的文件范围是由它的搜索路径来确定的。2）程序文件夹设定在“搜索路径”下的方法。第59页/共82页2023/2/1460第四章 数值计算第60页/共82页2023/2/14614.1多项式运算多项式的表示法多项式：P=an,an-1,.,a1,a0多项式的显示：poly2sym(p)；poly2str(p,x)；多项式的乘法和除法(conv()、deconv()c=conv(a,b);

29、q,r=deconv(c,b);多项式求根(roots()已知多项式的根，创建多项式(poly()r=roots(p);%返回向量r为一列向量p=poly(r1);第61页/共82页2023/2/1462求多项式的值(polyval()y=polyval(p,x);%x、y为矩阵、向量、标量。y=polyvalm(p,x);%x必须为方阵。多项式展开(residue()将分式多项式转换成部分分式。在电路的复频域分析中将复杂的传递函数转换成简单的传递函数之和。r,p,k=residue(b,a)b,a=residue(r,p,k)多项式求导(polyder()y=polyder(p);第62页/

30、共82页2023/2/1463多项式拟合(polyfit()p=polyfit(x,y,n);%x、y为数据向量，p为多项式系数向量，n为多项式阶次。一元插值(interp1()yi=interp1(x,y,xi,method);%计算xi对应的函数值yi，x、y为数据向量。method插值方法有：linear：线性插值；cubic:三次多项式插值；nearest：最近点插值；spline：三次样条插值。4.2数值分析第63页/共82页2023/2/1464数值积分(trapz()、quad()、quad8()梯形法：z=trapz(y)；%近似求解y的积分近似值。y为被积函数值，点距为1；z

31、=trapz(x,y)；%近似求解y对x的积分值，x、y必须是相同维数的向量。辛普生法：z=quad(fun,a,b)；%计算函数fun在a,b区间的相对误差为1.e-6的积分近似值，函数fun必须返回一个在x上的函数值向量。z=quad(fun,a,b,tol)；%与上相同，tol用以指定误差容度。牛顿-科特斯法：quad8调用格式和quad方式相同。第64页/共82页2023/2/1465数值微分(利用求差分的函数diff()实现)差分计算Y=diff(X);%用以计算数组中的差分。X=X(1),X(2),.,X(n)，Y=X(2)-X(1),X(3)-X(2),.,X(n)-X(n-1)

32、微分计算：设函数y=f(x)，可近似表示为：因此，可以用diff(y)./diff(x)近似计算函数的微分。第65页/共82页2023/2/1466常微分方程(ODE)的数值解高阶常微分方程的表示形式：求解方法：将高阶常微分方程转化为一阶常微分方程组。，第66页/共82页2023/2/1467l函数调用格式（7个函数调用格式完全一样）T,Y=ode45(odefun,tspan,y0)%odefun代表显式ODE方程组y=f(t,y)，t为标量，y和y为同长度的列向量。Tspan可以是区间t0,tf，也可以是t0,t1,tf，函数返回指定时间t0,t1,t2,上的解。y0是初始条件，是与y同长

33、度的列向量。解矩阵Y中的每一行对应于返回的时间列向量T中的一个时间点。T,Y=ode45(odefun,tspan,y0,options)%options用于设定微分方程解法器的参数（用函数odeset()生成）。T,Y=ode45(odefun,tspan,y0,options,p1,p2)第67页/共82页2023/2/1468函数最小值(fminbnd()、fminsearch()fminbnd：单变量非线性函数局部最小值函数 x=fminbnd(FUN,x1,x2)x=fminbnd(FUN,x1,x2,options)x=fminbnd(FUN,x1,x2,options,p1,p2

34、,)%返回值x是FUN函数在(x1,x2)上的局部最小值，FUN为单值非线性函数。options用来控制算法的参数设置。fminsearch为多变量无束缚非线性局部最小值函数，调用格式 与上一致。第68页/共82页2023/2/1469函数零点(fzero()xzero=fzero(FUN,x0)xzero=fzero(FUN,x0,tol,trace)%FUN：待求零点的函数；x0：搜索零点的初始值；tol：搜索相对精度；trace：指定迭代信息是否显示，默认为0，表示不显示迭代信息。第69页/共82页2023/2/1470第五章 符号运算第70页/共82页2023/2/14715.1 符号

35、对象和符号表达式符号对象定义基本的符号对象(常数、变量、表达式)f=sym(arg);%把数值、字符串或表达式arg转换为符号对象f=sym(arg,flagn);%把数值或数值表达式arg转换为flagn格式的符号对象。falgn可设定为：d：采用最接近的十进制精确表示；d：采用最接近的有理表示，缺省设置。syms(arg1,arg2);%把字符arg1，arg2定义为基本符号对象symsarg1,arg2;第71页/共82页2023/2/1472符号表达式和符号函数的操作符号表达式的化简collect()：合并同幂项simplify()：利用代数上的函数规则对表达式进行化简simple()

36、：以尽可能的办法将表达式再做化简，目的是使表达式以最少的字表示出来。horner()：对多项式进行嵌套型分解。第72页/共82页2023/2/1473符号函数的求反和复合g=finervse(f)：符号函数f和g满足：g(f(x)=xg=finervse(f,v)：符 号 函 数 g的 自 变 量 为 v，满 足g(f(x)=vg1=compose(f,g)：返回f=f(x)和g=g(y)时的复合函数f=f(g(y)g1=compose(f,g,z)：返回复合函数f=f(g(z)g1=compose(f,g,x,z)：返回复合函数f=f(g(z)，x是f的独立变量g1=compose(f,g,

37、x,y,z)：返回复合函数f=f(g(z)，x是f的独立变量，y是g的独立变量第73页/共82页2023/2/1474符号微积分g=int(f)：返回f的不定积分。g=int(f,v)：对f中指定的符号变量v计算不定积分。g=int(f,v,a,b)：对表达式f中指定的符号变量v计算从a到b的定积分。符号代数方程组的求解一般代数方程组的求解一般代数方程组包括线性、非线性和超越方程等，求解用函数solve()第74页/共82页2023/2/1475第六章 MATLAB计算结果可视化第75页/共82页2023/2/1476连续函数和离散函数的可视化离散函数的可视化对于离散实函数yn=f(xn)，其

38、自变量用向量形式表示为x=x1,x2,xn，由函数关系可得函数值向量为y=y1,y2,yn。该组向量可用平面上的一组点表示。连续函数的可视化连续函数的可视化也必须由一组离散自变量计算相应的函数值。离散的点不能表现函数的连续性，为了表示离散点之间的函数情况，有两种常用处理方法：1、把两点用直线连接，近似表现两点间的函数变化规律。2、对区间进行更细的分割，计算更多的点，去近似表现函数的连续变化。第76页/共82页2023/2/1477绘制二维图形的基本操作一般步骤1、数据准备：1）产生自变量采样向量e.g.t=0:0.001:1;2）计算相应的函数值向量e.g.y=sin(t).*sin(3*t)

39、+5*cos(12*t);2、选定图形串口及子图位置：1）用figure()函数创建或指定图形窗口e.g.figure(n);%指定n号图形窗2）指定子图号 e.g.Subplot(2,2,4);%指定4号子图 第77页/共82页2023/2/14783、调用绘图指令：e.g.plot(t,y,r-);%用红色虚线画曲线4、设置坐标轴的范围、刻度和网格线：e.g.axis(0,1,-6,6);%设置轴的范围gridon;5、图形注释：图名、坐标名、图例、文字说明e.g.title(波形)；%图名xlabel(t);ylabel(y);%轴名legend(曲线1);%图例text(2,1,y=s

40、in(t)sin(3t)+5cos(12t);%文字说明 第78页/共82页2023/2/14796、图形的精细修饰（图柄操作）1）利用对象属性值设置2）利用图形窗工具条进行说明：步骤一和三是最基本的绘图步骤，其他步骤并不完全必须。由这两步所画出的图形已具备足够的表现力。步骤四、五、六的前后次序不是一定的，可根据需要或习惯调整。步骤二一般在图形较多的情况下使用，以指定图形窗和子图。第79页/共82页2023/2/1480二维绘图指令plot()是最重要、最基本的指令。1、plot(x,s);x是实向量时，以该向量元素的下标为横坐标、元素值为纵坐标画出一条连续曲线；x是实矩阵时，按列绘制每列元素

41、值相对其下标的曲线。曲线数等于矩阵列数；x是复数矩阵时，则按列分别以元素实部和虚部为横、纵坐标绘制多条曲线。第80页/共82页2023/2/14812、plot(x,y,s);x、y是同维向量时，绘制以x,y元素为横坐标、纵坐标的曲线；x是向量，y是一维与x同维的实矩阵时，绘出多根不同色彩的曲线，曲线数等于y矩阵的另一维数；x是矩阵，y是向量时，绘出的多根曲线均以y为共同纵坐标。x、y是同维矩阵时，则以x，y对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。3、plot(x1,y1,s,x2,y2,s,x3,y3,s,);每三个绘图参数xi,yi,s构成一组，其结构和作用与plot(x,y,s)相同，不同参数组之间没有约束关系。第81页/共82页2023/2/1482感谢您的观看！第82页/共82页