短時間フーリエ変換

ウィキペディア

短時間フーリエ変換

出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2025/02/16 05:23 UTC 版)

短時間フーリエ変換 （たんじかんフーリエへんかん、（英: short-time Fourier transform, STFT）とは、関数に窓関数をずらしながら掛けて、それにフーリエ変換すること。

概要

短時間フーリエ変換（以下、STFT）は連続あるいは離散の信号に時間シフトした窓関数を掛けたうえでフーリエ変換する操作である（⇒ #定義）。これは信号の一部区間を取り出して周波数解析しそれを各時刻でおこなうと解釈できるため（⇒ #解釈）、音声など非定常信号の周波数と位相の変化を解析するためによく使われる。

理論上フーリエ係数を求めるには無限の区間に渡って積分を行わなければならないが、実験値等からフーリエ係数を求めるには範囲を区切らなければならない。そのために、ある範囲の実験値のフーリエ係数を求めるには、このある範囲の実験値が周期的に無限に繰り返されていると仮定して計算するのが一般的である。だがここで問題なのは、ある範囲の最初の値と最後の値を無理やりつなげることによって発生する不連続な要素である。これを解決するため、中央が 1 付近の値でその範囲外で 0 に収束する関数を掛けて、不連続な要素を極力排除することが行われる。これが短時間フーリエ変換である。このとき、この掛け合わせる関数を窓関数と言う。

定義

定義に先立ち、以下のように連続/離散での記号を定める：

対象	記号
	連続	離散
時刻	$t$ STFTの解像度の比較。左は時間分解能が良く、右は周波数分解能が良い この事実はウェーブレット変換を作る原因にもなった。ウェーブレット変換ではSTFTと異なり時間分解能と周波数分解能が両立することが出来る。 量子力学における運動量と位置に関するハイゼンベルクの不確定性原理とは普通区別されるが、実はフーリエ変換の不確定性原理に基因するものである。シュレディンガー方程式によれば、定常な場合、 $\psi (x,t)=\psi (x)e^{-{iEt \over {\hbar }}}=\psi (x)e^{-i\omega t}$ 25ミリ秒の窓 125ミリ秒の窓 375ミリ秒の窓 1000ミリ秒の窓 逆短時間フーリエ変換 逆短時間フーリエ変換 （ぎゃくたんじかんフーリエへんかん、（英: inverse short-time Fourier transform, ISTFT）は時間周波数表現を時間表現へ逆変換する操作である。各フレームのフーリエ係数を逆変換して時間表現へ戻し、合成窓 $w_{s}$ を掛け、フレーム群を位置合わせしながら重畳加算することで1つの時間表現を合成する。 時間的に重なった信号を重畳加算するため、元信号の完全再構成は一般に可能でない。 $w_{s}$ に矩形窓を採用して複数のフレームを重ねると元信号より大きくなることからこれは明らかである。「信号 →（STFT）→ 係数 →（ISTFT）→再構成信号」で完全な再構成を可能にする条件を 完全再構成条件 （かんぜんさいこうせいじょうけん）という。 派生 STFTの絶対値を2乗することで、パワースペクトルの時間変化が得られる： $\mathrm {power} _{x,w}(t,\omega )=|\mathrm {STFT} _{x,w}(t,\omega )|^{2}\,$ また、位相スペクトルの時間変化は、STFTの偏角で得られる： $\mathrm {phase} _{x,w}(t,\omega )=\arg \mathrm {STFT} _{x,w}(t,\omega )$ 参考文献 小野, 順貴 (2016). “短時間フーリエ変換の基礎と応用”. 日本音響学会誌. 72 (12): 764–769. doi:10.20697/jasj.72.12_764. 矢田部, 浩平 (2021). “第三回：短時間フーリエ変換”. 日本音響学会誌. 77 (6): 396–403. doi:10.20697/jasj.77.6_396. 関連項目 フーリエ変換 窓関数 不確定性原理（ハイゼンベルクの不確定性原理） ウェーブレット変換