《非线性控制系统分析 (2)2精选文档.ppt》由会员分享，可在线阅读，更多相关《非线性控制系统分析 (2)2精选文档.ppt（44页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、非线性控制系统分析MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答本讲稿第一页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答主要内容n原理要点n非线性系统概述n相平面法n描述平面法本讲稿第二页，共四十四页n原理要点原理要点n非线性系统的研究方法由于系统的复杂性和多非线性系统的研究方法由于系统的复杂性和多样性而成为控制界的研究热点，从而产生了很样性而成为控制界的研究热点，从而产生了很多理论方法。比较基本的有李雅普诺夫第二法，多理论方法。比较基本的有李雅普诺夫第二法，小范围线性近似法，描述函数法，相平面法
2、，小范围线性近似法，描述函数法，相平面法，计算机仿真等等。计算机仿真等等。本讲稿第三页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答1.典型的非线性特性典型的非线性特性典型的非线性特性有死区非线性、饱和非线性、间隙非线性、继典型的非线性特性有死区非线性、饱和非线性、间隙非线性、继电非线性等。电非线性等。Simulink给出了部分非线性特性模块。用户也给出了部分非线性特性模块。用户也可以自行构建非线性特性模块。可以自行构建非线性特性模块。2.非线性控制系统非线性控制系统含有非线性元件或环节的控制系统称为非线性控制系统。含有非线性元件或环节的控
3、制系统称为非线性控制系统。非线性系统输出暂态响应曲线的形状与输入信号的大小和初始非线性系统输出暂态响应曲线的形状与输入信号的大小和初始状态有关，非线性系统的稳定性亦与输入信号的大小和初状态有关，非线性系统的稳定性亦与输入信号的大小和初始状态有关。非线性系统常会产生持续振荡。始状态有关。非线性系统常会产生持续振荡。本讲稿第四页，共四十四页3.描述函数法描述函数法非线性特性的描述函数法是线性部件频率特性在非线性特性中非线性特性的描述函数法是线性部件频率特性在非线性特性中的推广。它是对非线性特性在正弦信号作用下的输出进行的推广。它是对非线性特性在正弦信号作用下的输出进行谐波线性化处理之后得到的，是非
4、线性特性的一种近似描谐波线性化处理之后得到的，是非线性特性的一种近似描述。述。4.用描述函数研究系统的稳定点的方法用描述函数研究系统的稳定点的方法用描述函数研究系统的稳定点的方法，是建立在线性系统用描述函数研究系统的稳定点的方法，是建立在线性系统Nyquist稳定判据基础上的一种工程近似方法。其基本思想稳定判据基础上的一种工程近似方法。其基本思想是把非线性特性用描述函数来表示，将复平面上的整个非是把非线性特性用描述函数来表示，将复平面上的整个非线性曲线线性曲线()理解为线性系统分析中的临界点，再将线性系理解为线性系统分析中的临界点，再将线性系统有关稳定性分析的结论用于非线性系统。统有关稳定性分
5、析的结论用于非线性系统。本讲稿第五页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答16.1非线性系统概述n含有非线性元件或环节的控制系统称为非线性控制系含有非线性元件或环节的控制系统称为非线性控制系统。统。n一般非线性系统的数学模型可表示为一般非线性系统的数学模型可表示为：写成多变量的形式为：写成多变量的形式为：本讲稿第六页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n在在F与与f函数中，如果相应的算子为线性，函数中，如果相应的算子为线性，则称为线性系统，否则称为非线性系统。则称为线性
6、系统，否则称为非线性系统。如果不显含如果不显含t，则为时不变系统，若显含，则为时不变系统，若显含t，则称为时变系统。，则称为时变系统。n非线性系统输出暂态响应曲线的形状与输非线性系统输出暂态响应曲线的形状与输入信号的大小和初始状态有关，非线性系入信号的大小和初始状态有关，非线性系统的稳定性亦与输入信号的大小和初始状统的稳定性亦与输入信号的大小和初始状态有关。非线性系统常会产生持续振荡。态有关。非线性系统常会产生持续振荡。本讲稿第七页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答16.2 非线性特性模块的构建及示例n典型的非线性特性有死区非线
7、性、饱和非线性、典型的非线性特性有死区非线性、饱和非线性、间隙非线性、继电非线性等。间隙非线性、继电非线性等。Simulink给出了给出了部分非线性特性模块。这在部分非线性特性模块。这在Simulink一章中已一章中已列出。在系统仿真中可以直接使用。但对于没有列出。在系统仿真中可以直接使用。但对于没有提供的模块则需要我们自己构建。那么如何根据提供的模块则需要我们自己构建。那么如何根据需要构建任意的非线性模块呢？事实上，任意的需要构建任意的非线性模块呢？事实上，任意的静态非线性模块，无论其是单值非线性，还是多静态非线性模块，无论其是单值非线性，还是多值非线性，都可以由值非线性，都可以由Simul
8、ink构建，并直接用构建，并直接用于仿真。于仿真。本讲稿第八页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n例例1：构建如图：构建如图16.1分段线性的非线性特性模块。分段线性的非线性特性模块。图图16.1 例例1非线性特性非线性特性本讲稿第九页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答1.新建一个空白模型。在模型中添加子模块集新建一个空白模型。在模型中添加子模块集Lookup Tables中的中的Lookup Table模块。模块。2.设置模块属性。双击设置模块属性。双击Look
9、up Table模块进入其属性设置窗口，模块进入其属性设置窗口，如图如图16.2，并添加非线性特性值。其中，并添加非线性特性值。其中，Vector of input values栏为横坐标向量，而栏为横坐标向量，而Table data栏为纵坐标向量。需要栏为纵坐标向量。需要注意的是，如果仅添加了图中的所有转折点坐标，则位于最左侧注意的是，如果仅添加了图中的所有转折点坐标，则位于最左侧与最右侧外边的特性将无法表现。因此还应该在特性曲线的两侧与最右侧外边的特性将无法表现。因此还应该在特性曲线的两侧再找两点，从而完整地表现非线性特性。根据非线性函数，位于再找两点，从而完整地表现非线性特性。根据非线性
10、函数，位于最左侧转折点最左侧转折点(-3，-1)之外的点取为之外的点取为(-4，-2)，位于最右侧转，位于最右侧转折点折点(4，1)之外的点取为之外的点取为(5，2)。本讲稿第十页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图图16.2 非线性特性属性设置窗口非线性特性属性设置窗口本讲稿第十一页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 例2：构建如图16.3的回环非线性特性模块。图图16.3 例例2非线性特性非线性特性本讲稿第十二页，共四十四页MATLAB与控制系统仿真实践，北京
11、航空航天大学出版社，2009.8.在线交流，有问必答n分析：该特性在输入信号增加时走一条折线，而在输入信号减小时走分析：该特性在输入信号增加时走一条折线，而在输入信号减小时走另一条折线。可以将特性分解为两个单值函数。如图另一条折线。可以将特性分解为两个单值函数。如图16.4。n根据例根据例1的结果，这两个单值函数都可以用查表模块实现。这里有两个问题的结果，这两个单值函数都可以用查表模块实现。这里有两个问题需要解决。一是如何判断输入是增加还是减小？在判断输入信号是否为增加需要解决。一是如何判断输入是增加还是减小？在判断输入信号是否为增加时，可通过比较输入信号的当前值和它的上一步值进行判断。而时，
12、可通过比较输入信号的当前值和它的上一步值进行判断。而Simulink离离散模块组中提供的散模块组中提供的Memory模块，可以用来记忆上一个计算步长的信号值，模块，可以用来记忆上一个计算步长的信号值，这样将输入信号的当前值和它的上一步值分别作为比较模块这样将输入信号的当前值和它的上一步值分别作为比较模块(Relational Operator)的输入，即可输出代表上升还是下降的逻辑值的输入，即可输出代表上升还是下降的逻辑值1 和和0。二是如何。二是如何控制特性曲线走不同折线？控制特性曲线走不同折线？Simulink中的中的Signal Routing子模块组中子模块组中Switch模块，使用比
13、较模块的输出作为输入控制，即可使模块对输入信号的模块，使用比较模块的输出作为输入控制，即可使模块对输入信号的不同变化走不同折线。具体实现如图不同变化走不同折线。具体实现如图16.5：本讲稿第十三页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答(a)输入上升分支输入上升分支 (b)输入下降分支输入下降分支图图16.4 特性分解后的两个单值函数特性分解后的两个单值函数本讲稿第十四页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图图16.5 例例2非线性特性的非线性特性的Simulink模型模
14、型本讲稿第十五页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答1.新建一个空白模型。添加所需的不同模块。新建一个空白模型。添加所需的不同模块。2.进行不同模块的连接并进行属性设置。进行不同模块的连接并进行属性设置。图图16.5中，输入上升分支和输入下降分支都是调用了查表模块。其设置中，输入上升分支和输入下降分支都是调用了查表模块。其设置见图见图16.6。(a)输入上升分支设置窗口输入上升分支设置窗口 本讲稿第十六页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答(b)输入下降分支设置窗口输
15、入下降分支设置窗口图图16.6 例例2非线性特性设置窗口非线性特性设置窗口本讲稿第十七页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答对输入信号当前值和其上一步的值比较，如果当前值大于等于前一步值，则对输入信号当前值和其上一步的值比较，如果当前值大于等于前一步值，则模块表现为上升分支的特性；反之，则表现为下降分支的特性。模块表现为上升分支的特性；反之，则表现为下降分支的特性。Relational Operator(比较模块比较模块)默认值为默认值为=。Switch(开关模块开关模块)的控制阀值的控制阀值(Threshold)可以设置。这里设
16、为可以设置。这里设为0.5。即控制端输入。即控制端输入=0.5时，按上升分支特性输出，否则按上时，按上升分支特性输出，否则按上降分支特性输出。降分支特性输出。3.给定输入，观察非线性模块的特性。给定输入，观察非线性模块的特性。本讲稿第十八页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图图16.7 例例2非线性特性在正弦输入的非线性特性在正弦输入的Simulink模型模型本讲稿第十九页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n本例给定输入为正弦信号，其幅值分别设为2、4，其输出可以
17、用示波器模块直接观察，也可以输出到工作空间后，使用plot函数绘制。其Simulink模型如图16.7。本例输出到工作空间变量名设为simout，其保存格式设为Array，在命令窗口使用plot函数绘制，运行结果如图。plot(tout,simout(:,1),tout,simout(:,2)本讲稿第二十页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答(a)输入信号幅值为输入信号幅值为2时的仿真输出时的仿真输出 (b)输入信号幅值为输入信号幅值为4时的仿真输出时的仿真输出 图图16.8 例例2不同输入下的仿真输出不同输入下的仿真输出本讲稿第
18、二十一页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答16.3相平面法n16.3.1 相平面法概述相平面法概述n相平面法是一种求解二阶以下线性或非线性微分方程的图解方法。n对于形如下式的二阶系统本讲稿第二十二页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n涉及的概念有：涉及的概念有：n1.相平面：以为横坐标，为纵坐标的直角坐标平面构成相平相平面：以为横坐标，为纵坐标的直角坐标平面构成相平面。面。n2.相轨迹：以时间为参变量，由表示运动状态的分别作为横相轨迹：以时间为参变量，由表示运动状
19、态的分别作为横坐标和纵坐标而绘制的曲线称为相轨迹，每根相轨迹与起始坐标和纵坐标而绘制的曲线称为相轨迹，每根相轨迹与起始条件有关。表示了质点在时刻的位置和速度。条件有关。表示了质点在时刻的位置和速度。n3.相平面图：同一系统，不同初始条件下的相轨迹是不同的。相平面图：同一系统，不同初始条件下的相轨迹是不同的。由所有相轨迹组成的曲线族所构成的图称为相平面图。由所有相轨迹组成的曲线族所构成的图称为相平面图。本讲稿第二十三页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答16.3.2 基于MATLAB的相轨迹图绘制实例n例例3：绘制如下系统的单位阶跃
20、输入时的相轨迹。其中，非：绘制如下系统的单位阶跃输入时的相轨迹。其中，非线性部分为饱和非线性，线性部分为饱和非线性，线性部分为线性部分为系统初始状态为系统初始状态为0。本讲稿第二十四页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n1.新建一个空白模型。将所需的不同模块添加到空白模型中。新建一个空白模型。将所需的不同模块添加到空白模型中。n2.连接各模块并设置各模块参数。这里将饱和非线性模块连接各模块并设置各模块参数。这里将饱和非线性模块upper limit设为设为0.3，lower limit设为设为-0.3。其它模块的设置不再赘述，模
21、型如图。其它模块的设置不再赘述，模型如图16.9。图图16.9 例例3的的Simulink模型模型 本讲稿第二十五页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n3.设置仿真参数。如图设置仿真参数。如图16.10，将，将Solver options下的下的Type项选为项选为Fixed-step，Solver项选项选ode5(Dormand-Prince)，Fixed-step size 设为设为0.01。图图16.10 仿真参数设置窗口仿真参数设置窗口本讲稿第二十六页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2
22、009.8.在线交流，有问必答n4.开始仿真。开始仿真。n相轨迹可以直接观察相轨迹可以直接观察XYGraph输出，也可使用输出到工作空间的参数绘输出，也可使用输出到工作空间的参数绘制，如图制，如图16.11所示。所示。plot(simout(:,1),simout1(:,1)grid图图16.11 例例3输出的相轨迹输出的相轨迹本讲稿第二十七页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n系统阶跃响应输出如图系统阶跃响应输出如图16.12所示。所示。图图16.12 系统阶跃响应输出系统阶跃响应输出由由16.11分析可知，系统的稳定点在分析
23、可知，系统的稳定点在(1,0)点，即稳态值为点，即稳态值为1。本讲稿第二十八页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答16.4描述函数法16.4.1描述函数法概述描述函数法概述P.J.Daniel于于1940年首先提出了描述函数法。非线性特性年首先提出了描述函数法。非线性特性的描述函数法是线性部件频率特性在非线性特性中的推的描述函数法是线性部件频率特性在非线性特性中的推广。它是对非线性特性在正弦信号作用下的输出进行谐广。它是对非线性特性在正弦信号作用下的输出进行谐波线性化处理之后得到的，是非线性特性的一种近似描波线性化处理之后得到的，
24、是非线性特性的一种近似描述。述。本讲稿第二十九页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答1.描述函数法的定义：描述函数法的定义：设非线性环节的输入输出关系为设非线性环节的输入输出关系为非线性环节输入正弦信号非线性环节输入正弦信号非线性环节的输出通常也为周期信号，可以分解为傅立叶级数非线性环节的输出通常也为周期信号，可以分解为傅立叶级数 其中，为直流分量，和是第其中，为直流分量，和是第n次谐波的幅值和相角，且次谐波的幅值和相角，且有有本讲稿第三十页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流
25、，有问必答若若，且，且 时时 很小，则非线性环节的输出近似为很小，则非线性环节的输出近似为本讲稿第三十一页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n可见，其近似结果和线性环节频率响应形式相似，依照线可见，其近似结果和线性环节频率响应形式相似，依照线性环节的频率特性的定义，非线性环节的输入输出特性可性环节的频率特性的定义，非线性环节的输入输出特性可由描述函数表示：由描述函数表示：对于非线性控制系统的描述函数分析方法，常用的负倒描述函数对于非线性控制系统的描述函数分析方法，常用的负倒描述函数为：为：本讲稿第三十二页，共四十四页MATLAB
26、与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答对于如图对于如图16.13 的等效非线性系统，且的等效非线性系统，且 在开环幅相平面上无右半平面的极点，稳定性判据为：在开环幅相平面上无右半平面的极点，稳定性判据为：如果如果 不被不被 包围，则系统是稳定的，如果包围，则系统是稳定的，如果 被被 包围，则系统是不稳定的系统。包围，则系统是不稳定的系统。包围的区域称为不稳定区域，不包围的区域称为稳定区域。包围的区域称为不稳定区域，不包围的区域称为稳定区域。如果如果 与与，则在交点处，若，则在交点处，若 沿着沿着 A 值增加的方向由不稳定区域进入稳定区域，则自激振荡是稳定的
27、，否值增加的方向由不稳定区域进入稳定区域，则自激振荡是稳定的，否则，自激振荡是不稳定的。则，自激振荡是不稳定的。在交点处有：在交点处有：本讲稿第三十三页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n2.用描述函数研究系统的稳定点的方法用描述函数研究系统的稳定点的方法n用描述函数研究系统的稳定点的方法，是建立在线性系统用描述函数研究系统的稳定点的方法，是建立在线性系统Nyquist稳定判据基础上的一种工程近似方法。其基本思想是稳定判据基础上的一种工程近似方法。其基本思想是把非线性特性用描述函数来表示，将复平面上的整个非线性把非线性特性用描述
28、函数来表示，将复平面上的整个非线性曲线曲线统有关稳定性分析的结论用于非线性系统。统有关稳定性分析的结论用于非线性系统。理解为线系理解为线系统分析中的临界点统分析中的临界点，再将线性系再将线性系本讲稿第三十四页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图图16.13 等效非线性系统等效非线性系统本讲稿第三十五页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答对于如图对于如图16.13 的等效非线性系统，且的等效非线性系统，且 在开环幅相平面上无右半平面的极点，稳定性判据为：在开环幅相平面
29、上无右半平面的极点，稳定性判据为：如果如果 不被不被 包围，则系统是稳定的，如果包围，则系统是稳定的，如果 被被 包围，则系统是不稳定的系统。包围，则系统是不稳定的系统。包围的区域称为不稳定区域，不包围的区域称为稳定区域。包围的区域称为不稳定区域，不包围的区域称为稳定区域。如果如果 与与，则在交点处，若，则在交点处，若 沿着沿着 A 值增加的方向由不稳定区域进入稳定区域，则自激振荡是稳定的，否则，值增加的方向由不稳定区域进入稳定区域，则自激振荡是稳定的，否则，自激振荡是不稳定的。自激振荡是不稳定的。在交点处有：在交点处有：由此可求出自激振荡的振幅由此可求出自激振荡的振幅 和振荡频率和振荡频率。
30、本讲稿第三十六页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答16.4.2 基于MATLAB的描述函数法 非线性系统分析实例 n例例4：考虑如图：考虑如图16.14的非线性系统，图的非线性系统，图中的继电器非线性模块。试判断系统是否中的继电器非线性模块。试判断系统是否存在自振；若有自振，求出自振的振幅和存在自振；若有自振，求出自振的振幅和频率。频率。本讲稿第三十七页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图图16.14 例例4系统框图系统框图本讲稿第三十八页，共四十四页MATLAB
31、与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n1.绘制非线性部分和线性部分的幅相图，判断系统稳绘制非线性部分和线性部分的幅相图，判断系统稳定情况定情况n程序如下：程序如下：x=1:0.1:20;disN=40/pi./x.*sqrt(1-x.(-2)-j*40/pi./x.2;%描述函数disN2=-1./disN;%负倒描述函数w=1:0.01:200;num=12;%线性部分分子den=conv(1 1,1 6 13);%线性部分分母rem,img,w=Nyquist(num,den,w);%线性部分Nyquist曲线参数plot(real(disN2),i
32、mag(disN2),rem,img)%同时绘制非线性部分和线性部分的极坐标图grid;%加网格本讲稿第三十九页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图图16.15 程序运行结果图程序运行结果图 本讲稿第四十页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 图图16.16 程序运行结果局部放大图程序运行结果局部放大图由图由图16.15可见，两曲线相交，系统存在自激振荡。可见，两曲线相交，系统存在自激振荡。本讲稿第四十一页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学
33、出版社，2009.8.在线交流，有问必答n2.利用交点坐标值求取振荡幅值和频率利用交点坐标值求取振荡幅值和频率n%读出线性部分和非线性部分交点的坐标值，并利用坐标值求出振荡幅值和频率nw0=spline(img,w,-0.0785)%当img=-0.0785时，所对应的w值nx0=spline(real(disN2),x,-0.166)%当disN2的实部为-0.166时，所对应的x值n由图16.15可见，两条曲线有交点，存在自激振荡。经局部放大，如图16.16，可得到交点坐标为(-0.0785，-0.166)。nw0=n 3.2087nx0=n2.3382n则系统中有的自激振荡。本讲稿第四十
34、二页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答n3.建立建立Simulink模型，如图模型，如图16.17，进行仿真。，进行仿真。图图16.17 例例4系统的系统的Simulink仿真模型仿真模型本讲稿第四十三页，共四十四页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图图16.18 系统的系统的Simulink仿真输出结果仿真输出结果由图由图16.18所示的仿真输出可见，系统中确实存在自所示的仿真输出可见，系统中确实存在自激振荡，进一步证实了前面的分析。激振荡，进一步证实了前面的分析。本讲稿第四十四页，共四十四页