Laboratori Om at Lab 1
Laboratori Om at Lab 1
Laboratori Om at Lab 1
Para obtener la respuesta de un sistema en el tiempo ante una entrada estándar, debe primero
definirse el sistema. Para ello puede definirse en MatLab la función de transferencia propia
del sistema o las ecuaciones de estado.
1. Función de Transferencia:
La función de transferencia de un sistema es una relación formada por un numerador y un
denominador:
En MatLab debe definirse el numerador Y(s) y el denominador X(s) como vectores, cuyos
elementos son los coeficientes de los polinomios del numerador y del denominador en
potencias decrecientes de S. Por ejemplo, para definir la función de transferencia:
Solucion:
Para el primer punto, se tiene que definor las variables tanto del numerador como del
denominador como se se muestra en la siguiente figura tomada de Matlab:
Figura N° 1, definición de variables.
2. Respuesta al escalón:
Para determinar la respuesta en el tiempo para una entrada escalón unitario de este sistema se
usa el comandos step indicando el vector del numerador y del denominador entre paréntesis:
step(num,den)
>>y=[1];
>>x=[1 0.25 1];
>>step(y,x)
Solucion:
Después de definir y escribir las variables, se tiene que, el paso siguiente es imprimir la
función, como resultado se obtiene la siguiente figura:
Figura N° 2, Impresión de la función de transferencia anteriormente definida.
3. Respuesta al impulso:
Para obtener la respuesta en el tiempo para una entrada impulso unitario se usa el comando
impulse, con sintaxis idéntica a la utilizada con el comando step:
Si se define el sistema en MatLab por los polinomios del numerador y denominador de la
función de transferencia tenemos:
» y=[1 5 4];
» x=[1 6 11 6];
» impulse(y,x)
MatLab permite, además de obtener la respuesta en el tiempo para una entrada escalón o
impulso, también obtener respuesta para otras entradas tal como rampas o sinusoides. El
comando lsim permite obtener la respuesta en el tiempo para un sistema con una entrada x,
donde x se define como una función del tiempo. La sintaxis de este comando es:
usando la función de transferencia.
Para obtener la respuesta en el tiempo para una función rampa, se define X de la siguiente
forma:
>>T=0:0.1:10
>>X=T;
>>NUM=[1];
>>DEN=[1 0.25 1];
>>[Y,X]=lsim(NUM,DEN,X,T);
>>PLOT(T,Y,T,X)
Al hacer X=T se está definiendo la función rampa. T es el vector de tiempo variando desde 0
hasta 10 seg. NUM y DEN son los vectores de los coeficientes decrecientes en potencia de S
de los polinomios del numerador y del denominador respectivamente. En la variable Y se
almacena la salida del sistema en función del tiempo T. El comando plot permite desplegar en
la ventana de figuras la variable Y (salida) y la entrada X (rampa) en función del tiempo.
80
60
40
20
-20
-10 -8 -6 -4 -2 0 2 4 6 8 10
t