16.06.2013 Views

Elaborazione Numerica dei Segnali

Elaborazione Numerica dei Segnali

Elaborazione Numerica dei Segnali

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

11.1. Formulazione nel Dominio delle Frequenze 217<br />

dove SYY(ω) = E[G(ω)G ∗ (ω)] è lo spettro di potenza di Y (t) e SXY(ω) = E[F (ω)G ∗ (ω)]<br />

è lo spettro di potenza incrociato di X(t) e Y (t).<br />

Il filtro di Wiener W (ω) che minimizza la funzione dist[G(ω)W (ω), F (ω)] verifica allora:<br />

Possiamo quindi concludere:<br />

W (ω)SYY(ω) − SXY(ω) = 0.<br />

Fatto 11.1 Il filtro di Wiener, che minimizza la distanza nel senso del quadrato dell’errore<br />

tra la risposta a un processo Y (t) e un processo X(t), ha la funzione di trasferimento W (ω):<br />

W (ω) = SXY(ω)<br />

SYY(ω) .<br />

dove SXX(ω) è lo spettro di potenza di X(t) e SXY(ω) lo spettro di potenza incrociato tra<br />

X(t) e Y (t).<br />

Allo scopo di dimostrare la (11.2), richiamiamo brevemente la nozione di derivata<br />

complessa. Se z = a + ib e Z(a, b) = A(a, b) + iB(a, b), allora la derivata complessa di Z<br />

rispetto a z è:<br />

In particolare vale che:<br />

∂z<br />

∂z<br />

= 0,<br />

∂Z<br />

∂z<br />

= ∂Z<br />

∂a<br />

∂z ∗<br />

∂z<br />

= 2,<br />

+ i∂Z<br />

∂b .<br />

∂zz ∗<br />

∂z<br />

= 2z.<br />

∂<br />

Applicando le regole sopra descritte al termine ∂W (ω) dist[G(ω)W (ω), F (ω)], si ottiene:<br />

∂E[e(ω)e ∗ (ω)]<br />

∂W (ω)<br />

=<br />

=<br />

Esempio 11.1.1<br />

∂<br />

∂W (ω) E[(G(ω)W (ω) − F (ω))(G∗ (ω)W ∗ (w) − F ∗ (ω))]<br />

∂<br />

∂W (ω) (W (ω)W ∗ (ω)E[G(w)G ∗ (ω)] − W (ω)E[G(ω)F ∗ (ω)])<br />

= 2W (ω)E[G(ω)G ∗ (ω)] − 2E[F (ω)G ∗ (ω)]<br />

= 2(W (ω)SYY − SXY).<br />

−<br />

∂<br />

∂W (ω) W ∗ (ω)E[G ∗ (ω)F (ω)]<br />

Filtro per la riduzione di rumore additivo. Consideriamo un processo Y (t) ottenuto<br />

contaminando un processo X(t) con rumore additivo N(t), a media 0 e scorrelato da<br />

X(t):<br />

Y (t) = X(t) + N(t).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!