標本化定理

提供: miniwiki
2018/8/19/ (日) 17:27時点におけるAdmin (トーク | 投稿記録)による版 (1版 をインポートしました)
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
移動先:案内検索

標本化定理(ひょうほんかていり、: sampling theorem: サンプリング定理とも)はアナログ信号デジタル信号へと変換する際に、どの程度の間隔で標本化(サンプリング)すればよいかを定量的に示す定理情報理論の分野において非常に重要な定理の一つである。

標本化定理は1928年にハリー・ナイキストによって予想され、1949年にクロード・E・シャノンと日本の染谷勲によってそれぞれ独立に証明された。そのためナイキスト定理ナイキスト・シャノンの定理シャノン・染谷の定理[1]とも呼ばれる。

概要

波形の最大周波数の2倍以上の周波数で標本化すれば完全再構成

標本化とは、数学的には連続関数の値からある点の値だけを標本として取り出して離散関数に変換する操作であり、与えられた連続関数と標本化関数の積を求めることと等しい。標本化関数とは、ある離散値(連続でない、飛び飛びの値)xに対してのみg(x)=1となり、その他のxに対してはg(x)=0となるような関数である。対象となる関数f(x)と標本化関数g(x)の積を取ると、関数[math]h(x)=g(x)f(x)[/math]が得られる。g(x)=1となるxに対してのみ[math]h(x)=f(x)[/math]となり、その他のxに対してはh(x)=0となる。

標本化定理とは、ある関数f(x)フーリエ変換した関数F(s)の成分(スペクトル)が、[math]|s| \ge W[/math]の範囲でF(s)=0であるような関数f(x)に対して、[math]s=2W[/math]に相当する周期より小さい周期をもつ標本化関数で標本化したときに得られる関数は、そのスペクトルのうち[math]|s|\lt W[/math]が原関数のスペクトルに一致するというものである。

工学的には、原信号に含まれる最大周波数成分を f とすると、2f よりも高い周波数 [math]f_s[/math] で標本化した信号は、低域通過(ローパス)フィルターで高域成分を除去することによって原信号を完全に復元することができるということを示している。たとえば原信号に含まれる周波数が最高で f=22.05kHz だった場合、[math]f_s[/math]=44.1kHz よりも高い[2]周波数で標本化(1秒間に44100回超、値を取得)すれば、原信号を完全に復元することができる。原信号が復元可能な最大周波数 [math]{f_s}/2[/math] を「ナイキスト周波数」と言い、ナイキスト周波数の逆数を「ナイキスト周期」と言う。

標本化周波数が 2f 以下であった場合、原信号にはない偽の周波数 [math]f_s - f[/math]エイリアス信号として、復元信号に現れる。よって、連続信号の標本化においては、ナイキスト周波数 2f よりも高い周波数で、標本化を行わなくてはならない。

なお、アナログ信号からデジタル信号への変換については、標本化のほかに量子化が必要である。

標本化定理の証明

標本化定理は、フーリエ級数を用いると簡単に証明することができる。

理想的な標本化パルス列s(t)は、Tをサンプリング周期とし、デルタ関数[math]\delta(t)[/math]を用いて、

[math] s(t) = \sum_{n=-\infty}^{\infty} \delta(t-nT) [/math]

と表される。標本化入力信号をg(t)とすると、出力信号p(t)

[math] p(t) = g(t)s(t) [/math]

であるから、

[math] p(t) = g(t)\sum_{n=-\infty}^{\infty} \delta(t-nT)=\sum_{n=-\infty}^{\infty} g(nT)\delta(t-nT) [/math]

となり、明らかにg(nT)の系列となる。

ここで、出力信号p(t)の周波数成分を計算するためにs(t)をフーリエ級数展開すると、

[math] s(t)=\frac{1}{T}\sum_{n=-\infty}^{\infty} e^{jn\omega_0 t} [/math]

となる。ただし、[math] \omega_0 = 2\pi f_0 = \frac{2\pi}{T} [/math]である。

扱いを容易にするために入力信号g(t)は振幅A、周波数[math]f_a = \frac{\omega_a}{2\pi} [/math]の単一正弦波として次のように置く。

[math] g(t) = A\cos(\omega_a t + \theta_a)=\frac{A}{2}e^{j(\omega_a t + \theta_a)} + \frac{A}{2}e^{-j(\omega_a t + \theta_a)} [/math]

これに対する出力信号p(t)は、上の式より

[math] p(t) = \frac{A}{2T}\sum_{n=-\infty}^{\infty} e^{j\{(n\omega_0+\omega_a)t+\theta_a\}} + \frac{A}{2T}\sum_{n=-\infty}^{\infty} e^{j\{(n\omega_0-\omega_a)t-\theta_a\}} [/math]

となる。この式から周波数スペクトルの図を描き検討すると証明ができる。

抵抗と電圧のゆらぎについてのナイキストの定理

抵抗[math]R[/math]と電圧のゆらぎとの比例関係。導体が温度[math]T[/math]にあるとき、その両端には電位差[math]V(t)[/math]が生じる。このとき

[math]\langle V(t)V(t') \rangle = 2Rkt\delta(t-t')[/math]

の関係をナイキストの定理という。この関係式は、角振動数[math]\omega[/math]に対する電気伝導度[math]\sigma(\omega)[/math][math]\omega[/math]によらず[math]\sigma(0)[/math]に等しい領域で成立する。これは一般の線形応答理論から基礎づけられる。これも歴史的には1つの揺動散逸定理の発見の例になっている[3]

脚注

  1. ただし、実際にこの呼び方を日本国外はもとより、日本国内でも、論文や文献などで目にする機会はほぼ皆無である。テブナンの定理を鳳-テブナンの定理とごく僅かに呼ばれるのと同様である。
  2. (注)“以上” ではない。
  3. 『物理学辞典』 培風館、1984年

関連項目

テンプレート:データ圧縮 テンプレート:デジタル信号処理