%% RC pasabajos

%defino escala horizontal lineal
x=linspace(0,20,100);

% defino las funciones transferencia, desfasaje y atenuacion
T=@(x) 1./sqrt(1+(x).^2);
phi=@(x) -atan(x)*180/pi;
A=@(x) 20*log10(T(x));

%grafico transferencia y desfasaje vs frecuencia
figure(1);clf
subplot(2,1,1)
plot(x,T(x))
xlabel('\omega/\omega_0')
ylabel('T')
grid on
subplot(2,1,2)
plot(x,phi(x))
xlabel('\omega/\omega_0')
ylabel('\phi [grados]')
grid on

% diagrama de bode, en atenuacion y desfasaje en escala logaritmica
%defino escala horizontal logaritmica
x=logspace(-2,4,100);

figure(2);clf
subplot(2,1,1)
semilogx(x,A(x))
xlabel('\omega/\omega_0')
ylabel('A')
grid on
subplot(2,1,2)
semilogx(x,phi(x))
xlabel('\omega/\omega_0')
ylabel('\phi [grados]')
grid on

% valores en la frecuencia de corte
fprintf('T(wc)=%f\n',T(1))
fprintf('phi(wc)=%f\n',phi(1))
fprintf('A(wc)=%f\n',A(1))