MoDelización y Análisis Temporal De Un Sistema Físico

Practica 2

Modelización y Análisis Temporal de un Sistema Físico



Laro de la Fuente Lastra
Oscar López Ruiz

Datos obtenidos teóricamente mediante Matlab:
4- Simular la transformadainversa calculada analíticamente y comparar dicha salida con la simulada por Matlab cuando la entrada al sistema es una entrada impulso. Obtener la velocidad y aceleración del sistema también conmatlab.
Apartado 4
t=0:0.01:20;
y=0.1.*(exp(-1/2.*t).*cos((3^(1/2)/2).*t));
figure(1), plot(t,y)
title('transformadainv')
grid;

ahora con una entrada impulso
y2=tf([0.1 0.1] , [1 1 1]);figure(2), impulse(y2,t)
title('entradaimpulso')
grid;

ahora calculo la velocidad
V=diff(y)./diff(t);
figure(3), plot(t(2:2001),V), title('velocidad'),grid;

ahora la aceleracionA=diff(diff(y)./diff(t));
figure(4), plot(t(3:2001), A)
title('aceleracion')
grid;

5- Repetir todo la anterior con una masa M=10Kg.
Apartado 5 m=10
syms s;;
t=0:0.01:100;
y3=(s+0.1)/(10*s^2+s+1);tiempo=ilaplace(y3);
y4=1/390.*exp(-1/20.*t).*(39.*cos(1/20.*39^(1/2).*t)+19.*39^(1/2).*sin(1/20.*39^(1/2).*t));
figure(5), plot(t,y4)

ahora con una entrada impulso
y5=tf([1 0.1] , [10 1 1]);
figure(6),impulse(y5,t),title('entradaimpulso'),grid;

ahora calculo la velocidad
V2=diff(y4)./diff(t);
figure(7), plot(t(2:10001),V2)
title('velocidad')
grid;

ahora la aceleracionA2=diff(diff(y4)./diff(t));
figure(8), plot(t(2:10000), A2)
title('aceleracion')
grid;

6- Repetir otra vez todo el procedimiento para M=1Kg, cuando la masa está en una posición inicial y0=0.05m y tiene una velocidadhacia arriba de 1m/s y la señal de excitación es cero.

Apartado 6, M=1kg, y0=0.05, va=1m/s
syms s;

t=0:0.01:100;
yss6=(0.05*s+1+0.05)/(s^2+s+1)
ys6=tf([(1/20) (21/20)],[1 1 1])yt6=ilaplace(yss6) % yt6 toma valores en función de t

figure(9)
impulse(ys6)
figure(10)
yt6val=subs(yt6,'t',t);
plot(t,yt6val)

V3=diff(yt6val)./diff(t);
figure(11)
size(t)
size(V3) %sale una...