MATLAB：经典控制系统的分析和设计 9.3.1 基本操作及命令1. 访问和退出MATLAB在大多数系统中，一旦安装了MATLAB，在调用时，应执行命令MATLAB。退出MATLAB应执行命令exit或quit。2. 如何应用MATLAB通常以命令驱动方式应用MATLAB。当输入单个命令时，MATLAB会立即对其进行处理，并且显示处理结果。MATLAB也能够执行存储在文件中的命令序列。通过键盘输入的命令，应用向上箭头键可以被存取。通过输入某个最新命令和调用特定的命令行，可以使屏幕内容向上滚动。3. MATLAB的变量MATLAB的一个特点是变量在应用之前不必是维数确定的。在MATLAB中，变量一旦被采用，会自动产生（如果必要，变量的维数以后还可以改变）。在命令exit或quit输入之前，这些变量将保留在存储器中。为了得到工作空间内的变量清单，可以通过键盘输入命令who，当前存放在工作空间内的所有变量便会显示在屏幕上。命令clear能从工作空间中清除所有非永久性变量。如果只需要从工作空间中清除某个特定变量，比如“x”，则应输入命令clear x。4. 以“%”开始的程序行在本书中，许多MATLAB程序在编写时附有注解和说明，这些注解和说明阐明了发生在程序中的具体进程。在MATLAB中以“%”开始的程序行，表示注解和说明。符号“%”类似于BASIC中“REM”。以“%”开始的行，用来存储程序的注解或说明，这些注解和说明是不执行的。这就是说，在MATLAB程序行中，出现在“%”以后的一切内容都是可以忽略的。如果注解或说明需要一行以上程序行，则每一行均需以“%”为起始。5. 应用分号操作符分号用来取消打印。如果语句的最后一个符号是分号，则打印被取消，但是命令仍在执行，而结果不再显示。这是一个有益的特性，因为打印中间结果可能不必要。此外，在输入矩阵时，除非最后一行，分号用来指示一行的结束。6. 应用冒号操作符冒号操作符在MATLAB中起着重要作用。该操作符用来建立向量，赋予矩阵下标和规定叠代。例如，j:k表示j j+1 k，A(:，j)表示矩阵A的第j列，A(i，:)表示矩阵A的第i行。7. 输入超过一行的长语句一个语句通常以回车键或输入键终结。如果输入的语句太长，超出了一行，则回车键后面应跟随由3个或3个以上圆点组成的省略号()，以表明语句将延续到下一行。下面是一个例子：x = 1.234+2.345+3.456+4.567+5.678+6.789+7.890+8.901-9.012符号=、+和-前后的空白间隔可以任选。这种间隔通常可以起到改善语句清晰度的效果。8. 在一行内输入数个语句如果在一行内可以把数个语句用逗号或分号隔开，则可以把这数个语句放在一行内。例如plot(x,y,o), text(1,20,System 1), text(1,15,System 2)和plot(x,y,o)；text(1,20,System 1)；text(1,15,System 2)9. 选择输出格式MATLAB中的所有计算都是以双精度方式完成的，但是显示输出可以是具有4个小数位的定点输出。例如，对于向量x = 1/3 0.00002MATLAB有下列输出：x=0.3333 0.0000如果在矩阵中至少有一个元素不是严格的整数，则有4种可能的输出格式。显示的输出量可以利用下列命令加以控制：format shortformat longformat short eformat long e一旦调用了某种格式，则这种被选用的格式将保持，直到对格式进行了改变为止。在控制系统分析中，format short和format long是经常采用的格式。一旦调用了MATLAB，即使没有输入格式命令，MATLAB也将以format short格式显示数值结果。如果矩阵或向量的所有元素都是严格的整数，则format short和format long的结果是相同的。10. 退出MATLAB时如何保存变量当键入“exit”或“quit”时，MATLAB中的所有变量将消失。如果在退出以前输入命令save，则所有的变量被保存在磁盘文件matlab.mat中。当再次进入MATLAB时，命令load将使工作空间恢复到以前的状态。9.3.2 MATLAB函数在当前的MATLAB5.0及以上版本中，MATLAB函数的调用格式和其它编程语言是不同的，其典型的调用格式为返回变量列表 = func_name(输入变量列表)其中等号左边的变量为返回变量，等号右边的变量为输入变量。MATLAB允许在函数调用时同时返回多个变量。而一个函数又可以由多种格式进行调用，例如bode( )函数可以由下面的格式调用mag, phase = bode(num, den, w)其中bode( )函数用来求取或绘制系统的Bode图，而系统在这里由传递函数分子num和分母den表示，还可以用下面的格式调用此函数mag, phase = bode(A，B，C，D，w)其中(A，B，C，D)为系统的状态方程模型。尽管两种调用格式是完全不同的，MATLAB函数还是会自动识别到底是用哪种格式调用该函数的，从而得出正确的结论。在当前版本的控制系统工具箱中还允许用线性模型对象变量G来描述系统，例如mag, phase = bode(G, w)MATLAB函数在返回变量的格式上可以不同，例如若上面的语句中若没有返回变量，则将自动地绘制系统的Bode图，否则将返回计算结果数据。9.3.3 绘制响应曲线MATLAB具有丰富的获取图形输出的程序集。命令plot可以产生线性x-y图形（用命令loglog、semilogx、semilogy或 polar取代plot，可以产生对数坐标图和极坐标图）。所有这些命令的应用方式都是相同的，它们只对如何对坐标轴进行分度和如何显示数据产生影响。1. x-y图如果x和y是同一长度的向量，则命令plot(x, y)将画出y值对于x值的关系图。2. 画多条曲线为了在一幅图上画出多条曲线，采用具有多个自变量的plot命令：plot(X1, Y1, X2, Y2, , Xn, Yn)变量X1, Y1, X2, Y2等等是一些向量对。每一个x-y对都可以图解表示出来，因而在一幅图上形成多条曲线。多重变量的优点是它允许不同长度的向量在同一幅图上显示出来。每一对向量采用不同的线型。在一幅图上画一条以上的曲线时，也可以利用命令hold。hold命令可以保持当前的图形，并且防止删除和修改比例尺。因此，随后的一条曲线将会重叠地画在原曲线图上。再次输入命令hold，会使当前的图形复原。3. 加进网格线、图形标题、x轴标记和y轴标记一旦在屏幕上显示出图形，就可以画出网格线，定出图形标题，并且标定x轴标记和y轴标记。MATLAB中关于网格线、标题、x轴标记和y轴标记的命令如下：grid(网格线)title(图形标题)xlabel(x轴标记)ylabel(y轴标记)应当指出，一旦恢复命令display，通过依次输入相应的命令，就可以将网格线、图形标题、x轴标记和y轴标记叠加在图形上。4. 在图形屏幕上书写文本为了在图形屏幕的点(x,y)上书写文本，采用命令：text(X, Y, text)例如，利用语句text(3, 0.45, sin t)将从点(3, 0.45)开始，水平地写出sin t。另外，下列语句：plot(x1, y1, x2, y2), text(x1, y1, 1), text(x2, y2, 2)标记出两条曲线，使它们很容易地区分开来。5. 图形类型下列语句：plot(X, Y, x)将利用标记符号x画出一个点状图，而语句：plot(X1, Y1, :, X2, Y2, +)将用虚线画出第一曲线，用加法符号（+）画出第二条曲线。MATLAB能够提供的线和点的类型如下：线的类型点的类型实线-圆点.短划线-加号+虚线：星号*点划线-.圆圈o号6. 颜色下列语句plot(X, Y, r)plot(X, Y, +g)表明，第一幅图采用红线，第二幅图采用绿色“+”号标记。MATLAB提供的颜色如下：红色r绿色g蓝色b白色w无色i7. 自动绘图算法在MATLAB 中，图形是自动定标的。在另一幅图形画出之前，这幅图形作为现行图将保持不变，但是在另一幅图形画出后，原图形将被删除，坐标轴自动地重新定标。关于暂态响应曲线、根轨迹、伯德图、奈魁斯特图等的自动绘图算法已经设计出来，它们对于各类系统具有广泛的适用性，但是并非总是理想的。因此，在某些情况下，可能需要放弃绘图命令中的自动坐标轴定标特性，改用手工选择绘图范围。8. 手工坐标轴定标如果需要在下列语句指定的范围内绘制曲线：v = x-min x-max y-min y-max则应输入命令axis(v)，式中v是一个四元向量。axis(v)把坐标轴定标建立在规定的范围内。对于对数坐标图，v的元素应为最小值和最大值的常用对数。执行axis(v)会把当前的坐标轴定标保持到后面的图中，再次键入axis恢复自动定标。axis(square)把图形的范围设定在方形范围内。对于方形长宽比，斜率为1的直线恰位于45o上，它不会因屏幕的不规则形状而变形。axis(normal)将使长宽比恢复到正常状态。9.3.4 MATLAB语言的联机帮助功能本书中将广泛应用MATLAB语言和其控制系统工具箱作为辅助工具，并将介绍相关函数的使用。但是以当前的篇幅想全面介绍MATLAB和控制系统工具箱是不可能的，不妨使用MATLAB的联机帮助功能。MATLAB的联机帮助既可以由help命令来直接获得，又可以由MATLAB图形界面下的Help菜单来查询。例如，可以由help lyap命令得出lyap()函数的联机帮助信息如下：help lyapLYAP Lyapunov equation.X = LYAP(A,C) solves the special form of the Lyapunov matrix equation:A*X + X*A = -CX = LYAP(A,B,C) solves the general form of the Lyapunov matrix equation:A*X + X*B = -CSee also DLYAP.可见，由这种方法得出的帮助信息包括该函数的解释，函数的调用格式和相关函数名等等，进一步的帮助内容可以查阅MATLAB或相应工具箱手册。还可以用lookfor命令在MATLAB路径下查询有关的关键词，例如若想查询关键词Hankel，则可以由下面的命令完成：lookfor hankelHANKEL Hankel matrix.BHRDEMO Demo of model reduction techniques (Hankel, Balanced, BST).HKSV Hankel singular values and grammians P, Q.OHKAPP Optimal Hankel norm approximation (stable plant).OHKDEMO Demo of optimal Hankel mod