Potreste testare per me questa funzione Matlab?

[anton89]

Ciao a tutti,
ho scritto una funziona che stima lo spettro di potenza di un segnale. Quando dico "spettro di potenza", intendo che sull'asse delle ordinate ho la potenza e NON la densità di potenza. E' importante questa cosa!

Inserisco il file .m come allegato a questo post.

Mi piacerebbe che la testaste perché io non sono convinto che le stime siano accurate. Ho notato che se genero delle sinusoidi di potenza determinata, lo spettro me le sbaglia di un paio di dBm massimo. Genero del rumore a -10 dBm e lui me lo visualizza a -30 dBm
Un'altra cosa che ho notato è che la stima sembra più accurata se non uso una potenza del 2 per il calcolo della fft e nessuna finestratura.

Insomma, non mi convince tanto. Spero che possiate aiutarmi a capire.
Grazie.

P.S.: la funzione non visualizza un grafico. Utilizzate la funzione plot(10*log10(X)+30) per visualizzare le potenze correttamente