ESTABILIDAD POR JURY Ejercicios PDF
ESTABILIDAD POR JURY Ejercicios PDF
ESTABILIDAD POR JURY Ejercicios PDF
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
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
%Entrada Escalon
Ez = z/(z-1);
%SALIDAS Y(z)
Yaz = Haz*Ez;
Ybz = Hbz*Ez;
Ycz = Hcz*Ez;
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))
b)
VI = 0
c)
VI = 0
VF = 33441179994024359/88643061634600359 = 0.37
Graficas de salidas y(k) para los casos a) b) c):
a)
b)
c)