matplotlibとscipyのフーリエ変換(FFT)を用いてパーシバルの定理を確認する方法 - Qiita

python でパワースペクトルを計算する方法 matplotlib と scipy で FFT を計算する方法 matplotlib と scipy のどちらもパワースペクトルを計算する方法があるが、両者でデフォルトで持っている関数と定義が違う。大概の使い分けは次の通り。 単純にパワースペクトルを計算する場合、mlab.psd が一番簡単。 scipy.fftpack を使うと実部と虚部が得られる。位相の計算など細かいことが可能だか、規格化は自分で行う必要あり。 どちらでやっても、パーセバルの定理から全パワーは同じになるが、両者の規格化の違いを考えないとパワースペクトルの縦軸が完璧に同じにならない。 ここでは単純な三角関数を両方の方法でフーリエ変換して、完璧にmatplotlibとscipyのFFTで結果が同じになる例を示す。おまけで、パワースペクトルから三角波の振幅を求めて、それがもと