Ingeniero Electronico
PRIMER SISTEMA
sys =
75
---------------
s^2 + 34 s + 75
Primero se definioel tiempo para el impulso, el cual va de 0 a 18, mediante el comando impulse
Figura 1. Respuesta al impulso
Luego se graficó la respuesta al escalon, con el comando step, paraun tiempo de 0 a 14.
Figura 2. Respuesta al escalón.
Posteriormente se hizo lo mismo para una entrada senoidal, con el comando lsim
Figura 3. Respuesta entrada senoidalLuego con simulink se procedió a comparar los datos obtenidos
Figura 4. Respuesta al impulso en simulink
Figura 5. Respuesta al escalón en simulink
Figura 6. Respuesta conentrada senoidal en simulink
Figura 7. Respuesta al impulso con datos traídos mediante workspace
Figura 8. Respuesta al escalón mediante el workspace
Figura 9. Respuestaentrada senoidal mediante workspace
Luego de ver los datos, y graficas obtenidas, se concluye que el comportamiento del sistema es bastante similar, exceptuando la respuesta alimpulso, ya que la gráfica dada mediante simulink tiene una amplitud de 15, mientras que la obtenida por matlab y la gráfica del workspace llegan hasta 1.
CODIGO
clc
close allclear all
num=75;
den=[1, 34,75];
sys=tf(num,den)
t1=[0:.1:18];
impulse(sys);
A=impulse(sys,t1); grid on
figure, plot(t1,A)
title('Respuesta al impulso')
xlabel('tiempo(segundos.)')
ylabel('Amplitud')
step(sys);
t2=[0:.1:14];
B=step(sys,t2);
figure, plot(t2,B); grid on
title('Respuesta al escalon')
xlabel('tiempo (segundos.)')ylabel('Amplitud')
t3=[0:0.2:30];
u=sin(2*pi*t3/3);
C=lsim(sys,u,t3);
figure, plot(t3,C); grid on
title('Respuesta a una entrada seno')
xlabel('tiempo (segundos.)')
ylabel('Amplitud')
Regístrate para leer el documento completo.