一、系统稳定及最小相位系统判据 本文来自www.eadianqi.com 对于连续时间系统,如果闭环极点全部在S平面左半平面,则系统是稳定的。 本文来自www.eadianqi.com 对于离散时间系统,如果系统全部极点都位于Z平面的单位圆内,则系统是稳定的。 自动控制网www.eadianqi.com版权所有 若连续时间系统的全部零极点都位于S左半平面;或若离散时间系统的全部零极点都位于Z平面单位圆内,则系统是最小相位系统。 自动控制网www.eadianqi.com版权所有 二、系统稳定及最小相位系统的判别方法 本文来自www.eadianqi.com (一)间接判别(工程方法) 本文来自www.eadianqi.com 劳斯判据:劳斯表中第一列各值严格为正,则系统稳定,如果劳斯表第一列中出现小于零的数值,系统不稳定。 自动控制网www.eadianqi.com版权所有
胡尔维茨判据:当且仅当由系统分母多项式构成的胡尔维茨矩阵为正定矩阵时,系统稳定。 本文来自www.eadianqi.com (二)直接判别 本文来自www.eadianqi.com MATLAB提供了直接求取系统所有零极点的函数,因此可以直接根据零极点的分布情况对系统的稳定性及是否为最小相位系统进行判断。 自动控制网www.eadianqi.com版权所有
例exp4_1.m 自动控制网www.eadianqi.com版权所有 已知某系统的模型如右所示: 自动控制网www.eadianqi.com版权所有 要求判断系统的稳定性及系统是否为最小相位系统。 本文来自www.eadianqi.com %exp4_1.m 自动控制网www.eadianqi.com版权所有 clear 本文来自www.eadianqi.com clc 自动控制网www.eadianqi.com版权所有 closeall 本文来自www.eadianqi.com %系统描述 自动控制网www.eadianqi.com版权所有 a=[12-12;2630;47-8-5;7216]; 本文来自www.eadianqi.com b=[-1001]'; 本文来自www.eadianqi.com c=[-2561];d=7; 自动控制网www.eadianqi.com版权所有 %求系统的零极点 本文来自www.eadianqi.com [z,p,k]=ss2zp(a,b,c,d) 本文来自www.eadianqi.com
%检验零点的实部;求取零点实部大于零的个数 自动控制网www.eadianqi.com版权所有 ii=find(real(z)>0) 本文来自www.eadianqi.com n1=length(ii); 本文来自www.eadianqi.com
%检验极点的实部;求取极点实部大于零的个数 本文来自www.eadianqi.com jj=find(real(p)>0) 本文来自www.eadianqi.com n2=length(jj); 本文来自www.eadianqi.com %判断系统是否稳定 本文来自www.eadianqi.com
if(n2>0) 本文来自www.eadianqi.com disp('thesystemisunstable') 本文来自www.eadianqi.com
disp('theunstablepoleare:') 本文来自www.eadianqi.com disp(p(jj)) 自动控制网www.eadianqi.com版权所有
else 自动控制网www.eadianqi.com版权所有 disp('thesystemisstable') 自动控制网www.eadianqi.com版权所有 end 自动控制网www.eadianqi.com版权所有 %判断系统是否为最小相位系统 本文来自www.eadianqi.com if(n1>0) 自动控制网www.eadianqi.com版权所有 disp('thesystemisanonminimalphaseone') 自动控制网www.eadianqi.com版权所有
else 本文来自www.eadianqi.com disp('thesyetemisaminimalphaseone') 本文来自www.eadianqi.com end 自动控制网www.eadianqi.com版权所有 %绘制零极点图 自动控制网www.eadianqi.com版权所有 pzmap(p,z) 自动控制网www.eadianqi.com版权所有
例exp4_2.m 自动控制网www.eadianqi.com版权所有 系统模型如下所示,判断系统的稳定性,以及系统是否为最小相位系统。 自动控制网www.eadianqi.com版权所有
本文来自www.eadianqi.com %exp4_2.m 自动控制网www.eadianqi.com版权所有 clear 自动控制网www.eadianqi.com版权所有 clc 自动控制网www.eadianqi.com版权所有
closeall 本文来自www.eadianqi.com %系统描述 自动控制网www.eadianqi.com版权所有 num=[3164128]; 自动控制网www.eadianqi.com版权所有 den=[11411052814942117112]; 自动控制网www.eadianqi.com版权所有 %求系统的零极点 自动控制网www.eadianqi.com版权所有 [z,p,k]=tf2zp(num,den) 自动控制网www.eadianqi.com版权所有 %检验零点的实部;求取零点实部大于零的个数 自动控制网www.eadianqi.com版权所有
ii=find(real(z)>0) 本文来自www.eadianqi.com n1=length(ii); 本文来自www.eadianqi.com %检验极点的实部;求取极点实部大于零的个数 本文来自www.eadianqi.com jj=find(real(p)>0) 自动控制网www.eadianqi.com版权所有 n2=length(jj); 本文来自www.eadianqi.com %判断系统是否稳定 自动控制网www.eadianqi.com版权所有 if(n2>0) 自动控制网www.eadianqi.com版权所有 disp('thesystemisunstable') 自动控制网www.eadianqi.com版权所有 disp('theunstablepoleare:') 本文来自www.eadianqi.com
disp(p(jj)) 自动控制网www.eadianqi.com版权所有 else 本文来自www.eadianqi.com disp('thesystemisstable') 自动控制网www.eadianqi.com版权所有
end 本文来自www.eadianqi.com %判断系统是否为最小相位系统 自动控制网www.eadianqi.com版权所有 if(n1>0) 本文来自www.eadianqi.com disp('thesystemisanonminimalphaseone') 本文来自www.eadianqi.com else 自动控制网www.eadianqi.com版权所有 disp('thesyetemisaminimalphaseone') 自动控制网www.eadianqi.com版权所有 end 自动控制网www.eadianqi.com版权所有 %绘制零极点图 自动控制网www.eadianqi.com版权所有 pzmap(p,z) 自动控制网www.eadianqi.com版权所有 p 本文来自www.eadianqi.com
z 本文来自www.eadianqi.com ii=find(条件式) 自动控制网www.eadianqi.com版权所有
用来求取满足条件的向量的下标向量,以列向量表示。 本文来自www.eadianqi.com 例如exp4_1.m中的条件式为real(p>0),其含义就是找出极点向量p中满足实部的值大于0的所有元素下标,并将结果返回到ii向量中去。这样如果找到了实部大于0的极点,则会将该极点的序号返回到ii下。如果最终的结果里ii的元素个数大于0,则认为找到了不稳定极点,因而给出系统不稳定的提示,若产生的ii向量的元素个数为0,则认为没有找到不稳定的极点,因而得出系统稳定的结论。 自动控制网www.eadianqi.com版权所有 pzmap(p,z) 本文来自www.eadianqi.com 根据系统已知的零极点p和z绘制出系统的零极点图 自动控制网www.eadianqi.com版权所有 |