⚡ electrical engineering

Fractional Order Filters with √f Characteristic

Este artículo deriva filtros de orden fraccionario con funciones de transferencia de $\sqrt{f}$ y $1/\sqrt{f}$ mediante la discretización espacial de una ecuación diferencial parcial de las ecuaciones de Maxwell para modelar el efecto pelicular, y presenta una implementación digital correspondiente vía transformación bilinear que se compara con el filtro recursivo de Oustaloup.

Autores originales: Hans Georg Brachtendorf

Publicado 2026-09-14

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Hans Georg Brachtendorf

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

En el mundo oculto de la electricidad que fluye a través del metal, existe un fenómeno peculiar que cambia la forma en que se comportan las señales dependiendo de su velocidad. Cuando una corriente eléctrica viaja a través de un conductor, no siempre llena todo el cable de manera uniforme. A velocidades muy altas, la corriente es empujada hacia la superficie exterior, un comportamiento conocido como efecto de piel (skin effect). Este efecto no es solo una curiosidad; altera fundamentalmente cómo se filtran y procesan las señales eléctricas. Los ingenieros han buscado durante mucho tiempo formas de crear componentes electrónicos que imiten una relación matemática específica donde la fuerza de la señal cambia en proporción a la raíz cuadrada de su frecuencia. Tales componentes son vitales para limpiar el ruido en mediciones sensibles, controlar sistemas mecánicos complejos y modelar cómo los tejidos biológicos interactúan con las ondas electromagnéticas. Durante décadas, las mejores herramientas para lograr esto han sido aproximaciones complejas que funcionan bien solo dentro de rangos limitados, requiriendo a menudo arreglos intrincados de polos y ceros para aplanar la respuesta de fase.

Un investigador de la Universidad de Ciencias Aplicadas de Alta Austria ha propuesto ahora un camino diferente, uno que no comienza con aproximaciones de circuitos, sino con las leyes fundamentales de la física misma. Al tomar las ecuaciones que describen cómo interactúan los campos eléctricos y magnéticos dentro de una barra de metal y desglosarlas en pasos pequeños y manejables, el autor ha derivado una nueva forma de construir estos filtros especiales. En lugar de forzar a un circuito estándar a actuar como un filtro de raíz cuadrada, este método construye un filtro que se comporta naturalmente de esa manera porque es una representación digital directa del efecto de piel físico. El trabajo demuestra que, al espaciar cuidadosamente estos pasos de cálculo —colocándolos más cerca unos de otros cerca de la superficie de la barra donde ocurre la acción—, se puede crear un modelo altamente preciso que funciona en un rango de frecuencias mucho más amplio que los métodos anteriores.

El viaje comienza con una imagen física simple: una barra de metal larga y rectangular. El investigador asume que el material es uniforme y que los campos eléctricos y magnéticos dentro de él siguen las reglas estándar del electromagnetismo, específicamente las ecuaciones de Maxwell. En el régimen donde el efecto de piel es fuerte, la profundidad a la que la corriente penetra en el metal es muy pequeña en comparación con el grosor de la barra. Esta realidad física dicta una relación específica entre la frecuencia y la pérdida de la señal. Para convertir esta ley física continua en algo que una computadora pueda usar, el investigador divide la barra en una serie de pequeñas rebanadas. Este proceso, conocido como discretización espacial, transforma las ecuaciones fluidas de la física en una cadena de ecuaciones diferenciales ordinarias. Estas ecuaciones pueden entonces mapearse directamente en un circuito eléctrico compuesto por resistencias y capacitores, creando un modelo de comportamiento que imita la respuesta de la barra a una señal entrante.

Un conocimiento crucial en este trabajo es cómo se disponen las rebanadas. Si las rebanadas fueran todas del mismo tamaño, el modelo requeriría un número impráctico de ellas para capturar los cambios rápidos que ocurren justo en la superficie del metal. En su lugar, el investigador utiliza una rejilla no uniforme, donde las rebanadas se vuelven progresivamente más pequeñas a medida que se acercan a la superficie. Este agrupamiento de puntos de cálculo refleja la realidad física de que los cambios más significativos ocurren en la capa más delgada del conductor. El estudio encuentra que la relación entre los tamaños de estas rebanadas adyacentes es crítica. Si las rebanadas crecen demasiado rápido en tamaño a medida que se alejan de la superficie, el modelo se vuelve inestable y produce efectos de oscilación (ringing) no deseados. Las simulaciones sugieren que la relación de los espaciamientos adyacentes debe mantenerse por debajo de una constante matemática específica, y un valor alrededor de 2.5 parece ofrecer el mejor equilibrio entre la precisión y el rango de frecuencias cubierto.

Una vez establecido el modelo analógico, el siguiente paso es traducir este a un formato digital que pueda ejecutarse en procesadores modernos. El investigador probó varios métodos estándar para convertir sistemas de tiempo continuo en sistemas de tiempo discreto. Un método popular, conocido como transformación bilineal, surgió como el claro ganador. A diferencia de otras técnicas que introducen un amortiguamiento artificial o requieren cálculos de orden superior que complican la implementación del hardware, la transformación bilineal preserva la estabilidad y la precisión del modelo físico original. El filtro digital resultante se comporta exactamente como el modelo de efecto de piel analógico, manteniendo la relación de raíz cuadrada deseada a través de un espectro amplio. En las simulaciones, este nuevo enfoque demostró ser superior al filtro recursivo de Oustaloup, una técnica estándar para aproximar sistemas de orden fraccionario. Mientras que el método de Oustaloup coloca polos y ceros en una escala logarítmica para aplanar la respuesta de fase, el modelo de efecto de piel derivado aquí logra naturalmente una mayor precisión de fase y sigue siendo válido en un rango de frecuencias que es mayor que el del modelo de Oustaloup.

Los resultados de las simulaciones son sorprendentes. Para un modelo que utiliza solo ocho puntos de rejilla, el nuevo filtro mantiene su comportamiento característico sobre un lapso de aproximadamente 6,500 veces la frecuencia inicial. Esta es una mejora significativa respecto a los métodos existentes, que a menudo luchan por mantener la precisión a través de bandas tan amplias. El estudio confirma que la intuición física detrás del efecto de piel proporciona una base robusta para diseñar estos filtros, evitando la necesidad de algoritmos de optimización complejos. El autor señala que, si bien el presente trabajo se centra en el caso específico de la raíz cuadrada de la frecuencia, el método subyacente de discretización de ecuaciones diferenciales parciales podría potencialmente extenderse a otros órdenes fraccionarios en el futuro. Por ahora, el trabajo ofrece una alternativa clara y físicamente fundamentada para los ingenieros que necesitan un control preciso sobre el procesamiento de señales, demostrando que, a veces, la mejor manera de resolver un problema digital es observar de cerca la física del mundo real.

Resumen Técnico: Filtros de Orden Fraccionario con Característica $\sqrt{f}$

Planteamiento del Problema

El artículo aborda el diseño y análisis de filtros de orden fraccionario caracterizados por una función de transferencia proporcional a $s^{\alpha}$ donde $\alpha = \pm 1/2$ (correspondiente a los comportamientos de frecuencia $\sqrt{f}$ y $1/\sqrt{f}$ ). Estos filtros son esenciales para aplicaciones tales como la preénfasis/desénfasis, el modelado de los espectros de potencia de ruido $1/f$ (ruido flicker), el control de orden fraccionario y la simulación de las propiedades dispersivas de medios como los tejidos biológicos.

Si bien las técnicas de aproximación existentes, como el filtro recursivo de Oustaloup y el método de Charef, son ampliamente utilizadas, estas se basan en la colocación de polos y ceros en una escala logarítmica para aplanar las respuestas de fase o magnitud. El artículo busca una derivación alternativa basada en leyes físicas en lugar de una pura aproximación matemática, con el objetivo de derivar un modelo de comportamiento que exhiba naturalmente las características de orden fraccionario deseadas mediante la discretización de las ecuaciones de Maxwell.

Metodología

1. Derivación Física a partir de las Ecuaciones de Maxwell

La metodología central comienza con las ecuaciones de Maxwell para una barra conductora lineal, isotrópica y homogénea. Bajo el supuesto de que la corriente de desplazamiento es despreciable ( $J \gg \partial D/\partial t$ ) y utilizando un enfoque cuasiestacionario para una barra rectangular donde la longitud $l \gg$ ancho $b$ y espesor $d$ , las ecuaciones se reducen a un sistema de ecuaciones diferenciales parciales (PDE) que describen los campos eléctrico ( $E_z$ ) y magnético ( $H_x$ ).

El artículo aprovecha el fenómeno físico conocido de que la profundidad de piel ( $\Delta x$ ) de un material conductor sigue una ley de $1/\sqrt{f}$ en el régimen de efecto de piel fuerte (donde la profundidad de piel es pequeña en comparación con la geometría). El objetivo es discretizar estas PDE para crear un filtro que imite este comportamiento.

2. Discretización Espacial y Modelado de Circuitos

Las PDE se discretizan utilizando un esquema de diferencias finitas (FD).

Grilla Uniforme: Inicialmente, se considera una grilla equidistante, lo que conduce a un sistema de ecuaciones diferenciales ordinarias algebraicas (DAE). Este sistema se mapea a modelos de circuitos de comportamiento pasivo (redes en escalera) que consisten en inductores y resistencias (LR) o capacitores y conductores (GC).
Grilla No Uniforme: Reconociendo que una grilla uniforme requiere un número impráctico de puntos para capturar el efecto de piel a altas frecuencias, el artículo propone una grilla no uniforme. El espaciamiento de la grilla es más pequeño en la superficie ( $y = d/2$ ) donde el efecto de piel es más pronunciado. La relación de los espaciamientos adyacentes está definida por un parámetro $c$ .
Realización de Circuitos: Las ecuaciones discretizadas se realizan como circuitos pasivos equivalentes. El artículo señala que, si bien el modelo LR es natural, el modelo GC (usando capacitores y conductores) es más adecuado para simuladores tipo SPICE, ya que reduce el número de incógnitas en el Análisis de Nodos Modificados (MNA).

3. Transformación Digital

Para obtener un equivalente digital, el filtro analógico se transforma utilizando la transformación bilineal (regla trapezoidal) y se compara contra las Fórmulas de Diferenciación hacia Atrás (BDF).

Transformación Bilineal: Se prefiere porque mapea el eje imaginario al círculo unitario sin introducir amortiguamiento numérico, un efecto secundario observado en los métodos BDF. También resulta en un orden de filtro menor en comparación con los métodos BDF de orden superior ( $p \ge 2$ ).
Implementación: El artículo discute estrategias de implementación eficientes para las ecuaciones de diferencia resultantes, incluyendo la introducción de vectores auxiliares para evitar multiplicadores en pasos recursivos específicos.

Contribuciones Clave

Derivación Física de Filtros Fraccionarios: El artículo deriva filtros de $\sqrt{f}$ y $1/\sqrt{f}$ directamente de la discretización de las ecuaciones de Maxwell para una barra conductora, proporcionando una base física para el comportamiento de orden fraccionario en lugar de depender únicamente de heurísticas de colocación de polos y ceros.
Optimización de Grilla No Uniforme: Demuestra que el uso de una grilla no uniforme con un parámetro de relación de espaciamiento específico ( $c$ ) reduce significativamente el orden del sistema preservando la precisión. El artículo identifica que la relación de los espaciamientos de grilla adyacentes debe ser menor a la constante de Euler para evitar efectos de oscilación (ringing), encontrando que $c=2$ a $c=2.5$ es óptimo.
Modelos de Circuitos de Comportamiento: La derivación produce modelos de circuitos en escalera pasivos (tanto LR como GC) que son inherentemente estables, ya que las pérdidas óhmicas aseguran que los polos permanezcan en el semiplano izquierdo.
Implementación Digital: El artículo proporciona un camino claro para la realización digital mediante la transformación bilineal, ofreciendo una estructura de filtro digital estable y eficiente.

Resultados

Comparación con Oustaloup: Los resultados de simulación para un orden de filtro de $K=11$ muestran que el modelo de profundidad de piel derivado en este trabajo ofrece un rango de validez de frecuencia aproximada mayor y una mayor precisión de fase en comparación con el filtro recursivo de Oustaloup.
Rango de Frecuencia: Para una grilla no uniforme con $K=8$ puntos y $c=2$ , el modelo exhibe el comportamiento esperado de $\sqrt{f}$ (pendiente lineal de 0.5 en un gráfico log-log de impedancia) durante aproximadamente cuatro décadas de frecuencia. Mientras que las estimaciones teóricas sugerían un rango más amplio, las simulaciones indican un rango de validez práctica de aproximadamente $f_1/10$ .
Sensibilidad de Parámetros: Las pruebas confirman que los valores del parámetro de relación de espaciamiento de la grilla $c$ entre 2 y la constante de Euler ( $\approx 2.718$ ) ofrecen el mejor equilibrio entre el orden del sistema, el rango de validez y la exactitud. Los valores de $c$ que exceden la constante de Euler resultan en efectos de oscilación indeseables.
Desempeño Digital: Los contrapartes digitales, implementados mediante la transformación bilineal, mantienen el comportamiento de raíz cuadrada durante aproximadamente cuatro décadas de frecuencia angular normalizada, confirmando la efectividad de la transformación.

Significado y Reivindicaciones

El artículo afirma que el método propuesto ofrece una alternativa superior a las aproximaciones recursivas existentes como la de Oustaloup, particularmente en cuanto al rango de validez y la precisión de fase. Al fundamentar el diseño del filtro en el efecto de piel físico de los conductores, los modelos resultantes no solo son matemáticamente sólidos, sino también físicamente interpretables como redes pasivas.

El autor enfatiza que la estrategia de discretización no uniforme permite modelar con precisión el efecto de piel con significativamente menos puntos de grilla (y, por tanto, un menor costo computacional) de lo que requeriría una discretización uniforme. La preferencia por la transformación bilineal se justifica por su capacidad para preservar la estabilidad y evitar el amortiguamiento numérico, lo que la hace adecuada para la implementación en hardware.

Finalmente, el artículo posiciona este trabajo como un caso específico de filtrado de orden fraccionario ( $\alpha = \pm 1/2$ ) y sugiere que la metodología subyacente de discretización de ecuaciones diferenciales parciales podría generalizarse en trabajos futuros para filtros de orden fraccionario con órdenes arbitrarios $\alpha$ más allá de $\pm 1/2$ .

Resumen Técnico: Filtros de Orden Fraccionario con Característica f\sqrt{f}f​