% Este comentario se puede ver con help
% H(s) = define la funcion como quieras

% Este no lo vas a poder ver, jaja

w = logspace(-1,4,100);  % define un espacio logaritmico

% variable de laplace [j es sqrt(-1)]
s = j * w;

H1 = (s + 10) .* (s + 1) ./ (s + 100) ./ (2*s + 2000); % ./ divide por elemento

H = H1;
Hdb = 20 * log10(abs(H));

clg; % borrar figura

%subplot(211)   % figuras 2 una arriba, otra abajo
figure(1)
semilogx(w,Hdb)
%legend('Hdb','Hdb+10');
grid on
xlabel('Frecuencia (rad/s)')
ylabel('H_2(jw) (dB)')
title('Diagrama de Bode')

%subplot(212);
figure(2)
Hang = angle(H)*180/pi;
semilogx(w,Hang)
grid on
%text(10,10,'o - aqui');
xlabel('Frecuencia (rad/s)')
ylabel('Angulo de H(j\omega) (^o)')