sto il porting del codice MATLAB per Python usando SciPy e rimasto bloccato con la seguente riga:Periodogramma in Octave/Matlab vs SciPy
Matlab/Octave codice
[Pxx, f] = periodogram(x, [], 512, 5)
codice Python
f, Pxx = signal.periodogram(x, 5, nfft=512)
Il problema è che ottengo output diversi sugli stessi dati. Più specificamente, i vettori Pxx sono diversi. Ho provato diverse finestre per signal.periodogram, eppure non ho avuto fortuna (e sembra che la finestra box car di default scypy sia la stessa della finestra rettangolare di matlab di default) Un altro strano comportamento è che in python, il primo elemento di Pxx è sempre 0, non importa quali dati l'input è.
Mi manca qualcosa? Qualsiasi consiglio sarebbe molto apprezzato!
semplice codice/Octave Matlab con dati reali: http://pastebin.com/czNeyUjs
Python Semplice + codice SciPy con i dati reali: http://pastebin.com/zPLGBTpn
Una trama sarebbe molto utile. –