⚡ electrical engineering

Fractional Order Filters with √f Characteristic

Cet article dérive des filtres d'ordre fractionnaire avec des fonctions de transfert en $\sqrt{f}$ et $1/\sqrt{f}$ en discrétisant spatialement une équation aux dérivées partielles issue des équations de Maxwell pour modéliser l'effet de peau, et présente une implémentation numérique correspondante via une transformation bilinéaire qui est évaluée par rapport au filtre récursif d'Oustaloup.

Auteurs originaux : Hans Georg Brachtendorf

Publié 2026-09-14

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Hans Georg Brachtendorf

Article original sous licence CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Dans le monde caché de l'électricité circulant à travers le métal, il existe un phénomène particulier qui modifie le comportement des signaux en fonction de leur vitesse. Lorsqu'un courant électrique voyage à travers un conducteur, il ne remplit pas toujours le fil de manière uniforme. À des vitesses très élevées, le courant est repoussé vers la surface extérieure, un comportement connu sous le nom d'effet de peau. Cet effet n'est pas seulement une curiosité ; il modifie fondamentalement la manière dont les signaux électriques sont filtrés et traités. Les ingénieurs cherchent depuis longtemps des moyens de créer des composants électroniques qui imitent une relation mathématique spécifique où l'intensité du signal change proportionnellement à la racine carrée de sa fréquence. De tels composants sont essentiels pour nettoyer le bruit dans les mesures sensibles, contrôler des systèmes mécaniques complexes et modéliser l'interaction des tissus biologiques avec les ondes électromagnétiques. Pendant des décennies, les meilleurs outils pour y parvenir ont été des approximations complexes qui ne fonctionnent bien que dans des plages limitées, nécessitant souvent des arrangements complexes de pôles et de zéros pour aplatir la réponse de phase.

Un chercheur de l'Université des Sciences Appliquées de Haute-Autriche a maintenant proposé une voie différente, qui ne part pas d'approximations de circuits, mais des lois fondamentales de la physique elle-même. En prenant les équations qui décrivent comment les champs électriques et magnétiques interagissent à l'intérieur d'une tige métallique et en les décomposant en petites étapes gérables, l'auteur a dérivé une nouvelle façon de construire ces filtres spéciaux. Au lieu de forcer un circuit standard à agir comme un filtre racine carrée, cette méthode construit un filtre qui se comporte naturellement de cette façon car il est une représentation numérique directe de l'effet de peau physique. Le travail démonte qu'en espaçant soigneusement ces étapes de calcul — en les plaçant plus près les unes des autres près de la surface de la tige là où l'action se produit — on peut créer un modèle hautement précis qui fonctionne sur une plage de fréquences beaucoup plus large que les méthodes précédentes.

Le voyage commence par une image physique simple : une longue tige métallique rectangulaire. Le chercheur suppose que le matériau est uniforme et que les champs électriques et magnétiques à l'intérieur suivent les règles standard de l'électromagnétisme, spécifiquement les équations de Maxwell. Dans le régime où l'effet de peau est fort, la profondeur à laquelle le courant pénètre le métal est très petite par rapport à l'épaisseur de la tige. Cette réalité physique dicte une relation spécifique entre la fréquence et la perte de signal. Pour transformer cette loi physique continue en quelque chose qu'un ordinateur peut utiliser, le chercheer divise la tige en une série de fines tranches. Ce processus, connu sous le nom de discrétisation spatiale, transforme les équations fluides et continues de la physique en une chaîne d'équations différentielles ordinaires. Ces équations peuvent ensuite être mappées directement sur un circuit électrique composé de résistances et de condensateurs, créant un modèle comportemental qui imite la réponse de la tige à un signal entrant.

Une intuition cruciale dans ce travail réside dans la manière dont les tranches sont disposées. Si les tranches étaient toutes de même taille, le modèle nécessiterait un nombre impraticable d'entre elles pour capturer les changements rapides se produisant juste à la surface du métal. Au lieu de cela, le chercheur utilise une grille non uniforme, où les trages deviennent progressivement plus petites à mesure qu'elles approchent de la surface. Ce regroupement des points de calcul reflète la réalité physique selon laquelle les changements les plus significatifs se produisent dans la couche la plus mince du conducteur. L'étude démontre que le rapport entre les tailles de ces tranches adjacentes est critique. Si les tranches augmentent trop rapidement en taille à mesure que l'on s'éloigne de la surface, le modèle devient instable et produit des effets de résonance indésirables. Les simulations suggèrent que le rapport des espacements adjacents doit rester inférieur à une constante mathématique spécifique, une valeur d'environ 2,5 semblant offrir le meilleur équilibre entre précision et plage de fréquences couverte.

Une fois le modèle analogique établi, l'étape suivante consiste à le traduire dans un format numérique capable de fonctionner sur des processeurs modernes. Le chercheur a testé plusieurs méthodes standards pour convertir des systèmes à temps continu en systèmes à temps discret. Une méthode populaire, connue sous le nom de transformation bilinéaire, est apparue comme la grande gagnante. Contra fait aux autres techniques qui introduisent un amortissement artificiel ou nécessitent des calculs d'ordre supérieur qui compliquent la mise en œuvre matérielle, la transformation bilinéaire préserve la stabilité et la précision du modèle physique original. Le filtre numérique résultant se comporte exactement comme le modèle analogique de l'effet de peau, maintenant la relation racine carrée souhaitée sur un large spectre. Dans les simulations, cette nouvelle approche s'est révélée supérieure au filtre récursif d'Oustaloup, une technique standard pour approximer les systèmes d'ordre fractionnaire. Alors que la méthode d'Oustaloup place des pôles et des zéros sur une échelle logarithmique pour aplatir la réponse de phase, le modèle d'effet de peau dérivé ici atteint naturellement une précision de phase plus élevée et reste valide sur une plage de fréquences plus large que celle du modèle d'Oustaloup.

Les résultats des simulations sont frappants. Pour un modèle utilisant seulement huit points de grille, le nouveau filtre maintient son comportement caractéristique sur une étendue d'environ 6 500 fois la fréquence de départ. C'est une amélioration significative par rapport aux méthodes existantes, qui peinent souvent à maintenir la précision sur des bandes aussi larges. L'étude confirme que l'intuition physique derrière l'effet de peau fournit une base robuste pour la conception de ces filtres, évitant ainsi le recours à des algorithmes d'optimisation complexes. L'auteur note que si le présent travail se concentre sur le cas spécifique de la racine carrée de la fréquence, la méthode sous-jacente de discrétisation des équations aux dérivées partielles pourrait potentiellement être étendue à d'autres ordres fractionnaires à l'avenir. Pour l'instant, ce travail offre une alternative claire et physiquement ancrée pour les ingénieurs ayant besoin d'un contrôle précis sur le traitement du signal, prouvant que parfois, la meilleure façon de résoudre un problème numérique est de regarder attentivement la physique du monde réel.

Résumé technique : Filtres d'ordre fractionnaire avec caractéristique en $\sqrt{f}$

Énoncé du problème

L'article traite de la conception et de l'analyse de filtres d'ordre fractionnaire caractérisés par une fonction de transfert proportionnelle à $s^{\alpha}$ où $\alpha = \pm 1/2$ (correspondant aux comportements en fréquence $\sqrt{f}$ et $1/\sqrt{f}$ ). Ces filtres sont essentiels pour des applications telles que la préaccentuation/déaccentuation, la modélisation des spectres de puissance de bruit $1/f$ (flicker), le contrôle d'ordre fractionnaire, et la simulation des propriétés dispersives de milieux comme les tissus biologiques.

Alors que les techniques d'approximation existantes, telles que le filtre récursif d'Oustaloup et la méthode de Charef, sont largement utilisées, elles reposent sur le placement de pôles et de zéros sur une échelle logarithmique pour aplatir les réponses de phase ou d'amplitude. L'article cherche une dérivation alternative ancrée dans les lois physiques plutôt que dans une pure approximation mathématique, visant à dériver un modèle comportemental qui présente naturellement les caractéristiques d'ordre fractionnaire souhaitées grâce à la discrétisation des équations de Maxwell.

Méthodologie

1. Dérivation physique à partir des équations de Maxwell

La méthodologie centrale commence par les équations de Maxwell pour une tige conductrice linéaire, isotrope et homogène. Sous l'hypothèse que le courant de déplacement est négligeable ( $J \gg \partial D/\partial t$ ) et en utilisant une approche quasi-stationnaire pour une tige rectangulaire où la longueur $l \gg$ la largeur $b$ et l'épaisseur $d$ , les équations se réduisent à un système d'équations aux dérivées partielles (EDP) décrivant les champs électrique ( $E_z$ ) et magnétique ( $H_x$ ).

L'article exploite le phénomène physique connu selon lequel la profondeur de peau ( $\Delta x$ ) d'un matériau conducteur suit une loi en $1/\sqrt{f}$ dans le régime d'effet de peau fort (où la profondeur de peau est petite par rapport à la géométrie). L'objectif est de discrétiser ces EDP pour créer un filtre qui imite ce comportement.

2. Discrétisation spatiale et modélisation de circuits

Les EDP sont discrétisées à l'aide d'un schéma de différences finies (FD).

Grille uniforme : Initialement, une grille équidistante est considérée, menant à un système d'équations différentielles-algébriques ordinaires (DAE). Ce système est mappé vers des modèles de circuits comportementaux passifs (réseaux en échelle) composés d'inductances et de résistances (LR) ou de condensateurs et de conducteurs (GC).
Grille non uniforme : Reconnaissant qu'une grille uniforme nécessite un nombre impraticable de points pour capturer l'effet de peau à haute fréquence, l'article propose une grille non uniforme. L'espacement de la grille est minimal à la surface ( $y = d/2$ ) là où l'effet de peau est le plus prononcé. Le rapport des espacements adjacents est défini par un paramètre $c$ .
Réalisation de circuit : Les équations discrétisées sont réalisées sous forme de circuits passifs équivalents. L'article note que bien que le modèle LR soit naturel, le modèle GC (utilisant des condensateurs et des conducteurs) est plus adapté aux simulateurs de type SPICE car il réduit le nombre d'inconnues dans l'analyse nodale modifiée (MNA).

3. Transformation numérique

Pour obtenir un équivalent numérique, le filtre analogique est transformé à l'aide de la transformation bilinéaire (règle du trapèze) et comparé aux formules de différenciation vers l'arrière (BDF).

Transformation bilinéaire : Préférée car elle projette l'axe imaginaire sur le cercle unité sans introduire d'amortissement numérique, un effet secondaire observé avec les méthodes BDF. Elle permet également d'obtenir un ordre de filtre inférieur par rapport aux méthodes BDF d'ordre supérieur ( $p \ge 2$ ).
Implémentation : L'article traite des stratégies d'implémentation efficaces des équations aux différences résultantes, incluant l'introduction de vecteurs auxiliaires pour éviter les multiplicateurs dans certaines étapes récursives spécifiques.

Principales contributions

Dérivation physique de filtres fractionnaires : L'article dérive des filtres $\sqrt{f}$ et $1/\sqrt{f}$ directement de la discrétisation des équations de Maxwell pour une tige conductrice, fournissant une base physique au comportement d'ordre fractionnaire plutôt que de s'appuyer uniquement sur des heuristiques de placement de pôles et de zéros.
Optimisation de la grille non uniforme : Il démontre que l'utilisation d'une grille non uniforme avec un rapport d'espacement spécifique ( $c$ ) réduit considérablement l'ordre du système tout en préservant la précision. L'article identifie que le rapport des espacements de grille adjacents doit être inférieur à la constante d'Euler pour éviter les effets de pompage (ringing), un ratio de $c=2$ à $c=2.5$ étant jugé optimal.
Modèles de circuits comportementaux : La dérivation produit des modèles de circuits passifs en échelle (tant LR que GC) qui sont intrinsèquement stables, car les pertes ohmiques garantissent que les pôles restent dans le demi-plan gauche.
Implémentation numérique : L'article fournit une voie claire vers la réalisation numérique via la transformation bilinéaire, offrant une structure de filtre numérique stable et efficace.

Résultats

Comparaison avec Oustaloup : Les résultats de simulation pour un ordre de filtre $K=11$ montrent que le modèle de profondeur de peau dérivé dans cet article offre une plage de validité approximative et une précision de phase plus larges par rapport au filtre récursif d'Oustaloup.
Plage de fréquences : Pour une grille non uniforme avec $K=8$ points et $c=2$ , le modèle présente le comportement $\sqrt{f}$ attendu (pente linéaire de 0,5 dans un graphique log-log de l'impédance) sur environ quatre décades de fréquence. Bien que les estimations théoriques suggéraient une plage plus large, les simulations indiquent une plage de validité pratique d'environ $f_1/10$ .
Sensibilité des paramètres : Des tests confirment que les valeurs de $c$ entre 2 et la constante d'Euler ( $\approx 2,718$ ) offrent le meilleur compromis entre l'ordre du système, la plage de validité et l'exactitude. Des valeurs de $c$ dépassant la constante d'Euler entraînent des effets de pomping indésirables.
Performance numérique : Les homologues numériques, implémentés via la transformation bilinéaire, maintiennent le comportement de racine carrée pendant environ quatre décades de fréquence angulaire normalisée, confirmant l'efficacité de la transformation.

Signification et affirmations

L'article affirme que la méthode proposée offre une alternative supérieure aux approximations récursives existantes comme celle d'Oustaloup, particulièrement en ce qui concerne la plage de validité et la précision de phase. En ancrant la conception du filtre dans l'effet de peau physique des conducteurs, les modèles résultants ne sont pas seulement mathématiquement rigoureux, mais aussi physiquement interprétables comme des réseaux passifs.

L'auteur souligne que la stratégie de discrétisation non uniforme permet de modéliser précisément l'effet de peau avec nettement moins de points de grille (et donc un coût de calcul moindre) que ce que nécessiterait une discrétisation uniforme. La préférence pour la transformation bilinéaire est justifiée par sa capacité à préserver la stabilité et à éviter l'amortissement numérique, ce qui la rend adaptée à une implémentation matérielle.

Enfin, l'article positionne ce travail comme un cas spécifique de filtrage d'ordre fractionnaire ( $\alpha = \pm 1/2$ ) et suggère que la méthodologie sous-jacente de discrétisation des équations aux dérivées partielles pourrait être généralisée dans des travaux futurs pour des filtres d'ordre fractionnaire avec des ordres $\alpha$ arbitraires au-delà de $\pm 1/2$ .

Résumé technique : Filtres d'ordre fractionnaire avec caractéristique en f\sqrt{f}f​