質問: signal.periodogram(y, fs) plt.figure() plt.plot(freq1, 10*np.log10(P1), "b", label="periodogram") plt.ylim(-60, 0) plt.legend(loc="upper right") plt.xlabel("Frequency[Hz]") plt.ylabel("Power/frequency[dB/Hz]") plt.show() ```
回答: > > freq1とP1は同じものであるという事でしょうか? > > > いいえ、違います。 これはsignal.periodogramが返す2つの戻り値をそれぞれfreq1とP1に代入している記述です。 関数の仕様を見るとf(周波数[ndarray型]), Pxx(パワースペクトル[ndarray型])の2つの戻り値を返すようです。 なので、freq1には周波数、P1にはパワースペクトルが入ります。 scipy.signal.periodogramの仕様 https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.signal.periodogram.html