Tarea 4: Curva de Respuesta

Para obtener la salida del sistema se usa la función lsim, describida a continuación.

Function File: y= lsim (sys, u, t)

Parámetros

• sys: Sistema invariante en el tiempo (LTI).
• u: Vector de señal de entrada. Debe ser del mismo largo de filas que el vector de tiempos y el mismo largo que entradas de columnas.
• t: Vector de tiempos espaciado uniformemente.

Salida

y: Vector de respuesta de salida. Tiene tantas filas como el largo de t y tantas columnas como salidas.

Proceso:

Creamos un vector con los digitos del numerador, y un vector con los digitos del denominador y creamos un sistema con ellos.
 

	num = [0 0 10 10];
	den = [1 4 0 0];
	sys = tf(num, den);

view raw sys.m hosted with ❤ by GitHub

Después creamos un vector de tiempos uniformes(0.1 de espacio) y obtenemos el vector entrada de aceleración unitaria en base a éste:

	t = 0:0.1:5;
	r = 0.5*t.^2;

view raw vect.m hosted with ❤ by GitHub

Ahora podemos utilizar la función lsim para obtener la salida del sistema, dando como parametros el sistema, el vector de tiempos y vector de entrada.

y = lsim(sys, r, t);

view raw lsim.m hosted with ❤ by GitHub

Después graficamos los resultados:

	plot(t, r, "-",t, y, "o");
	grid;
	title('Respuesta a la Aceleracion unitaria del sistema');
	xlabel('t segundos');
	ylabel('Aceleracion Unitaria');
	text(0.7, 4.5, 'Curva de Respuesta');

view raw plot.m hosted with ❤ by GitHub

Resultado: