Respuesta escalón unitario con MATLAB Ejemplo

Tenemos dos sistemas con las siguientes funciones de transferencia:

sys1: y sys2:.

Realizar con MATLAB una gráfica donde veamos la respuesta de los dos sistemas ante un escalón unitario con un tiempo de simulación de 30s. También representar en la misma gráfica, la función escalón unitario

¿Como se hace en matlab?

%----------------------------------------------------------------------
%REG.AUTOMATICA Y MATLAB
%En este ejemplo veremos el uso del comando step
%----------------------------------------------------------------------

%Definicion de los sistemas:
sys1=tf([1],[1 0.5 1]);
sys2=tf([1],[1 0.5 4]);
%Representacion de la respuesta:
t=0:0.01:30; %Respuesta hasta los 30 s.
step(sys1,'r', sys2,'g',t);%Representacion en la misma grafica
%Aplicaremos rejilla y pondremos un titulo con text:
grid
text(5, 1.4,'Respuesta de dos sistemas','FontSize',13);
%Representacion de la entrada
t0 = -2.0:0.01:-0.01; % definicion u(t)=0, -2<=t<=-.01
u0 = zeros(size(t0));
t1 = 0:0.01:30; % definicion u(t)=1, 0<=t<=25
u1 = ones(size(t1));
t = [t0 t1]; % creamos t and u(t)
u = [u0 u1];
hold on
plot(t,u);
legend('Sistema 1','Sistema 2','Escalon unitario');