Sistema de segundo orden Ejemplo

Dado un sistema de segundo orden, realizar unas gráficas en 3D y 2D donde se vea claramente la diferencia entre un sistema sin amortiguamiento, subamortiguado, críticamente amortiguado y sobreamortiguado

¿Como se hace en matlab?

%----------------------------------------------------------------------------
%REG.AUTOMATICA Y MATLAB
%En este ejemplo variaremos el parametro zeta en un sist.de 2ºorden
%----------------------------------------------------------------------------
clear all;

%Declaracion de los parametros:
wn=1;%la frecuencia sera fija
zeta=[0 0.5 1 2];% Variaremos zeta
t=0:15/200:(15-15/200);


%Por cada parametro zeta almacenaremos la respuesta
%del sistema en una matriz y:
for n=1:length(zeta)
    num=wn^2;
    den=[1 2*zeta(n)*wn wn^2];
    W=tf(num,den);
    y(:,n)=step(W,t);
end;

%Representacion en 3D de la respuesta
figure(1);
mesh(zeta,t,y);
xlabel('\zeta');
ylabel('t');
zlabel('y(t)');
title('Respuesta de un sistema de segundo orden variando \zeta');
%Representacion en 2D de la respuesta
figure(2)
plot(t,y);xlabel('t');ylabel('y(t)');
title('Respuesta de un sistema de segundo orden variando \zeta');
legend('\zeta=0 (sin amortiguamiento)',...
'\zeta=0.5 (Subamortiguado)',...
'\zeta=1 (Criticamente Amortiguado)','\zeta=2 (Sobreamortiguado)');