非线性系统分析 自动控制原理件.pptx

《非线性系统分析 自动控制原理件.pptx》由会员分享，可在线阅读，更多相关《非线性系统分析 自动控制原理件.pptx（70页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1第 8 章 非线性控制系统8.1 概述8.2 非线性系统的特点8.3 相平面法8.4 描述函数法8.5 MATLAB在非线性控制系统分析中的应用第1页/共70页28.1 概述 非线性系统与线性系统有着很大的差别，诸如非线性系统的响应取决于输入信号的幅值和形式，不能应用叠加原理，目前还没有统一的且普遍适用的处理方法。由于非线性系统的复杂性和特殊性，受数学工具限制，一般情况下难以求得非线性微分方程的解析解，通常采用工程上适用的近似方法。（1）相平面法 （2）描述函数法 （3）逆系统法第2页/共70页3 1.相平面法:一种图解分析方法，适用于具有严重非线性特性的一阶、二阶系统，该方法通过在相平面绘

2、制相轨迹曲线，确定非线性微分方程在不同初始条件下解的运动形式。2.描述函数法:一种等效线性化的图解分析方法，该方法对于满足结构要求的非线性系统，通过谐波线性化，将非线性特性近似为复变增益环节，然后推广应用频率法，分析非线性系统的稳定性或自激振荡。3.逆系统法:运用内环非线性反馈控制，构造伪线性系统，以此为基础，设计外环控制网络，该方法直接应用数学工具研究非线性控制问题，是非线性系统研究的一个发展方向。但是，这些方法主要是解决非线性系统的“分析”问题，且以稳定性问题为主展开的。非线性系统的“综合”方法的研究成果远不如稳定性问题研究所取得的成果。第3页/共70页41.饱和特性 2.死区特性 3.间

3、隙特性 图8.1 饱和非线性特性图8.2 死区非线性特性 图8.3 间隙非线性特性 典型非线性特性第4页/共70页54.继电器特性 图8.4 继电器型非线性特性第5页/共70页6非线性系统的运动特点 由于描述非线性系统运动的数学模型为非线性微分方程，叠加原理不再适用，因此非线性系统的运动表现出以下特点：1.稳定性分析复杂 2.系统的零输入响应形式 3.自激振荡(极限环)4.频率响应 第6页/共70页7 1.稳定性分析复杂:在研究非线性系统的稳定性问题时，必须要明确两点：a.指明给定系统的初始状态或输入信号 b.指明相对于哪一个平衡状态来分析系统的 稳定性。2.系统的零输入响应形式:线性系统的零

4、输入响应形式与系统初始状态的幅值无关。某些非线性系统的零输入响应形式与系统的初始状态有关。当初始状态不同时，同一个非线性系统可能有不同的响应形式，如单调收敛、振荡收敛或振荡发散等。第7页/共70页8例如：具有分段线性特性的非线性增益控制系统，当初始状态(初始误差)。E时，系统的零输入响应形式为振荡收敛形式，如图8.5所示第8页/共70页9 3.自激振荡(极限环):线性定常系统：例如典型二阶线性系统，如果阻尼比=0，在初始状态的激励下，系统的零输入响应为等幅周期振荡,其角频率 取决于系统的参数，其振幅A与初始状态有关。但是，实际的线性系统要维持振幅A和角频率 不变的等幅周期振荡是不可能的。一是系

5、统的参数会发生变化，即使很微小的变化，也将导致 0；二是假定系统的参数不变，0，然而，系统不可避免地会受到扰动，将使响应的振幅A发生变化，因此，原来的等幅周期振荡不复存在。第9页/共70页10 4.频率响应:在正弦输入信号作用下，非线性系统呈现出一些在线性系统中见不到的特殊现象，诸如跳跃谐振和多值响应、倍频振荡和分频振荡、频率捕捉(跟踪)现象等。图8.6为一机械系统，由重物、阻尼器、非线性弹簧组成。其动态特性的微分方程为 (8-1)式中 重物的位移；M质量；B阻尼器的粘性摩擦系数；非线性弹簧力。第10页/共70页11 图8.7描绘了大于零、等于零、小于零三种情况下频率和振幅的关系。在对图8.6

6、所示的系统进行强迫振荡实验 时，系统的微分方程为式中为外作用函数。式中 为外作用函数 图8.6 机械系统 图8.7 频率和振幅的关系曲线第11页/共70页12 因此，响应曲线实际上是不连续的，并且对于频率增加和减小的两种情况,响应曲线上的点沿着不同的路线移动。点2与点5之间曲线对应的振荡是不稳定的振荡，在实验中是观测不到的。(a)具有硬弹簧的机械系统 (b)具有软弹簧的机械系统 图8.8 机械系统的频率响应第12页/共70页138.3相平面法 相平面法是庞卡莱(H.Poincare)提出来的一种用图解法求解一阶、二阶微分方程的方法，它实质上属于状态空间分析法在二维空间中的应用，该方法适合于研究

7、给定初始状态的二阶自由运动系统和给定初始状态及非周期输入信号(如阶跃、斜坡或脉冲信号等)的二阶系统 相平面的基本概念 相平面图的绘制方法 奇点和极限环 相平面分析举例第13页/共70页14相平面的基本概念 考虑二阶线性系统 (8-2)式中 与 是阻尼比和无阻尼自然振荡频率。设系统仅由初始条件激励。这一系统的状态可以用两个变量，和 来描述。若令，则方程(8-2)可化为 (8-3)(8-4)只要给定初始条件 、或 、，由这两个一阶联立微分方程便可唯一地确定系统的状态。如此定义的变量和称为相变量(或状态变量)。图8.9(a)绘出了初始条件为及时，和在不同阻尼下的时间响应曲线。第14页/共70页15（

8、a）(b)图8.10 相平面图图8.9 时间响应与相轨迹第15页/共70页16相平面图的绘制方法 设描述二阶系统的微分方程为 (8-5)是的线性函数或非线性函数。若令为相变量，并将式(8-5)化为两个一阶微分方程 (8-6)(8-7)用式(8-6)去除式(8-7)，于是得到一个以x为自变量，为因变量，不显含时间t的一阶微分方程 (8-8)式(8-8)给出了相轨迹通过点的斜率。根据此式，用解析法或图解法即可绘出相平面图。第16页/共70页17 1.相平面图的特点:相平面图的对称性 相平面图往往是关于原点或坐标轴对称的，故绘制时可只画其中的一部分，而另一部分可根据对称原理添补上。相平面图的对称性可

9、以从相轨迹的斜率来判断。若相平面图关于轴对称，则相轨迹曲线在 和 点上的斜率相等，符号相反。由式(8-8)，应有 即 是关于x的奇函数。若相平面图关于x轴对称，则相轨迹曲线和的斜率相等，符号相反，应有 即 是关于 的偶函数。第17页/共70页18 若相平面图关于原点对称，则相轨迹曲线在 和 点上的斜率相等，符号相同，应有 即有 。第18页/共70页19 1.相平面图的特点:相平面图上的奇点和普通点 相平面上任一点 ，只要不同时满足 和 ，则由式(8-8)确定的斜率是唯一的，通过该点的相轨迹有且仅有一条，这样的点称为普通点。在相平面上，同时满足 和 的点，由于 相轨迹的斜率不是一个确定的值，说明

10、通过该点的相轨迹曲线有一条以上，这样的点称为奇点，显然奇点只分布在相平面的x轴上。第19页/共70页20（3）相轨迹通过x轴的斜率在x轴上，所有点都满足 。除奇点外相轨迹在x轴上的斜率为所以，除了奇点外，相轨迹和x轴垂直相交。（4）相轨迹移动的方向 在相平面的上半平面，由于，则x随着参变量时间t的加而增大，所以系统状态沿相轨迹由左向右运动；反之，下半平面，由于，则x随着时间t的增加而减小，所以系统态沿相轨迹由右向左运动。系统状态沿相轨迹的移动方向相轨迹上的箭头表示。第20页/共70页21 2.绘制相平面图的解析法 若系统微分方程比较简单，则对式(8-8)直接积分即可求得相轨迹方程 (8-9)式

11、中包含了初始条件。给定不同初始条件，由 式(8-9)可直接绘出系统的相平面图。若系统微分方程不能直接积分求解，则可先求得时间解 ，然后消去变量t求得相轨迹方程。若消去t有困难或过于烦杂，则可求得不同t时的x，值，据此数值关系画出相轨迹曲线。第21页/共70页22例8-1 二阶系统的微分方程为，试绘制系统的相平面图。解 系统方程可改写为 (8-10)方程(8-10)可用分离变量法进行积分，求得相轨迹方程为 (8-11)式中C为常量，由初始条件确定。设初始状态为 ，则C=。由方程(8-11)可知，系统相轨迹为一组以坐标原点为中心的椭圆轨迹簇，如图8.11所示。其中粗实线是初始条件为 的相轨迹。图8

12、.11 例8-1的相平面图第22页/共70页233.绘制相平面图的图解法等倾线法例8-2 试用等倾线法求下列方程的相平面图。(8-17)解 式(8-17)是非线性微分方程,但可分解为两个线性微分方程 ，(8-18),(8-19)由方程(8-17)可知 ,而 。因此相平面图对称于x轴。由式(8-18)可得上半平面的等倾线方程：第23页/共70页24奇点和极限环 1.奇点 对于二阶系统 (8-21)相轨迹的斜率可表示为 (8-22)在奇点处，相轨迹的斜率不确定，即同时满足 (8-23)如果把相变量x视为位移，于是和可以理解为速度和加速度。在奇点处，由于系统的速度和加速度均为零，因此奇点就是系统的平

13、衡点。第24页/共70页25系统奇点的分类(1)焦点(2)节点(3)中心点(4)鞍点第25页/共70页262极限环（1）稳定极限环（2）不稳定极限环（3）半稳定极限环 图8.16 极限环示意图第26页/共70页27 解 由 求得系统的奇点为 根据式(8.25)在奇点处进行线性化来确定奇点的性质。在(xi,0)奇点附近，系统的线性化方程为 在奇点(0，0)处，xi=x1=0，则系统的线性化方程为式中阻尼比0 1，因此奇点(0，0)为稳定焦点。在奇点(-2，0)处，xi=x2=-2，代入前式得线性化方程为由奇点类型可知，奇点(-2，0)为鞍点，是不稳定奇点。例8-3 某系统方程如下，试分析系统的稳

14、定性第27页/共70页28 利用等倾线法可求得相平面图，如图8.17所示。可以看到通过鞍点的一条分界线，把相平面分为两个区域。在阴影区域内，所有相轨迹都收敛于稳定焦点(0，0)，是稳定区域。在此范围外，则所有相轨迹都将趋于无穷，是不稳定区域。这证实了非线性系统的重要特点：系统的稳定性与初始条件有关。图8.17 例8-3的相平面图第28页/共70页29相平面分析举例1.继电型控制系统的分析:（1）理想继电器特性 图8.18 理想继电器型非线性系统 0 设继电型控制系统如图8.18(a)所示,试分析在阶跃信号作用下系统的性能。继电型特性为：当e0时，m=M；当e0时，m=-M。因此分界线为直线e=

15、0。它把相平面分成两个线性区域区、区。第29页/共70页30 在区域内，e0，m=M，系统方程为 (8-26)由(8-26)式可得等倾线方程 等倾线是平行于e轴的直线，其中有一条特殊的等倾线，即当a=0时的等倾线 ，此时，相轨迹的斜率与相应的等倾线斜率相等，全部相轨迹曲线都趋近于该直线 。相轨迹曲线簇如图8.18(b)右半平面所示。在区域内，e0，m=-M，系统方程为 (8-27)比较方程(8-26)、(8-27)可知，其相平面图对称于原点。利用对称性求得相轨迹曲线簇如图8.18(b)左半面所示。第30页/共70页31（2）滞环继电特性图8.18(a)所示的非线性系统中，若继电器元件换成如图8

16、.19(a)所示的滞环特性，则该非线性特性可用以下方程描述：在0时的平面内，分界线为e=+。在0的平面内，分界线为e=-。它们把相平面分为两部分。其右半平面，系统在+M信号作用下，系统方程为式(8-26)，相轨迹为曲线簇。其左半平面，系统在-M信号作用下，系统方程为式(8-27)，相轨迹为曲线簇。相平面如图8.19(b)所示。在阶跃输入作用下，系统的运动轨迹如图中粗实线所示。相轨迹收敛于稳定极限环，极限环随 的增大而增大。第31页/共70页32mMe-M(a)-mM-M(a)e图8.19 滞环继电型非线性系统 图8.20 死区继电型非线性系统第32页/共70页33(3)死区继电特性:图:8.1

17、8所示的非线性系统中，若继电元件具有如图8.20(a)所示的死区特性，则可用以下方程描述 当 e ，m=+M 当 e-，m=-M 当-e ，m=0 分界线为e=+和e=-，它们将相平面分为三个区域，如图8.20(b)所示。在区域、中，系统方程分别用式(8-26)、(8-27)描述，相轨迹分别为曲线族、。在区域中，m=0，系统的误差方程为 可求得相轨迹的斜率 为常数，即其相轨迹是一组 斜率为 的直线。由上式还可得到：当 时,必有 。因此在区域内，直线 上所有点都是奇点(又称奇线或平衡线)。系统的相平面图如8.20(b)所示。由图可知系统可能稳定在奇线上任一点。第33页/共70页34 为了缩短调节

18、时间，减少振荡次数，继电控制系统可采用速度反馈校正，如图8.21(a)所示。继电元件的输入信号为 ，当系统在阶跃信号r(t)=1(t)作用下，由e=r-c可得继电元件输入信号 ，因此 当 则 当 则 分界线由方程 确定，这是一条通过原点，斜率为-1/Kt的直线。它将相平面分为、两个区域，分别由方程(8-26)、(8-27)描述。图8.21(b)中给出了分界线及其相轨迹曲线、。第34页/共70页35在阶跃输入作用下，系统状态的运动如图中实线所示。相轨迹由初始点A0开始，沿相轨迹移动到达分界线上的衔接点A1；进入线性区后，沿相轨迹移动到下一个衔接点A2，。当衔接点位于分界线B1B2线段内时，相轨迹

19、将沿分界线向原点滑动，最后趋近于原点，这就是非线性系统的“滑动”现象，该现象可以缩短系统的调节时间。比较图8.19(b)及图8.21，可以明显看到速度反馈校正的效果：超调量减小，调节时间缩短，振荡次数减少。图8.21 继电型非线性系统的速度反馈校正第35页/共70页362.具有非线性增益控制系统的分析在线性系统中，增益的选择需要兼顾调节时间、超调量及振荡次数等性能指标，当增益K值取得较大时，系统快速性较好，但超调量大，振荡次数多,如图曲线1所示。若K值较小，超调量、振荡次数将减小，但系统快速性较差，如图曲线2所示。在线性系统中只能选取折中方案。若采用非线性校正，则可能得到较好效果如图曲线3所示

20、。第36页/共70页37（1）阶跃响应分析 设系统输入信号为阶跃函数r(t)=R，当t0时，故在区内系统方程为 (8-34)系统奇点位于原点(0，0)，为稳定节点。在区内，系统方程为 (8-35)系统奇点位于原点(0，0)，为稳定焦点，但是由于该奇点不在区内，系统状态实际上不能到达该点，故称该点为虚奇点。在阶跃输入作用下，系统 状态的运动轨迹如图8.24 所示。相轨迹的起点A由 初始条件e(0)=R，所决定，它经过BCDEF 最终趋向于相平面的原点，虽然响应曲线是振荡的，但 超调量，振荡次数都减小很多。第37页/共70页38（2）斜坡响应分析 设输入信号r(t)=R+V t，当t0时，由式(8

21、.32)和(8.33)可得 (8-36)(8-37)与前面阶跃输入情况相比，相平面的分界线没有变化，但奇点的位置不同。式(8-36)对应的奇点P1位于(V/kK，0)，式(8-37)对应的奇点P2位于(V/K，0)，因为k1，所以P1总在P2的右边。图8.25 图8.26 VkKe0，Re0时的相轨迹 kKe0VKe0，R=0时的相轨迹 第38页/共70页39 3.时间最优控制系统的分析及综合 以图8.27所示直流电动机随动系统为例，电动机输出力矩极限为，其中为电动机力矩系数。若负载为纯惯性的，则加速度的上限为，J为负载转动惯量。图8.27 直流电动机随动控制系统 图8.28系统的响应曲线 据

22、上述限制，为使系统在阶跃输入作用下，响应时间最短，系统的响应过程应如图8.28所示。当t0时，控制器应给出信号，使系统以加速，此时电动机的速度以直线上升。为使当系统输出c=R时,也同时为零而结束响应过程，要求在c到达R值之前的适当时刻，控制器就给出信号，使系统以减速，以直线下降。当c=R时，恰好使=0，这时控制器给出零信号，阶跃响应过程结束。0tRrci第39页/共70页401.1.谐波线性化 系统中常见的非线性特性，当输入为正弦函数时，其输出一般为同周期的非正弦函数。例如对理想继电特性加以正弦输入信号，则输出y y(t t)为与输入同周期的方波，见图8.308.30。图8.30 理想继电特性

23、在正弦输入时的输出波形 描述函数的基本概念8.4描述函数分析法描述函数分析法 描述函数法的基本原理描述函数法的基本原理：当系统满足一定条件时，系统中非线性环节在正弦信号作用下的输出可用一次谐波分量来近似，由此导出非线性环节的近似等效频率特性。第40页/共70页41将方波信号用傅里叶级数表示，即为(8-42)设非线性环节描述为非线性特性的输入信号为输出信号可以表示为傅氏级数形式式中(8-43)若非线性特性具有奇对称特性，则A0=0，如果略去输出高次谐波分量，仅以基波分量近似地代替整个输出，则有(8-44)第41页/共70页42式中2描述函数 非线性特性在进行谐波线性化后，参照幅相频率特性的定义，

24、建立非线性特性的等效幅相特性，即描述函数。把非线性元件输出信号y(t)中的一次谐波分量y1(t)与正弦输入信号x(t)的复数比，称为非线性元件的描述函数，其数学表达式为(8-45)式中 A为非线性元件正弦输入信号的振幅 为非线性元件输出信号中一次谐波分量的振幅；为非线性元件输出信号中一次谐波分量的相位移。第42页/共70页43典型非线性特性的描述函数1.饱和特性饱和特性以及它对正弦输入的输出波形如图8.31(8-46)图8.31 饱和非线性及其输入、输出波形第43页/共70页44饱和特性的描述函数为(8-48)2.死区特性死区特性以及它对正弦输入的输出波形如图8.32所示于是死区特性的描述函数

25、为(8-49)(8-51)第44页/共70页45图8.32 死区非线性及其输入、输出波形第45页/共70页46于是间隙特性描述函数为(8-53)(8-54)3.间隙特性 间隙特性以及它对正弦输入的输出波形如图8.33所示 第46页/共70页47图8.33 间隙非线性及其输入、输出波形 第47页/共70页48(8-55)式中4.继电器特性 （1）具有死区和滞环的继电器 其输出量y(t)的方程为 具有死区和滞环的继电器特性以及它对正弦输入的输出波形如图8.34所示。第48页/共70页49图8.34 具有死区和滞环的继电特性及其输入、输出波形第49页/共70页50于是，具有死区和滞环继电器的描述函数

26、为(8-56)（2）双位继电器 双位继电器非线性的描述函数(8-57)图8.35 双位继电器非线性 第50页/共70页51（4）具有滞环的继电器 三位继电器特性的描述函数 (8-58)图8.36 三位继电器非线性 具有滞环继电器非线性的描述函数(8-59)图8.37 滞环继电器非线性（3）三位继电器第51页/共70页52用描述函数法分析非线性系统1.描述函数法的应用条件 （1）非线性系统能简化成一个非线性环节和一个线性部分且闭环连接 的典型结构形式，如图8.38所示，其中G(s)代表系统的线性部分。图8.38 非线性控制系统（2）非线性环节输入输出特性y(x)应是x的奇函数，即以保证非线性环节

27、的正弦响应不含有常值分量，即 （3）系统的线性部分应具有较好的低通滤波性能。当非线性环节的输入为正弦信号时，实际输出必定含有高次谐波分量，但经线性部分传递之后，由于低通滤波的作用，高次谐波分量将被大大削弱，从而保证描述函数法所分析的结果比较准确。第52页/共70页53 当非线性特性采用描述函数近似等效时，闭环系统的特征方程为1+N(A)G(jw)=0 或 N(A)G(jw)=-1 (8-60)即G(jw)=-1/N(A)(8-61)-1/N(A)称为非线性环节的负倒描述函数。由线性控制系统理论知，线性系统的特征方程为G(jw)=-1(8-62)根据复平面内系统的开环频率特性G(jw)曲线与临界

28、点(-1，j0)的相对位置，应用奈魁斯特(Nyquist)稳定判据，可以分析线性控制系统的稳定性。将方程(8-61)与(8-62)对照，显然可以把奈魁斯特稳定判据，推广应用于谐波线性化的非线性系统，需要修改的仅仅是将复平面内的临界点(-1，j0)扩展为临界曲线，即-1/N(A)曲线。根据奈魁斯特稳定性判据，如果-1/N(A)曲线不被G(jw)曲线包围(8.39(a)则系统是稳定的。2.非线性系统的稳定性分析 第53页/共70页54 如果-1/N(A)曲线被G(jw)曲线全部包围(图8.39(b)，则系统状态在干扰作用下，不能回到平衡状态，所以系统是不稳定的。如果-1/N(A)曲线与线性部分频率

29、特性G(jw)曲线相交(图8.39(c)，交点处的参数，即振幅Ai和频率wi使方程(8-60)或(8-61)成立，非线性系统可能产生的自激振荡.图8.39 非线性系统零平衡状态的稳定性3描述函数分析举例 例8-4 双位继电器非线性系统（图8.40）线性部分的传递函数为 系统的参考输入r(t)=0，系统开始处于静止状态。(1)分析非线性系统零平衡状态的稳定性和自激振荡的稳定性；第54页/共70页55(2)如果系统产生自激振荡，确定自激振荡的参数A和w。解 由式(8-57)求得-1/N(A)曲线在复平面内与负实轴重合。线性部分的频率特性G(jw)为(8-63)在复平面内画出双位继电器的-1/N(A

30、)曲线和线性部分的G(jw)曲线如图8.41所示。图8.40 双位继电器非线性系统 第55页/共70页56 图8.41 例8-4系统的-1/N(A)曲线和G(jw)曲线 G(jw)曲线与负实轴(即双位继电器的-1/N(A)曲线)相交时，G(jw)的虚部为零。令G(jw)的虚部为零，求得G(jw)曲线与-1/N(A)曲线交点处的频率w 于是 (弧度/秒)(弃之)。将 代入式(8-63)得到交点处的参数应满足系统的特征方程第56页/共70页57即有：求得交点处自激振荡的振幅A为A=8M/p。双位继电器的输入信号x(t)为线性部分G(s)的输出为 双位继电器的输出信号y(t)是振幅为M的方波，根据式

31、(8-42)其基波分量为经线性部分后相应的基波分量为即线性部分输出中所含基波分量的振幅为8M/p。由式(8-42)得方波y(t)中的三次谐波分量为 第57页/共70页58而 于是y(t)中的三次谐波分量y3(t)经过线性部分滤波后，相应的振幅成为 4非线性系统的简化（1）非线性特性的并联图8.46 两个非线性部件并联的系统第58页/共70页59（2）非线性特性的串联 当两个非线性环节串联环节时，则先将两个环节的特性等效为一个特性，然后求总描述函数N(A)。图8.47表示了死区特性与带死区的继电特性相串联的等效非线性图形。应当注意，调换串联环节的前后次序，等效特性将会不同（3）存在局部反馈情况

32、图8.48给出了非线性系统的另一种结构，其中非线性部分被线性部分局部反馈所包围。对这种结构，可视 、为并联联结，合并为一个线 G1(s)N1G2(s)N2第59页/共70页1.利用Simulink绘制非线性系统相轨迹图 已知系统的微分方程为求系统的奇点，并通过奇点附近的相轨迹确定奇点的类型 解：由 ，求出系统的奇点为8.5 MATLAB在非线性控制系统分析中的应用第60页/共70页在Simulink下建立如下模型模型代表微分方程 第61页/共70页 附近相轨迹 相轨迹的起点从（2，0）开始，收敛到（0，0）点。可见奇点（0，0）是稳定焦点。第62页/共70页 附近相轨迹 相轨迹的起点从（5，0

33、）开始顺时针移动，在（-4，0）点附近转向后趋于无穷，可见奇点（-4，0）是鞍点，是不稳定奇点。第63页/共70页.利用 MATLAB判断非线性系统的稳定性及自激振荡 具有滞环非线性的系统如下图所示，其中非线性特性特性参数M=4，=1。分析非线性系统是否有自激振荡存在，若有，求自激振荡的参数。第64页/共70页在命令窗口执行以下程序，在复平面内画出-1/N(A)和G(jw)曲线 A=1:0.1:20;%A的范围和步长M=4;delt=1;%定义非线性特征值disN1=4*M/pi./A.*sqrt(1-(delt./A).2)-j*4*M*delt/pi./A.2;%非线性部分的描述函数dis

34、N=-1./disN1;%负倒描述函数w=1:0.01:200;%频率范围和步长num=4;%线性部分传递函数分子den=conv(1 1 0,1 2);%线性部分传递函数分母r,i,w=nyquist(num,den,w);plot(real(disN),imag(disN),r,i);grid;第65页/共70页结果见下图 随着A的增加，-1/N(A)曲线由右向左移动，随着w的增加，G(jw)曲线由下向上移动，G(jw)曲线和-1/N(A)曲线相交，系统出现稳定的自激振荡。第66页/共70页将交点处的图像局部放大 交点处的坐标为（-0.978，-0.1963）第67页/共70页用下面的程序求出自激振荡的幅值和频率A0=spline(real(disN),A,-0.978);%用样条函数插值出-0.978处的幅值w0=spline(i,w,-0.1963);%用样条函数插值出-0.1963处的频率值求出自激振荡的频率w0=1.14，自激振荡的振幅A0=5.08 第68页/共70页69本章结束，谢谢。第69页/共70页自动控制原理70感谢您的观看！第70页/共70页