modulazione_claudia.py

# istruzioni per claudia (e il gatto):
# 1. salva questo file come "sdr.py" e mettilo nella stessa cartella dove hai
#    estratto i file di input (dati_re.txt e dati_im.txt).
# 2. Apri un terminale (Ctrl-Alt-T per fare prima).
# 3. Spostati nella cartella che contiene i file con:
#    $ cd percorso/cartella
#  ( ^ il carattere dollaro $ serve a far capire che devi eseguire il comando
#    nel terminale. NON devi scriverlo. )
#    Ad esempio, se hai messo tutto nella cartella "Downloads":
#    $ cd Downloads
#  ( Suggerimento: se premi "tab", una o due volte, mentre scrivi il percorso
#    della cartella, ti vengono forniti dei suggerimenti.)
# 4. Verifica che ci sono tutti i file necessari (sdr.py, dati_re.txt,
#    dati_im.txt) usando "ls":
#    $ ls
# 5. Aggiusta il file dati_im.txt che è in formato sbagliato, sostituendo
#    tutti i caratteri "tab" (\t) con dei caratteri "newline" (\n):
#    $ sed 's/\t/\n/g' dati_im.txt > dati_im_corretto.txt
#  ( Suggerimento: be patient.)
# 6. Sostituisci il file corretto con quello sbagliato.
#    $ mv dati_im.txt dati_im_cacca.txt
#    $ mv dati_im_corretto.txt dati_im.txt
#  ( Suggerimento: esegui "ls" per verificare che hai combinato.)
#  ( Suggerimento: volendo puoi eliminare dati_im_cacca.txt con
#    $ rm dati_im_cacca.txt
#    Altrimenti lascialo pure lì.)
# 7. Genera i campioni \o/
#    $ python sdr.py > out.txt
#    Se è andato tutto bene i campioni sono in out.txt
#    Puoi vederli eseguendo "less out.txt" oppure puoi vedere gli ultimi campioni con
#    "tail out.txt" oppure fai un po' quello che ti aggrada.
# THE END.
#
# Appendice stupida: potevi estrarre i file anche con "unzip dati.zip".

from itertools import izip
import math

def dati():
    for r, i in izip(open('dati_re.txt'), open('dati_im.txt')):
        yield complex(float(r), float(i))

def proc():
    input = dati()
    j = next(input)
    for i in input:
        x = i * j.conjugate()
        yield math.atan2(x.imag, x.real)
        j = i

if __name__ == '__main__':
    for x in proc():
        print(x)