Visualizza wikitesto di Elaborazione numerica dei segnali (comdig) T1/2006-2007

[[Categoria:Corsi 2006-2007]]
<!-- non rimuovere la categorizzazione precedente -->

==INFORMAZIONI GENERALI==
*Docente: '''Paola Campadelli'''
*Luogo ed orario: '''LUN 8.45-10.30 aula 405, GIO 10.30-12.30 aula G24'''
*Sito del professore: [http://homes.dsi.unimi.it/~campadel/DSP/index.htm]
*Info sul professore (webcen): [http://www.dico.unimi.it/persona.php?z=0;id_persona=28]
----




----

== Lezione 1 - 2 ottobre 2006 ==
*'''Segnale'''= veicolo per le informazioni, strumento. Viaggia su un canale di comunicazione. Non ha bisogno di alcuna ''semantica'' (grammatica) poichè è solo un'onda pressoria (con probabile disturbo-rumore). Inteso anche come variazione di una grandezza fisica.

Il segnale può essere rappresentato tramite una funzione (con relativi dominio e codominio).
Un'ipotetica funzione-segnale può essere la classica f:R<math>\rightarrow</math>R (funzione:dominio<math>\rightarrow</math>codominio), cioè una funzione che ha come dominio e come codominio una retta.
(NB: un'immagine ha come dominio il piano reale vale a dire RxR, cioè <math>R^2</math> e come codominio la retta reale cioè R).
Il '''suono''' è un tipo di segnale che varia in una sola dimensione: il tempo.

Considerando f:R<math>\rightarrow</math>R ed x variabile, t variabile: f(x) e f(t) hanno la stessa '''semantica''', cioè la funzione è sempre la stessa, ha lo stesso comportamento sia su x che su t, nonostante queste ultime siano differenti.

*'''Costante'''= un numero fisso (anche se potrebbe non essere espresso il valore)
*'''Variabile'''= può denotare un numero, assume valori in un determinato dominio.
'''Una variabile NON è una costante!'''

*'''Segnale (2)'''= funzione da un dominio a un codominio. ''Se un segnale è definito R<math>\rightarrow</math>R è un segnale continuo.''

*'''Segnale continuo o analogico'''= è un segnale a valori continui ed a tempo continuo. Questo genere di segnale ha però valori '''infiniti'''

Come rendere FINITO il segnale?
*'''Campionamento'''= divisione del tempo in un tot di '''campioni''', rendendo il tempo non più una serie di valori numerici infinita, ma una serie di valori finita. Il tempo diventa '''discreto''' (finito).
''Il segnale campionato diventa una funzione f:Z<math>\rightarrow</math>R.'' Poniamo di prendere un campione <math>\tau</math>(lettera greca minuscola tau): tutti gli altri campioni saranno multipli di <math>\tau</math> (cioè <math>2\tau</math>, <math>3\tau</math>, ecc...).
Il '''segnale campionato''' si dice a '''tempo discreto''' e a '''valori continui''': valori continui poichè, pur dividendo il segnale in campioni, fra un campione e l'altro sussistono ancora valori di segnale (infiniti).
In un'immagine (bidimensionale) dopo un campionamento avrei f:<math>Z^2</math><math>\rightarrow</math>R.

*'''Quantizzazione'''= divisione della funzione (dei valori della funzione) in valori finiti: avrò, per esempio, da una serie infinita di numeri, solo "pochi" numeri (ad esempio la serie 1 2 4 6) che possono ripetersi (ad esempio la mia funzione potrà avere valori 1 2 4 6 6 6 4 1...).

Il segnale (solamente) quantizzato si dice a '''valori discreti''' e a '''tempo continuo'''.
La funzione (solamente) quantizzata è definita come f:R<math>\rightarrow</math>Z.

Se noi eseguiamo la seguante (con f:R<math>\rightarrow</math>R, quindi segnale continuo):
'''Segnale<math>\rightarrow</math>Campionamento<math>\rightarrow</math>Quantizzazione<math>\rightarrow</math>Segnale  Digitale'''
otteniamo, per l'appunto, il '''segnale digitale''', una funzione definita in Z<math>\rightarrow</math>Z, cioè a '''tempo discreto''' e a '''valori discreti'''.
----
'''ESERCIZIO'''
Quanti bit sono necessari a memorizzare N valori interi?
(Si prenderanno in esempio i livelli di grigio di un'immagine, cioè 256 valori)

N= 256= {0, 1, ..., 255}

Dato che   <math>log_2</math>256= 8, serve 1 Byte (8 bit) per memorizzare 256 valori.

Infatti   <math>2^8</math>= 256.
----
'''ESERCIZIO'''
Dato un certo intervallo di campionamento tau <math>\tau</math>, sapendo che il reciproco di <math>\tau</math> è <math>1\over\tau</math> =<math>\nu</math> (nu): calcolare la frequenza, dato <math>\tau</math>=1000.

Troviamo   <math>1\over\tau</math>=  <math>1\over1000</math>=  1 KHz
----
'''ESERCIZIO'''
Con un intervallo di durata T, un intervallo tra campioni <math>\tau</math> ed N valori possibili calcolare i bit necessari per rappresentare il segnale con i valori dati.

La formula è   <math>T\over\tau</math><math>log_2</math>N (in bit).
----




----

==Lezione 2 - 5 ottobre 2006==
----
'''ESERCIZIO'''
Abbiamo 3 minuti di musica campionata a 40'000 Hz. Prendo 256 valori per campione. Quanti bit ho?

3x60 sec x40'000x<math>log_2</math>256= 57'600'000

Considerazione: se nel logaritmo mi dovesse venire un numero indeciso (ad es 7,3) prendo comunque il numero arrotondato per eccesso (ad es 8).
----
*'''Sistema'''= qualsiasi dispositivo che prenda un segnale in ingresso e che produca in uscita un segnale differente.
*'''Sistema deterministico'''= un sistema nel quale non venga mai cambiato il comportamento, nonostante l'input differente.
*'''DSP'''= Digital Signal Processor, strumento usato per realizzare nel modo più efficiente possibile un sistema. Prende in input un segnale continuo e lo immette in un blocco chiamato '''ADC''' (Analog to Digital Converter) trasformando il segnale da analogico a digitale.
*'''ADC'''= costituito da un campionatore, da un quantizzatore e da un filtro iniziale ('''filtro di input''').
*'''DAC'''= Digital to Analog Converter.
Se anche noi ipotizzassimo di inviare un ipotetico segnale f(t) all'interno di un ADC<math>\rightarrow</math>DSP<math>\rightarrow</math>DAC, in uscita non avremo mai di nuovo f(t), bensì un segnale g(t) simile a f ma mai uguale. Infatti, ADC fa perdere informazioni del segnale e crea quindi un segnale distorto.
*'''Filtro'''= elimina alcune informazioni.

===Operazioni con le funzioni - il caso analogico===
*'''Addizione'''= f(t)+g(t)
*'''Sottrazione'''= f(t)-g(t)
*'''Prodotto'''= f(t)g(t)
*'''Moltiplicazione per un reale'''= af(t), :a<math>\in</math>R. Fa in modo di cambiare tutti i valori della funzione.
*'''Composizione'''= f(g(t)) oppure f<math>\circ</math>g. Viene applicata prima la funzione g(t) (per semplicità poniamolo =y) e dopodichè viene applicata f(y).
*'''Variazione di scala'''=f(at). La funzione si dilata o si comprime.
*'''Traslazione'''= f(t-k), :k<math>\in</math>R. ''Contrariamente a quello che si potrebbe pensare, quando la traslazione è a -k, la funzione si sposta (sul piano delle x) '''verso destra''', mentre se è a +k allora si sposta verso sinistra.''

NB. Purtroppo, ma solo per il momento, non sono disponibili esempi grafici.

*'''La funzione rettangolo rect(t)'''= E' una funzione che '''vale 1 se x<math>\leq</math>1/2, e 0 altrimenti'''. Dal punto di vista grafico, è un rettangolo di altezza 1 che ha metà del suo "corpo" a destra dell'asse delle y e metà a sinistra, ed è poggiato (ha la base poggiata) sull'asse x. La base di questo rettangolo è lunga a ed ha i due lati di altezza 1 nelle posizioni -a/2 ed a/2 (sull'asse x ovviamente).
La Rect è importante per quanto riguarda i trattamenti sui segnali perchè posso, facendo f(t)rect(t) otteniamo un "pezzo" della funzione f(t) compreso nella regione di piano della rect(t) (ritagiando quindi il segnale).

'''''A rect((t-b)/a)''' è una funzione rettangolo di ampiezza A, centrato in b (traslato) con base a.''
----




----

==Lezione 3 - 9 ottobre 2006==
Parte grafica riguardante le funzioni trigonometriche seno e coseno (circonferenza in radianti o gradi).

*'''Funzione periodica'''= f(t) è periodica di un periodo T se f(t)=f(t+T), cioè se si ripete uguale dopo il periodo T.

Se f(t) è una funzione periodica di periodo T la funzione g(t)=f(at) (cioè variata di scala) ha periodo T'=T/a.
DIMOSTRAZIONE: g(t)=f(at) per definizioe. g(t+T')=g(t). f(a(t+T'))=f(at) applico la definizione. f(at+aT')=f(at) f è periodica di t quindi aT'=T quindi T'=T/a.

Per trovare la frequenza espressa in cicli al secondo o in radianti:

'''Funzione'''       -   '''Periodo'''        -  '''Freq in cicli/sec'''    -  '''Freq in radianti'''

sen (t)      -    2<math>\pi</math>     -     1/2<math>\pi</math>   -   1

sen (2<math>\pi</math>t)         - 2<math>\pi</math>/2<math>\pi</math>= 1      -    1      -2<math>\pi</math>

sen (10<math>\pi</math>t)     -     1/5     -   5    -    10<math>\pi</math>= 5x2<math>\pi</math>

sen (1/2t)     -     4<math>\pi</math>      -  1/4<math>\pi</math>    -  1/2


NB: 1 Hz= 2<math>\pi</math> rad/sec. 1 rad/sec= 1/2<math>\pi</math> Hz.

'''f(t)= A sen (2<math>\pi\nu</math>t + <math>\phi</math>)''', con <math>\phi</math> angolo ('''fase istantanea''') ed f periodica.

Esercizi ed esempi grafici non disponibili.

----
'''ESERCIZIO'''
f(t)= 2sen(t) se 0<math>\leq</math>t<math>\leq</math>3<math>\pi</math> e 0 altrimenti.
CAMPIONARLA ad un intervallo di campionamento <math>\pi</math>/2 (il nostro <math>\tau</math>).
QUANTIZZARLA con 3 valori per campione.

Campionando f(n<math>\tau</math>) = [f(0), f(1x<math>\pi</math>/2), f(2x<math>\pi</math>/2), f(3x<math>\pi</math>/2), f(4x<math>\pi</math>/2), f(5x<math>\pi</math>/2), f(6<math>\pi</math>/2)] Vettore

Quindi campionando ho 3<math>\pi</math>/(<math>\pi</math>/2)= 6 campioni più lo 0.
La funzione campionata a <math>\pi</math>/4 sarebbe stato ancora meglio (avrei preso il doppio dei campioni e perciò avrei un segnale più preciso!).

Ora quantiziamo con 3 valori per campione. Fortunatamente con 3 valori non ha problemi: f(0)= 0, f(<math>\pi</math>/2)= 2, f(<math>\pi</math>)= 0, f(3/2<math>\pi</math>)= -2, f(2<math>\pi</math>)= 0, f(5/2<math>\pi</math>)= 2, f(3<math>\pi</math>)= 0.
Quindi posso utilizzare 0, -2 e 2 e non avere perdita. (Poniamo di aver quantizzato con i 3 valori ma di avere campionato a <math>\pi</math>/4. Avrei avuto perdita perchè avrei dovuto approssimare maggiormente e il segnale sarebbe cambiato molto).
----




----

==Lezione 4 - 12 ottobre 2006==
*'''Funzione sinc(t)'''= sinc (t)= sen(<math>\pi</math>t)/<math>\pi</math>t, è una funzione che oscilla perchè ha sen al nominatore ma quando t aumenta il valore della funzione diminuisce. Per ogni valore di t intero ho 0. Invece, in t=0 ho il valore della funzione pari a 1. Inoltre la funzione è simmetrica, cioè sinc(t)=sinc(-t). E' inoltre limitata e continua.

*'''Impulso o Delta di Dirac <math>\delta</math>(t)'''= '''<math>\lim_{\Delta \to 0}P_\Delta (t)</math>''' dove '''<math>P_\Delta(t)</math>=1/<math>\Delta</math>rect(t/<math>\Delta</math>)'''. La Delta di Dirac '''ha ''sempre'' area 1''', poichè, dal punto di vista grafico, è un rettangolo con una base che tende a 0 (piccolissima) ed un'altezza sempre più alta (proprio perchè mantiene l'area uguale a 1). Infatti, '''<math>\int_{-\infty}^{\infty} \delta(t) dt = 1</math>'''. Spesso l'impulso viene rappresentato con una freccia verso l'alto sulle x (in posizione dell'eventuale traslazione).

===Operazioni sui segnali===
Le operazioni sui segnali vengono fatte attraverso i '''sistemi''', che possono effettuare tutte le operazioni che abbiamo visto: amplificazione (Af(t)), traslazione (f(t-k) o f(t+k)), campionamento (f(n<math>\tau</math>)) e la quantizzazione (argmin|x-y|).

I sistemi possono essere composti come le funzioni, e con più segnali posso generalizzare il comportamento di segnali composti tramite una 
*'''combinazione lineare''':   <math>\sum_{i=1}^n a_if_i(t)</math>.

L'operazione di composizione non è sempre commutativa, mentre traslazione ed amplificazione lo sono. Scalatura e traslazione non sono commutative.

Considereremo i segnali come dei vettori.
----



----

==Lezione 5 - 16 ottobre 2006==
'''Proprietà di uno spazio vettoriale''': ( , ): VxV<math>\rightarrow</math>R
*'''(v1,v)<math>\geq</math>0'''
*'''(v1,v2)=(v2,v1)'''
*'''(av1+bv2, v3)=a(v1,v3)+b(v2,v3)

(Parte mancante...)

----




----

==Lezione 6 - 19 ottobre 2006==
Due vettori v1 e v2 sono *'''Ortogonali''' se (v1,v2)=0, cioè se il loro '''prodotto interno''' è 0.
*'''Base'''= è un insieme, dato uno spazio di grandezza n, di '''vettori linearmente indipendenti''' {v1, v2,... vn}.
*'''Base ortogonale'''= quando tutti i vettori della base sono tra loro '''ortogonali a coppie''' cioè ho (vi, vj)=0 oppure norma pari a 1.

----
'''ESEMPIO'''
Il prodotto interno di {(1,0), (0,1)}= 1x0 + 0x1= 0.
La norma è <math>\sqrt{1x1 + 0x0}</math>= 1 cioè del prodotto internop del vettore con se stesso.
----
'''ESERCIZIO'''
Consideriamo lo spazio vettoriale di <math>R^4</math> e il vettore (2, 5, 3, 7).
Voglio trovare la norma e l'angolo rispetto al vettore (1, 0, 0, 0).

Norma |(2, 5, 3, 7)|= <math>\sqrt{(2, 5, 3, 7) + (2, 5, 3, 7)}</math>= <math>\sqrt{2^2+5^2+3^2+7^2}</math>
L'angolo lo calcolo dalla formula (v1, v2)= |v1||v2|cos<math>\alpha</math> da cui cos<math>\alpha</math>= (v1, v2)/|v1||v2|.
----

Se V è uno spazio vettoriale a dimensione finita e {v1,...vn} è una base ortogonale per V allora per ogni vettore v che appartiene a V posso scrivere v come '''combinazione lineare''' degli elementi della base v= <math>\sum_{i=1}^n a_i v_i(t)</math> dove vi è l'elemento della base. L'elemento ai lo posso scrivere come (v, vi)/(vi, vi).

continua.......

==Lezione 7 - 23 ottobre 2006==
*Concetto di supporto finito del segnale
*Combinazione lineare di segnali
*Complementi sulla Delta di Dirac e dimostrazione
*Sistemi lineari ed operazioni dei sistemi lineari

==Lezione 8 - 30 ottobre 2006==
*Spazi vettoriali a dimensione finita con base ortogonale
*Spiegazioni più approfondite sui Sistemi Lineari
*Esercizi vari
*Prodotto di convoluzione
*Sistemi Lineari Tempo Invarianti (LTI)
*Il caso discreto (LTI a tempo discreto)
*Risposta all'impulso

==Lezione 9 - 2 novembre 2006==
*Il prodotto di convoluzione 
*Proprietà del prodotto di convoluzione nel caso discreto e nel caso continuo: commutativo, associativo e distributivo rispetto alla somma
*Prodotto di convoluzione e risposta all'impulso
*Esercizi d'esempio

==Lezione 10 - 7 novembre 2006==
*Sistemi LTI causali e non causali (caso discreto e caso continuo)
*Filtro FIR
*Funzione gradino unitario U(n)
*Esercizi su filtro FIR causale
*Compito da fare su filtro FIR non causale
*Quali tra le funzioni che abbiamo visto sono LTI e quali non lo sono
*Sistema stabile (dimostrazioni)