《《非线性系统分析》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《非线性系统分析》PPT课件.ppt（103页珍藏版）》请在金锄头文库上搜索。
1、第第7 7章章 非线性系统的分析非线性系统的分析7.1 7.1 非线性系统基本概念非线性系统基本概念7.2 7.2 二阶线性和非线性系统的相平面分析二阶线性和非线性系统的相平面分析7.3 7.3 线性系统的相轨迹线性系统的相轨迹7.4 7.4 非线性系统的相平面分析非线性系统的相平面分析7.5 7.5 描述函数法描述函数法7.7.6 6 分析分析非线性系统的谐波平衡分析法非线性系统的谐波平衡分析法7.7.7 7 非线性系统的非线性系统的计算机仿真计算机仿真小小 结结Nonlinear Science 研究什么？n客观世界是非线性的、非平衡的复杂世界nPrigogine：“自然用一千种声音说话，
2、我们刚开始去听”。n自古希腊：人们笃信和向往世界的稳定性、规则性、和谐性、有序性、因果性、稳定性、规则性、和谐性、有序性、因果性、本质简单性、周期性、对称性、本质简单性、周期性、对称性、n现在：人们越来越认识到：我们所处的大千世界是以不稳定动力系统不稳定动力系统为特征的，充满了：非平衡、非线性、非平衡、非线性、非稳定、非均匀、非结构、非确定、非可积、非可逆、非晶非稳定、非均匀、非结构、非确定、非可积、非可逆、非晶态、非规则、非连续、非光滑、非周期、非对称、非标准分态、非规则、非连续、非光滑、非周期、非对称、非标准分析、非析、非von Neumannvon Neumann计算机、计算机、 人类理
3、智夸入“想入非非想入非非”时代非线性科学的四个发展阶段非线性科学的四个发展阶段n40年代：组织理论：年代：组织理论：控制论，信息论，一般系统论n60年代：自组织理论年代：自组织理论（系统如何从无序有序）：Catastrophic Theory （Thom, Arnold)，超循环论（Eigen），Dissipative Structure（Prigogine），Synergetics （Haken）n70年代：非线性科学年代：非线性科学 （系统如何从有序 混沌和无序 更高层次的有序）Chaotic Dynamics（Feigenbaum， Ford， Kadanoff）， Integrable
4、 SystemSoliton Theory（Scott，扎哈罗夫）， Fractals （Mandelbrot）n90年代：复杂性科学年代：复杂性科学（复杂性的定义及量度，复杂系统的行为及模型）Neural Network （Hoppfield）， Cellular Automaton （Wolfram），人工生命什么是非线性系统？什么是非线性系统？n平衡与非平衡平衡与非平衡：物理概念物理概念线性与非线性：线性与非线性：数学概念数学概念n线性系统线性系统：整体的行为或性质是部分之和整体的行为或性质是部分之和1 1 复杂性不因叠加产生复杂性不因叠加产生2 2 只要知道初始条件，即可了解过去，预测
5、未来只要知道初始条件，即可了解过去，预测未来n非线性系统：非线性系统： 叠加原理失效叠加原理失效整体的行为和性质整体的行为和性质各部分的行为与性质（本质区别）各部分的行为与性质（本质区别）系统行为对初始条件极端敏感依赖系统行为对初始条件极端敏感依赖ChaosChaos长期行为不可预测，决定论性混沌、内在随机性长期行为不可预测，决定论性混沌、内在随机性非非线线性性系系统统：只只要要系系统统中中包包含含一一个个或或一一个个以以上上具具有有非非线线性性静静特特性性的的元元件件，即即称称为为非非线线性性系系统统(P267)(P267)。线性控制系统：线性控制系统：由线性元件组成，输入输出具有叠加性和均
6、匀性由线性元件组成，输入输出具有叠加性和均匀性由线性元件组成，输入输出具有叠加性和均匀性由线性元件组成，输入输出具有叠加性和均匀性. . . .非线性控制系统：非线性控制系统： 系统中有非线性元件，输入输出不具有叠加性和均匀性系统中有非线性元件，输入输出不具有叠加性和均匀性系统中有非线性元件，输入输出不具有叠加性和均匀性系统中有非线性元件，输入输出不具有叠加性和均匀性. . . . 非本质非线性非本质非线性: :能够用小偏差线性化方法进行线性化处理的非线性能够用小偏差线性化方法进行线性化处理的非线性能够用小偏差线性化方法进行线性化处理的非线性能够用小偏差线性化方法进行线性化处理的非线性. .
7、. .本质非线性：本质非线性：用小偏差线性化方法不能解决的非线性用小偏差线性化方法不能解决的非线性用小偏差线性化方法不能解决的非线性用小偏差线性化方法不能解决的非线性. . . .例：例：对于一由非线性微分方程对于一由非线性微分方程x=-x(1x)描述的非线性系统，显然有两个平衡点，即描述的非线性系统，显然有两个平衡点，即x1=0和和x2=1。将上式改写为将上式改写为非线性系统的稳定性不仅取决于控制系统的固非线性系统的稳定性不仅取决于控制系统的固有结构和参数，而且与系统的初始条件以及外有结构和参数，而且与系统的初始条件以及外加输入有关系。加输入有关系。设设t0时，系统的初态为时，系统的初态为x
8、0。积分上式可得积分上式可得x(t)t10图图1 1 一阶非线性系统一阶非线性系统系统的稳定性除与结构参数有关外，还与起始偏系统的稳定性除与结构参数有关外，还与起始偏差的大小有关差的大小有关。统的响应形式与输入信号的大小和初始条件有关。统的响应形式与输入信号的大小和初始条件有关。在没有外界周期变化信号输入时，非线性系统完全在没有外界周期变化信号输入时，非线性系统完全可能产生具有固定周期和幅值的稳定振荡过程。可能产生具有固定周期和幅值的稳定振荡过程。(自自激振荡激振荡P272)非线性系统的主要特征：非线性系统的主要特征：1.研究非线性系统的意义研究非线性系统的意义1）实实际际的的控控制制系系统统
9、，存存在在着着大大量量的的非非线线性性因因素素。这这些些非非线线性性因因素素的的存存在在，使使得得我我们们用用线线性性系系统统理理论论进进行行分分析析时时所所得得出出的的结结论论，与与实实际际系系统统的的控控制制效效果果不不一一致致。线线性性系系统理论无法解释非线性因素所产生的影响统理论无法解释非线性因素所产生的影响.2）非线性特性的存在，并不总是对系统产生不良影响非线性特性的存在，并不总是对系统产生不良影响.2.研究非线性系统的方法研究非线性系统的方法 1）相平面法）相平面法是用图解的方法分析一阶，二阶非线是用图解的方法分析一阶，二阶非线性系统的方法。通过绘制控制系统相轨迹，达到分性系统的方
10、法。通过绘制控制系统相轨迹，达到分析非线性系统特性的方法。析非线性系统特性的方法。 2）描述函数法描述函数法是近似分析法是近似分析法, ,受线性系统频率法受线性系统频率法启发，而发展出的一种分析非线性系统的方法。它启发，而发展出的一种分析非线性系统的方法。它是一种谐波线性化的分析方法，是频率法在非线性是一种谐波线性化的分析方法，是频率法在非线性系统分析中的推广。系统分析中的推广。3）计算机求解法）计算机求解法是利用计算机运算能力和高速度是利用计算机运算能力和高速度对非线性微分方程的一种数值解法。对非线性微分方程的一种数值解法。 非线性系统与线性系统的区别，相平面的基本概念，相轨迹，极限环，描述
11、函数的基本思想，描述函数的定义和求取，描述函数法分析非线性系统的自持振荡，非线性系统的计算机仿真。 知 识 要 点 静态非线性特性中，死区特性、饱和特性、继电特性、间隙特性是最常见的，也最简单。一、数学描述 一个单输入单输出静态非线性特性的数学描述为：7.1 7.1 非线性系统基本概念非线性系统基本概念二 非线性特性分类：(P269)1 1、死区特性、死区特性常常是由放大器、传感器、执行机构的不灵敏区造成的。实际的死区特性一般如图7-1中的点划线所示，为了分析的方便，我们将它用图7-1中的三段直线（实线）来近似，并称之为理想死区特性。理想型死区特性的的数学描述为： 图 71 死区特性(7-2)
12、 死区特性可能给控制系统带来不利影响，它会使控制的灵敏度下降，稳态误差加大；死区特性也可能给控制系统带来有利的影响，有些系统人为引入死区以提高抗干扰能力。2 2、饱和特性、饱和特性 可以说，任何实际装置都存在饱和特性，因为它们的输出不可能无限增大，磁饱和就是一种饱和特性。实际的饱和特性一般如图7-2中的点划线所示，为了分析的方便，我们将它用图7-2 中的三段直线来近似，并称之为理想饱和特性。 理想饱和特性的数学描述为： (7-2)图7-2 饱和特性 继电特性顾名思义就是继电器所具有的特性, 继电特性有双位特性如图7-3（a）和（b），三位特性如图7-3（c）等，图7-3（b）（c）的继电特性还
13、带有滞环。当然，不限于继电器，其它装置如果具有类似的非线性特性，我们也称之为继电特性，比如：电磁阀、斯密特触发器等。 分析继电特性有十分重要的意义，因为采用继电器、电磁阀等元件的的控制系统比比皆是，例如大多数家用电冰箱、空调就是继电器控制系统。 3 3、继电特性、继电特性图7-3 几种典型的继电特性图7-3（a）所示继电特性的数学描述为：图（c）所示继电特性的数学描述为： 图（b）所示继电特性的数学描述由读者自行导出。 传传动动机机构构的的间间隙隙也也是是控控制制系系统统中中常常见见的的非非线线性性特特性性，齿齿轮轮传传动动是是典典型型的的间间隙隙特特性性，图图7-4（a）表表示示齿齿轮轮传传
14、动动原原理理，图图7-4（b）表表示示主主动动轮轮位位移移与与从从动动轮轮位位移移的的关关系系。设设主主动动轮轮与与从从动动轮轮间间的的最最大大间间隙隙为为2b，那那么么当当主主动动轮轮改改变变方方向向时时，主主动动轮轮最最大大要要运运动动2b从从动动轮轮才才能能跟跟随随运运动动。间间隙隙特特性性类类似似于于线线性性系系统统的的滞滞后后环环节节，但但不不完完全全等等价价，它它对对控控制制系系统统的的动动态态、稳稳态态特特性性都都不不利利。设设齿齿轮轮传传动动速比为，则图速比为，则图7-4间隙特性的数学描述为：间隙特性的数学描述为：4 4、间隙特性、间隙特性图7-4 间隙特性返回式中，为常数，它
15、等于主动轮改变方向时的值。 相相平平面面法法是是庞庞加加莱莱（PoincarePoincare）18851885年年首首先先提提出出的的，本本来来它它是是一一种种求求解解二二元元一一阶阶非非线线性性微微分分方方程程组组的的图图解解法法, ,两两个个变变量量构构成成的的直直角角坐坐标标系系称称为为相相平平面面，方方程程组组的的解解在在相相平平面面上上的的图图象象称称为为相相轨轨迹。迹。 这里是将相平面法用于分析一阶尤其是二阶这里是将相平面法用于分析一阶尤其是二阶非线性控制系统，并形成了一种特定的相平面法，非线性控制系统，并形成了一种特定的相平面法，它对弄清非线性系统的稳定性、稳定域等基本属它对弄
16、清非线性系统的稳定性、稳定域等基本属性性, ,解释极限环等特殊现象，起到了直观形象的作解释极限环等特殊现象，起到了直观形象的作用。用。 7.2 7.2 二阶系统的相平面分析法二阶系统的相平面分析法 因为绘制两维以上的相轨迹是十分困难的，所以相平面法对于二阶以上的系统几乎无能为力。但是，如果我们将相平面概念推广到到抽象空间，就得到n维状态空间以后再专门介绍。下面讨论相平面和相轨迹的基本概念。 考察二阶非线性时不变微分方程:7.2.1 7.2.1 相平面的基本概念相平面的基本概念 为了引入相平面法，将二阶微分方程改写成二元一阶微分方程组: 微分方程组(7-6)有两个变量: x可以看作广义位移， 可
17、以看作广义速度。 一般，直接对微分方程（7-5）求解，可以得到该系统的时间解 x(t),还可以作出x(t)与t的关系图时间响应曲线。 如果我们对微分方程组(7-6)求解，可以得到解x(t)和 ,如果我们取 x 和 为坐标，以时间 t 为参变量,则系统的每一时刻的状态均对应于该平面上的一点，此平面即为相平面。当 t 变化时，这一点在 平面上绘出的曲线，表征了系统的运动过程，这个曲线就是相轨迹。我们用一个二阶线性时不变系统来体验一下相平面和相轨迹。 例7-1 考虑二阶系统：将它写成微分方程组：两式相除得到：即：两边积分得: 在相平面上绘出的相轨迹如图7-5(a)所示椭圆，如果取遍所有的初始值，就会
18、得到无数一环套一环的椭圆称为相轨迹场，相轨迹场布满了整个相平面，相轨迹场从全局上展示了动态系统的运动过程，图（a）只绘出了相轨迹场中的2根相轨迹。当 xo=0 时，响应曲线如图（b）。图7-5 例7-1的相轨迹与时间响应 7.2.2 7.2.2 相轨迹图相轨迹图绘制相轨迹图有多种办法，概括起来有如下几类：第一类：手工绘制概略图。概略图就象相轨迹的素描，它是根据相轨迹的基本特征、特殊点、特殊线等信息而随手画出的草图，它虽然在具体细节上缺乏精度，但却能提供许多重要的定性结论。第二类：手工图解绘制近似图。在计算机未得到广泛应用的年代，人们研究出好几种手工近似作图法，如等倾线法、法等。这些手工作图法要
19、绘出有一定精度的相轨迹图是十分繁琐的，如今已没有多大实用价值。第三类：计算机绘制精确图。借助计算机数值解法以及SIMULINK等软件绘制相轨迹图。 相轨迹的基本特征有：1 1）奇点）奇点 对于二阶系统，相平面上满足 且 的点叫做奇点，记作 。对照方程（7-6）知，奇点座标 是代数方程 的解，显然奇点一定在轴上。 对于二阶系统，对于二阶系统， 和和 就是就是速度速度和和加速度加速度均为零，也就意味着不均为零，也就意味着不再运动，所以，奇点又称再运动，所以，奇点又称平衡点平衡点。相平面上任。相平面上任何其它点，都叫普通点。奇点又分何其它点，都叫普通点。奇点又分稳定奇点稳定奇点和和不稳定奇点不稳定奇
20、点，稍后将讨论。，稍后将讨论。 2 2）相轨迹切线斜率）相轨迹切线斜率 由方程（7-6）知，相轨迹上任一点的切线斜率为：某点的切线斜率就是相轨迹通过该点的运动方向，前面提到的等倾线就是相轨迹场上所有切线斜率等于某一常数的点的连线。 3 3）相轨迹图形特征）相轨迹图形特征 如果微分方程（7-6）满足解的存在性和唯一性条件，那么，相轨迹（场）图一定有如下基本特征： 1)任一普通点有且只有一条相轨迹通过（解的存在性和唯一性）； 2)相轨迹必垂直通过横轴。 3)横轴上方的相轨迹从左向右运动，横轴下方的相轨迹从右向左运动。例7-2作出下列二阶系统的相轨迹 将它写成微分方程组：容易求出奇点为（0，0）。图
21、7-6 例7-2的根轨迹ABCDO对应初始条件为EFO对应初始条件为。 从相轨迹图可以直观地看到：所有的相轨迹都最终收敛到奇点（0，0），这说明系统是渐近稳定的；可以证明，每一条相轨迹都是向心螺旋线，这说明系统的运动过程是衰减振荡的。 返回 研究二阶线性系统相轨迹的意义主要在两个研究二阶线性系统相轨迹的意义主要在两个方面：方面：一是一是许多非线性特性可以近似为分段线性许多非线性特性可以近似为分段线性的，如死区特性、饱和特性、继电特性等，而分的，如死区特性、饱和特性、继电特性等，而分段线性系统的相轨迹可以由几段线性系统相轨迹段线性系统的相轨迹可以由几段线性系统相轨迹连接而成；连接而成；二是二是大
22、多数非线性系统在奇点附近的大多数非线性系统在奇点附近的相轨迹，与其在奇点附近的线性化系统的相轨迹相轨迹，与其在奇点附近的线性化系统的相轨迹十分接近。十分接近。 7.3 7.3 线性系统的相轨迹线性系统的相轨迹二阶线性系统的微分方程是：令：则可写成：即得：再可写为：由此式解出x1与x2的关系，即为二阶系统相轨迹 的方程。7.3.1 7.3.1 二阶线性系统的相轨迹二阶线性系统的相轨迹另外，系统的特征方程为： 极点为：下面由系统闭环极点的位置分析系统相轨迹性质：1.：共轭虚根。相轨迹方程：可见其为一族同心的椭圆，坐标原点有一孤立奇点(中心点)，每一椭圆对应一个简谐运动。(图7-2-2 a) 2.
23、为一对负实部共轭复根。相轨迹是收敛于相平面原点的对数螺线,奇点为稳定的焦点.(图7-2-2b) 均为负实根。相轨迹是一族趋向相平面原点的抛物线,奇点为稳定的节点.（图7-2-2 c）3. 4. 当 为实根，且符号相反（一左一右）时，奇点为鞍点。相轨迹如（图7-2-2d） 为一对正实部共轭复根。相轨迹是从相平面原点出发的对数螺线,奇点为不稳定的焦点.(图7-2-2e） 均为正实根。相轨迹是由相平面原点出发的发散型抛物线族,奇点为不稳定的节点.（图7-2-2f）5. 6. 图7-2-2 二阶线性系统的特征根与奇点(P275)n由上可见，二阶线性系统的相轨迹和奇点的性质由系统的特征根决定，与初始状态
24、无关。初始状态不同，相轨迹性质不变，形状相似，且对应于不同初始状态的相轨迹不会相交，只可能部分重合，而在奇点处则相交。另外，由于相轨迹的性质与系统初始状态无关，相平面中局部范围内相轨迹的性质就有决定性意义，从局部范围内相轨迹的性质可推知全局。 显然，当取不同值时，特征根在根平面上的分布也不同，响应曲线和相轨迹的形态也不同。（见P277表） 返回7.4 7.4 非线性系统的相平面分析非线性系统的相平面分析 非线性系统相平面分析的关键是绘出相轨迹图，有了相轨迹图，我们就可以得到系统稳定性、稳定域、振型、稳态误差等方面的结论。我们也可以绘出不同初始条件、不同输入、不同系统参数所对应的相轨迹图，研究其
25、中的规律。解析法图解法等倾线法法绘制方法 为了介绍相平面法基本原理，本节主要讨论一类最为简单的动态非线性系统。系统的结构如图7-8 所示，其显著特点是：系统具有静态非线性环节和动态线性环节的分离结构，且静态非线性环节是分段线性的，动态线性环节为一阶或二阶。图7-8 具有分离结构的非线性系统 分析图7-8 所示非线性系统的具体方法是：先将静态非线性特性分为几个线性段，划分出每段对应的相平面分区；然后在每个分区按线性系统绘出相轨迹；最后将各分区的相轨迹进行衔接就得到整个非线性系统的相轨迹。 （一）、解析法：（一）、解析法：首先要将相变量方程改写为：然后对其积分，得到x1与x2的关系式，即x2 =
26、f(x1)，这就是相轨迹方程，可按其绘图.讲解例7-3-1(P280).（二）、等倾线法（二）、等倾线法式表示了相轨迹的斜率，对于相平面上任意一点(x1,x2),为相轨迹通过该点时的切线的斜率。若取定q为常数，则可在相平面上找到一些点，在它们上面相轨迹的斜率都为q。可在它们上面作斜率均为q切线等倾线。给定不同的q，可在相平面上画出许多等倾线。等倾线描述了相轨迹运动的方向。当给定了初始状态后，便可沿着给定的相轨迹切线方向画出系统的相轨迹.(注意！等倾线不同于相轨迹！).讲解例7-3-2(P281).（三）、（三）、法：法：系统的相变量方程为：在上式中的后一式中等号右边加、减 项，得式中：从而得到
27、：在点(x1，x2)附近的小邻域内,可将(x1，x2)视为常量.n如果选取新坐标为(0x1，x2)，则以上方程是在新坐标系中的一个圆，其圆心在 ，半径R为从圆心到所取点 的距离。在所取点附近以此画弧就近似地表示了所选取点附近的相轨迹。因此，相轨迹可用许多段小圆弧连接而成。n例7-3-4(P283)自学. 上图中，如果非线性特性为式(7-2)所描述的饱和特性，线性部分的方程为：联列式（7-2）和（7-13）得到整个系统的数学模型为：7.4.1 阶跃函数 现在研究不同的输入函数作用下系统的相轨迹与性能： （1）阶跃函数 ，因为 ,系统 的分区域方程简化为： 区域2的奇点为(0,0),当 时是稳定焦
28、点，当 时是稳定节点。 设 ,式（7-15） 具体化为：（2）斜坡函数r(t)=Vt ，因为 ,这时系统的分区域方程为：7.4.2 斜坡函数 区域2的奇点为(V/Kk，0),当4KTk1时是稳定焦点，当4KTk1时是稳定节点。 固定T=1，K=4，k=1，g=0.2具体化为： 接下来研究相轨迹与参数的关系，分三种情况来分析： 情况一：情况一：VKkg(0.8),比如取V=1.2，这时区域2的相轨迹应收敛到稳定焦点（0.3，0）, 但（0.3，0）不在区域2范围内，故称为之虚焦点。绘出系统的相轨迹场如图7-9（b）, 显然，每一条相轨迹的座标都趋向无穷大，即稳态误差为无穷大。 情情况况二二：V
29、a时，饱和特性输出时，饱和特性输出x(t)为为式中，式中，由于输出波形为奇函数，由于输出波形为奇函数，A10,饱和特性描述函数求得如下：饱和特性描述函数求得如下：3）继电特性：）继电特性： 图7-13 求继电特性的描述函数 考虑图7- 13（a）带滞环的继电特性，当输入为时，输出y(t)如图7-13（b）所示，并且4)间隙特性：间隙特性： 当输入当输入 时时，由间隙由间隙的数学描述可知，间隙输出的数学描述可知，间隙输出x(tx(t) ) 为为X(t)=式中式中B1于是，可求得间隙的描述函数于是，可求得间隙的描述函数N（A）为为5）一般非线性）一般非线性： 描述函数不仅适合于分段线性系统，也适合
30、于一般非线性系统，只要能求出非线性环节的描述函数。我们举一个例子:因为它是单值、奇对称的， ，先求出y(t)： 所以 概括起来，求描述函数的过程是：先根据已知的输入x(t)=Xsint和非线性特性y=f(x)求出输出 ,然后由积分式求出 ，求出N(X) 。主要工作量和技巧主要在积分。 此外，描述函数也可以由实验近似获得。当系统具有良好的低通特性时，给系统施加正弦信号，其输出也近似为正弦信号。改变输入正弦信号的幅值，记录输出信号的幅值和相位，即可近似求出 。 系统开环部分可分离为：系统开环部分可分离为： 非线性环节非线性环节N(A) N(A) 线性部分线性部分G(s)G(s)假定：假定：非线性环
31、节非线性，即不是时间的函数；非线性环节非线性，即不是时间的函数；非线性环节特性是斜对称的；非线性环节特性是斜对称的； 系统的线性部分具有较好的低通滤波性能。系统的线性部分具有较好的低通滤波性能。类似传递函数类似传递函数谐波线性化方法谐波线性化方法非线性系统的频率特性法非线性系统的频率特性法7.7.6 6 分析分析非线性系统的谐波平衡分析法非线性系统的谐波平衡分析法1)1)非线性闭环系统的稳定性非线性闭环系统的稳定性( (乃奎斯特判据乃奎斯特判据) )若开环稳定，则闭环稳定若开环稳定，则闭环稳定的充要条件是的充要条件是G(jG(j ) ) 轨迹轨迹不包围不包围G G平面的平面的(-1,j0)(-
32、1,j0)。负倒描述函数（描述函数负倒特性负倒描述函数（描述函数负倒特性) )线性系统：线性系统：(-1,j0)? 将判断线性系统稳定的结论推广到N(A)为非线性函数的情况。因为X连续变化时N(A)是复平面上的一根曲线，所以闭环系统是否稳定取决于曲线G(j)是否包围1N(A)曲线。具体讲就是(P310)：1)、在复平面上，如果曲线G(j)不包围1N(A)曲线，那么闭环系统稳定；2)、如果G(j)曲线包围1N(A)曲线，那么闭环系统不稳定；3)、如果曲线G(j)与曲线1N(A)相交，那么闭环系统出现自振荡（极限环）。图示： G(j G(j ) ) 与负倒描述函数相交与负倒描述函数相交 闭环系统出
33、现自持振荡闭环系统出现自持振荡( (极限环振荡极限环振荡) ) 稳定稳定 ？不稳定？不稳定？ 振幅（振幅（A A）？）？ 频率频率( ( ) )？设：系统开环的线性部分设：系统开环的线性部分G(jG(j ) )稳定稳定 G(j G(j ) )不包围负倒描述不包围负倒描述函数函数 闭环系统稳定闭环系统稳定 G(j G(j ) )包围负倒描述函包围负倒描述函数数 闭环系统不稳定闭环系统不稳定2）极限环的稳定性）极限环的稳定性 正如相平面法中所讨论的，极限环本身存在一个稳定性问题，极限环的稳定性也可以用描述函数来分析。参见图7-16 : 7-16 极限环的稳定性 图中A、B两点都出现极限环，先看A点
34、：如果因某种干扰使振荡幅值略有减小，比如工作点移到D，D点不被G(j)曲线包围，这时闭环系统应趋向稳定振荡幅值应逐渐减小到零（停振）；反之，如果因某种干扰使振荡幅值略有增大，比如工作点移到C，C点被G(j)曲线包围，这时闭环系统应趋向不稳定振荡幅值应逐渐增大，工作点移到、B；可见A点属不稳定极限环。 再看B点：如果因某种干扰使振荡幅值略有减小，比如工作点移到F，F点被 曲线包围，这时闭环系统应趋向不稳定振荡幅值应增大，增大后又回到B电；反之，如果因某种干扰使振荡幅值略有增大，比如工作点移到E，E点不被 曲线包围，这时闭环系统应趋向稳定振荡幅值应减小，减小后又回到B电；可见A点属稳定极限环。返回
35、7.7.7 7 基于基于SIMULINK SIMULINK 的非线性系统分析的非线性系统分析 SIMULINK 对于非线性系统的分析与设计是很有用的，SIMULINK提供了死区、饱和、继电等多种非线性模块，也能构成很复杂的非线性函数。注：本节内容将上机操作练习！ 考虑图7-17所示非线性系统，图7-17 非线性系统先考虑非线性环节为死区特性，用SIMULINK构成彷真模型如图7-18。图7-18 SIMULINK仿真模型 在Matlab6.5的MATLAB窗中双击SIMULINK图标就打开Simulink Library Browser窗，在此窗口进入FileNewModel，就会打开一个un
36、titled窗(可以用Save as保存此窗口并改名)。在Simulink library Browser 窗口下有DISCONTINUTIES、MATH OPERATIONS、SINK、SOURC等子目录，每个子目录下都包含若干可利用的模块，可直接拖到untitled窗口。图7-18中，积分环节(Integrator)、死区特性(Dead Zone)来自SIMULINKDISCONTINUTIES，求和、增益(Gain)、乘积(Dot Product)来自SIMULINKMATH OPERATIONS，XY Graph 来自SIMULINKSINK，阶跃函数来自SIMULINKSOURCE。
37、 在SIMULATION SIMULATION PARAMETERSSOLVER中设置SolverType为Fixed Step，Solver(步长)为0.05，Stop Time 为50。运行SIMULATIONSTART，绘出的相轨迹如图7-19（a）。可以看到由于死区的存在,稳态误差不为0。 在XY Graph上绘出相轨迹，关键是得到x和信号 ,图7-18中它们分别来自Integrator2的输出和输入。双击某一个模块，就会出现该模块的设置窗口，我们设置:Integrator1的初始条件为1（相当于 ),Integrator2的初始条件为0（相当于 ），阶跃函数的幅值为0(相当于输入为0
38、)死区特性的死区为-0.3，0.3。图7-19 SIMULINK绘出的相轨迹图 然后，我们将非线性模块依次换为饱和特性、继电特性和带滞环的继电特性.图7-18（b）是饱和特性对应的相轨迹，饱和特性的线性区为-0.5，0.5。图7-18（c）是继电特性对应的相轨迹，这个继电特性没有滞环。图7-18（d）是带滞环的继电特性对应的 相 轨 迹 ， 滞 环 宽 度 为 -0.1， 0.1。 图 7-18（c）（d）都出现了极限环。 如果要将生成的相轨迹图保存为图形文件，可执行下列程序：subplot(1,1,1),plot(xout(:,1),xout(:,2)grid on,axis(-1,1,-1
39、,1)xlabel(X),ylabel(dX) 执行后进入FileExport，即可保存为一定格式的图形文件。注：例7-10-1、7-10-2上机操作。小小 结结 本章介绍了非线性系统分析的基础知识，主要包括相平面法和描述函数法。本章的内容还揭示了非线性系统的一些特殊现象，诸如：非线性系统存在全局稳定和局部稳定问题，非线性系统的稳定性与初始条件有关，非线性系统的振型与初始条件和输入幅值都有关，非线性系统的稳定极限环等，这些在线性系统中是不会出现的。 返回nP326: 7- HomeworkThanks for AttentionThanks for AttentionThanks for AttentionThanks for Attention！
