Codifica numerica del segnale audio/Appendici

Appendici

../Capitolo 7 ../Bibliografia IncludiIntestazione 17 marzo 2021 75% Da definire

Capitolo 7 Bibliografia
[p. 298 modifica]

Appendice A


PROCESSI AUTOREGRESSIVI

A.1 SISTEMI LINEARI TEMPO INVARIANTI A TEMPO DISCRETO

Un sistema discreto è individuato dalla trasformazione T che lega l’ingresso x con l’uscita y

  (A.1)

Affinché un sistema discreto sia lineare deve risultare che la sua uscita ad una combinazione lineare di ingressi sia la stessa combinazione lineare alle singole uscite, cioè

  (A.2)

dove a e b sono delle costanti. Invece, un sistema è detto invariante alle traslazioni se, nota l’uscita y ad un certo ingresso x, la sua uscita ad un ingresso costituito dal segnale x traslato nel tempo coincida con la versione traslata della y

 
  (A.3)

Rappresentando un segnale discreto x(n) come una serie di impulsi δ(n)

  (A.3)
[p. 299 modifica]dove
  (A.5)

per sistemi lineari tempo-invarianti, la relazione

  (A.6)

può essere espressa, per la linearità, come

  (A.7)

Indicando con

  (A.8)

la risposta impulsiva del sistema, per l’invarianza alla traslazione si ottiene

  (A.9)

da cui

  (A.10)

Si può quindi concludere che un sistema discreto lineare tempo invariante è completamente caratterizzato dalla sua risposta impulsiva o, equivalentemente, dalla funzione di trasferimento ottenuta come trasformata z della stessa

  (A.11)

Lo svantaggio fondamentale della sommatoria di convoluzione introdotta per l’analisi dei sistemi nel dominio del tempo risiede nel fatto che gli estremi della sommatoria risultano pari a ±∞. Per evitare tale inconveniente, il legame ingresso-uscita può essere approssimato tramite equazioni alle [p. 300 modifica]differenze. Una equazione lineare alle differenze a coefficienti costanti di ordine n è un legame tra tra gli ultimi M campioni dell’ingresso e gli ultimi N dell’uscita del tipo

  (A.12)

A seconda della struttura dell’equazione alle differenze, si possono distinguere tre differenti tipi di sistemi. Se N = 0 (sistema Moving Average: MA), l’equazione alle differenze si riduce a

  (A.13)

Tale equazione non è più ricorsiva, in quanto dipende esclusivamente dall'ingresso. Confrontando tale equazione con quella di definizione della risposta impulsiva di un sistema, si nota come i coefficienti bk/a0 coincidano con i campioni della risposta impulsiva h(k). Posto x(n) = δ(n) si ha che la risposta impulsiva è data da

  (A.14)

ed è di durata finita: il sistema è definito FTR (Finite Impulse Response). Se M = 0 (sistema Auto Regressive: AR), l’equazione alle differenze si semplifica come

  (A.15)

A causa della retrazione dell’uscita, la risposta impulsiva del filtro è di durata finita (sistema IIR: Infinite Impulse Response). Il caso più generale è quello in cui sia M che N non sono nulli. Il sistema è quindi governato dall'equazione alle differenze

  (A.16)
[p. 301 modifica]ed il sistema (di tipo IIR) è detto ARMA. Esprimendo l’equazione alle differenze
  (A.17)

tramite la sua trasformata Z, si ottiene

  (A.18)

con funzione di trasferimento

  (A.19)

data da un rapporto di polinomi in z con a0 = 1. Se b0 ≠ 0, è possibile eliminare le potenze negative di z fattorizzando i termini b0z-M e a0z-N ed ottenendo

  (A.20)

dove G = b0/a0 = b0. Per sistemi FIR il denominatore non è presente e la funzione di trasferimento, è composta da soli zeri. Per sistemi AR la funzione di trasferimento è composta da soli poli.

Essendo la funzione di trasferimento dipendente dal contributo dei singoli poli e zeri, è possibile prevedere le caratteristiche della risposta impulsiva di un sistema a tempo discreto (e, quindi, il suo comportamento nel dominio del tempo) in funzione della loro localizzazione. Mentre il contributo del numeratore è poco intuitivo (ma sarà analizzato nel seguito), l’influenza della localizzazione dei poli è, invece, molto evidente. [p. 302 modifica]Imponendo il vincolo che l'uscita del sistema discreto sia una serie di valori reali, anche l’equazione alle differenze deve risultare a coefficienti reali. Nel caso di funzione di trasferimento con un solo polo, ciò comporta che il polo stesso si trovi sull'asse reale. Per quanto riguarda la risposta impulsiva si ha

  (A.21)

Nel caso in cui il polo si trovi sul semiasse positivo (r > 0), la risposta impulsiva risulta decadere esponenzialmente, essere costante o crescere esponenzialmente a secondo che il suo modulo risulti minore, uguale o maggiore di 1. Il valore di r determina la rapidità del decadimento. Nel caso in cui il polo si trovi sul semiasse negativo (r < 0), l’inviluppo della risposta impulsiva ha Io stesso andamento precedentemente descritto in funzione del modulo, solo che i campioni risultano a segni alternati. Il comportamento non cambia sostanzialmente nel caso di poli reali coincidenti

  (A.22)

Per quanto riguarda poli complessi coniugati

  (A.23)

la risposta impulsiva è oscillante con inviluppo decrescente, costante o crescente, secondo che i poli si trovino all'interno, sulla o all'esterno della circonferenza di raggio unitario.

Se un sistema è stabile, l’uscita prodotta a seguito di un ingresso limitato deve risultare limitata. Essendo l’uscita funzione della risposta impulsiva, dall'analisi precedente segue che un sistema è stabile se ha poli esclusivamente all'interno della circonferenza unitaria.

Per quanto riguarda la risposta in frequenza del sistema discreto LTI, definita come trasformata della sua risposta impulsiva, si nota come essa coincida con la H(z) calcolata sulla circonferenza di raggio unitario [p. 303 modifica]

 
  (A.24)

Interpretando i termini nelle produttorie, e rappresenta il vettore uscente dall'origine con vertice sulla circonferenza unitaria nel punto di pulsazione CO, mentre zk e pk rappresentano i vettori che congiungono l’origine con gli zeri ed i poli (fig. A.l). I termini tra parentesi tonde, quindi, sono i vettori congiungenti i poli e gli zeri con la circonferenza, che in forma polare sono esprimibili come

  (A.25)

con

  (A.26)

Noti i moduli Vk e Uk dei vettori congiungenti i poli e gli zeri con la circonferenza unitaria e le loro rotazioni Θk e Φk, il modulo della funzione di trasferimento si ottiene come prodotto dei moduli

  (A.27)

Analizzando il contributo dei singoli poli e zeri, si vede che il modulo della funzione di trasferimento si riduce per valori di © corrispondenti a punti della circonferenza unitaria prossimi ad uno zero, in quanto si riduce il relativo termine Vk. Per quanto riguarda i poli, invece, la riduzione del temine Uk si [p. 304 modifica]Fig. A.1 - Interpretazione geometrica del contributo di un polo al modulo ed alla fase della funzione di trasferimento.

traduce in un’esaltazione della risposta in frequenza, dato che tale termine si trova al denominatore della H(ω). È importante sottolineare come, a parità di modulo di due poli complessi coniugati, l’ampiezza del picco di risonanza dipenda dalla loro frequenza (fig. A.2). Viceversa, a parità di frequenza e interpretando il contributo di un polo come un filtraggio passa banda nell'intorno della sua frequenza centrale, la banda del filtro si riduca all'aumentare del modulo del polo. Tale miglioramento dal punto di vista della selettività in frequenza si paga, a causa della riduzione dello smorzamento, dal punto di vista della selettività temporale.

Per quanto riguarda la fase, essa dipende dalla differenza tra la somma delle fasi di ciascuno zero e quella dei poli

  (A.28)
[p. 305 modifica]Nel caso di sistemi discreti LTI posti tra loro in cascata, si può verificare che la funzione di trasferimento complessiva è data dal prodotto delle funzioni di trasferimento dei singoli sistemi. Infatti
  (A.29)

Ipotizzando di voler realizzare il sistema inverso ad uno dato. Ciò vuol dire che si vuole realizzare un sistema che messo in cascata ad un secondo produce una funzione di trasferimento complessiva unitaria. A tal fine è immediato verificare che è sufficiente realizzare un sistema che abbia come poli gli zeri del sistema da invertire e viceversa. Affinché anche il sistema inverso sia stabile, è necessario che i suoi poli siano all'interno della circonferenza unitaria, il che si traduce con il vincolo che il sistema originario, oltre ad avere i propri poli all'interno della circonferenza unitaria per la sua stabilità, vi abbia anche gli zeri. Sistemi con sia poli che zeri all'interno della circonferenza unitaria sono detti a fase minima.

A.2 FUNZIONE DI AUTOCORRELAZIONE DI UN PROCESSO AUTOREGRESSIVO

L’equazione caratteristica di un processo AR, soddisfatta sia dal processo u(n)

  (A.30)

che dall’autocorrelazione R(n)

  (A.31)

è pari a

  (A.32)
[p. 306 modifica]Fig. A.2a - Legame tra posizione dei poli, risposta impulsiva e risposta in frequenza. [p. 307 modifica]Fig. A.2b - Legame tra posizione dei poli, risposta impulsiva e risposta in frequenza. [p. 308 modifica]Per la stabilità, le radici p dell’equazione caratteristica debbo risiedere all’interno del circolo di raggio unitario. Si vuole mostrare come, variando i coefficienti del sistema AR, si riesca a variare le caratteristiche della funzione di autocorrelazione del processo generato.

Per un processo del secondo ordine, la funzione caratteristica è pari a

  (A.33)

con poli

  (A.34)

Imponendone la stabilità, si ricava la zona di esistenza dei coefficienti (fig. A.3)

  (A.35)

Fig. A.3 - Regione di esistenza per i coefficienti del processo AR. [p. 309 modifica]All'interno di tale zona, cambiando i coefficienti del processo, la funzione di autocorrelazione assume le seguenti caratteristiche

  (A.36)
L’andamento della funzione di autocorrelazione è riportata in figura A.4 per i tre casi. [p. 310 modifica]Fig. A.4a - Funzioni di autocorrelazione di processi autoregressivi. [p. 311 modifica]Fig. A.4b - Funzioni di autocorrelazione di processi autoregressivi. [p. 312 modifica]

Appendice B


METODO DEI MÌNIMI QUADRATI RICORSIVO

Il vantaggio principale dell’LMS è la sua semplicità. Per contro, il mantenere fisso il passo di aggiornamento è svantaggioso in quanto esso risulta essere solitamente troppo piccolo nelle fasi iniziali, ritardando la convergenza, ed eccessivo nelle fasi finali, deteriorando le prestazioni. Inoltre, dipendendo l’LMS solamente dall'errore istantaneo, non ha memoria dei precedenti passi di elaborazione. È possibile introdurre un differente algoritmo, derivato dall’LS, in cui l'aggiornamento dei coefficienti del predittore avviene indipendentemente e ricorsivamente per ciascuno di essi in funzione dell’andamento della stima (Recursive Least-Square: RLS). Ciò introduce memoria nell'algoritmo, come mostrato nel seguito, ma questo non rappresenta uno svantaggio nel caso di fenomeni stazionari. Nel caso di processi variabili, invece, è necessario attenuare la memoria sui campioni più remoti man mano che si procede nella stima. Ciò si ottiene con l’adozione di un nuovo criterio di ottimizzazione che introduce un fattore di pesatura w < 1 nella definizione dell’errore quadratico

  (B.1)

È poi necessario introdurre la funzione di errore a priori, definita come

  (B.2)
[p. 313 modifica]nella quale si stima il campione corrente in funzione dei pesi calcolati nel passo precedente, che risulta solitamente differente dalla funzione di errore a posteriori
  (B.3)

nella quale si utilizza il vettore dei coefficienti aggiornato. Ripetendo i passi già presentati nel capitolo 5, la m inim izzazione dell’errore in funzione dei coefficienti del predittore porta al passo n-esimo all’equazione normale

  (B.4)

che risulta formalmente identica a quella utilizzata nella derivazione delle equazioni di Wiener-Hopf, ma nelle quale è implicitamente presente il fattore di pesatura w. La matrice di autocorrelazione dell’ingresso ed il vettore di cross-correlazione tra l’ingresso e l’uscita desiderata risultano, infatti, pari a

 
  (B.5)

La soluzione di tale sistema è ancora dato da

  (B.6)

anche se, in questo caso, ne ricerca una sua soluzione ricorsiva. Per quanto riguarda la matrice di autocorrelazione, è possibile dame una definizione ricorsiva isolandone dalla definizione il termine per i = n

 
  (B.7)
[p. 314 modifica]Analogamente si ottiene
  (B.8)

L'inversione della matrice di autocorrelazione può essere evitata, in quanto risulta [Hay86]

  (B.9)

Per semplificare la notazione si introduce la matrice di autocorrelazione inversa

  (B.10)

ed il vettore dei guadagni

  (B.11)

che ha le stesse dimensioni dell'ordine del predittore. In tal modo, per l'equazione dell'inversa della matrice di autorrelazione si ottiene la relazione ricorsiva

  (B.12)

Riorganizzando la definizione del vettore dei guadagni si ottiene, invece

  (B.13)

con la quale si ottengono le seguenti relazioni ricorsive per i coefficienti del predittore [p. 315 modifica]

 
 
 
  (B.14)

Sostituendo l’espressione dell’errore

  (B.15)

si ottiene, infine


  (B.16)

Da tale equazione si nota come il vettore dei guadagni k(n) ha le stesse funzioni del parametro p nell’LMS. A differenza di quest’ultimo, però, esso risulta variabile e, dipendendo dalla matrice di correlazione del segnale, ha memoria del processo di stima.

In tal modo sono state ottenute equazioni ricorsive per il calcolo di k, p, a e P. Al fine di ridurre la complessità computazionale, però, è opportuno definire la grandezza

  (B.17)

con la quale i passi per la predizione sono ridotti ai seguenti: - si calcola la grandezza

  (B.18)

- da questa si calcola il vettore dei guadagni

  (B.19)


errore quadratico medio [p. 316 modifica]Fig. B.1 - Curve di apprendimento per LMS e RLS.

  • si calcola l’errore
  (B.20)
  • si aggiornano i parametri del predatore
  (B.21)
  • si calcola, infine, il valore aggiornato della matrice di autocorrelazione
  (B.22)

Per quanto riguarda le condizioni iniziali, è necessario scegliere la matrice P(0) in modo tale che risulti non singolare. La scelta più immediata è quella di scegliere una costante δ positiva “piccola” (es.: 0.1) e fissare [p. 317 modifica]

  (B.23)

Per i coefficienti del predatore il loro valore iniziale è tipicamente nullo, mentre il coefficiente di smorzamento w deve risultare eventualmente minore, ma comunque molto prossimo ad 1 (es.: 1-10-3).

Dal punto di vista delle prestazioni, la variabilità del vettore dei guadagni fa si che sia i tempi di convergenza che le fluttuazioni nell'intorno del vettore ottimo siano ridotti rispetto all’LMS (fig. B.l). In compenso, però, la complessità computazionale è aumentata, anche considerando versioni più sofisticate dell’algoritmo [Hay86]. [p. 318 modifica]

Appendice C


ALGORITMI A BLOCCHI

PER IL FILTRAGGIO ADATTATIVO

Gli algoritmi per la soluzione dell’equazione normale della predizione a blocchi possono essere ricondotti, fondamentalmente, alle seguenti due classi (fig.C.1)

  • metodo della covarianza;
  • metodo dell’autocorrelazione;

La loro differenza risiede nel differente modo con il quale vengono stimati gli elementi della matrice di autocorrelazione Φ. La stima della funzione di auto correlazione viene calcolata come sommatoria dei prodotti ottenuti utilizzando due finestre di campioni di pari lunghezza progressivamente traslate. Nel secondo caso si utilizza un’unica finestra di campioni. Le traslazioni richieste nel calcolo dei prodotti avvengono estendendo tale finestra base con dei campioni nulli. Nel primo caso, invece, la finestra base di campioni viene estesa utilizzando effettivi campioni del segnale. Tale differenza porta poi ad una differente struttura della matrice Φ e quindi a differenti algoritmi per la soluzione del sistema lineare di equazioni associato.

Il metodo dell autocorrelazione è quello attualmente più utilizzato sia per l’esistenza di algoritmi più robusti nei riguardi di errori di quantizzazione (adottando strutture a traliccio), sia per resistenza di più efficienti algoritmi di calcolo (ricorsione di Schür). [p. 319 modifica]Fig. C.1 -Differenti finestre utilizzate nel metodo della covarianza e dell’autocorrelazione.

C.l METODO DELLA COVARIANZA

Nel metodo della covarianza gli elementi della matrice Φ vengono stimati utilizzando sequenze di “n” campioni del segnale, progressivamente sfasate. Ciascun elemento della matrice viene cioè calcolato come

 
  (C.1)
[p. 320 modifica]Considerare blocchi di “n” campioni del segnale vuol dire moltiplicare il segnale stesso per una finestra di pari lunghezza. La finestra utilizzata nel seguito è quella rettangolare, tralasciando considerazioni sull'opportunità di utilizzare finestre differenti.

Per il calcolo degli elementi Φ(i, k) da utilizzare nella predizione, è necessario considerare il prodotto di sequenze traslate nel tempo al massimo di “p” campioni. Se la finestra di interesse è quella relativa alla sequenza x(n), è necessario considerare fino a “p” campioni esterni a tale intervallo. Con tale approssimazione, la matrice Φ dell’equazione normale deterministica viene indicata come matrice di covarianza, da cui il nome dell’algoritmo. Si mette in evidenza come tale nome non indica che i termini della matrice rappresentino la covarianza del segnale, ma sono sempre termini di autocorrelazione. La soluzione dell’equazione si ottiene poi come

  (C.2)

Es.: si consideri il vettore

  (C.3)

ottenuto considerando la funzione f = cos(t) + cos(2t) e prelevando 5 campioni nel suo periodo. Considerando un predittore del secondo ordine, l’equazione normale è

  (C.4)

Per il calcolo degli elementi della matrice di covarianza è necessario estendere l’intervallo di osservazione, includendo anche due campioni del precedente periodo. Nel nostro caso traslazioni verso destra o verso sinistra si equivalgono.

 
  (C.5)

È possibile, quindi, calcolare gli elementi della matrice di covarianza [p. 321 modifica]

 
 
 
 
  (C.6)

e calcolare il vettore dei coefficienti ottimi come

  (C.7)

Il calcolo della Φ può avvenire anche evitando di utilizzare campioni esterni al blocco corrente. A tal fine è necessario considerare sequenze di n-p campioni opportunamente sfasate ritagliate all'interno del blocco. Il calcolo della Φ può poi avvenire costruendo la matrice delle osservazioni

  (C.8)

dalla quale sì ricavano la matrice di auto-correlazione ed il vettore di cross-correlazione

 
  (C.9)


La matrice di covarianza risulta essere definita positiva e, per costruzione, simmetrica. La soluzione dell’equazione normale può essere quindi trovata evitando l’inversione della stessa, ma scomponendola tramite l’algoritmo di Cholesky nel prodotto [p. 322 modifica]

  (C.10)

dove V è una matrice triangolare inferiore con la diagonale principale unitaria, mentre D è una matrice diagonale

  (C.11)

Per eseguire la fattorizzazione della matrice di covarianza si procede come segue: dato

  (C.12)

per gli elementi sulla diagonale segue

  (C.13)

a partire dalla condizione iniziale

  (C.14)

mentre gli elementi della matrice V si ricavano riscrivendo l’espressione del generico elemento della matrice di covarianza isolando dalla sommatori il termine per k=j

 
  (C.15)
[p. 323 modifica]In tal modo è possibile calcolare ricorsivamente gli elementi Vij e dj. Una volta fattorizzata la matrice di covarianza, anche la soluzione dell’equazione normale può essere ottenuta ricorsivamente. Infatti, definendo il vettore y come
  (C.16)

risulta

 
  (C.17)

Essendo V una matrice triangolare, gli elementi della y possono essere calcolati ricorsivamente come

 
  (C.18)

Una volta calcolati i componenti del vettore y è possibile ricavare i coefficienti ottimi del predittore in quanto, risultando

  (C.19)

Si ottengono le relazioni ricorsive

  (C.20)

con condizioni iniziali

  (C.21)
[p. 324 modifica]Es. si consideri lo stesso vettore di campioni precedente. Applicando la

fattorizzazione di Cholesky risulta

 
 
  (C.22)

Fattorizzata la matrice di covarianza come

 
  (C.23)

si calcolano dapprima le variabili ausiliarie y

 
  (C.24)

tramite la relazione

 
 
  (C.25)
[p. 325 modifica]ed infine i coefficienti del predittore dalla relazione
 
 
  (C.26)
C.2 METODO DELL’AUTOCORRELAZIONE

II metodo dell'autocorrelazione stima gli elementi della matrice F utilizzando esclusivamente campioni all'interno di una finestra di lunghezza “n”. Ciò equivale ad assumere implicitamente che il segnale sia nullo al di fuori di tale intervallo e quindi la traslazione di sequenze di campioni si ottengono con il suo completamento tramite campioni nulli. In tal modo il calcolo dell’autocorrelazione è eseguito su di un numero di campioni non nulli via via inferiore man mano che si calcolano valori di autocorrelazione di indice crescente. Gli elementi della matrice di autocorrelazione andrebbero quindi calcolati come

  (C.27)

In tal modo, però, al diminuire del denominatore, verrebbe esaltato lo scostamento della stima dal valore effettivo dell'autocorrelazione, scostamento dovuto al ridursi del numero di campioni utili. Si preferisce, quindi, non eseguire il calcolo degli elementi della matrice della autocorrelazione (non polarizzata) secondo tale definizione, ma calcolare gli elementi della matrice della autocorrelazione (polarizzata) come

  (C.28)
[p. 326 modifica]Ciò permette di attenuare il contributo degli elementi esterni in quanto,

al ridursi delle dimensioni delle sequenze, il valore del denominatore rimane costante. L’equazione normale deterministica in forma matriciale diventa

 
  (C.29)

con soluzione

  (C.30)

Es.: ripetiamo il calcolo dei coefficienti del predittore del secondo ordine

  (C.31)

per lo stesso segnale precedentemente considerato. Gli elementi della matrice di autocorrelazione sono dati da

 
 
  (C.32)

I coefficienti del predittore sono, quindi, calcolati come

  (C.33)
[p. 327 modifica]Anche in questo caso è possibile utilizzare la matrice delle osservazioni,

avente dimensioni [n + m -1, m], che assume la forma

  (C.34)

Dalla H si ricavano la matrice di autocorrelazione ed il vettore di cross-correlazione

 
  (C.35)

La soluzione diretta dell’equazione di Yule-Walker tramite l’inversione della matrice di auto correlazione può essere anche in questo caso evitata, dato che essa risulta essere una matrice di Toeplitz, cioè simmetrica con tutti gli elementi di una diagonale identici. Ciò permette di ricavare ricorsivamente i coefficienti di un predittore di ordine p in funzione dei coefficienti di predittori di ordine inferiore tramite una relazione ricorsiva del tipo

  (C.36)

dove αk(m) indica il coefficiente k-esimo del predittore di ordine “m” e quindi al passo m-esimo dell’algoritmo. Le condizioni iniziali sono quelle di un predittore del primo ordine, che possono essere ricavate direttamente dall'equazione normale

 
  (C.37)
[p. 328 modifica]Dall’espressione del minimo dell’MSE, per un predatore del primo ordine si ricava
  (C.38)

Considerando poi il problema di calcolare i coefficienti di un predittore di ordine m in funzione dei coefficienti di un predittore di ordine (m-1) è necessario poter risolvere il sistema


  (C.39)

con δ(m-1) e k(m) da determinare. In tal modo si ottiene l’equazione aumentata del predittore

  (C.40)

dove rb(m-1) rappresenta il vettore di cross-correlazione al passo m-1 r(m-1) con gli elementi in ordine invertito. Da questa relazione matriciale derivano le due relazioni

 
  (C.41)

Per l’equazione del predittore di ordine (m-1)

  (C.42)

dalla prima equazione si ottiene

 
  (C.43)
[p. 329 modifica]dove αb(m-1) rappresenta il vettore dei coefficienti del predittore di ordine m-1 α(m-1) con gli elementi in ordine invertito. Sostituendo l’espressione di δ(m-1) nell’equazione iniziale è possibile ricavare ricorsivamente i coefficienti del predittore come
 
  (C.44)

In tale espressione la costante k(m) risulta ancora incognita. Per poterla ricavare si sostituisce l’espressione di nella seconda componente dell’equazione normale aumentata, ottenendo

 
  (C.45)

Riconoscendo nel denominatore l’espressione dell’errore minimo al passo m-1, si ha

  (C.46)

Anche l’errore può essere calcolato ricorsivamente, dato che

  (C.47)

Sostituendo l’espressione ricorsiva degli a k si ottiene

  (C.48)
[p. 330 modifica]Confrontando l’espressione in parentesi quadre con il numeratore dell’espressione che fornisce k( m ), si ottiene
  (C.49)

ed essendo , si ha, infine

  (C.50)

Riepilogando, è possibile calcolare ricorsivamente tramite l’algoritmo di Levinson-Durbin i coefficienti di un predatore di ordine p in funzione dei coefficienti di predittori di ordine inferiore. Ciò si ottiene tramite i seguenti passi:

  • si determina il valore dell’errore dalla relazione
 
  (C.51)
  • si ricava il valore del coefficiente k(m) dalla relazione
 
  (C.52)
  • si determinano i coefficienti del predittore
 
  (C.53)
  • si determina il valore aggiornato dell’errore
  (C.54)
[p. 331 modifica]Tale procedura si ripete per un numero di passi pari all'ordine del predittore.

Es.: ripetiamo il calcolo dei coefficienti del predittore per lo stesso segnale precedentemente considerato. L’errore al passo iniziale è

  (C.55)

Il ciclo per m = 1 porta al calcolo delle grandezze di un predittore del primo ordine

 
 
  (C.56)

Per m = 2 si ottengono le grandezze del predittore del secondo ordine

 
 
  (C.57)
C.3 METODI UTILIZZANTI STRUTTURE A TRALICCIO

Sia il metodo della covarianza che il metodo dell’autocorrelazione precedentemente esposti sono concettualmente organizzati su due fasi: una stima preliminare della matrice di autocorrelazione del segnale ed il successivo calcolo dei coefficienti del predittore. Il metodo che si vuole presentare in questo paragrafo (lattice) può essere visto come un raffinamento del metodo dell’autocorrelazione, nel quale tali due fasi sono fuse in un’unica procedura ricorsiva. [p. 332 modifica]Per descrivere l’algoritmo è necessario riprendere l’equazione del predittore

  (C.58)

e della sua funzione di trasferimento

  (C.59)

Definita la funzione d’errore del predittore di ordine “p” come

  (C.60)

è possibile definire la funzione di trasferimento del sistema d’errore. Si considera, cioè, il sistema che, dato come ingresso il segnale, produce in uscita la funzione d’errore. Tale funzione di trasferimento è pari a

  (C.61)

In tal modo è possibile esprimere l’errore di predizione in z come

  (C.62)

Ora ci si può porre il problema di esprimere la A(z) di un predittore di ordine “m” in funzione della A(z) di un predittore di ordine Per tale problema si può utilizzare la trattazione svolta per l’algoritmo di Levinson-Durbin che fornisce un’espressione ricorsiva per i coefficienti del predittore a k. Sostituendo tale relazione nella definizione di A(z) si ottiene

 
[p. 333 modifica]
 
 
  (C.63)

Sostituendo nell'espressione dell’errore si ottiene

  (C.64)

Il primo termine del secondo membro rappresenta la trasformata z della funzione d’errore di un predittore di ordine “m - 1”. Tale predittore stima il campione x(n) utilizzando i campioni [ x(n - m)... x(n -1) ] e verrà nel seguito indicato come predittore in avanti. Per il secondo termine si può dare un’interpretazione simile se si definisce

  (C.65)

che ha un’antitrasformata del tipo

  (C.66)

Questo è interpretabile come l’errore di un predittore che tenta di stimare x(n - m) in funzione degli “m-1” campioni [x(n m 1)... x(n) ]. Tale predittore verrà nel seguito indicato come predittore all'indietro. Sostituendo la B(z) nella E(z) si ottiene

  (C.67)

che, antitrasformata fornisce

  (C.68)
[p. 334 modifica]Tale espressione rappresenta una relazione ricorsiva della funzione d’errore della predizione in avanti. Sostituendo l’espressione ricorsiva della A(z) nella definizione della B(z) si ottiene
 
  (C.69)

che, antitrasformata, fornisce la relazione ricorsiva della funzione d’errore della predizione all'indietro

  (C.70)

Se si traccia il diagramma di flusso di tale algoritmo si ottiene una struttura a traliccio (lattice), da cui il nome dell’algoritmo (fig. C.2). È quindi possibile il calcolo ricorsivo delle due funzioni d’errore. Le condizioni iniziali sono quelle per un predittore di ordine zero, cioè senza alcuna predizione e quindi

  (C.71)

In forma matriciale le relazioni ricorsive in z sono esprimibili come

  (C.72)

Per il calcolo dei coefficienti k( m ) si impone la simultanea minimizzazione degli errori quadratici di predizione in avanti e all'indietro

  (C.73)

La minimizzazione, ottenuta annullando le derivate rispetto ai k [p. 335 modifica]

 
  (C.74)

porta alla seguente espressione ricorsiva per il calcolo dei coefficienti

  (C.75)

Dato che tale espressione ha la forma di una cross-correlazione tra le funzioni di errore di predizione in avanti e all’indietro, i coefficienti k( m ) vengono detti PARCOR (PARtial CORrelation). Noti i coefficienti k( m ) è possibile determinare i coefficienti del predittore come già visto nel Levinson-Durbin (dal quale il presente algoritmo deriva)

 
  (C.76)

Per chiarire ulteriormente la struttura degli errori di predizione è opportuno analizzare con un esempio la loro struttura. Per un predittore di ordine zero risulta

  (C.77)
[p. 336 modifica]Fig. C.2 - Predittore di ordine p implementato tramite una struttura a traliccio.

Per un predittore di ordine uno la predizione in avanti utilizza solamente il campione che precede il corrente, mentre la predizione all'indietro utilizza il campione corrente per predire il precedente. Infatti risulta

  (C.78)

Analogamente, per un predittore di ordine due risulta

  (C.79)

e così via. [p. 337 modifica]Riassumendo, i passi per il metodo lattice sono i seguenti:

  • si impostano le condizioni iniziali
  (C.80)
  • si calcolano i coefficienti PARCOR
  (C.81)


  • si calcolano i coefficienti del predittore
 
  (C.82)
  • si aggiornano le funzioni di errore
 
  (C.83)

con vettori che, inizialmente di dimensioni pari ad “n”, assumono dimensioni crescenti fino a “n + p”.

Es.: si ripete il calcolo del predittore del secondo ordine con il segnale precedentemente utilizzato. Le condizioni iniziali sono

 
 
 
  (C.84)
[p. 338 modifica]Il coefficiente PARCOR del primo ordine è
 
  (C.85)

che porta alle nuove funzioni di errore

 
 
 
  (C.86)

ed al coefficiente

  (C.87)

Per il secondo ordine si ha

  (C.88)

Si possono quindi calcolare i coefficienti del precettore del secondo ordine

 
  (C.89)
[p. 339 modifica]Si noti come i valori dei coefficienti coincidono con quelli calcolati con

il metodo dell’auto-correlazione

C.4 METODO RICORSIVO DI SCHÜR

Si osserva che l’algoritmo di Levinson-Durbin fa uso di una relazione ricorsiva per Terrore al passo m-esimo in funzione dell’errore e dei coefficienti PARCOR al passo (m-1) del tipo

  (C.90)

È possibile esprimere tale errore il funzione dell’autocorrelazione esplicitandola come

 
 
 
  (C.91)

Nell'algoritmo di Schür si esprimono tutte le relazioni al generico passo m-esimo in funzione degli elementi della matrice di autocorrelazione R(i) (i=0,l,...,m) e dei parametri k(i) (i=l,2,...,m). Per l’algoritmo di Schür si definiscono innanzitutto gli elementi

 
 
 
  C.92
[p. 340 modifica]dove le dimensioni di G0 sono pari a 2 * (p+1), con p ordine del predittore. Si esegue poi il prodotto V1G1 e si definisce G2 come il risultato della precedente operazione con la seconda riga traslata di un posto verso destra. Si ripete il procedimento per p passi, ottenendo al generico passo m-esimo le seguenti relazioni ricorsive
 
  C.93

Seguendo passo-passo Io svolgimento dell’algoritmo di Schür, si può facilmente verificare che al passo m-esimo l’elemento Gm(2,m+1), ovvero il denominatore della relazione che definisce K(m), corrisponde al termine e(m-1) nell'algoritmo di Levinson-Durbin. Ad esempio, se implementiamo la stima dei parametri di un predittore di ordine 3 la matrice G 0 è definita come

  C.94

Si ottiene poi:

 
 
  C.95

dove, al generico passo m-esimo, si è ricavato R(m) invertendo la relazione che definisce k(m) nell'algoritmo di Levinson Durbin.

Si nota anche che al generico passo m-esimo dell’algoritmo di Schür, la matrice Gm presenta le prime m colonne nulle, il che può essere sfruttato per ridurre sia la complessità di calcolo che l’occupazione di memoria.

Infatti, la versione dell’algoritmo di Schifi proposta dallo standard GSM organizza le informazioni contenute nelle matrici G m in due vettori (K e P) ed evita, ad ogni passo, di aggiornarne quegli elementi che non apportano [p. 341 modifica]contributo ai passi successivi. In questo modo non si utilizzano relazioni matriciali; tutte le operazioni sono ridefinite sugli elementi di K e di P, vettori rispettivamente di dimensione (p-1) e (p+1). Al generico passo m, la determinazione di r(m) è sempre effettuata come rapporto tra il secondo elemento di P ed il primo; vengono aggiornati tutti gli elementi di P e di K tranne gli ultimi m. [p. 342 modifica]

Appendice D


RICHIAMI SU FILTRI NUMERICI

D.1 STRUTTURE PER MULTI RATE DSP

Per Multirate Digital Signal Processing si intende la realizzazione di sistemi discreti tali che il periodo T’dei segnali in uscita sia differente da quello T dei segnali in ingresso. Nel caso in cui la frequenza di campionamento del segnale in uscita sia maggiore di quella del segnale d’ingresso, è necessaria un’operazione di interpolazione per ricostruire il segnale x(t) mancante. Viceversa, nel caso in cui la frequenza del segnale in uscita sia minore di quella d’ingresso, è necessario decimare il segnale, eliminando campioni. È, però, necessario far precedere alla decimazione un’operazione di filtraggio che limiti la banda del segnale da ricampionare, al fine di evitare aliasing. Indicando con h(n) la risposta impulsiva dei filtri coinvolti nel cambiamento di frequenza (anti aliasing o di interpolazione), l’uscita è prodotta come

  (D.1)

In un sistema discreto è conveniente che il rapporto tra i periodi di campionamento in ingresso ed in uscita sia riconducibile ad un rapporto tra interi del tipo

  (D.2)
[p. 343 modifica]Inoltre, data la linearità del filtraggio, è possibile considerare separatamente i due casi di interpolazione ( M = 1, L > 1 ) e di decimazione (M>1,L=1). Il caso più generale { M > 1,L> 1 } è poi risolvibile come cascata di una interpolazione seguita da una decimazione.

Nel caso della decimazione, il ricampionamento si ottiene semplicemente eliminando M-l campioni del segnale originario ogni M. Al fine di evitare aliasing, il segnale dovrebbe essere preventivamente limitato in banda tramite un filtro passa basso ideale con funzione di trasferimento

  (D.3)

Indicando con h(n) la corrispondente risposta impulsiva, l’uscita si ottiene come

  (D.4)

Con una trasformazione di variabili, l’equazione precedente può essere riscritta come

  (D.5)

dalla quale si nota come i campioni derivanti dalla decimazione si ottengano traslando progressivamente la risposta impulsiva del filtro passa basso con un passo pari ad M campioni. Dal punto di vista della rappresentazione in frequenza, l’operazione di decimazione riespande lo spettro del segnale, limitato dal filtraggio nell’intervallo [-π/M, π/M], ad occupare tutto l’intervallo [-π, π].

Per eseguire, invece, un’interpolazione, è necessario introdurre L-l campioni ogni due campioni adiacenti del segnale originale. Un modo per ottenere questo risultato è quello di considerare una versione sovracampionata [p. 344 modifica]di un fattore L del segnale PAM corrispondente alla x(n), ottenuto introducendo dei campioni nulli nel segnale

  (D.6)

Dal punto di vista della rappresentazione in frequenza, ciò non altera lo spettro del segnale, coincidente con quello del corrispondente segnale PAM, il quale presenta repliche a partire dalla pulsazione ω = π/L. Per passare allo spettro della x(t) sovracampionata con un fattore L è necessario eliminare tali repliche tramite un filtraggio passa basso, con frequenza di taglio corrispondente a ω = π/L. Nel tempo, tale filtraggio si traduce in un’interpolazione che genera i campioni mancanti nella x(n). Per mantenere nel segnale sovracampionato l’ampiezza del segnale originario, il filtro di interpolazione deve avere un guadagno non unitario, ma pari ad L. Dato che, nel filtraggio, la risposta impulsiva h(n) si trova ad essere sistematicamente combinata con campioni nulli della w(n), la sommatoria di convoluzione può essere semplificata come

  (D.7)

Analizzando quali campioni della h(n) vengono interessati dalla sommatoria, si nota che vengono utilizzati periodicamente L suoi sottoinsiemi gm, tra loro scalati di un campione (fig. D.l). Tale sfasamento corrisponde ad un offset temporale che che è una frazione del periodo T pari a

  (D.8)

E quindi possibile riscrivere l’equazione che fornisce la y(m) come l’uscita di un filtro tempo variante con risposta impulsiva g m, tale che

  (D.9)
[p. 345 modifica]dove come indice della x deve considerarsi, ovviamente, un intero.

Analogamente, è possibile considerare i campioni della y(m) come prodotti da un banco di L filtri in parallelo, ciascuno con una differente risposta impulsiva gm(n). L’uscita dell’interpolatore è poi ottenuto selezionando periodicamente l’uscita di uno degli L filtri che operano tra loro sfasati di un colpo di clock.

Fig. D.1 -Insieme di campioni della risposta impulsiva coinvolti in un’operazione di inteipolazione. [p. 346 modifica]Passando al problema dell’implementazione, è necessario considerare che i filtri utilizzati nella conversione di frequenza di campionamento sono tipicamente FIR a lunghezza finita N (fig. D.2), per cui l’uscita è ottenuta come

  (D.10)

Fig. D.2 - Strutture di filtri FIR per decimazione. [p. 347 modifica]Considerando per il momento il problema della decimazione, un inconveniente che si incontra implementando direttamente tale filtro è che esso lavora alla frequenza del segnale sovracampionato, mentre sarebbe auspicabile che lavorasse alla frequenza del segnale decimato. Per ottenere ciò è possibile anticipare l’operazione di decimazione prima del prodotto con i coefficienti del filtro come

  (D.11)

ottenendo un filtro con clock pari a quello dell’uscita. Imponendo poi che i filtri siano a fase lineare e quindi con risposta impulsiva risulta simmetrica, la lunghezza del filtro può essere dimezzata in quanto

  (D.12)

Strutture analoghe a quelle presentate nel caso della decimazione possono essere realizzate nel caso di interpolazione. In ogni caso, volendo avere filtri sufficientemente selettivi, le loro dimensioni risultano essere notevoli. L’efficienza computazionale può essere incrementata riducendo la lunghezza del filtro tramite filtri polifase (fig. D.3). È stato mostrato come i campioni prodotti da un interpolatore possano essere pensati come l’uscita di uno di L filtri in parallelo, ciascuno con risposta impulsiva g m (n). Dato che la risposta impulsiva di ciascun ramo del filtro coincide con una versione decimata della h(n), che è diversa da zero solo negli istanti L x T ed è nulla altrimenti, è possibile eliminare i contributi di tali campioni nulli considerando un insieme di L filtri di dimensioni N/L tali che

  (D.13)

Ovviamente è opportuno che la lunghezza N del filtro prototipo sia multipla di L. Con una struttura polifase, per ogni campione dell’ingresso vengono prodotti L campioni d’uscita, uno per ogni ramo del filtro. L-l di tali [p. 348 modifica]campioni costituiscono l’interpolazione tra due campioni adiacenti dell’ingresso, necessari per il richiesto incremento della frequenza di campionamento.

Fig. D.3 - Strutture di filtri FIR polifase.

Per quanto riguarda la funzione di trasferimento di ciascuno degli L rami del filtro polifase, si ricorda che essi sono stati ottenuti come decimazione di un fattore L di un filtro passa basso ideale con frequenza di taglio corrispondente a ω = π/L. Dato che la decimazione riespande lo spettro di un fattore L, [p. 349 modifica]ciascun ramo del filtro polifase risulterà con funzione di trasferimento piatta nell'intervallo [-π, πJ. L’unica differenza tra i differenti rami è quindi data dalla fase, da cui il nome del filtro.

Tipicamente la struttura di un filtro polifase non può essere semplificata ulteriormente sfruttando eventuali simmetrie. Infatti, pur utilizzando una h(n) simmetrica, le L risposte impulsive dei rami del filtro polifase non lo sono. D’altra parte, l’indipendenza tra i vari rami del filtro, ne rende possibile un'implementazione hardware parallela. Per quanto riguarda la decimazione, la relazione

  (D.14)

può essere riscritta come

  (D.15)

Definendo i coefficienti del filtro polifase di decimazione come

  (D.16)

si ottiene

  (D.17)

dove xk(n) = x(nM - k) è l’ingresso del ramo k-esimo del filtro. Di tali filtri è possibile dare un’interpretazione tramite un commutatore che preleva campioni nel caso della interpolazione e distribuisce campioni nel caso della decimazione ruotando in senso antiorario. Ciò vuol dire, ad esempio, che nella decimazione il primo campione ingresso viene assegnato al ramo di ordine maggiore e così via. È possibile realizzare strutture con senso di rotazione del commutatore orario definendo i filtri polifase come

  (D.18)
[p. 350 modifica]Un’applicazione importante di multirate DSP si ha nell'elaborazione disegnali passa banda, cioè di segnali con una banda limitata ωΔ, ma centrati in ω0 ≠ 0 (fig. D.4). Tale esigenza si incontra, ad esempio, in un banco di filtri: il fine è quello di ricavare il corrispondente segnale in banda base X(m), cioè il segnale che, a parità di banda, risulta centrato nell'origine. Ciò è concettualmente ottenibile innanzitutto traslando il segnale passa banda nell'origine tramite una modulazione e quindi decimandolo dopo il necessario filtraggio antialiasing tramite un passa basso con frequenza di taglio corrispondente a ωΔ/2. Indicando con h(n) la risposta impulsiva di tale filtro, si ottiene
  D.19

dove M è il fattore di decimazione che si ricava dal rapporto tra la frequenza centrale del segnale e la sua banda. Per ricostruire il segnale originale dalla sua versione decimata, invece, è innanzitutto necessario interpolarlo per poi modularlo, riportandolo su ω0. Indicandolo con f(n) la risposta impulsiva del filtro di interpolazione, si ottiene

  D.20

Fig. D.4 - Differenti modulazioni per segnali passabanda. [p. 351 modifica]Dato che il segnale d’uscita risulta essere complessa, a causa del prodotto con gli esponenziali, tale tecnica è detta di modulazione complessa (fig. D.5). La struttura del sistema di decimazione/interpolazione può essere anche differente da quella presentata, portando le modulazioni come ultima operazione nella decimazione e come prima della interpolazione. Con tale soluzione, però, è necessario l’uso di filtri passa banda centrati sul segnale stesso, al posto dei passa basso del caso precedente.

Fig. D.5 - Modulazione complessa ed a banda laterale unica.

Analiticamente l’espressione del decimatore può essere ottenuta dalla precedente moltiplicandola per e riarrangiando i termini come

  (D.21)
[p. 352 modifica]dove il termine tra parentesi quadre corrisponde alla risposta impulsiva del filtro passa banda di ampiezza ωΔ centrato su ω0. Analogamente, per l’interpolazione
  (D.22)

Le relazioni presentate forniscono un mappaggio del segnale passa banda in un segnale complesso in banda base nell'intervallo [- ωΔ/2, ωΔ/2] e viceversa. Per applicazioni di codifica del segnale, invece, è più interessante considerare un mappaggio nell'intervallo [0, ωΔ ], in grado di fornire un segnale reale. Per tale mappaggio, noto come modulazione SSB per l’analogia con le analoghe tecniche trasmissive, al segnale ottenuto dalla modulazione complessa viene applicata una traslazione in frequenza di metà larghezza di banda, tramite il prodotto con e ejnωΔ/2.

D.2 BANCHI DI FILTRI POLIFASE

La funzione di un banco di filtri può essere sia quella di separare un segnale nelle componenti che fanno capo a differenti bande spettrali (filtro di analisi), sia quella di ricombinare i contributi di bande differenti in un unico segnale (filtro di sintesi) (fig. D.6). Per tali operazioni è necessario considerare i differenti segnali passa banda relativi a ciascuna sottobanda, applicando ad essi i risultati ottenuti in Appendice D.l.

Per quanto ci riguarda, il caso di maggiore interesse è quello di un banco di filtri in grado di generare K bande uniformemente spaziate. In tal caso, la sottobanda relativa alle frequenze inferiori risulta modulata SSB. Tale disposizione delle sottobande è normalmente indicata come even-stacking, a differenza del caso in cui, pur avendo bande uniformi, la prima sottobanda è centrata nell'origine, nel qual caso si parla di struttura odd-stacking del banco di filtri.

Considerando modulazioni complesse, il segnale relativo alla banda k-esima X k (m) è quindi ottenibile come

  (D.23)
[p. 353 modifica]mentre per il filtro di sintesi, indicando con il risultato dell’elaborazione eseguita sul segnale relativo alla banda k-esixna, si ricava
  (D.24)

dove h(m) è la risposta impulsiva del filtro passa basso di analisi e f(m) quella utilizzata nella sintesi. Riscrivendo il filtro di sintesi come

  (D.25)

si vede come le uscite dei banchi di analisi e sintesi siano ottenibili tramite trasformata discreta di Fourier, da cui il nome di DFT filter banks utilizzati nel caso di banchi con bande uniformemente spaziate. Il massimo valore ammissibile per il fattore di decimazione M coincide con il numero di sottobande K, nel qual caso si parla si bande critically-sampled.

Fig. D.6 - Struttura del sistema di analisi e sintesi. [p. 354 modifica]Utilizzando filtri passa banda complessi la struttura cambia, con le modulazioni che diventano l’ultima operazione nell’analisi e la prima della sintesi ed ottenendo

 
  (D.26)

Nel caso di K bande uniformi, quindi, le funzioni di trasferimento di ciascun filtro passa banda sono derivate per traslazione di un unica funzione prototipo passa basso H(ffl), in modo tale che per la banda k-esima si abbia

  (D.27)

Un miglioramento dal punto di vista dell’efficienza computazionale si ha utilizzando filtri polifase (fig. D7). In tal caso, la risposta impulsiva del ramo di ordine p per la banda k-esima ed il suo ingresso sono pari a

  (D.28)

Secondo tale impostazione sembrerebbe che i filtri polifase da utilizzare sarebbero in numero pari al prodotto del numero di sottobande per il fattore di decimazione. Fortunatamente, però, sfrattando la relazione

  (D.29)

si ricava

  (D.30)
[p. 355 modifica]Se si definisce il p-esimo ramo dell’equivalente filtro polifase passa basso come
  (D.31)

allora è possibile riscrivere il filtro polifase della banda k-esima come

  (D.32)

Sostituendo nelle equazioni precedenti ed indicando con qρ(m) il ρ-esimo ramo del polifase passa basso equivalente del filtro di sintesi, si ottiene

  (D.33)

Fig. D.7 - Banco di filtri polifase.

Dall'analisi di tali relazioni, si nota come l’efficienza computazionale nasca, innanzitutto, dal poter condividere il filtro pp(m) tra le differenti bande, Inoltre, si nota come l’uscita del filtro di analisi sia ottenibile come DFT delle uscite dei singoli filtri polifase e che il filtro di sintesi si ottiene applicando i filtri polifase alla DFT inversa degli Xk. Un ulteriore incremento di efficienza computazionale nasce, quindi, dall'uso di FFT. [p. 356 modifica]Un problema che si incontra nel filtraggio passa banda è la presenza di aliasing, introdotto dalla decimazione, a seguito dell’uso di filtri non ideali. In tal modo l’energia in ciascuna sottobanda è influenzata dal contributo delle bande adiacenti. Questo fenomeno è tanto più sensibile per quelle bande a bassa energia, per le quali il contributo dovuto alla non idealità dei filtri può essere predominante. La realizzazione ottimale del banco dei filtri è, quindi, un punto chiave e la scelta dei filtri h(n) ed f(n) deve essere tale da minimizzare l’aliasing. È possibile analizzare tale fenomeno sia nel dominio del tempo o in frequenza.

Analizzando l’aliasing nel dominio del tempo, si vuole vedere sotto quali condizioni ne è possibile 1*eliminazione. Imporre che x(n) = x(n) è possibile sostituendo l’espressione del filtro di analisi in quello di sintesi, ottenendo

 
  (D.34)

Dato che l’esponenziale è pari ad 1 per r = n - sK, con s intero, l’espressione precedente può essere riscritta come

  (D.35)

affinché sia valida l’uguaglianza x(n) = x(n), la seconda sommatoria deve annullarsi per s * 0 e deve valere uno per s = 0. Indicando le versioni traslate delle risposte impulsive come

  (D.36)

la relazione precedente si riscrive come

  (D.37)
[p. 357 modifica]Questa equazione ci permette di definire quali siano le condizioni di perfetta

ricostruzione del segnale. Affinché sia valida l'uguaglianza x(n) = x(n),dalla prima sommatoria si ricava che

  (D.38)

il che indica che la somma del prodotto delle risposte impulsive traslate deve essere unitaria per qualsiasi valore di n. Dalla seconda sommatoria, che è necessario annullare, si comprende come sia necessario annullare le infinite repliche poste ad intervalli di K campioni, introdotte dalla decimazione. Ciò è ottenibile imponendo

  (D.39)


il che si raggiunge solo con filtri di lunghezza inferiore a K. Dato che la lunghezza dei filtri risultante da tale vincolo sarebbe estremamente modesta, essi non risulterebbero soddisfacenti dal punto di vista dell’analisi in frequenza in quanto poco selettivi. Passando, infatti, all’analisi in frequenza, è conveniente considerare un modello a filtri passa banda complessi, semplificato eliminando le modulazioni terminali (che si elidono) e sostituendo la decimazione/interpolazione con un modulatore che moltiplica il segnale per la sequenza

  (D.40)

La rappresentazione in frequenza del legame ingresso/uscita che si ottiene è pari a

  (D.41)
[p. 358 modifica]dove
  (D.42)

sono i filtri di analisi e sintesi della banda k-esima. Mentre il primo termine di tale equazione implica che la funzione di trasferimento complessiva sia unitaria, il secondo impone funzioni passa basso ideali, approssimabili solo con filtri di dimensione elevata. In conclusione, non è possibile raggiungere compietamente la cancellazione dell’aliasing né nel dominio del tempo né in frequenza.

Tutte le relazioni precedenti sono state ricavate per una struttura even-stacking del banco di filtri. Per un’organizzazione odd-stacking deve risultare

  (D.43)

Si è dimostrato come sia il filtro di analisi che quello di sintesi sono legati ad un mappaggio tramite DFT dal dominio del tempo alla frequenza e viceversa. Il mappaggio tramite DFT, però, è sempre tale da legate l’origine temporale n = 0 a quella in frequenza co = 0. Volendo permettere il mappaggio tra due qualsiasi origini nei due domini, è possibile definire una trasformata discreta di Fourier generalizzata come

  (D.44)

dove no e ko sono le nuove origini nel tempo ed in frequenza. Applicando tale generalizzazione a banchi di filtri di analisi e sintési che, nel caso di struttura odd-stacking, sono caratterizzati ko = 1/2 , si ottiene

  (D.45)
[p. 359 modifica]Anche per tale banco, la struttura può evolvere verso strutture con filtri

passa banda complessi o polifase. Inoltre, al fine di ottenere segnali reali, è possibile trasformarlo applicando una modulazione SSB

  (D.46)

da cui, posto M = K/2, deriva

 
  (D.47)

Tale relazione può essere formalmente semplificata in

  (D.48)
D.3 BANCHI DI FILTRI QMF

Una tecnica differente per la realizzazione del banco è quella dei Quadrature-Mirror Filters (QMF) (fig. D.8). Un banco QMF permette la divisione dello spettro del segnale in due sotto-bande. Ciò è ottenuto tramite un filtraggio passa-alto abbinato ad un filtraggio complementare passa-basso. Banchi con un numero maggiore di due sottobande si ottengono ripetendo il filtraggio su entrambi (per sottobande uniformi) o uno solo dei segnali d’uscita (fig. D.9). La struttura risultante è ad albero, bilanciato o meno. Al fine di minimizzare l’effetto dell’aliasing, la somma delle due funzioni di trasferimento deve essere piatta. Se si indica con H1(ejωT) e H1(ejωT) la [p. 360 modifica]Fig. D.8 - Diagramma di un banco QMF.

funzione di trasferimento del filtro passa alto e passa basso di analisi, i segnali da loro prodotti sono pari a

 
  (D.49)

Indicando con F0(ω) e F1(ω) la funzione di trasferimento del filtro passa basso e passa alto di sintesi, in fase di ricostruzione si ottiene

  (D.50)

Sostituendo in tale espressione quella dei segnali prodotti da ciascuna sottobanda si ottiene

 
  (D.51)

Imponendo che si annulli l’effetto dell’aliasing è necessario annullare il secondo termine della precedente equazione, ottenendo

  (D.52)
[p. 361 modifica]Fig. D.9 - Filtraggio passa banda tramite cascata di filtri complementari.

Per trovare le condizioni tali da soddisfare tale vincolo è necessario ridurre le variabili del problema. A tal fine si imporre innanzitutto che il filtro passa altro sia una replica traslata del passa basso di analisi

  (D.53)

da cui deriva

  (D.54)

dove il fattore di scala 2 serve per rendere unitaria la funzione di trasferimento complessiva. Infatti

  (D.55)
[p. 362 modifica]L’esser riusciti ad eliminare il secondo termine dell’equazione precedente vuol dire che l’aliasing introdotto nel filtro di sintesi da parte della decimazione è cancellato dalle repliche introdotte dall’interpolazione nel filtro di sintesi. Nel dominio del tempo, le risposte impulsive dei filtri sono date da
  (D.56)

Rimane una sola incognita, data da H(co), per la determinazione della quale si possono seguire differenti criteri [Cro83]. Una classe di filtri particolarmente diffusa è quella dei filtri FIR a fase lineare. Il vantaggio dell’assenza di distorsioni di fase risulta particolarmente utile nel caso di connessioni in cascata, eliminando la necessità di equalizzazioni ad ogni stadio. La risposta impulsiva di un filtro FIR a fase lineare risulta essere simmetrica, per cui

  (D.57)

e la funzione di trasferimento è data da una funzione reale ed un termine di fase lineare

  (D.58)

La funzione di trasferimento complessiva nel caso di N pari e dispari risulta

 
  (D.59)

Dato che con N dispari, la funzione di trasferimento si annulla per ω = π/2 in quanto H(π/2) = H(3π/2), si utilizzano solamente filtri con N pari. Fissato l’ordine del filtro, l’unico filtro che rende unitaria la funzione di trasferimento complessiva è il filtro

  (D.60)
[p. 363 modifica]le cui caratteristiche selettive, però, sono povere. Ogni altro filtro introduce delle distorsioni di ampiezza. Per la definizione di differenti funzioni di trasferimento si ricorre a procedure iterative per la contemporanea minimizzazione della distorsione introdotta dal filtro e dell’energia in banda soppressa, cioè dell’errore quadratico
  (D.61)
il cui primo termine rappresenta l’energia in banda soppressa, il secondo la distorsione in banda e w un coefficiente di pesatura. [p. 364 modifica]

Appendice E


RICHIAMI SU TRASFORMATE NUMERICHE

E.1 TRASFORMATA DI FOURIER DISCRETA

Data una funzione x(t) continua in -∞ < t < ∞, la sua trasformata ed antitrasformata di Fourier è data da

 
  (E.1)

La trasformata di Fourier può essere estesa dall’analisi dei segnali continui ai segnali a tempo-discreto (Discrete Time Fourier Transform: DTFT) tramite la sostituzione della variabile continua t con l’indice n, ottenendo

  (E.2)

La trasformata di Fourier rappresenta la sequenza x(n) come sovrapposizione di esponenziali complessi ejωn di ampiezza

  (E.3)

ed è quindi una funzione complessa continua. Inoltre, data la periodicità dell’esponenziale complesso, la DTFT risulta essere una funzione periodica in [p. 365 modifica]frequenza. Ciò si ripercuote nell’antitrasformata, nella quale gli estremi di integrazione dell’integrale di definizione sono limitati in un intervallo di ampiezza pari a 2p

  (E.4)

Es.: si consideri la DTFT di una rect(n) discreta

  (E.5)

la sua DTFT è data da

  (E.6)

II modulo della DTFT è dato dal solo rapporto tra sinusoidi, che è una funzione periodica, il cui primo zero è in ω = 2π/N. L’esponenziale è relativo alla sola traslazione della rect(n), che non risulta centrata nell'origine. È interessante confrontare questa DTFT con la trasformata di Fourier di un’analoga rect(t) continua, data da

  (E.7)

cioè una sinc(ω) il cui primo zero è in ω = 2π/T. Diagrammando le due funzioni in modo che gli zeri coincidano (che simula l’ottenimento della rect(n) dalla rect(t) campionata con N punti) si nota l’effetto dell’aliasing che fa discostare la DTFT dalla FT tanto più quanto ci si allontana dall'origine (fig. E.1).

Come la DTFT è l’analoga della trasformata, data una sequenza di N elementi periodica con periodo T, o una sequenza di durata finita periodicizzata, è possibile introdurre nel caso discreto l’equivalente della serie di Fourier (Discrete Fourier Series: DFS) ripetendo la sostituzione della variabile continua t con l’indice n. A causa della periodicità delle sequenze sinusoidali, [p. 366 modifica]Fig. E.1 - FT e DTFT della rect.

l’analisi in frequenza ha senso solamente nell'intervallo limitato [0 ≤ ω < 2π]. Per sequenze con pulsazione ω0 rappresentate da N punti, inoltre, si possono considerare solamente le prime N armoniche, per cui gli indici delle sommatorie nelle definizioni della serie vengono limitati ad N.

 
  (E.8)

dove T0=T/N e WN=e-j(2π/N) è radice n-esima dell’unità interpretabile come campionamento di funzioni trigonometriche. La limitazione in frequenza comporta un effetto di aliasing che, come per la DTFT, fa si che i coefficienti della serie siano differenti dagli analoghi dei sistemi continui. Il calcolo numerico [p. 367 modifica]della serie di Fourier è ricondotto alla soluzione del sistema di N equazioni espresse dalla definizione, una per ogni valore di m.

Esistendo una relazione diretta tra coefficienti della serie e valori dello spettro, le equazioni precedenti stabiliscono una corrispondenza tra N campioni di una sequenza x(n) e N campioni di una trasformata X(k), per cui la serie di Fourier per segnali tempo discreto è anche interpretabile come trasformata discreta di Fourier (DFT), definita come: