%  Exemplo 16.1
clear all
close all
clc

tau = 1;
H = tf([1],[1 1]); % Função de transferência
bode(H)     % Plotando os diagramas de Bode de módulo (dB) e fase (graus)
grid on

M = [];     % Inicializando a variável de módulo
Phi = [];   % Inicializando a variável de fase
for w = 0:0.05:10
    M = [M, 1/sqrt(1+w^2)];         % Calculando o módulo
    Phi = [Phi, -atan(w)*180/pi];   % Calculando a fase
end

w = 0:0.05:10;

figure
plot(w,M)   % Plotando o módulo
grid on
title('Módulo da Resp. em Freq.')
xlabel('\omega')
ylabel('M(\omega)')

figure
plot(w,Phi) % Plotando a fase
grid on
title('Fase da Resp. em Freq.')
xlabel('\omega')
ylabel('\Phi(\omega)')