⚡ electrical engineering

Fractional Order Filters with √f Characteristic

Questo articolo deriva filtri di ordine frazionario con funzioni di trasferimento $\sqrt{f}$ e $1/\sqrt{f}$ discretizzando spazialmente un'equazione differenziale parziale derivata dalle equazioni di Maxwell per modellare l'effetto pelle, e presenta una corrispondente implementazione digitale tramite trasformata bilineare che viene confrontata con il filtro ricorsivo di Oustaloup.

Autori originali: Hans Georg Brachtendorf

Pubblicato 2026-09-14

📖 5 min di lettura🧠 Approfondimento

CC BY 4.0

Autori originali: Hans Georg Brachtendorf

Articolo originale sotto licenza CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Questa è una spiegazione generata dall'IA dell'articolo qui sotto. Non è stata scritta né approvata dagli autori. Per precisione tecnica, consulta l'articolo originale. Leggi il disclaimer completo

Nel mondo nascosto dell'elettricità che scorre attraverso il metallo, esiste un fenomeno peculiare che cambia il modo in cui i segnali si comportano a seconda della loro velocità. Quando una corrente elettrica viaggia attraverso un conduttore, non sempre riempie l'intero filo in modo uniforme. Ad altissime velocità, la corrente viene spinta verso la superficie esterna, un comportamento noto come effetto pelle (skin effect). Questo effetto non è solo una curiosità; altera fondamentalmente il modo in cui i segnali elettrici vengono filtrati ed elaborati. Gli ingegneri cercano da tempo modi per creare componenti elettronici che imitino una specifica relazione matematica in cui l'intensità del segnale cambia in proporzione alla radice quadrata della sua frequenza. Tali componenti sono vitali per pulire il rumore nelle misurazioni sensibili, controllare sistemi meccanici complessi e modellare l'interazione dei tessuti biologici con le onde elettromagnetiche. Per decenni, i migliori strumenti per ottenere ciò sono stati complesse approssimazioni che funzionano bene solo entro intervalli limitati, richiedendo spesso intricati arrangiamenti di poli e zeri per appiattire la risposta di fase.

Un ricercatore presso l'Università delle Scienze Applicate dell'Alta Austria ha ora proposto una strada diversa, una che non parte dalle approssimazioni dei circuiti, ma dalle leggi fondamentali della fisica stessa. Prendendo le equazioni che descrivono come i campi elettrici e magnetici interagiscono all'interno di una barra metallica e scomponendole in piccoli passaggi gestibili, l'autore ha derivato un nuovo modo per costruire questi speciali filtri. Invece di costringere un circuito standard ad agire come un filtro a radice quadrata, questo metodo costruisce un filtro che si comporta naturalmente in quel modo perché è una rappresentazione digitale diretta dell'effetto pelle fisico. Il lavoro dimostra che, disponendo attentamente questi passaggi di calcolo — posizionandoli più vicini tra loro vicino alla superficie della barra dove avviene l'azione — è possibile creare un modello altamente accurato che funzioni su un intervallo di frequenze molto più ampio rispetto ai metodi precedenti.

Il viaggio inizia con un semplice quadro fisico: una lunga barra metallica rettangolare. Il ricercatore assume che il materiale sia uniforme e che i campi elettrici e magnetici al suo interno seguano le regole standard dell'elettromagnetismo, specificamente le equazioni di Maxwell. Nel regime in cui l'effetto pelle è forte, la profondità a cui la corrente penetra nel metallo è molto piccola rispetto allo spessore della barra. Questa realtà fisica detta una specifica relazione tra frequenza e perdita di segnale. Per trasformare questa legge fisica continua in qualcosa che un computer possa utilizzare, il ricercatore divide la barra in una serie di sottili fette. Questo processo, noto come discretizzazione spaziale, trasforma le equazioni fluide della fisica in una catena di equazioni differenziali ordinarie. Queste equazioni possono poi essere mappate direttamente su un circuito elettrico composto da resistori e condensatori, creando un modello comportamentale che imita la risposta della barra a un segnale in entrata.

Un'intuizione cruciale in questo lavoro è il modo in cui le fette sono disposte. Se le fette fossero tutte della stessa dimensione, il modello richiederebbe un numero impraticabile di esse per catturare i rapidi cambiamenti che avvengono proprio sulla superficie del metallo. Inveve, il ricercatore utilizza una griglia non uniforme, dove le fette diventano progressivamente più piccole man mano che ci si avvicina alla superficie. Questo raggruppamento dei punti di calcolo rispecchia la realtà fisica secondo cui i cambiamenti più significativi avvengono nello strato più sottile del conduttore. Lo studio trova che il rapporto tra le dimensioni di queste fette adiacenti è critico. Se le fette crescono troppo rapidamente di dimensione mentre ci si allontana dalla superficie, il modello diventa instabile e produce effetti di oscillazione (ringing) indesiderati. Le simulazioni suggeriscono che il rapporto tra le spaziature adiacenti debba essere mantenuto al di sotto di una specifica costante matematica, con un valore intorno a 2,5 che sembra offrire il miglior equilibrio tra accuratezza e intervallo di frequenze coperto.

Una volta stabilito il modello analogico, il passo successivo è tradurlo in un formato digitale che possa girare su moderni processori. Il ricercatore ha testato diversi metodi standard per convertire sistemi a tempo continuo in sistemi a tempo discreto. Un metodo popolare, noto come trasformata bilineare, è emerso come il vincitore chiaro. A differenza di altre tecniche che introducono uno smorzamento artificiale o richiedono calcoli di ordine superiore che complicano l'implementazione hardware, la trasformata bilineare preserva la stabilità e l'accuratezza del modello fisico originale. Il filtro digitale risultante si comporta esattamente come il modello analogico dell'effetto pelle, mantenendo il desiderato rapporto di radice quadrata su uno spettro ampio. Nelle simulazioni, questo nuovo approccio si è dimostrato superiore al filtro ricorsivo di Oustaloup, ampiamente utilizzato, una tecnica standard per approssimare sistemi a ordine frazionario. Mentre il metodo di Oustaloup posiziona poli e zeri su una scala logaritmica per appiattire la risposta di fase, il modello dell'effetto pelle derivato qui raggiunge naturalmente una maggiore accuratezza di fase e rimane valido su un intervallo di frequenze più grande rispetto al modello di Oustaloup.

I risultati delle simulazioni sono sorprendenti. Per un modello che utilizza solo otto punti di griglia, il nuovo filtro mantiene il suo comportamento caratteristico su un intervallo di circa 6.500 volte la frequenza iniziale. Questo è un miglioramento significativo rispetto ai metodi esistenti, che spesso faticano a mantenere l'accuratezza su bande così ampie. Lo studio conferma che l'intuizione fisica dietro l'effetto pelle fornisce una base robusta per la progettazione di questi filtri, evitando la necessità di complessi algoritmi di ottimizzazione. L'autore osserva che, sebbene l'attuale lavoro si concentri sul caso specifico della radice quadrata della frequenza, il metodo sottostante di discretizzazione delle equazioni differenziali alle derivate parziali potrebbe potenzialmente essere esteso ad altri ordini frazionari in futuro. Per ora, il lavoro offre un'alternativa chiara e fisicamente fondata per gli ingegneri che necessitano di un controllo preciso sull'elaborazione dei segnali, dimostrando che a volte il modo migliore per risolvere un problema digitale è guardare attentamente alla fisica del mondo reale.

Sintesi Tecnica: Filtri di Ordine Frazionario con Caratteristica $\sqrt{f}$

Definizione del Problema

Il documento affronta la progettazione e l'analisi di filtri di ordine frazionario caratterizzati da una funzione di trasferimento proporzionale a $s^{\alpha}$ dove $\alpha = \pm 1/2$ (corrispondenti ai comportamenti in frequenza $\sqrt{f}$ e $1/\sqrt{f}$ ). Questi filtri sono essenziali per applicazioni quali preemphasis/deemphasis, modellazione degli spettri di potenza del rumore $1/f$ (flicker), controllo di ordine frazionario e simulazione delle proprietà dispersive di mezzi come i tessuti biologici.

Sebbene le tecniche di approssimazione esistenti, come il filtro ricorsivo di Oustaloup e il metodo Charef, siano ampiamente utilizzate, esse si basano sul posizionamento di poli e zeri su una scala logaritmica per appiattire la risposta di fase o di ampiezza. Il documento cerca un'alternativa di derivazione radicata nelle leggi fisiche piuttosto che in una pura approssimazione matematica, con l'obiettivo di derivare un modello comportamentale che esibisca naturalmente le caratteristiche di ordine frazionario desiderate attraverso la discretizzazione delle equazioni di Maxwell.

Metodologia

1. Derivazione Fisica dalle Equazioni di Maxwell

La metodologia centrale inizia con le equazioni di Maxwell per una barra conduttrice lineare, isotropa e omogenea. Sotto l'assunto che la corrente di spostamento sia trascurabile ( $J \gg \partial D/\partial t$ ) e utilizzando un approccio quasi-stazionario per una barra rettangolare dove la lunghezza $l \gg$ larghezza $b$ e spessore $d$ , le equazioni si riducono a un sistema di equazioni differenziali alle derivate parziali (PDE) che descrivono i campi elettrico ( $E_z$ ) e magnetico ( $H_x$ ).

Il documento sfrutta il noto fenomeno fisico per cui la profondità di pelle (skin depth, $\Delta x$ ) di un materiale conduttore segue una legge $1/\sqrt{f}$ nel regime di forte effetto pelle (dove la profondità di pelle è piccola rispetto alla geometria). L'obiettivo è discretizzare queste PDE per creare un filtro che imiti questo comportamento.

2. Discretizzazione Spaziale e Modellazione Circuitale

Le PDE vengono discretizzate utilizzando uno schema a differenze finite (FD).

Griglia Uniforme: Inizialmente viene considerata una griglia equidistante, che porta a un sistema di equazioni differenziali ordinarie-algebriche (DAE). Questo sistema viene mappato su modelli circuitali comportamentali passivi (reti a scala) composti da induttori e resistori (LR) o condensatori e conduttori (GC).
Griglia Non Uniforme: Riconoscendo che una griglia uniforme richiederebbe un numero impraticabile di punti per catturare l'effetto pelle ad alte frequenze, il documento propone una griglia non uniforme. La spaziatura della griglia è minima in superficie ( $y = d/2$ ), dove l'effetto pelle è più pronunciato. Il rapporto tra le spaziature adiacenti è definito da un parametro $c$ .
Realizzazione Circuitale: Le equazioni discrete vengono realizzate come circuiti passivi equivalenti. Il documento nota che, sebbene il modello LR sia naturale, il modello GC (utilizzando condensatori e conduttori) è più adatto per simulatori di tipo SPICE poiché riduce il numero di incognite nell'Analisi Modale Nodale Modificata (MNA).

3. Trasformazione Digitale

Per ottenere un equivalente digitale, il filtro analogico viene trasformato utilizzando la trasformazione bilineare (regola trapezoidale) e confrontato con le Formule di Differenziazione Posteriore (BDF).

Trasformazione Bilineare: Preferita perché mappa l'asse immaginario sul cerchio unitario senza introdurre smorzamento numerico, un effetto collaterale osservato nei metodi BDF. Inoltre, risulta in un ordine di filtro inferiore rispetto ai metodi BDF di ordine superiore ( $p \ge 2$ ).
Implementazione: Il documento discute strategie di implementazione efficienti per le equazioni alle differenze risultanti, inclusa l'introduzione di vettori ausiliari per evitare moltiplicatori in specifiche fasi ricorsive.

Contributi Chiave

Derivazione Fisica dei Filtri Frazionari: Il documento deriva i filtri $\sqrt{f}$ e $1/\sqrt{f}$ direttamente dalla discretizzazione delle equazioni di Maxwell per una barra conduttrice, fornendo una base fisica per il comportamento di ordine frazionario anziché affidarsi esclusivamente a euristiche di posizionamento di poli e zeri.
Ottimizzazione della Griglia Non Uniforme: Dimostra che l'uso di una griglia non uniforme con un rapporto di spaziatura specifico ( $c$ ) riduce significativamente l'ordine del sistema preservando l'accuratezza. Il documento identifica che il rapporto tra le spaziature adiacenti dovrebbe essere inferiore alla costante di Eulero per evitare effetti di ringing, con un valore ottimale compreso tra $c=2$ e $c=2.5$ .
Modelli Circuitali Comportamentali: La derivazione produce modelli circuitali a scala passivi (sia LR che GC) che sono intrinsecamente stabili, poiché le perdite ohmiche garantiscono che i poli rimangano nel semipiano sinistro.
Implementazione Digitale: Il documento fornisce un percorso chiaro per la realizzazione digitale tramite la trasformazione bilineare, offrendo una struttura di filtro digitale stabile ed efficiente.

Risultati

Confronto con Oustaloup: I risultati delle simulazioni per un ordine di filtro $K=11$ mostrano che il modello di profondità di pelle derivato in questo lavoro offre un intervallo di validità approssimativa più ampio e una maggiore accuratezza di fase rispetto al filtro ricorsivo di Oustaloup.
Intervallo di Frequenza: Per una griglia non uniforme con $K=8$ punti e $c=2$ , il modello esibisce il comportamento $\sqrt{f}$ atteso (pendenza lineare di 0.5 in un grafico log-log dell'impedenza) su circa quattro decadi di frequenza. Sebbene le stime teoriche suggerissero un intervallo più ampio, le simulazioni indicano un intervallo di validità pratica di circa $f_1/10$ .
Sensibilità ai Parametri: I test confermano che i valori del rapporto di spaziatura della griglia $c$ tra 2 e la costante di Eulero ( $\approx 2.718$ ) offrono il miglior compromesso tra ordine del sistema, intervallo di validità ed esattezza. Valori di $c$ superiori alla costante di Eulero causano effetti di ringing indesiderati.
Prestazioni Digitali: I corrispettivi digitali, implementati tramite la trasformazione bilineare, mantengono il comportamento di radice quadrata per circa quattro decadi di frequenza angolare normalizzata, confermando l'efficacia della trasformazione.

Significato e Rivendicazioni

Il documento sostiene che il metodo proposto offre un'alternativa superiore alle attuali approssimazioni ricorsive come quella di Oustaloup, in particolare per quanto riguarda l'intervallo di validità e l'accuratezza di fase. Fondando la progettazione del filtro sull'effetto pelle fisico dei conduttori, i modelli risultanti non sono solo matematicamente solidi, ma anche fisicamente interpretabili come reti passive.

L'autore sottolinea che la strategia di discretizzazione non uniforme consente di modellare accuratamente l'effetto pelle con significativamente meno punti di griglia (e quindi un costo computazionale inferiore) rispetto a quanto richiederebbe una discretizzazione uniforme. La preferenza per la trasformazione bilineare è giustificata dalla sua capacità di preservare la stabilità ed evitare lo smorzamento numerico, rendendola adatta all'implementazione hardware.

Infine, il documento pone questo lavoro come un caso specifico di filtraggio di ordine frazionario ( $\alpha = \pm 1/2$ ) e suggerisce che la metodologia sottostante di discretizzazione delle equazioni differenziali alle derivate parziali potrebbe essere generalizzata in lavori futuri per filtri di ordine frazionario con ordini $\alpha$ arbitrari oltre $\pm 1/2$ .

Sintesi Tecnica: Filtri di Ordine Frazionario con Caratteristica f\sqrt{f}f​