%Estabilidad segun método de Jury clc disp('Estabilidad disp('Estabilidad segun metodo de Jury'); Jury'); disp(' disp(' '); '); disp('Seleccione disp('Seleccione presionando'); presionando'); disp('1:Para disp('1:Para función de la planta en TF'); TF'); disp('2:Para disp('2:Para función de la planta en ZPK'); ZPK'); disp('3:Ecuación disp('3:Ecuación caracteristica del sistema'); sistema'); n=input('Seleccione n=input('Seleccione la opción'); opción'); disp(' disp(' '); '); switch n switch n case 1 case 1 num=input('Entre num=input('Entre coeficientes del numerador: num='); den=input('Entre den=input('Entre coeficientes del denominador:den='); denominador:den='); disp('La disp('La función de transferencia de la planta es: Gp(s)='); Gp=tf(num,den) clc disp(' disp(' '); '); T=input('Enter T=input('Enter tiempo de muestreo:T='); muestreo:T='); disp(' disp(' '); '); disp('La disp('La funcion de transferencia discreta de la planta: G(z)='); Gz = c2d(Gp,T,'zoh' c2d(Gp,T,'zoh') ) disp(' disp(' '); '); disp('La disp('La función de transferencia de lazo cerrado es: Glc(z)='); Glc=feedback(Gz,1) [numz,denz]=tfdata(Glc,'v' [numz,denz]=tfdata(Glc, 'v'); ); case 2 case 2 Z=input('Entre Z=input('Entre vector de ceros:Z='); ceros:Z='); P=input('Entre P=input('Entre vector de polos:P='); polos:P='); K=input('Ganancia K=input('Ganancia es igual a:K='); a:K='); disp('La disp('La función de transferencia de la planta es:Gp(s)='); es:Gp(s)='); Gp=zpk(Z,P,K) clc disp(' disp(' '); '); T=input('Enter T=input('Enter tiempo de muestreo:T='); muestreo:T='); disp(' disp(' '); '); disp('LA disp('LA función de transferencia discreta de la planta :G(z)='); Gz=c2d(Gp,T,'zoh' Gz=c2d(Gp,T,'zoh') ) disp(' disp(' '); '); disp('La disp('La funcion de transferencia lazo cerrado es: Glc(z)='); Glc(z)='); Glc=feedback(Gz,1) [numz,denz]=tfdata(Glc,'v' [numz,denz]=tfdata(Glc, 'v'); ); case 3 case 3 denz=input('coeficientes denz=input('coeficientes de la ecuación caracteristica:denz='); end disp(' disp(' '); '); clc disp('El disp('El orden de la ecuación caracteristica es:n=') es:n=') n=length(denz)-1 %Obtención de matriz de Jury disp(' disp(' '); '); m=n+1; for i=1:m for i=1:m
a(i)=denz(i); end pz11=polyval(denz,1); if abs(pz11)<1.0e-5 pz1=0; else pz1=pz11; end pz2=polyval(denz,-1); x=n-1; switch x %Sistema de segundo orden case 1 disp(''); disp('La matriz de Jury es:a2 a1 a0='); [a(3) a(2) a(1)] %prueba de estabilidad if(pz1==0)&(abs(a(m))
a(1) disp('el sistema es inestable'); elseif pz1<=0 disp('El sistema es inestable'); elseif pz2>0 disp('El sistema es estable'); else disp('El sistema es inestable'); end %Sistema de 3 orden case 2 disp(' '); disp('La matriz de jury es:'); disp(''); disp('a3 a2 a1 a0'); disp('a0 a1 a2 a3'); disp('b2 b1 b0'); for k=1:m-1 b(k)=det([a(m) a(m-k);a(1)a(k+1)]); end A1=[a(4) a(3) a(29 a(1)]; A2=[a(1) a(2) a(3) a(4)]; B=[b(3) b(2) b(1) 0]; [A1;A2;B] pause clc disp(' '); disp('Prueba de estabilidad:'); disp('a(n)='); a(m) disp('a(0)='); a(1) disp('Para z=1,P(z)='); [pz1] if(pz1==0)&(abs(a(m))
disp('El sistema es críticamente estable') end disp('Para z=-1,P(z)='); [pz2] pause clc disp(' '); if abs(a(m))>=a(1) disp('El sistema es inestable'); elseif pz1<0 disp('El sistema es inestable'); elseif pz2>=0 disp('El sistema es inestable'); elseif(abs(b(m-1))>abs(b(1))) disp('b(n-1)='); b(m-1) disp('b(0)'); b(1) disp('El sistema es estable'); else disp('El sistema es inestable'); end %Sistema de 4 orden
case 3 disp(' '); disp('La matriz de Jury es:'); disp(' '); disp('a4 a3 a2 a1 a0'); disp('a0 a1 a2 a3 a4'); disp('b3 b2 b1 b0'); disp('b0 b1 b2 b3'); disp('c2 c1 c0'); for k=1:m-1 b(k)=det([a(m) a(m-k);a(1) a(k+1)]); end for k=1:m-2 c(k)=det([b(m-1)b(m-1-k);b(1) b(k+1)]); end A1=[a(5) a(4) a(3) a(2) a(1)]; A2=[a(1) a(2) a(3) a(4) a(5)]; B1=[b(4) b(3) b(2) b(1) 0]; B2=[b(1) b(2) b(3) b(4)]; C=[c(3) c(2) c(1) 0 0]; [A1;A2;B1;B2;C] %Prueba de estabilida if abs(a(m))>a(1) disp('El sistema es criticamente estable'); elseif pz1<0 disp('El sistema es inestable'); elseif pz2<=0 disp('El sistema es inestable'); elseif (abs(b(m-1))<=abs(b(1))) disp('El sistema es inestable'); elseif(abs(c(m-2))>abs(c(1))) disp('El sistema es estable'); else disp('El sistema es inestable'); end end
