%Datos acontrol=35
datos=load('datosPos35.txt');
figure
AccionControl=(datos(100,2));
sim('ModeloPosicion');
tiempo=0:0.05:4.95;
plot(tiempo,datos(:,1))
hold
plot(tiempo,RespuestaModelo())
plot(tiempo,RespuestaModeloSI())
title('Datos posición con acción de control 35')
xlabel('Tiempo(s)')
ylabel('Posición(rad)')
legend('Posición real','Modelo analítico','Modelo Sys. Ident.')

%Datos acontrol=50
datos=load('datosPos50.txt');
figure
AccionControl=(datos(100,2));
sim('ModeloPosicion');
tiempo=0:0.05:4.95;
plot(tiempo,datos(:,1))
hold
plot(tiempo,RespuestaModelo())
plot(tiempo,RespuestaModeloSI())
title('Datos posición con acción de control 50')
xlabel('Tiempo(s)')
ylabel('Posición(rad)')
legend('Posición real','Modelo analítico','Modelo Sys. Ident.')

%Datos acontrol=80
datos=load('datosPos80.txt');
figure
AccionControl=(datos(100,2));
sim('ModeloPosicion');
tiempo=0:0.05:4.95;
plot(tiempo,datos(:,1))
hold
plot(tiempo,RespuestaModelo())
plot(tiempo,RespuestaModeloSI())
title('Datos posición con acción de control 80')
xlabel('Tiempo(s)')
ylabel('Posición(rad)')
legend('Posición real','Modelo analítico','Modelo Sys. Ident.')

%Datos acontrol=-40
datos=load('datosPos-40.txt');
figure
AccionControl=(datos(100,2));
sim('ModeloPosicion');
tiempo=0:0.05:4.95;
plot(tiempo,datos(:,1))
hold
plot(tiempo,RespuestaModelo())
plot(tiempo,RespuestaModeloSI())
title('Datos posición con acción de control -40')
xlabel('Tiempo(s)')
ylabel('Posición(rad)')
legend('Posición real','Modelo analítico','Modelo Sys. Ident.')