控制系统的稳定性分析精.ppt
控制系统的稳定性分析控制系统的稳定性分析MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答第1页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答本章主要内容本章主要内容n n原理要点n n11.1 系统稳定性的MATLAB直接判定 n n11.1.1 MATLAB直接判定的相关函数n n11.1.2 MATLAB直接判定实例n n11.2 系统稳定性的MATLAB图解判定n n11.2.1 MATLAB图解判定的相关函数n n11.2.2 MATLAB图解判定实例n n11.3 MATLAB LTI Viewer稳定性判定实例第2页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答原理要点系统稳定的概念 经典控制分析中，关于线性定常系统稳定性的概念是：若控制系统在初始条件和扰动作用下，其瞬态响应随时间的推移而逐渐衰减并趋于原点(原平衡工作点)，则称该系统是稳定的。反之，如果控制系统受到扰动作用后，其瞬态响应随时间的推移而发散，输出呈持续振荡过程，或者输出无限制地偏离平衡状态，则称该系统是不稳定的。第3页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答原理要点系统稳定的意义 系统稳定性是系统设计与运行的首要条件。只有稳定的系统，才有价值分析与研究系统自动控制的其它问题，例如只有稳定的系统，才会进一步计算稳态误差。所以控制系统的稳定性分析是系统时域分析、稳态误差分析、根轨迹分析与频率分析的前提。对一个稳定的系统，还可以用相对稳定性进一步衡量系统的稳定程度。第4页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 系统的相对稳定性越低，系统的灵敏性和快速性越强，系统的振荡也越激烈。第5页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答原理要点系统特征多项式 设线性定常系统闭环传递为：式中，称为系统特征多项式。第6页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答为系统特征方程。第7页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答原理要点系统稳定的判定 对于线性连续系统，如果系统的所有特征根(极点)的实部为负，则系统是稳定的；如果有实部为零的根，则系统是临界稳定的（在实际工程中视临界稳定系统为不稳定系统）；反之，如有正实部的根，则系统不稳定。线性连续系统稳定的充分必要条件是：描述该系统的微分方程的特征方程的根全具有负实部，即第8页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 全部根在左半复平面内。或者说系统的闭环传递函数的极点均位于左半s平面内。线性离散系统稳定的充分必要条件是：如果闭环线性离散系统的特征方程根或者闭环脉冲传递函数的极点为则当所有特征根的模都小于1时，即 该线性离散系统是稳定的：如果模的值大于1时，则该线性离散系统是不稳定的。第9页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答原理要点其它稳定性判据 除上述判据之外，还有很多其它判据（其它分析方法中，后面各章将阐述）从各个不同的角度对系统的稳定性加以判别，说明系统稳定性是系统能够成立与运行的首要条件。第10页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答11.1 系统稳定性的 MATLAB直接判定第12页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答本节主要内容本节主要内容n n11.1.1 MATLAB直接判定的相关函数n n11.1.2 MATLAB直接判定实例第13页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 由系统的稳定判据可知，实际上是判定系统闭环特征方程的根的位置。其前提需要求出特征方程的根。MATLAB提供了与之相关的函数，见表11.1：第14页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答表11.1 判定系统稳定的MATLAB函数p=eig(G)p=eig(G)求取矩阵特征根。系统的模型求取矩阵特征根。系统的模型GG可以是传递函数、状态方程可以是传递函数、状态方程和零极点模型，可以是连续或和零极点模型，可以是连续或离散的离散的 P=pole(G)/Z=zero(G)P=pole(G)/Z=zero(G)分别用来求系统的极点和零点。分别用来求系统的极点和零点。GG是已经定义的系统数学模型是已经定义的系统数学模型 p,z=pzmap(sys)p,z=pzmap(sys)求系统的极点和零点。求系统的极点和零点。syssys是是定义好的系统数学模型定义好的系统数学模型 r=roots(P)r=roots(P)求特征方程的根。求特征方程的根。P P是系统闭是系统闭环特征多项式降幂排列的系数环特征多项式降幂排列的系数向量向量 第15页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答11.1.2MATLAB直接判定实例直接判定实例例1：已知系统闭环传递函数为用MATLAB判定 稳定性。num=1 0 2 1;num=1 0 2 1;den=1 2 8 12 20 16 16;den=1 2 8 12 20 16 16;G=tf(num,den)G=tf(num,den)%得到系统模型得到系统模型第16页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答Transfer function:Transfer function:s3+2 s+1 s3+2 s+1-s6+2 s5+8 s4+12 s3+20 s2+16 s+16s6+2 s5+8 s4+12 s3+20 s2+16 s+16 p=eig(G)p=eig(G)%求系统的特征根求系统的特征根p=p=0.0000+2.0000i 0.0000+2.0000i 0.0000-2.0000i 0.0000-2.0000i -1.0000+1.0000i -1.0000+1.0000i-1.0000-1.0000i-1.0000-1.0000i 0.0000+1.4142i 0.0000+1.4142i 0.0000-1.4142i 0.0000-1.4142i p1=pole(G)p1=pole(G)%求系统的极点求系统的极点第17页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答p1=p1=0.0000+2.0000i 0.0000+2.0000i 0.0000-2.0000i 0.0000-2.0000i-1.0000+1.0000i-1.0000+1.0000i -1.0000-1.0000i -1.0000-1.0000i 0.0000+1.4142i 0.0000+1.4142i 0.0000-1.4142i 0.0000-1.4142i r=roots(den)r=roots(den)%求系统特征方程的根求系统特征方程的根r=r=0.0000+2.0000i 0.0000+2.0000i 0.0000-2.0000i 0.0000-2.0000i -1.0000+1.0000i -1.0000+1.0000i -1.0000-1.0000i -1.0000-1.0000i 0.0000+1.4142i 0.0000+1.4142i 0.0000-1.4142i 0.0000-1.4142i第18页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 系统特征根有2个是位于s左半平面的，而4个位于虚轴上。由于有位于虚轴的根，系统是临界稳定的。在实际工程应用上看，系统可认为是不稳定的。分析：由不同MATLAB函数求得的系统特征方程根是一致的。在需要时根据情况选择使用。第19页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答例2：给定系统如图11.1，给出MATLAB程序判 定系统是否稳定，要求程序给出适当提 示。第20页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答um0=1 3;um0=1 3;den0=2 4 5 8 10;den0=2 4 5 8 10;G=tf(num,den);G=tf(num,den);Gc=feedback(G,1);Gc=feedback(G,1);num,den=tfdata(Gc,v);num,den=tfdata(Gc,v);r=roots(den);r=roots(den);disp(disp(系统闭环极点：系统闭环极点：););disp(r)disp(r)a=find(real(r)0);a=find(real(r)0);b=length(a);b=length(a);if b0if b0 disp(disp(系统不稳定系统不稳定.);.);else disp(else disp(系统稳定系统稳定.);.);endend第21页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答程序运行结果程序运行结果：系统闭环极点：系统闭环极点：0.0682+2.1811i 0.0682+2.1811i 0.0682-2.1811i 0.0682-2.1811i-1.1469+0.7535i-1.1469+0.7535i-1.1469-0.7535i-1.1469-0.7535i 0.0786+1.4147i 0.0786+1.4147i 0.0786-1.4147i 0.0786-1.4147i 系统不稳定系统不稳定.例3：某控制系统的方框图如图11.2所示。试用MATLAB确定当系统稳定时，参数K的取值范围（假设 ）。第22页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图11.2 例3系统框图由题，闭环系统的特征方程为：整理得：第23页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 当特征方程的根均为负实根或实部为负的共轭复根时，系统稳定。先假设K的大致范围，利用roots()函数计算这些K值下特征方程的根，然后判断根的位置以确定系统稳定时K的取值范围。程序如下：k=0:0.01:100;k=0:0.01:100;for index=1:100for index=1:100 p=2 15 27 k(index)+12 k(index)+1;p=2 15 27 k(index)+12 k(index)+1;r=roots(p);r=roots(p);if max(real(r)0 if max(real(r)0 第24页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答break;break;end endendendsprintf(sprintf(系统临界稳定时系统临界稳定时K K值为值为:K=%7.4fn,k(index):K=%7.4fn,k(index)程序运行结果为：ans=系统临界稳定时K值为:K=0.9900第25页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答11.2系统稳定性的 MATLAB图解判定第27页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答本节主要内容本节主要内容n n 11.2.1MATLAB图解判定的相关函数n n 11.2.2MATLAB图解判定实例第28页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答11.2.1MATLAB图解判定的相关函数图解判定的相关函数 对于给定系统G，pzmap(G)函数在无返回参数列表使用时，直接以图形化的方式绘制出系统所有特征根在S复平面上的位置，所以判定系统是否稳定只需看一下系统所有极点在S复平面上是否均位于虚轴左侧即可。这种图形化的方式更直观。第29页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答11.2.2MATLAB图解判定实例图解判定实例例4：已知一控制系统框图，如图11.3所示，试 判断系统的稳定性。第30页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 G1=tf(1 1,2 1);G1=tf(1 1,2 1);G2=tf(5,2 3 1);G2=tf(5,2 3 1);H1=tf(1,2 1);H1=tf(1,2 1);Gc=feedback(G2*G1,H1)Gc=feedback(G2*G1,H1)%得到闭环系统传递函数得到闭环系统传递函数 Transfer function:Transfer function:10 s2+15 s+5 10 s2+15 s+5-8 s4+20 s3+18 s2+12 s+68 s4+20 s3+18 s2+12 s+6 pzmap(Gc)pzmap(Gc)第31页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答第32页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答分析：由于特征根全部在S-平面的左半平面，所以此负反馈系统是稳定的。例5：给定离散系统闭环传递函数分别为：和采样周期均为0.1秒。分别绘制系统零极点分布图，并判定各系统稳定性。第33页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 num=1 4.2 5.43;num=1 4.2 5.43;den=1-2.7 2.5 2.43-0.56;den=1-2.7 2.5 2.43-0.56;Gc=tf(num,den,0.1)Gc=tf(num,den,0.1)Transfer function:Transfer function:z2+4.2 z+5.43 z2+4.2 z+5.43-z4-2.7 z3+2.5 z2+2.43 z-0.56z4-2.7 z3+2.5 z2+2.43 z-0.56 Sampling time:unspecified Sampling time:unspecified pzmap(Gc)pzmap(Gc)第34页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图11.5 例5运行结果第35页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 由上图可知，系统G在单位圆外有极点存 在，系统是不稳定的。num=0.68 5.43;num=0.68 5.43;den=1-1.35 0.4 0.08 0.002;den=1-1.35 0.4 0.08 0.002;G2=tf(num,den,0.1)G2=tf(num,den,0.1)Transfer function:Transfer function:0.68 z+5.43 0.68 z+5.43-z4-1.35 z3+0.4 z2+0.08 z+0.002z4-1.35 z3+0.4 z2+0.08 z+0.002Sampling time:0.1Sampling time:0.1 pzmap(G2)pzmap(G2)第36页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答图11.6 例5运行结果第37页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 由图可知，系统G2闭环传递函数的所有极点都位于单位圆内部，据此可知此闭环系统是稳定的。第38页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答11.3MATLAB LTI Viewer稳定性判定实例稳定性判定实例 MATLAB LTI Viewer是MATLAB为LTI（Linear Time Invariant）系统的分析提供的一个图形化工具。用它来可以很直观简便地分析控制系统的时域和频域响应。用MATLAB LTI Viewer来观察闭环系统的零极点分布情况，需要首先在MATLAB中建立系统的闭环系统传递函数模型 第40页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 例6：已知单位负反馈控制系统的开环传递函数为 用MATLAB LTI Viewer观察闭环系统的零极点分布情况，并判断此闭环系统的稳定性。第41页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答1.建立系统模型。z=-3;z=-3;p=0-2-5;p=0-2-5;k=3;k=3;G=zpk(z,p,k)G=zpk(z,p,k)Zero/pole/gain:Zero/pole/gain:3(s+3)3(s+3)-s(s+2)(s+5)s(s+2)(s+5)Gc=feedback(G,1)Gc=feedback(G,1)Zero/pole/gain:Zero/pole/gain:3(s+3)3(s+3)-(s+4.599)(s2 +2.401s+1.957)(s+4.599)(s2 +2.401s+1.957)第42页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答2.打开LTI Viewer。在命令窗口输入：ltiview即进入LTI Viewer窗口，如图11.6。图图11.6 LTI Viewer11.6 LTI Viewer窗口图窗口图11.7 LTI Viewer11.7 LTI Viewer导入系统模型窗口导入系统模型窗口 第43页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 3.导入在MATLAB中建立好的系统模型。在LTI Viewer窗口上选FileImport，出现如图11.7的窗口。我们可以从Workspace项中选刚建立好的系统Gc。系统默认会给出系统的阶跃响应曲线。图图11.8 11.8 选择系统响应类型图图选择系统响应类型图图11.9 11.9 系统零极点分布图系统零极点分布图第44页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答 4.观察系统的零极点分布。在图11.8的窗口中点击右键，选Plot typesPole/Zero。即绘制出系统的零极点分布图。如图11.9所示。由图可知，系统的闭环极点全部位于S-平面左半平面。可判定系统是稳定的。第45页，本讲稿共46页MATLAB与控制系统仿真实践，北京航空航天大学出版社，2009.8.在线交流，有问必答结束结束第46页，本讲稿共46页
