ナイキスト周波数とは?求め方とサンプリング定理を図解でわかりやすく解説

アナログ信号をデジタルに変換するとき、必ず出てくるのが「どのくらいの細かさで測れば、元の波形を失わずに済むのか」という問いです。この問いに完全な答えを与えるのがサンプリング定理(Nyquist-Shannon sampling theorem)であり、その答えの中心にある数値がナイキスト周波数です。

ナイキスト周波数とは:サンプリング周波数の半分

結論から書きます。ナイキスト周波数 $f_N$ とは、サンプリング周波数 $f_s$ のちょうど半分の値で、

$$ f_N = \frac{f_s}{2} $$

で求まります。$f_s = 44100$ Hz なら $f_N = 22050$ Hz。この $f_N$ が「そのサンプリングで正しく記録できる周波数の上限」を意味し、これを超える成分は低い周波数に折り返されて化けます。逆に言えば、信号に含まれる最高周波数が $f_N$ より低ければ、飛び飛びのサンプル点から元の連続波形を完全に復元できる——これがサンプリング定理の主張です。

この定理を押さえておくと、次のような場面で判断に迷わなくなります。

  • A/D変換器を選ぶとき: 測りたい現象の最高周波数から必要なサンプリング周波数が決まる
  • 謎の低周波が観測されたとき: それが本物か、高周波の折り返しかを切り分けられる
  • オーディオの仕様を読むとき: 44.1 kHz や 96 kHz という数値が何から決まっているかが分かる
  • 振動診断・センサ計測: 歯車の噛み合い成分など、見たい高次成分を取りこぼさない設計ができる

本記事の内容

  • ナイキスト周波数の定義と求め方
  • 混同しやすい「ナイキストレート」との違い
  • サンプリング定理のフーリエ変換による証明
  • エイリアシングで何 Hz に化けるかの計算方法
  • 「ちょうど2倍」では足りない理由(実測)
  • sinc 補間による完全復元
  • アンチエイリアシングフィルタと、CD が 44.1 kHz である理由

前提知識

この記事を読む前に、以下の記事を読んでおくと理解が深まります。

未読でも本記事単独で読み進められます。

ナイキスト周波数とナイキストレートは別物

先に用語を整理しておきます。日本語の解説でもよく混同されているので、ここを分けておくと以降が一気に読みやすくなります。

ナイキスト周波数とナイキストレートの違い

ナイキスト周波数 $f_N = f_s/2$ は、サンプラー側で決まる量です。「この装置で測れる上限はいくらか」を表します。一方のナイキストレート $2f_{\max}$ は、信号側で決まる量です。「この信号を測るには最低いくらのサンプリング周波数が必要か」を表します。

どちらも単位は Hz なので紛らわしいのですが、主語が違うと覚えてください。そしてサンプリング定理は、この2つを比べているだけです。

$$ f_s > 2f_{\max} \quad \Longleftrightarrow \quad f_N > f_{\max} $$

左は「サンプリング周波数がナイキストレートを超えている」、右は「ナイキスト周波数が信号の最高周波数を超えている」。同じことを言っています。

具体例で確認しましょう。20 kHz までの音を録りたいとします。このときナイキストレートは $2 \times 20 = 40$ kHz。実際に $f_s = 44.1$ kHz を選べばナイキスト周波数は 22.05 kHz となり、20 kHz より高いので条件を満たします。

用語が整理できたところで、なぜこの条件で完全復元できるのかを見ていきます。

サンプリング定理の証明

ステップ1:サンプリングの数学的表現

連続信号 $x(t)$ をサンプリング周期 $T_s$ で標本化することは、デルタ関数の櫛(comb function)を掛けることに等しいです。

$$ x_s(t) = x(t) \cdot \sum_{n=-\infty}^{\infty} \delta(t – nT_s) = \sum_{n=-\infty}^{\infty} x(nT_s) \, \delta(t – nT_s) $$

デルタ関数を掛けるという操作は、「その瞬間の値だけを抜き出して、他は捨てる」ことの数学的な表現です。

ステップ2:サンプリング後のスペクトル

デルタ櫛のフーリエ変換は、周波数領域でも周期 $f_s = 1/T_s$ のデルタ櫛になります。

$$ \sum_{n=-\infty}^{\infty} \delta(t – nT_s) \quad \longleftrightarrow \quad f_s \sum_{k=-\infty}^{\infty} \delta(f – kf_s) $$

時間領域の乗算は周波数領域の畳み込みに対応するので、

$$ X_s(f) = X(f) * f_s \sum_{k=-\infty}^{\infty} \delta(f – kf_s) = f_s \sum_{k=-\infty}^{\infty} X(f – kf_s) $$

デルタ関数との畳み込みは「その位置に平行移動したコピーを置く」操作なので、この式は元のスペクトルが $f_s$ 間隔で無限に複製されることを意味します。

サンプリングはスペクトルをfs間隔で複製する操作

上段が元の信号のスペクトル($f_{\max} = 8$ Hz で帯域制限された三角形)、下段がサンプリング後です。まったく同じ形が $f_s = 40$ Hz 間隔でずらりと並んでいます。サンプリングという操作の正体は、周波数領域では「コピーを等間隔にばらまくこと」なのです。

そして中央のコピー($k = 0$)だけを切り出せれば、元のスペクトルがそのまま手に入ります。切り出せる条件は明白で、隣のコピーと重なっていないことです。

ステップ3:復元の条件

$X(f)$ の帯域が $[-f_{\max}, f_{\max}]$ のとき、$k = 0$ のコピーは $f_{\max}$ まで、$k = 1$ のコピーは $f_s – f_{\max}$ から始まります。重ならない条件は

$$ f_s – f_{\max} > f_{\max} \quad \Longrightarrow \quad \boxed{f_s > 2f_{\max}} $$

これでサンプリング定理が導かれました。条件を満たしていれば、カットオフ $f_s/2$、ゲイン $T_s$ の理想低域通過フィルタを掛けるだけで $X(f)$ が取り出せます。

証明はこれだけです。「コピーが重なるか重ならないか」という一点に集約されているところが、この定理の美しさでもあります。では、重なってしまうとどうなるのでしょうか。

エイリアシング — 重なると何が起きるか

サンプリング周波数が足りないとコピーが重なる

上段は $f_s = 20 > 2f_{\max} = 16$ で条件を満たしており、コピーどうしに隙間があります。下段は $f_s = 12 < 16$ で、隣のコピーが $k=0$ の領域に食い込んでいます(赤で塗った部分)。

いったん足し合わされてしまうと、どちらのコピー由来の成分なのか区別する手段はありません。これがエイリアシング(aliasing、日本語では折り返し歪み)です。ノイズのように後から取り除けるものではなく、情報そのものが失われています。

何 Hz に化けるのか

実務でよく必要になるのが「$f_N$ を超えた信号は結局いくつに見えるのか」という計算です。答えは

$$ f_{\text{観測}} = \left| f – k f_s \right| \quad (\text{結果が } 0 \leq f_{\text{観測}} \leq f_N \text{ に収まるような整数 } k) $$

で与えられます。図にすると三角波の形になります。

折り返しの地図:fNを超えると三角形に折り返す

$f_s = 8$ Hz の場合、$f_N = 4$ Hz までは入力周波数がそのまま観測されます(緑の直線)。そこを超えると折り返しが始まり、5 Hz は 3 Hz に、11 Hz も 3 Hz に、13 Hz も 3 Hz に見えます。まったく違う周波数が同じ値に化けるわけで、観測されたデータだけを見て元がどれだったかを言い当てることは原理的に不可能です。

時間領域で見ると、この「区別できなさ」がさらに実感できます。

5Hzと3Hzがサンプル点では完全に一致する

本物の 5 Hz(水色)と、見かけの 3 Hz(赤の破線)は、まったく違う波形です。ところが $f_s = 8$ Hz のサンプル点(オレンジ)では両者の値が完全に一致します。実測すると差は $1.2 \times 10^{-15}$、つまり倍精度の丸め誤差そのものでした。サンプリングした後にデータをいくら眺めても、この2つを見分ける情報は残っていません。

だからこそ、サンプリングする前に高周波を除去しておく必要があります。それが後述のアンチエイリアシングフィルタです。

身のまわりのエイリアシング

この現象は信号処理の中だけの話ではありません。日常でも同じことが起きています。

車のホイールが逆回転して見える。映画やテレビは毎秒24〜30コマの「サンプリング」をしています。ホイールのスポークがコマ間隔にちょうど近い周期で回ると、実際より遅く見えたり、止まって見えたり、逆回転して見えたりします。$f_s = 24$ fps に対してスポークの通過周波数が 25 Hz なら、$|25 – 24| = 1$ Hz——つまり「毎秒1回ゆっくり前進」して見えるわけです。26 Hz なら $|26 – 24| = 2$ Hz で少し速く、逆に 23 Hz なら折り返して「毎秒1回の逆回転」に見えます。

モアレ縞。細かいストライプ模様の服をカメラで撮ると、実物には無いうねった縞が現れます。これは空間方向のサンプリング(画素の格子)に対して、模様の空間周波数が高すぎるために起きる2次元のエイリアシングです。デジタルカメラに光学ローパスフィルタが入っていたのは、まさにこれを防ぐためでした。

扇風機やヘリコプターのローター。蛍光灯のちらつき(商用電源の2倍の周波数)が実質的なサンプリングとして働き、羽根が止まって見えることがあります。

どれも「速すぎる動きを、遅い間隔で覗いた」結果であり、原理はまったく同じです。エイリアシングは信号処理特有の落とし穴ではなく、離散的に観測することの本質的な帰結だと捉えると腑に落ちます。

「ちょうど2倍」では足りない

サンプリング定理の条件は $f_s > 2f_{\max}$ であって、$f_s \geq 2f_{\max}$ ではありません。等号が入らないこの違いは、単なる数学的な厳密さの話ではなく、実際に破綻します。

ちょうど2倍だと位相しだいで全サンプルがゼロになる

5 Hz の正弦波を $f_s = 10$ Hz(ちょうど2倍)でサンプリングした結果です。位相が 90° なら振幅 1.000 がきちんと取れます。ところが位相が 0° だと、すべてのサンプル点がゼロ交差に当たってしまい、振幅 0.000——信号が存在しないのと区別がつきません。位相 45° では 0.707 と、本来の振幅より小さく測れます。

つまり「ちょうど2倍」では、測れるかどうかが位相という制御できない要素に左右されるのです。実務で $f_s = 2f_{\max}$ ぴったりを選ぶことはなく、必ず余裕を取ります。

余裕をどれだけ取るべきかは、復元誤差を測ってみると見当がつきます。

サンプリング周波数と復元誤差の関係

5 Hz の正弦波を各サンプリング周波数で標本化し、sinc 補間で復元したときの RMSE です。$f_s/f_0 = 1.6$ では 1.065(波形が別物)、ちょうど 2.0 でも 0.707 とまだ大きく、2.2 で 0.145、4.0 で 0.014 と急速に落ちていきます。2倍を境に誤差が桁で改善するが、ちょうど2倍ではまだ足りないという傾向がはっきり出ています。

なお $f_s$ を上げても RMSE がゼロに漸近しないのは、有限個のサンプルで sinc の無限和を打ち切っているためです。理論上は無限個のサンプルがあれば完全復元されます。

sinc 補間による復元

条件を満たしていれば、離散サンプルから元の連続信号を次の式で復元できます。

$$ \boxed{x(t) = \sum_{n=-\infty}^{\infty} x(nT_s) \, \mathrm{sinc}\!\left(\frac{t – nT_s}{T_s}\right)}, \qquad \mathrm{sinc}(u) = \frac{\sin(\pi u)}{\pi u} $$

導出

理想低域通過フィルタの伝達関数は

$$ H(f) = \begin{cases} T_s & (|f| \leq f_s/2) \\ 0 & (|f| > f_s/2) \end{cases} $$

これを逆フーリエ変換したものがインパルス応答 $h(t)$ です。

$$ h(t) = \int_{-f_s/2}^{f_s/2} T_s \, e^{j2\pi ft} \, df = T_s f_s \, \mathrm{sinc}(f_s t) = \mathrm{sinc}\!\left(\frac{t}{T_s}\right) $$

サンプリング信号にこれを畳み込めば、上の復元式が得られます。

$$ x(t) = x_s(t) * h(t) = \sum_{n=-\infty}^{\infty} x(nT_s) \, \mathrm{sinc}\!\left(\frac{t – nT_s}{T_s}\right) $$

なぜこれで「完全」なのか

sinc 関数は $u = 0$ で 1、$u$ が 0 以外の整数で 0 になります。この性質から、サンプル点上では

$$ x(mT_s) = \sum_{n} x(nT_s) \, \mathrm{sinc}(m – n) = x(mT_s) $$

と、元の値がそのまま再現されます。サンプル点の間は、他のすべてのサンプルからの寄与が足し合わされて埋まります。

sinc補間による復元

上段が各サンプル点に立てた sinc 関数、下段がそれを全部足した結果です。元の波形(水色)と復元(赤の破線)がほぼ完全に重なりました(中央部の RMSE 0.0276)。

ここで理解しておきたいのは、sinc 補間は「点をなめらかにつなぐ」処理ではないということです。直線補間やスプライン補間は「それらしく見える曲線」を引くだけですが、sinc 補間は「そのサンプル点を通り、かつ $f_N$ 以下の成分しか含まない信号」——そういう信号は数学的にただ1つしか存在しません——を復元しています。だから「完全」と言えるのです。

実務での決め方

アンチエイリアシングフィルタ

エイリアシングは後から取り除けないので、A/D変換の手前にアナログの低域通過フィルタを入れて、$f_N$ より上の成分をあらかじめ落としておきます。これがアンチエイリアシングフィルタです。

CDが44.1kHzである理由:遷移帯域の余裕

ここで問題になるのが、現実のフィルタは垂直に落とせないということです。20 kHz までは通し、20 kHz を超えたら即座にゼロ、という理想フィルタは作れません。必ず「だんだん落ちていく」遷移帯域が要ります。

CD の 44.1 kHz は、まさにこの事情から決まっています。可聴上限 20 kHz に対してナイキストレートは 40 kHz ですが、それでは $f_N = 20$ kHz でフィルタの遷移帯域がゼロになってしまいます。$f_s = 44.1$ kHz にすれば $f_N = 22.05$ kHz となり、20 kHz から 22.05 kHz までの 2.05 kHz が遷移帯域として使えるわけです。

実用値の目安

用途別のサンプリング周波数とナイキスト周波数

分野ごとに慣例となっている値を並べました。どれも考え方は同じで、「測りたい最高周波数の2倍」を下限にして、フィルタの遷移帯域ぶんを上乗せするという手順で決まっています。

電話の 8 kHz は極端に低く見えますが、音声の明瞭度に必要な 3.4 kHz が $f_N = 4$ kHz に収まるので用途としては十分です。逆に振動診断では、歯車の噛み合い周波数など高次成分を見たいので「回転数の10倍以上」といった経験則が使われます。

意図的にエイリアシングを使う場合

最後に例外にも触れておきます。エイリアシングは常に敵というわけではなく、アンダーサンプリング(帯域通過サンプリング)として積極的に利用されることがあります。

たとえば 100 MHz を中心とする幅 1 MHz の狭帯域信号を扱うとき、素直にやれば 200 MHz 超のサンプリングが必要です。しかし信号が $[99.5, 100.5]$ MHz という狭い帯域にしか存在しないと分かっているなら、うまく $f_s$ を選ぶことで折り返し後に低い周波数へ降ろしてから低速な A/D で受けられます。ソフトウェア無線の受信機などで使われる手法です。

これが成立するのは「信号がどの帯域にあるかを事前に知っている」場合に限られます。サンプリング定理の条件 $f_s > 2f_{\max}$ は「$0$ から $f_{\max}$ まで何が入っているか分からない」という前提での要求であり、帯域の位置が既知ならもっと少ないサンプル数で足りる、というだけの話です。

よくある質問

検索されやすい疑問を、簡潔にまとめておきます。

Q. ナイキスト周波数の求め方は? A. サンプリング周波数を 2 で割るだけです。$f_N = f_s/2$。$f_s = 1000$ Hz なら $f_N = 500$ Hz、$f_s = 48$ kHz なら $f_N = 24$ kHz。信号の中身とは無関係に、サンプリング周波数だけで決まります。

Q. サンプリング周波数はナイキストレートの何倍にすべき? A. 理論上の下限は「わずかでも超えていればよい」ですが、実務では 2.2〜2.5倍が一つの目安です。アンチエイリアシングフィルタの遷移帯域が必要なためで、CD の 44.1 kHz は 20 kHz に対して 2.2 倍にあたります。波形そのものを目で見て確認したい用途(オシロスコープ的な使い方)なら 10 倍以上取ることも珍しくありません。

Q. ナイキスト周波数を超える成分が入っていたらどうなる? A. その成分は消えるのではなく、$|f – kf_s|$ の位置に移動して重なります。元からそこにあった正しい成分と混ざってしまうため、後処理では分離できません。サンプリング前にアナログフィルタで落とすのが唯一の対策です。

Q. サンプリング周波数を上げれば音質は良くなる? A. サンプリング周波数が決めるのは「記録できる周波数の上限」だけです。人の可聴上限が 20 kHz なので、44.1 kHz あれば聴こえる範囲は原理的にすべて記録できています。96 kHz や 192 kHz のハイレゾに意味があるとすれば、それはフィルタ設計に余裕ができて位相特性が改善するといった間接的な理由によるもので、「44.1 kHz では高音が欠けている」という説明は正確ではありません。

Q. ビット深度(16bit / 24bit)とはどう違う? A. 別の軸の話です。サンプリング周波数は時間方向の細かさ、ビット深度は振幅方向の細かさを決めます。サンプリング定理が扱うのは前者だけで、後者は量子化誤差(ノイズフロア)の大きさを左右します。16 bit で約 96 dB、24 bit で約 144 dB のダイナミックレンジに相当します。

まとめ

本記事では、ナイキスト周波数とサンプリング定理について解説しました。

  • ナイキスト周波数 $f_N = f_s/2$ は、サンプリング周波数の半分。正しく記録できる周波数の上限を表す
  • ナイキストレート $2f_{\max}$ は信号側の量で、$f_N$ とは主語が違う。サンプリング定理は $f_N > f_{\max}$ を要求している
  • 証明の核心は「サンプリングとはスペクトルを $f_s$ 間隔で複製する操作」であり、コピーが重ならない条件が $f_s > 2f_{\max}$
  • エイリアシングでは $|f – kf_s|$ に折り返される。5 Hz と 3 Hz がサンプル点で $10^{-15}$ の精度で一致するように、後から見分けることはできない
  • 「ちょうど2倍」では不十分。位相 0° では全サンプルがゼロになり、振幅が測れない
  • sinc 補間は「点をつなぐ」のではなく「唯一の帯域制限信号を復元する」操作
  • 実務ではアンチエイリアシングフィルタの遷移帯域ぶんの余裕が要る。CD の 44.1 kHz はここから決まっている

次のステップとして、以下の記事も参考にしてください。