% Exemplo 16.2 - Diagrama assintótico do circuito RC
clear all
close all
clc

w = logspace(-2, 2);
for i = 1:length(w)
    MdB(i) = mag2db((1/sqrt(1+w(i)^2))); 
end
figure
semilogx(w,MdB,'LineWidth',2);
axis([10^(-2) 10^2 -60 20])
grid on
hold on

MAb = [0 0];    % Assíntota da magnitude para baixa frequência
plot([10^(-2) 1],MAb,'r','LineWidth',2);

MAa = [0 -40];  % Assíntota da magnitude para alta frequência
plot([1 10^2],MAa,'r','LineWidth',2);

title('Diagrama assintótico de módulo')
xlabel('\omega')
ylabel('M_{dB}(\omega)')


w = logspace(-2, 2);
for i = 1:length(w)
    Phi(i) = -atan(w(i))*180/pi; 
end
figure
semilogx(w,Phi,'LineWidth',2);
axis([10^(-2) 10^2 -130 50])
grid on
hold on

PAb = [0 0]; % Assíntota da fase para baixa frequência
plot([10^(-2) 10^(-1)],PAb,'r','LineWidth',2);

PAwc = [0 -90]; % Assíntota da fase para frequencias próximas a wc
plot([10^(-1) 10^1],PAwc,'r','LineWidth',2);

PAa = [-90 -90]; % Assíntota da fase para baixa frequência
plot([10^1 10^2],PAa,'r','LineWidth',2);
title('Diagrama assintótico de fase')
xlabel('\omega')
ylabel('\Phi(\omega)')
