ESTABILIDAD POR JURY Ejercicios PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 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