% 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

w0 = 1;
xi = 1/(2*w0);

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

for i=1:length(w)
    Mw(i) = w0^2/sqrt((w0^2-w(i)^2)^2 + 4*xi^2*w0^2*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

Mwr = max(Mw);
ind = find(Mw==Mwr);
wr = w(ind);

fprintf('M(wr) = %1.4f\n',Mwr)
fprintf('wr = %1.4f\n',wr)

hold on
plot(wr,Mwr,'ro')




% Forma alternativa - Matlab (usando a função de transferência)

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

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

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


    