基于Matlab的控制系统稳定性判定

Matlab(400037):;Bode;MatlabBode,Bode,Matlab:;;Matlab;:TP273:B:1004373X(2005)2107302ApplicationinDecidingofControlSystemStabilityBasedonMatlabHUANGChongfu(ChongqingVocationalInstituteofEngineering,Chongqing,400037,China)Abstract:Thepaperintroducesthecontrolsystemstabilityidea1Itexpoundsthedifferenceandrelationofthealgebrastabilitycriterion,thefootlocuslawandBodeschartintheclassicscontroltheory1ItintroducesLyapunovslawapplicationinthecontemporarycontroltheoryindecidingofcontrolsystemstability,proposesinthewayofMatlabsoftwaretoapplyinthealgebrastabilitycriterion,footlocuslaw,BodeschartandLyapunovslawindecidingofcontrolsystemstability1Itdecidesthemeanssimplyastothedifferentindecidingofcontrolsystemstability1Algebrastabilitycriterion,footlocuslaw,BodeschartandLyapunovslaw,theapplicationgivesoutcorrelationinstance,orderandrelevantfigure1ExpoundingtodecidethestabilitymeansonthebasisofMatlabscomputationalsolution1Keywords:controlsystem;stability;Matlab;Lyapunovsfunction:200508051,,,,SISOBode,,2:,,,,,2Matlab211,,Matlabroots(P):G(S)=100(S+2)S(S+1)(S+20):n1=[100200];d1=conv([10],conv([11],[120]));sys1=tf(n1,d1);sys=feedback(sys1,1);roots(sys1den{1}):ans=37200521212þü-1218990-510000-311010,212,S()rlocus(sys)SISOLTIG(S),K(K0),SYS=G(S)1+KG(S),/[k,poles]=rlocfind(sys),,,K,poles,(),G(S)=K(S+2)S(S+1)(S+20):n1=[12];d1=conv([10],conv([11],[120]));G=tf(n1,d1);rlocus(G)[k,poles]=rlocfind(G)11Matlab,,Kpoles:k=11215352poles=-919772-813073-217155,,213BodeBodeC0,,Matlabmargin()BodeG(S)=100(S+2)S(S+1)(S+20):num=[100200];den=conv([10],conv([11],[120]));s1=tf(num,den);[Gm,Pm,Wcp,Wcg]=margin(s1)margin(s1)22:C=65140,214x=f(x,t),f(0,t)=0tt0V(x,t),:V(x,t);x=f(x,t)V(x,t)x=0V(x,t)x=Ax,x=0:Q(Q),P,ATP+PA=-QV(x)=xTPxMatlablyap(A,C)P,V(x)=xTPx(76)47:Matlab,,215Ls,(010V),AöD,C4C[5],,1024,,,1,,2min,,,1024,,,33:voidexint0(void)interrupt0using0{temp=(uint)(ADHI4)+(ADLO&0x0f);save_data[temp]=save_data[temp]+1;return;}save_data[1024],1024,,,54;,,,,[1],,,1[M]1:,19971[2],189C51[J]1,2003,(3):61621[3]1[M]1:,20011[4]1[M]1:,20021[5]1CMCS51[M]1:,20031,1977,,1949(74):G(S)=100õ(S+2)S(S+1)(S+20),:num=[100200]den=conv([10],conv([11],[120]))q=[100;010;001]sys1=tf(num,den)sys2=feedback(sys1,1)sys=ss(sys2):ans=010299-010156-010889-010156013557-010625-010889-0106252153043,,,Matlab[1],,1[M]1:,20031[2]1[M]1:,19981[3]1Matlab[M]1:,20011,1996,67: