Respuesta ante una señal cualquiera con MATLAB Ejemplo

Sea un sistema con la siguiente función de transferencia: realizar:
a) Representar la respuesta del sistema ante la entrada 
b) Representar la respuesta del sistema para una señal cuadrada con un periodo de 5s.,     durante 30s.
c) Representar la respuesta del sistema ante una entrada rampa unitaria. Representar     también la entrada rampa unitaria

¿Como se hace en matlab?

%-----------------------------------------------------------------------------
%REG.AUTOMATICA Y MATLAB
%En este ejemplo veremos el uso de los comandos lsim y gensig
%-----------------------------------------------------------------------------
%Primero definimos el sistema 
num=[1 4];
den=[1 4 25];
sys=tf(num,den);
% A)
t=-5:0.01:30;% Establecemos el tiempo de simulacion en 30 s
u1=5*exp(-3*t);% Introducimos la primera funcion de entrada
subplot(2,2,1),lsim(sys,u1,t);%Aplicamos lsim con la funcion
legend('respuesta ante u(t)');
% B)
[u,t]=gensig('square',5,30,0.1);%generamos la onda cuadrada 
subplot(2,2,2),lsim(sys,u,t);%y dibujamos la respuesta
legend('respuesta ante onda cuadrada');

% C)
%Esta es otra manera de conseguir la respuesta ante una rampa
u=t;
subplot(2,2,3),lsim(sys,u,t,'r')%Dibujamos la respuesta en rojo
hold
plot(t,t,'g');%Esta es la rampa unitaria, en verde
legend('respuesta del sistema','entrada rampa unitaria')