% Esboce o módulo M(w) e a fase Phi(w) da função de transferência
% H(s)=1/(s+1)

clear all
close all
clc

w = 0:0.5:5;        % vetor de frequências

for i=1:length(w)
    Mw(i) = 1/sqrt(1+w(i)^2);   % Calcula o módulo de H(jw) pela expressão obtida
    Phiw(i) = -atan(w(i))*180/pi;      % Calcula a fase de H(jw) pela expressão obtida
end

% Desenha os gráficos:

figure(1)
plot(w,Mw)
xlabel('\omega')
ylabel('|H(j\omega)|')
grid on

figure(2)
plot(w,Phiw)
xlabel('\omega')
ylabel('\Phi(\omega)')
grid on


% Forma alternativa - Matlab

H = tf([1],[1 1]);

figure
BODE(H,{1e-5,5})
grid on

[MAG,PHASE,W] = BODE(H);


    