短時間フーリエ変換
[Wikipedia|▼Menu]
.mw-parser-output .ambox{border:1px solid #a2a9b1;border-left:10px solid #36c;background-color:#fbfbfb;box-sizing:border-box}.mw-parser-output .ambox+link+.ambox,.mw-parser-output .ambox+link+style+.ambox,.mw-parser-output .ambox+link+link+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+style+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+link+.ambox{margin-top:-1px}html body.mediawiki .mw-parser-output .ambox.mbox-small-left{margin:4px 1em 4px 0;overflow:hidden;width:238px;border-collapse:collapse;font-size:88%;line-height:1.25em}.mw-parser-output .ambox-speedy{border-left:10px solid #b32424;background-color:#fee7e6}.mw-parser-output .ambox-delete{border-left:10px solid #b32424}.mw-parser-output .ambox-content{border-left:10px solid #f28500}.mw-parser-output .ambox-style{border-left:10px solid #fc3}.mw-parser-output .ambox-move{border-left:10px solid #9932cc}.mw-parser-output .ambox-protection{border-left:10px solid #a2a9b1}.mw-parser-output .ambox .mbox-text{border:none;padding:0.25em 0.5em;width:100%;font-size:90%}.mw-parser-output .ambox .mbox-image{border:none;padding:2px 0 2px 0.5em;text-align:center}.mw-parser-output .ambox .mbox-imageright{border:none;padding:2px 0.5em 2px 0;text-align:center}.mw-parser-output .ambox .mbox-empty-cell{border:none;padding:0;width:1px}.mw-parser-output .ambox .mbox-image-div{width:52px}html.client-js body.skin-minerva .mw-parser-output .mbox-text-span{margin-left:23px!important}@media(min-width:720px){.mw-parser-output .ambox{margin:0 10%}}

この記事は検証可能参考文献や出典が全く示されていないか、不十分です。出典を追加して記事の信頼性向上にご協力ください。(このテンプレートの使い方
出典検索?: "短時間フーリエ変換" ? ニュース ・ 書籍 ・ スカラー ・ CiNii ・ J-STAGE ・ NDL ・ dlib.jp ・ ジャパンサーチ ・ TWL(2023年1月)

短時間フーリエ変換(たんじかんフーリエへんかん、short-time Fourier transform、short-term Fourier transform、STFT)とは、関数に窓関数をずらしながら掛けて、それにフーリエ変換すること。音声など時間変化する信号の周波数位相(の変化)を解析するためによく使われる。

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

STFTは以下のように数式表現できる( i {\displaystyle i} は虚数単位):

S T F T x , w ( t , ω ) = ∫ − ∞ ∞ x ( τ ) w ( τ − t ) e − i ω τ d τ {\displaystyle \mathrm {STFT} _{x,w}(t,\omega )=\int _{-\infty }^{\infty }x(\tau )w(\tau -t)e^{-i\omega \tau }\,d\tau }

ここで w ( t ) {\displaystyle w(t)\,} は窓関数であり、普通 t = 0 {\displaystyle t=0\,} 付近に中心をもつ山の概形をしていて、 t = 0 {\displaystyle t=0\,} の付近から離れると 0 になる関数である。また x ( t ) {\displaystyle x(t)\,} は変換される関数(信号)である。そして、 S T F T x , w ( t , ω ) {\displaystyle \mathrm {STFT} _{x,w}(t,\omega )\,} は時刻 t {\displaystyle t\,} 角周波数 ω {\displaystyle \omega \,} のスペクトルを表現する複素数である。

離散時間に関するSTFTは次のようになる:

S T F T x , w [ n , ω ] = ∑ m = − ∞ ∞ x [ n + m ] w [ m ] e − i ω m {\displaystyle \mathrm {STFT} _{x,w}[n,\omega ]=\sum _{m=-\infty }^{\infty }x[n+m]w[m]e^{-i\omega m}}

連続時間での式と同様に w [ n ] {\displaystyle w[n]\,} は窓関数、 x [ n ] {\displaystyle x[n]\,} は変換される関数である。この式では n {\displaystyle n\,} が離散値であるが ω {\displaystyle \omega \,} が連続でもよい。しかし通常高速フーリエ変換を用いて計算機で計算されるので、 ω {\displaystyle \omega \,} も離散化される。

STFTの絶対値を2乗することで、パワースペクトルの時間変化が得られる:

p o w e r x , w ( t , ω ) = 。 S T F T x , w ( t , ω ) 。 2 {\displaystyle \mathrm {power} _{x,w}(t,\omega )=|\mathrm {STFT} _{x,w}(t,\omega )|^{2}\,}

また、位相スペクトルの時間変化は、STFTの偏角で得られる:

p h a s e x , w ( t , ω ) = arg ⁡ S T F T x , w ( t , ω ) {\displaystyle \mathrm {phase} _{x,w}(t,\omega )=\arg \mathrm {STFT} _{x,w}(t,\omega )}
STFTの問題点である不確定性原理

不確定性原理(フーリエ変換の不確定性原理)とは、時刻の不確定さと周波数の不確定さの間に

Δ x Δ ω ≥ 1 2 {\displaystyle \Delta x\Delta \omega \geq {\frac {1}{2}}}

の関係があることである。一般化された言い方では、フーリエ変換で結ばれた2つの変数の対に対して上のような関係がなりたつことを指す。

STFTの問題点の一つは解像度が限られてしまうことである。窓関数の窓の幅などの形状によって、周波数分解能を良くするか時間分解能を良くするかのトレードオフが決まってしまう。幅の広い窓は周波数分解能が良いが時間分解能は悪い。


次ページ
記事の検索
おまかせリスト
▼オプションを表示
ブックマーク登録
mixiチェック!
Twitterに投稿
オプション/リンク一覧
話題のニュース
列車運行情報
暇つぶしWikipedia

Size:19 KB
出典: フリー百科事典『ウィキペディア(Wikipedia)
担当:undef