Scribd
Cargar
217 vistas9 páginas

ESTABILIDAD POR JURY Ejercicios PDF

Este documento presenta el estudio de estabilidad de tres sistemas de control (SCLC) aplicando el método de Jury. Determina que el primer sistema es inestable, mientras que los otros dos son estables. Calcula las respuestas temporales y(k) de cada sistema ante una entrada escalón unitario usando Matlab y grafica los resultados. También determina los valores inicial y final de cada respuesta.

Cargado por

alder anthony anahue
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
217 vistas9 páginas

ESTABILIDAD POR JURY Ejercicios PDF

Este documento presenta el estudio de estabilidad de tres sistemas de control (SCLC) aplicando el método de Jury. Determina que el primer sistema es inestable, mientras que los otros dos son estables. Calcula las respuestas temporales y(k) de cada sistema ante una entrada escalón unitario usando Matlab y grafica los resultados. También determina los valores inicial y final de cada respuesta.

Cargado por

alder anthony anahue
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
Está en la página 1/ 9

UNIVERSIDAD NACIONAL DE SAN AGUSTIN

FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS


ESCUELA PROFESIONAL DE INGENIERÍA ELECTRÓNICA

ESTABILIDAD POR JURY


GRUPO: A

DOCENTE:
Ing. Oscar Salazar Alarcón

INTEGRANTES:
Anahue Catunta Alder 20170297
Condorvica Cáceres Carlos Daniel 20153563
Mamani Huisa Carlos Esteban 20171900
Rojas Sanchez Jose Luis 20173168
Ydme Revilla Paul Cristian 20051429
Torres Aranibar Jorge Luis 20170312

06 de julio del 2020


Arequipa – Perú
0.7445 0.7446

La condición 7 no cumple, entonces


La E.C. : F(z) No tiene todas sus raíces
dentro del circulo de radio 1 en el plano z.
Por lo tanto el SCLC es Inestable

0.7445 0.7446
Comprobación de las raices con matlab

El valor absoluto de la primera raiz es 1.00013 y esta fuera del circulo de radio
1 del plano z, Por lo tanto el SCLC es Inestable.
Problema 2
 Estudiar estabilidad en LC aplicando el metodo JURY
 Hallar la salida temporal y(k)=? Y graficar para entrada e(k) escalon
unitario.
Usando Matlab para hallar y(k) para los casos a) b) y c) se obtiene:
%Problema 2
clc
clear
syms z k

%FT LAZO ABIERTO


Gaz = (0.031794*z-0.023558)/(z^2-1.18788*z+0.203635);
Gbz = (0.024468*z-0.016794)/(z^2-1.05232*z+0.066139);
Gcz = ((0.000276/z)+(0.0069754/(z^2)))/(z^2-1.116099*z+0.128069);

%FT LAZO CERRADO


Haz = Gaz/(1+Gaz);
Hbz = Gbz/(1+Gbz);
Hcz = Gcz/(1+Gcz);

%Entrada Escalon
Ez = z/(z-1);

%SALIDAS Y(z)
Yaz = Haz*Ez;
Ybz = Hbz*Ez;
Ycz = Hcz*Ez;

%VALORES INICIALES Y FINALES


VIa = limit(Yaz,z,inf)
VFa = limit((z-1)*Yaz,z,1)

VIb = limit(Ybz,z,inf)
VFb = limit((z-1)*Ybz,z,1)

VIc = limit(Ycz,z,inf)
VFc = limit((z-1)*Ycz,z,1)

%SALIDAS y(k)
yak = iztrans(Yaz,z,k)
ybk = iztrans(Ybz,z,k)
yck = iztrans(Ycz,z,k)

%Graficas de y(k)
k = 0:1:150;
ya = double(subs(yak,sym('k'),k));
yb = double(subs(ybk,sym('k'),k));
yc = double(subs(yck,sym('k'),k));
figure()
stem(k,ya),xlabel('k'),ylabel('ya(k) para entrada Escalon')
figure()
stem(k,yb),xlabel('k'),ylabel('yb(k) para entrada Escalon')
figure()
stem(k,yc),xlabel('k'),ylabel('yc(k) para entrada Escalon')
Salidas y(k) para los casos a) b) c):
a) y(k) =(178468751233900995215525595444952980649519918170716787109375*(-1)^k*900719925474099200000^(1 -
k)*31246412913483362178593725262600912967401^(1/2)*(- 2*31246412913483362178593725262600912967401^(1/2) -
520654847880824722598)^(k -
1))/703101785264276923549675830135568765434772586514790855196682625118372638682845684 -
(3031502297670592806562897551759730006250*(-1)^k*162198942019599371875^k*cos(k*(pi -
acos((260327423940412361299*51903661446271799^(1/2))/43539946963748718326579200000))))/(562546001049032787
1996539153069167996021*(1677721600000*51903661446271799^(1/2))^k) -
(178468751233900995215525595444952980649519918170716787109375*(-1)^k*900719925474099200000^(1 -
k)*31246412913483362178593725262600912967401^(1/2)*(2*31246412913483362178593725262600912967401^(1/2) -
520654847880824722598)^(k -
1))/703101785264276923549675830135568765434772586514790855196682625118372638682845684 +
7418329306204684375/21609171732049126679
b) y(k) = (57770704271392376019916954483177175741460082850039404296875*(-1)^k*900719925474099200000^(1 -
k)*696984507737163355267566597465901160371361^(1/2)*(- 696984507737163355267566597465901160371361^(1/2)/2
- 925806776838403810769/2)^(k -
1))/3122987873771756174155116023653701621864602068505314866759282198946678778474196451 -
(2899617604549879847068734466048874143750*(-1)^k*44446024722519428125^k*cos(k*(pi -
acos((925806776838403810769*14222727911206217^(1/2))/47723555655107104630374400000))))/(4480713472256189 47
2163051559547930226691*(1677721600000*14222727911206217^(1/2))^k) -
(57770704271392376019916954483177175741460082850039404296875*(-1)^k*900719925474099200000^(1 -
k)*696984507737163355267566597465901160371361^(1/2)*(696984507737163355267566597465901160371361^(1/2)/2 -
925806776838403810769/2)^(k -
1))/3122987873771756174155116023653701621864602068505314866759282198946678778474196451 +
1728031177022059375/4839793339553704339
c) y(k) = 33441179994024359/88643061634600359 - symsum(-(1733188924791366804193480759524719*r3*r3^k -
1075747738060524024345360526485600*r3^k - 17904813579305664689881199269102592*r3^k*r3^2 +
154220222218155815776645265317953536*r3^k*r3^3)/(- 456254538858335276953130305926448128*r3^3 +
104707669368126192737232442365149184*r3^2 + 338481405479891664833838093777843*r3 +
11406005776736407326545716179542400), r3 in RootOf(z3^4 - (5026463040508583*z3^3)/4503599627370496 +
(36048218792357*z3^2)/281474976710656 + (1272825341085959*z3)/4611686018427387904 +
1005261082904325/144115188075855872, z3))

Valore Iniciales y Finales para los casos a) b) c):


a)
VI = 0

VF = 7418329306204684375 / 21609171732049126679 = 0.34

b)
VI = 0

VF = 17280.35031177022059375 / 4839793339553704339 = 0.35

c)
VI = 0

VF = 33441179994024359/88643061634600359 = 0.37
Graficas de salidas y(k) para los casos a) b) c):
a)

b)
c)

También podría gustarte