ウィグナー分布

ウィグナー分布 (テンプレート:Lang-en-short、 WDF) は、信号処理の分野でテンプレート:仮リンクに用いられる変換である。

ウィグナー分布はもともと、1932年にユージン・ウィグナーにより古典統計力学への量子補正として提案され、テンプレート:仮リンクにおいて重要である（ウィグナー関数、またはウィグナー・ビレ分布も比較参照のこと）。

代数的に、位置-運動量の関係は時間-周波数の関係と同様に正準共役関係にあるので、この変換は信号処理の分野において時間-周波数解析に用いられる。テンプレート:仮リンクなどの短時間フーリエ変換に比べて、ウィグナー分布はより明瞭な結果を与える場合がある。

数学的定義

ウィグナー分布の定義にはいくつかの異なる流儀がある。ウィーナー＝ヒンチンの定理を参照されたい。以下に示す定義は時間周波数分析特有のものである。時系列信号 $x [t]$ に対する自己相関関数は次のように定義される。

C_{x} (t_{1}, t_{2}) = ⟨ (x [t_{1}] - μ [t_{1}]) {(x [t_{2}] - μ [t_{2}])}^{*} ⟩,

ここで $⟨ ⟩$ は全ての可能なプロセスにわたっての平均を意味し、 $μ (t)$ は時間に依存する、もしくは依存しない平均値を意味する。ウィグナー分布 $W_{x} (t, f)$ はまずこの自己相関関数 $C_{x}$ を、平均時間 $t = (t_{1} + t_{2}) / 2$ と時間差 $τ = t_{1} - t_{2}$ の関数に直し、時間差 $τ$ についてフーリエ変換を施すことによって得られる。

W_{x} (t, f) = \int_{- \infty}^{\infty} C_{x} (t + \frac{τ}{2}, t - \frac{τ}{2}) e^{- 2 π i τ f} d τ .

よって、 single (mean-zero) time series に対しては、ウィグナー分布は次のように単純に与えられる。

W_{x} (t, f) = \int_{- \infty}^{\infty} x (t + \frac{τ}{2}) x {(t - \frac{τ}{2})}^{*} e^{- 2 π i τ f} d τ .

ウィグナー分布を用いる動機は、定常過程についてはそれが全ての時間 $t$ に対してスペクトル密度関数に帰着し、非定常過程については自己相関関数と完全に一致することである。そのため、ウィグナー分布により、スペクトル密度が時間の経過につれてどのように変化するかを（おおよそ）知ることができる。

時間周波数解析の例

ここでは、ウィグナー分布が時間周波数解析においてどのように用いられるかについての例を挙げる。

定常入力信号

入力信号が定数の場合、その時間周波数分布は時間軸に沿った水平線となる。たとえば、 x(t) = 1 の場合、以下のようになる。

W_{x} (t, f) = \int_{- \infty}^{\infty} e^{- i 2 π τ f} d τ = δ (f) .

正弦波信号

入力信号が正弦波関数の場合、その時間周波数分布は時間軸から正弦波信号の周波数分だけ平行移動した、水平線となる。たとえば、テンプレート:Math の場合、以下のようになる。

\begin{matrix} W_{x} (t, f) & = \int_{- \infty}^{\infty} e^{i 2 π k (t + \frac{τ}{2})} e^{- i 2 π k (t - \frac{τ}{2})} e^{- i 2 π τ f} d τ \\ = \int_{- \infty}^{\infty} e^{- i 2 π τ (f - k)} d τ \\ = δ (f - k) . \end{matrix}

チャープ信号

入力信号がチャープ信号の場合、その瞬間周波数は線形関数になる。つまり、時間周波数分布は直線となる。たとえば、

x (t) = e^{i 2 π k t^{2}}

の場合、瞬間周波数は次のようになり、

\frac{1}{2 π} \frac{d (2 π k t^{2})}{d t} = 2 k t,

そのウィグナー分布は

\begin{matrix} W_{x} (t, f) & = \int_{- \infty}^{\infty} e^{i 2 π k {(t + \frac{τ}{2})}^{2}} e^{- i 2 π k {(t - \frac{τ}{2})}^{2}} e^{- i 2 π τ f} d τ \\ = \int_{- \infty}^{\infty} e^{i 4 π k t τ} e^{- i 2 π τ f} d τ \\ = \int_{- \infty}^{\infty} e^{- i 2 π τ (f - 2 k t)} d τ \\ = δ (f - 2 k t) . \end{matrix}

デルタ関数

入力信号がデルタ関数の場合、 t=0 でのみ非零で、無限の周波数成分を含むため、その時間周波数分布は原点を通り時間軸に垂直な線となる。つまり、デルタ関数の時間周波数分布もまたデルタ関数となる。ウィグナー分布は以下のようになる。

\begin{matrix} W_{x} (t, f) & = \int_{- \infty}^{\infty} δ (t + \frac{τ}{2}) δ (t - \frac{τ}{2}) e^{- i 2 π τ f} d τ \\ = 4 \int_{- \infty}^{\infty} δ (2 t + τ) δ 2 t - τ) e^{- i 2 π τ f} d τ \\ = 4 δ (4 t) e^{i 4 π t f} \\ = δ (t) e^{i 4 π t f} \\ = δ (t) . \end{matrix}

ウィグナー分布は、入力信号の位相が二次以下の場合に最も時間周波数解析に適する。このような信号については、ウィグナー分布は入力信号の時間周波数分布に完全に一致する。

矩形関数

x (t) = {\begin{matrix} 1 & | t | < 1 / 2 \\ 0 & otherwise \end{matrix}

上のような矩形関数の場合、ウィグナー関数は以下のようになる。

W_{x} (t, f) = \frac{1}{π f} \sin (f [1 - 2 | t |])

交差項の性質

ウィグナー分布は線形変換ではない。複数の周波数成分が入力信号に存在する場合、各成分の干渉によるうなりに似た交差項（「時間うなり」）が生じる。もともとのウィグナー関数では、この項は期待値を正確に与えるために必要であり、物理的に重要である。対照的に、短時間フーリエ変換ではこの交差項は生じない。ウィグナー分布における交差項の性質のうちいくつかを以下に示す。

$x (t) = {\begin{matrix} \cos (2 π t) & t \leq - 2 \\ \cos (4 π t) & - 2 < t \leq 2 \\ \cos (3 π t) & t > 2 \end{matrix}$

$x (t) = e^{i t^{3}}$

交差項問題を緩和するため、様々な変換が提案されている。テンプレート:仮リンクやテンプレート:仮リンク、テンプレート:仮リンクなどが挙げられる。