matlab非线性控制系统分析优秀PPT.ppt

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