clc
%Lab.2
%Name: Chávez Choque Alvaro
%Adress: Pedro Ferrari #4, Toledo, Madrid
%Thelephone: 5264639
%I.N. 5775394 or.
%e-mail: alvaro.chavez_choque@hotmail.com
'1 Hallar la respuesta en el tiempo de la funcion dada'
% determinar G(s)
syms s G U1 U0
num=[-1]
den=[0.188 0.87 1]
G=tf(num,den)
'Respuesta grafica del sistema'
% ajuste de la ganancia
G1=-5*G
step(G1)
% para que simule 20 seg
subplot(1,1,1)
step(G1,25)
grid
title('respuesta del sistema')
'Hallar la respuesta analitica'
syms g ui uo s
g=(-1)/(0.188*s^2+0.87*s+1)
pretty(g)
vi=(-5)/s
v0=g*vi
pretty(v0)
pretty(simplify(ilaplace(v0)))
'2 Polos y ceros de las funciones de transferencia'
[ceros,polos,ganancia]=tf2zp(num,den)
'3 Expansión en fracciones parciales'
nume=-5
dene=[1 0]
Ui=tf(nume,dene)
U0=G*Ui
[nums,dens]=tfdata(U0,'v')
[residue,polos,directo]=RESIDUE(nums,dens)
Haciendo correr el programa vemos que el command window no muestra lo siguiente.
1 Hallar la respuesta en el tiempo de la funcion dada
num =
-1
den =
0.1880 0.8700 1.0000 Transfer function:
-1
----------------------
0.188 s^2 + 0.87 s + 1
ans =
Respuesta grafica del sistema
Transfer function:
5
----------------------
0.188 s^2 + 0.87 s + 1
ans =
Hallar la respuesta analitica
g =
-1/(47/250*s^2+87/100*s+1)
1
- ------------------
47 2 87
--- s + --- s + 1
250 100
vi =
-5/s
v0 =
5/(47/250*s^2+87/100*s+1)/s
5
----------------------
/47 2 87 \
|--- s + --- s + 1| s
\250 100 /
435 35 435 35
5 - 5 exp(- --- t) cosh(--- t) - 435/7 exp(- --- t) sinh(--- t)
188 188 188 188
ans =
No hay comentarios:
Publicar un comentario