{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Síntesis de funciones de excitación no disipativas\n",
"![](\"./img/logo_UTN.svg\")
\n",
"\n",
"#### Por Mariano Llamedo Soria"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Resumen \n",
"Se presentan los métodos de Foster y Cauer para la síntesis de redes no disipativas. Ambos métodos nos permitirán implementar cualquier función de excitación, real y positiva (FRP), en una red eléctrica canónica, es decir, con la menor cantidad posible de componentes.\n",
"\n",
"* Dibujo de redes: [dibujar_foster_derivacion](https://pytc2.readthedocs.io/en/latest/autoapi/pytc2/dibujar/index.html#pytc2.dibujar.dibujar_foster_derivacion), [dibujar_foster_serie](https://pytc2.readthedocs.io/en/latest/autoapi/pytc2/dibujar/index.html#pytc2.dibujar.dibujar_foster_serie)\n",
"* Funciones de síntesis de dipolos: [foster](https://pytc2.readthedocs.io/en/latest/autoapi/pytc2/sintesis_dipolo/index.html#pytc2.sintesis_dipolo.foster)\n",
"* Funciones para presentación de markdown y latex: [print_subtitle](https://pytc2.readthedocs.io/en/latest/autoapi/pytc2/general/index.html#pytc2.general.print_subtitle), [print_latex](https://pytc2.readthedocs.io/en/latest/autoapi/pytc2/general/index.html#pytc2.general.print_latex), [a_equal_b_latex_s](https://pytc2.readthedocs.io/en/latest/autoapi/pytc2/general/index.html#pytc2.general.a_equal_b_latex_s)\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Introducción a las funciones de excitación no disipativas\n",
"\n",
"Una función de excitación que representa una inmitancia (una impedancia ó admitancia) $F(s)$ será no disipativa si cumple que:\n",
"\n",
"* $F(s)=\\frac{P(s)}{Q(s)}$ sea una función real positiva (FRP)\n",
"* $\\Re \\{F\\}=0$\n",
"\n",
"La primer condición es necesaria para que la inmitancia sea realizable, ver Apéndice E de libro de Cogollos Borrás, S. para más detalles. La segunda restricción implica que \n",
"\n",
"$$ F = \\frac{m_P}{n_Q} \\lor F = \\frac{n_P}{m_Q} $$\n",
"\n",
"siendo $m$ la parte par y $n$ la parte impar de sendos polinomios numerador y denominador (P y Q respect.).\n",
"\n",
" [1] Cogollos Borrás, Santiago. [*FUNDAMENTOS DE LA TEORÍA DE FILTROS*](https://www.lalibreria.upv.es/portalEd/UpvGEStore/products/p_287-4-1). ISBN: 978-84-9048-443-2\n",
"\n",
"### Un teorema de reactancias (1924): la propuesta de Foster\n",
"\n",
"Queda definir la forma propuesta por Foster en [*A reactance theorem*](https://en.wikipedia.org/wiki/Foster%27s_reactance_theorem):\n",
"\n",
"$$ F(s)=\\frac{P(s)}{Q(s)} = \\frac{k_0}{s} + k_\\infty.s + \\sum_{i=1}^N\\frac{2.k_i.s}{s^2+\\omega_i^2} $$\n",
"\n",
"donde la inmitancia F puede representarse como una suma de términos, o fracciones simples, que individualmente están asociandas a cada elemento circuital de la red. Teniendo la función $F(s)$ N polos finitos, la red resultante de la expansión Foster tendrá $2.N + 2$ componentes. Es decir, los residuos $k$ serán los valores de los componentes, como se observa que $k_0$ y $k_\\infty$ serán sendos capacitor e inductor, mientras que cada polo finito cuyo residuo es $k_i$, se asocian a tanques resonantes L-C. Justamente la frecuencia de resonancia de cada tanque será $\\omega_i$. En los ejemplos siguientes quedará evidenciada esta asociación.\n",
"\n",
"Se comienza por cargar las librerías."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"\n",
"import sympy as sp\n",
"\n",
"# Ahora importamos las funciones de PyTC2\n",
"\n",
"from pytc2.sintesis_dipolo import foster\n",
"from pytc2.dibujar import dibujar_foster_serie, dibujar_foster_derivacion\n",
"from pytc2.general import print_latex, print_subtitle, a_equal_b_latex_s\n",
"\n",
"# Importante importar símbolos de variables \n",
"from pytc2.general import s\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Ahora se definen las funciones que usaremos en los ejemplos"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle F(s)=\\frac{2 s^{4} + 20 s^{2} + 18}{s^{3} + 4 s}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"# Sea la siguiente función de excitación\n",
"FF = (2*s**4 + 20*s**2 + 18)/(s**3 + 4*s)\n",
"\n",
"print_latex(a_equal_b_latex_s('F(s)', FF))\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Expansión en fracciones simples: Método de Foster\n",
"\n",
"Se procede a la expansión Foster"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle k_0=\\frac{9}{2}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/latex": [
"$\\displaystyle k_1 = \\left[ \\frac{1}{ \\frac{1}{s. \\frac{\\omega_i^2}{2.k_i} } + s . \\frac{1}{2.k_i} } \\right] = \\\n",
" \\left[ \\frac{1}{ \\frac{k_0}{s} + s . k_\\infty } \\right] = \\\n",
" \\left[ k_0, k_\\infty \\right] = \\\n",
" \\left[ \\\n",
" \\left[ \\frac{\\omega_1^2}{2k_1}, \\frac{1}{2k_1} \\right] \\\n",
" \\right]=\\left[ \\left[ \\frac{8}{15}, \\ \\frac{2}{15}\\right]\\right]$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/latex": [
"$\\displaystyle k_\\infty=2$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/latex": [
"$\\displaystyle F(s)=\\frac{2 s^{4} + 20 s^{2} + 18}{s^{3} + 4 s}=2 s + \\frac{15 s}{2 \\left(s^{2} + 4\\right)} + \\frac{9}{2 s}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Se expande FF a la Foster\n",
"k0, koo, ki_wi, _, FF_foster = foster(FF)\n",
"\n",
"print_latex(a_equal_b_latex_s('k_0', k0))\n",
"\n",
"print_latex(a_equal_b_latex_s(r'k_1 = \\left[ \\frac{1}{ \\frac{1}{s. \\frac{\\omega_i^2}{2.k_i} } + s . \\frac{1}{2.k_i} } \\right] = \\\n",
" \\left[ \\frac{1}{ \\frac{k_0}{s} + s . k_\\infty } \\right] = \\\n",
" \\left[ k_0, k_\\infty \\right] = \\\n",
" \\left[ \\\n",
" \\left[ \\frac{\\omega_1^2}{2k_1}, \\frac{1}{2k_1} \\right] \\\n",
" \\right]', ki_wi ))\n",
"\n",
"print_latex(a_equal_b_latex_s('k_\\infty', koo))\n",
"\n",
"\n",
"print_latex(a_equal_b_latex_s(a_equal_b_latex_s('F(s)', FF)[1:-1], FF_foster ))\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Luego queda solo implementar la inmitancia como impedancia o admitancia. Para ello se dispone de funciones de dibujo de redes:"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/markdown": [
"#### Foster serie"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/latex": [
"$\\displaystyle Z(s)=F(s)=\\frac{2 s^{4} + 20 s^{2} + 18}{s^{3} + 4 s}=2 s + \\frac{15 s}{2 \\left(s^{2} + 4\\right)} + \\frac{9}{2 s}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"print_subtitle('Foster serie')\n",
"\n",
"print_latex(a_equal_b_latex_s(a_equal_b_latex_s('Z(s)=F(s)', FF)[1:-1], FF_foster ))\n",
"\n",
"# Tratamos a nuestra función imitancia como una Z\n",
"dibujar_foster_serie(k0, koo, ki_wi, z_exc = FF)\n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Como se puede ver, la suma de términos simples, resulta en que cada término será una impedancia asociada en serie. Luego si se asume F(s) como admitancia Y(s), la red resultante estará asociada en paralelo."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/markdown": [
"#### Foster derivación"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/latex": [
"$\\displaystyle Y(s)=F(s)=\\frac{2 s^{4} + 20 s^{2} + 18}{s^{3} + 4 s}=2 s + \\frac{15 s}{2 \\left(s^{2} + 4\\right)} + \\frac{9}{2 s}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"print_subtitle('Foster derivación')\n",
"\n",
"print_latex(a_equal_b_latex_s(a_equal_b_latex_s('Y(s)=F(s)', FF)[1:-1], FF_foster ))\n",
"\n",
"# Tratamos a nuestra función imitancia como una Y\n",
"dibujar_foster_derivacion(k0, koo, ki_wi, y_exc = FF)\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Como puede verse, el método de Foster da lugar a redes con $2.1 + 2 = 4$ componentes reactivos en ambos casos. Se dice que es un método **canónico** dado que 4 es la menor cantidad de componentes con la que la función F(s) puede ser implementada con elementos pasivos.\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Expansión en fracciones continuas: Método de Cauer\n",
"\n",
"Del mismo modo que el método de Foster se sutenta en la expansión en fracciones simples, el de Cauer consiste en la expansión en [fracciones continuas](https://en.wikipedia.org/wiki/Network_synthesis#Cauer_synthesis). Con la salvedad que puede representarse la expansión mediante residuos $k_0$ y $k_\\infty$.\n",
"\n",
"$$ F(s)= \\frac{k_0}{s} + \\frac{1}{ \\frac{k_1}{s} + \\frac{1}{ \\frac{k_2}{s} + \\cdots } } = \n",
" k_{\\infty,1}.s + \\frac{1}{ k_{\\infty,2}.s + \\frac{1}{ k_{\\infty,3}.s + \\cdots } } $$\n",
"\n",
"en este caso cualquiera de ambas expansiones da lugar a redes escalera, por ejemplo si se tratara de una impedancia:\n",
"\n",
"$$ Z(s)= \\frac{1}{s.C_1} + \\frac{1}{ \\frac{1}{s.L_1} + \\frac{1}{ \\frac{1}{s.C_2} + \\cdots } } = \n",
" s.L_1 + \\frac{1}{ s.C_1 + \\frac{1}{ s.L_2 + \\cdots } } $$\n",
"\n",
"o de una admitancia:\n",
"\n",
"$$ Y(s)= \\frac{1}{s.L_1} + \\frac{1}{ \\frac{1}{s.C_1} + \\frac{1}{ \\frac{1}{s.L_2} + \\cdots } } = \n",
" s.C_1 + \\frac{1}{ s.L_1 + \\frac{1}{ s.C_2 + \\cdots } } $$\n",
"\n",
"se ve claramente como una función de inmitancia da lugar a 4 redes diferentes. Analizaremos la siguiente función de inmitancia\n"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle F(s)=\\frac{2 s^{4} + 20 s^{2} + 18}{s^{3} + 4 s}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from pytc2.sintesis_dipolo import cauer_LC\n",
"from pytc2.dibujar import dibujar_cauer_LC\n",
"\n",
"\n",
"# Sea la siguiente función de excitación\n",
"FF = (2*s**4 + 20*s**2 + 18)/(s**3 + 4*s)\n",
"\n",
"print_latex(a_equal_b_latex_s('F(s)', FF))\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\n",
"## Cauer 1: remociones en $\\infty$\n",
"\n",
"Se comienza con la primera forma del método, es decir cuando se realizan remociones de los residuos en infinito."
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle F(s)=\\frac{2 s^{4} + 20 s^{2} + 18}{s^{3} + 4 s}=2 s + \\frac{1}{\\frac{s}{12} + \\frac{1}{\\frac{24 s}{5} + \\frac{36}{5 s}}}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"# Implementaremos FF mediante Cauer 1 o remociones continuas en infinito\n",
"koo, F_cauer_oo, rem = cauer_LC(FF, remover_en_inf=True)\n",
"\n",
"print_latex(a_equal_b_latex_s(a_equal_b_latex_s('F(s)', FF)[1:-1], F_cauer_oo ))\n",
"\n",
"# Tratamos a nuestra función inmitancia como una Z\n",
"dibujar_cauer_LC(koo, z_exc = F_cauer_oo)\n",
"\n",
"# Tratamos a nuestra función inmitancia como una Y\n",
"dibujar_cauer_LC(koo, y_exc = F_cauer_oo)\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Se observa, como es de esperarse, que interpretar la inmitancia como impedancia o admitancia utilizando el mismo método, Cauer 1 en este caso, da lugar a redes **duales**.\n",
"\n",
"## Cauer 2: remociones en 0 o DC\n",
"\n",
"Repetimos ahora para remociones de los residuos en 0 Hz, obtendremos dos redes también duales."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$\\displaystyle F(s)=\\frac{2 s^{4} + 20 s^{2} + 18}{s^{3} + 4 s}=\\frac{1}{\\frac{1}{\\frac{62 s}{15} + \\frac{961}{30 s}} + \\frac{8}{31 s}} + \\frac{9}{2 s}$"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"# Implementaremos F mediante Cauer 2 o remociones continuas en cero\n",
"k0, F_cauer_0, rem = cauer_LC(FF, remover_en_inf=False)\n",
"\n",
"print_latex(a_equal_b_latex_s(a_equal_b_latex_s('F(s)', FF)[1:-1], F_cauer_0 ))\n",
"\n",
"# Tratamos a nuestra función inmitancia como una Z\n",
"dibujar_cauer_LC(k0, z_exc = F_cauer_0)\n",
"\n",
"# Tratamos a nuestra función inmitancia como una Y\n",
"dibujar_cauer_LC(k0, y_exc = F_cauer_0)\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Ambas redes duales, tienen características en común como ser que ambas tienen capacitores en serie e inductores en derivación. Esto contrasta con las redes halladas mediante Cauer 1, donde observamos lo contrario: inductores en serie y capacitores en derivación. Está claro que Cauer 1 da lugar a redes que, en caso **que imaginemos un puerto de salida en el extremo derecho de la red**, serían filtros **pasabajo**, mientras que sintetizar una red mediante Cauer 2 da lugar a una red escalera **pasa-altos**.\n",
"\n",
"Finalmente, decir que habrá redes que al tener la misma función de exitación, cargarán de forma idéntica al puerto de entrada, pero en caso de asumir una eventual transferencia respecto a un puerto en el otro extremo de la red, la transferencia **será complementaria** (pasabajo-pasaalto):"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/svg+xml": [
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Tratamos a nuestra función inmitancia como una Z\n",
"dibujar_cauer_LC(koo, z_exc = F_cauer_oo)\n",
"\n",
"# Tratamos a nuestra función inmitancia como una Z\n",
"dibujar_cauer_LC(k0, z_exc = F_cauer_0)\n"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.3"
},
"varInspector": {
"cols": {
"lenName": 16,
"lenType": 16,
"lenVar": 40
},
"kernels_config": {
"python": {
"delete_cmd_postfix": "",
"delete_cmd_prefix": "del ",
"library": "var_list.py",
"varRefreshCmd": "print(var_dic_list())"
},
"r": {
"delete_cmd_postfix": ") ",
"delete_cmd_prefix": "rm(",
"library": "var_list.r",
"varRefreshCmd": "cat(var_dic_list()) "
}
},
"types_to_exclude": [
"module",
"function",
"builtin_function_or_method",
"instance",
"_Feature"
],
"window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 4
}