スペクトル密度

**スペクトル密度**（スペクトルみつど、）とは、信号を構成する成分を周波数の関数として表したものであり、単にスペクトルと言うこともある。成分として信号の全エネルギーについて表したものを**エネルギースペクトル密度**（、ESD）という。また、仕事率（パワー）について表したものを**パワースペクトル密度**（）といい、仕事率として電力が対象となる場合、**電力スペクトル密度**とも呼ばれる。スペクトル密度は定常過程に関する周波数値の正実数の関数または時間に関する決定的な関数である。直観的には、確率過程の周波数要素を捉えるもので、周期性を識別するのを助ける。

## 概要
信号のエネルギーは振幅の二乗和でしばしば定義される。信号を定常波の和すなわちスペクトルとして見たとき（フーリエ変換）、信号全体のエネルギーは部分定常波エネルギーの総和になると考えられる。より正確には、連続値である各周波数にエネルギー密度が定義出来てその積分値が信号全体のエネルギーになると考えられる（パーセバルの定理）。各周波数におけるエネルギー密度を**エネルギースペクトル密度**という。

また、信号の仕事率（パワー）は時間当たりのエネルギーでしばしば定義される。全く同じ議論がパワーに関してもでき、各周波数におけるパワー密度を**パワースペクトル密度**という。

物理学の観点では、信号とは波動であり、代表的な波動には電磁波や音波がある。信号がどのような物理的次元を伝わるのかは問題ではないが、以下の議論では時間と共に変化する信号について解説する。次元解析の観点では、パワースペクトル密度の単位はワット毎ヘルツ (W/Hz) か、ワット毎ナノメートル (W/nm) で表される（後者は周波数の代わりに波長を用いる）。

## 定義

### エネルギースペクトル密度

#### 連続信号
連続信号 *f*(*t*) の**エネルギースペクトル密度**は次の式で定義される。
\int_{-\infty}^\infty f(t)e^{-i\omega t}\,dt\right|^2 = \frac{F(\omega)F^*(\omega)}{2\pi}}}
*ω は*角周波数*、F*(*ω*) は *f*(*t*) の連続フーリエ変換、*F**(*ω*) はその複素共役である。1/2\pi  という係数は絶対的なものではなく、フーリエ変換での正規化定数の定義に依存する。*f*(*t*) が有限エネルギー信号であるとき、その信号のスペクトル密度 *ESD*(*ω*) は、信号をフーリエ変換したときの大きさの2乗である。

すなわちESDは信号のエネルギーが周波数についてどのように分布するかを示す。

#### 離散信号
離散信号 *f<sub>n</sub>* = *f*(*n*,*dt*) が無限に続くとするなら**エネルギースペクトル密度**は次の式で定義される。
\sum_{n=-\infty}^\infty f_n e^{-i\omega n}\right|^2=\frac{dt^2}{2\pi}F_d(\omega)F_d^*(\omega)}}
ここで、*F*(*ω*) は *f<sub>n</sub>* の離散時間フーリエ変換である。数学ではサンプリング間隔 *dt* を 1 として扱うことが多い。しかしながら、正確な物理単位を維持するためと、*dt* → 0 とした場合に連続時間の関数へ逆変換できることを保証するためには *dt* が必要となる。

#### 次元解析
ここで、エネルギーは信号の2乗を積分したものであり、その信号を電圧として 1オームの負荷に加えたときの物理エネルギーに等しい。*f*(*t*) が伝送路を通って伝播する電気信号の電位を（ボルトで）表す場合、スペクトル密度 *ESD*(*ω*) の測定単位は volt2×seconds2 として現れるが、物理学のエネルギースペクトル密度としてはまだ次元的に正確ではない。しかしながら、（オームで表される）伝送路の特性インピーダンス *Z* によって除算すると、*ESD*(*ω*) の次元は1オーム当たり volt2×seconds2 になる。これは、物理学で定義されるエネルギースペクトル密度の国際単位であるジュール毎ヘルツ（J/Hz）と等価となる。

### パワースペクトル密度
上述のエネルギースペクトル密度の定義は、信号のフーリエ変換が存在するパルスのような信号に最も適している。たとえば定常物理過程を示す連続信号について、**パワースペクトル密度**あるいは**電力スペクトル密度** (PSD) を定義することは価値があり、信号や時系列のパワーが周波数についてどのように分布しているかを示す。抽象的な信号についても、信号の2乗と定義できる。このとき、信号 *f*(*t*) のある一瞬の力は次のように与えられる。
:P(t) = f(t)^2

平均（あるいは期待値）としての *P*(*t*) は、全周波数領域にわたるパワースペクトル密度の積分である。

正規化されたフーリエ変換:
: \mathcal{F}_T(\omega) = \frac{1}{\sqrt{T}} \int_0^T f(t) \exp(-i\omega t) dt
を使用して、次のようにパワースペクトル密度を定義できる。
: PSD(\omega) = \lim_{T \rightarrow \infty} \mathbf{E} \left[ | \mathcal{F}_T(\omega) | ^ 2 \right]
確率論的な信号については、フーリエ変換の二乗値は一般的に極限に近づけないが、期待は行う。（を参照。）

**見解**：取り扱う多くの信号が積分可能ではなく、その信号の *非正規化（＝通常の）* フーリエ変換は存在しない。何人かの著者（たとえば Risken）は、まだ非正規化フーリエ変換を使ってパワースペクトル密度の定義
: \langle F(\omega) F^\ast(\omega') \rangle = 2\pi\,PSD(\omega)\,\delta(\omega-\omega')
を公式化している。ここで、δ(*ω* &minus; *ω*') はディラックのデルタ関数である。このような公式の文献は直観を導くには有用であるが、十分な注意と共に使用されるべきである。

このような形式推論を用いると、定常ランダム過程とパワースペクトル密度 *PSD*(*ω*) およびこの信号の自己相関関数 *R*(τ) =  がフーリエ変換対でなければならないことに気づくだろう。このことは真実であり、ノーバート・ウィーナーおよびアレクサンドル・ヒンチンによって作り出された意味深い定理（ウィーナー・ヒンチンの定理）となる。
:PSD(f)=\int_{-\infty}^\infty \,R(\tau)\,e^{-i\omega\tau}\,d \tau=\mathcal{F}(R(\tau))
多くの著者が、実際にパワースペクトル密度を *定義する* ためにこの等式を使用している。そうする理由は「数学的曖昧さ」を回避するためであると、多くの書籍に記載されている。

ある周波数帯域 [ω1, ω2] における信号の力は、正の周波数と負の周波数について積分することで計算できる。
:P=\int_{\omega_1}^{\omega_2}\,PSD(\omega)+PSD(-\omega) \,d \omega
信号のパワースペクトル密度は、その信号が広義の定常過程であるときだけ存在する。信号が広義、もしくは、狭義の定常過程でない場合、その自己相関関数は2つの変数の関数となる。広義の  のような場合、PSD は存在する可能性がある。
より一般に、似たような技法で時と共に変化するスペクトル密度の近似を求めることができる。

パワースペクトル密度の定義は、全測定時間 *T* =*n dt* の間に離散時間 *f<sub>n</sub>* = *f*(*n*,*dt*) でサンプリングされた信号のような有限の時系列 *f<sub>n</sub>* = *f*(*n*,*dt*)（ただし 1 ≤ *n* ≤ *N*）を直接的に一般化する。
:PSD(\omega)=\frac{dt^2}{T}\left|\sum_{n=1}^N f_n e^{-i\omega n}\right|^2.
実世界の応用では、観察された物理過程の基礎となる実際の PSD のより正確な推定を行うために、一度の測定で得られる PSD の結果を複数回反復測定し平均化することが一般的である。このように計算された PSD はと呼ばれる。平均する時間間隔Tを無限に近づける場合、ピリオドグラムが真のパワースペクトル密度 (PSD) に近づくことを証明できる（ブラウンとホワン）。

2つの信号共に ***パワースペクトラ*** （正確な用語；spectrum の複数形は spectra）を有する場合、これらの相互相関関数を用いてクロスパワースペクトルを計算できる。

#### パワースペクトル密度の特性
PSD には次のような特性がある。
*実際に使われる過程のスペクトルは対称である: *S*(&minus; *f*) = *S*(*f*) 言い換えると、偶関数である。
*[&minus; 1/2, +1/2] の範囲で連続しており、微分可能である。
*PSD の微分は *f* = 0 で 0 となる。（このことはパワースペクトルが偶関数となるために必要である。）そうでない場合、微分は *f* = 0 で存在しない可能性がある。
*自己共分散関数はフーリエ逆変換を使うことにより再構成することができる。
*PSD は、時間軸上の分散の分布を示している。とりわけ、
*:\text{Var}(X_t) = \gamma_0 = 2 \int_0^{1/2} S(\omega) d\omega である。
*PSD は自己共分散関数の一次関数となる。
*:もし γ が2つの関数 γ(τ) = α1γ1(τ) + α1γ2(τ) に再構成される場合、
*:*S*(*f*) = α1*S*1(*f*) + α2*S*2(*f*) となる。
*:: ここで S_i(f) = \mathcal{F}\{ \gamma_i \}
**パワースペクトル** *G*(*f*) は次式で定義される。
:G(f)=  \int _{-\infty}^f S(f^\prime) \, df^\prime. 

## 推定
スペクトル密度推定の目的は、連続した時間サンプルからランダム信号のスペクトル密度を推定（estimate）することである。
信号から何が知られているかに依存するが、推定方法は  と非パラメトリック推定の2つの方法があり、時間領域または周波数領域の分析が基本となる。たとえば、 で共通の技術は自己回帰モデルに観測を適応させることを含んでいる。非パラメトリック推定で共通の技術はである。

スペクトル密度は通常フーリエ変換法を使用して推定されるが、や最大エントロピー法といった他の技術も使用することができる。

## 特性
**f*(*t*) のスペクトル密度と *f*(*t*) の自己相関は、フーリエ変換対を形成する（PSD と ESD とで、自己相関関数の異なる定義が使われる）。
*フーリエ解析の1つの結果としてパーセバルの定理がある。それによると、エネルギースペクトル密度\Phi(\omega)の曲線の面積は、信号の振幅f(t)の自乗すなわち全エネルギーの面積に等しい。

この定理は離散的な場合でも成り立つ。同様にパワースペクトル密度の積分したものは、それに対応する信号の全エネルギーの平均に等しい（それはまた、遅延ゼロでの自己相関関数である）。

## 関連する概念
*周波数分布を示すグラフは、ほとんどの場合スペクトル密度を表している。完全な周波数スペクトルを描く場合、振幅と周波数のグラフ（スペクトル密度に相当）と位相と周波数のグラフ（スペクトル密度以外の情報）で表される。信号 *f*(*t*) の波形は、完全な周波数スペクトルがあれば再現できる。信号 *f*(*t*) をスペクトル密度情報だけから再現することはできない。
*スペクトル密度関数の中点を、その信号のスペクトル重心と呼ぶ。すなわち、その周波数を分割点として、上と下でエネルギーが拮抗する。
*スペクトル密度は周波数の関数であって、時間の関数ではない。しかし、長い信号の非常に短い期間のスペクトル密度を計算することもでき、それらを時系列に並べることもできる。そのようなグラフをスペクトログラムと呼ぶ。これは、短時間フーリエ変換やウェーブレット変換などのスペクトル解析技法の基本である。
*スペクトル密度を信号とみなし、フーリエ変換して得られる信号をケプストラムと呼ぶ。すなわち、スペクトルのスペクトルである。

## 応用
### 電子工学
信号のパワースペクトル密度は電子工学の基本概念の1つであり、特に電子通信システム（無線、レーダーなど）で重要である。電気信号のパワースペクトルを測定して表示する機器としてスペクトラムアナライザがある。

スペクトラムアナライザは、入力信号の短時間フーリエ変換 (STFT) の絶対値を測るのが基本である。解析対象の信号が定常的ならば、STFT はパワースペクトル密度のよい近似となる。

### 測色法
right
光のスペクトルとは、色に対応した各周波数で運ばれる力を示したものである。光スペクトルは周波数よりも波長で表されることが多く、厳密にはスペクトル密度ではない。分光測色器によっては、1から2ナノメートル単位の分解能を持つ。値は他の用途に使われたり、光源のスペクトル属性を示すために図示されたりする。これを使って光源の色特性を解析する。

## 関連項目
*
*ノイズ
*スペクトル漏れ
*窓関数
*周波数領域
*周波数スペクトル
*
*スペクトル解析
*スペクトログラム
*スペクトラムアナライザ

## 脚注

## 外部リンク
*時系列データ解析におけるパワースペクトル密度関数について Cygnus Research International
*スペクトル解析の基礎知識