\documentclass[svgnames, border=2mm]{standalone}
\usepackage[spanish]{babel}
\usepackage[T1]{fontenc}
\usepackage[utf8]{inputenc}
\usepackage{pgfplots, tgpagella, verbatim}
\usepgfplotslibrary{units}
\pgfplotsset{compat=newest}
\begin{comment}
El presente es un gráfico doble de una Función de transferencia de un circuíto RLC en serie definida por la relación entre el voltaje en la resistencia ($V_R$) y el voltaje en la fuente ($V_S$). La ecuación obtenida es $H(j\omega)=\frac{j\omega RC}{LC(j\omega)^2+RCj\omega+1}$. En el listado, se definen los valores de R, L y C en función de los cuales se traza la gráfica de la ganancia (en rojo) y del desfase (en azul), para el rango de frecuencias considerado.
La función de transferencia fue obtenida de la web "Guía práctica para construir un diagrama de Bode" http://panamahitek.com/guia-practica-para-construir-un-diagrama-de-bode/ y se aprovecha el uso de GNUPLOT para graficar la función de transferencia, la cual depende de variables complejas.
\end{comment}
\begin{document}
\begin{tikzpicture}
\begin{semilogxaxis}[
% Configuración del eje de los gráficos
axis background/.style={fill=Linen}, % Fondo del gráfico
cycle list name=color list,
grid=major,
every tick label/.append style={font=\textbf\footnotesize},
legend style={font=\tiny, legend pos=north east, draw=none, fill=none}, % Estilo de la leyenda
grid=both,
line width=0.25mm,
minor tick num=9,
title style={text width=7cm},
title=\textbf{\centering \color{DarkCyan}Función de Transferencia \\ de un Filtro Pasabanda\\},
xlabel=\small{\color{DarkCyan} Frecuencia Angular $\omega$},
x unit=\small{\color{DarkCyan} rads/s},
ylabel=\small{\color{DarkCyan} $H(j\omega)$},
]
% Valores de los componentes pasivos del filtro
\def\R{2000}; % 2000 Ohmios
\def\C{1e-8}; % 10 nanoFaradios
\def\L{5.5e-3}; % 5.5 milihenrios
\foreach \i in {0, 1}
\addplot+[no marks, smooth, id=transfer] gnuplot [raw gnuplot]{
set xr [1e4:2.3e6];
h(x) = (x * \R * \C * {0, 1})/(-\L * \C * x^2 + \R * \C * x * {0, 1} + 1); % Función de Transferencia
f0(x) = abs(h(x)); % Ganancia
f1(x) = arg(h(x)); % Fase en radianes
plot f\i(x)
};
\legend{$|H(j\omega)| \,(Voltios)$, $\phi(H(j\omega)) \, (radianes)$}
\end{semilogxaxis}
\end{tikzpicture}
\end{document}
