標本化定理とエイリアシング
音声も画像もなぜ折り返し雑音で崩れるのか。標本化定理とナイキスト周波数、アンチエイリアスフィルタの要否を原理から押さえ、標本化レートを正しく選べるようになります。
- 1.帯域が B までの信号は、標本化周波数 fs が 2B を超えれば標本列から完全に復元できる(標本化定理)。fs/2 をナイキスト周波数と呼び、これが表現できる周波数の上限。
- 2.fs/2 を超える成分は、標本化でスペクトルが fs 間隔に複製される際に低域へ写り込み、別名(エイリアス)の偽周波数として現れる。一度折り返した成分は事後処理では分離できない。
- 3.対策は標本化の前段でアナログのアンチエイリアスフィルタをかけ、fs/2 以上を落とすこと。再構成は sinc 補間(理想ローパス)が理論解で、実機はゼロ次ホールド+補正で近似する。
連続信号を「点」で捉える
アナログの連続時間信号を計算機で扱うには、一定間隔で値を拾う標本化(サンプリング)が要ります。素朴には「細かく取れば元に戻せそう」ですが、では「どれだけ細かければ十分か」に厳密な答えを与えるのが標本化定理です。ここでは、標本化がスペクトルに何をするか・ナイキスト周波数の意味・折り返し雑音(エイリアシング)の発生機構・アンチエイリアスフィルタと再構成までを原理で追います。記号は LaTeX を使わず、プレーンな式で示します。
連続信号 x(t) を間隔 Ts ごとに拾った標本列を考えます。標本化周波数(サンプリングレート)は fs = 1/Ts です。数学的には、間隔 Ts のインパルス列(ディラックの櫛)を掛ける操作としてモデル化します。
x_s(t) = x(t) · Σ_k δ(t − k·Ts) k = …,−1,0,1,…
標本値 x[k] = x(k·Ts)
時間軸で「櫛を掛ける」この操作が、周波数軸で何を引き起こすかが核心です。
標本化はスペクトルを複製する
時間領域の積は周波数領域では畳み込みになります。インパルス列のフーリエ変換は、間隔 fs のインパルス列(やはり櫛)です。したがって標本化後のスペクトル Xs(f) は、元のスペクトル X(f) を fs おきに無限に並べた複製の重ね合わせになります。
Xs(f) = fs · Σ_n X(f − n·fs) n = …,−1,0,1,…
つまり標本化とは、元のスペクトルを fs 間隔でコピー&ペーストして足し合わせる操作に等しい、というのが周波数領域から見た正体です。元の帯域(ベースバンド)のコピーが ±fs, ±2fs, … の位置に「イメージ」として現れます。
連続信号のスペクトルは1つの塊ですが、標本列は「Ts ごとの値」しか持たないため、周波数軸で fs だけずれた成分どうしを区別できません。ある周波数 f の正弦波と、f + fs の正弦波は、標本点上でまったく同じ値を取ります。この見分けのつかなさが、スペクトル上では fs 周期の複製として現れます。
ナイキスト周波数と折り返し
複製されたスペクトルどうしが重ならなければ、ベースバンドだけをローパスフィルタで切り出して元に戻せます。重ならない条件は、元信号が帯域制限されていて最高周波数を B とするとき、隣の複製の下端 fs − B がベースバンド上端 B を上回ること、すなわち次式です。
fs − B > B ⇔ fs > 2·B (標本化定理/ナイキストの条件)
この下限 2B をナイキストレート、標本化周波数の半分 fN = fs/2 をナイキスト周波数と呼びます。fN が、その標本化レートで一意に表現できる周波数の上限です。逆に fs が 2B に満たないと、隣の複製がベースバンドに食い込みます。ベースバンドの上端付近にあった高周波成分が、あたかも低い周波数の成分であるかのように重なって現れる——これが折り返し(エイリアシング)です。
fs/2 を超える周波数 f(ただし f が fs 未満)は、ベースバンド内の見かけの周波数 f_alias に写ります。
f_alias = | f − round(f / fs) · fs | (0 〜 fs/2 に折り返す)
例えば fs = 8 kHz で 7 kHz(fN=4kHz を超える)の正弦波を標本化すると、|7 − 8| = 1 kHz の正弦波と区別がつかなくなります。高い音が低い音に化けて紛れ込むわけです。
一度エイリアシングが起きると、偽の低周波成分は本物の低周波成分と数値上まったく同じ標本列に混ざります。標本化してしまった後では、どんなデジタルフィルタを掛けても両者を分離できません。だからこそ対策は「標本化する前」に打つ必要があります。ここが初学者の誤解しやすい急所です。
アンチエイリアスフィルタ:標本化の前に削る
実世界の信号は理想的な帯域制限を満たさず、fs/2 以上の成分(雑音や高調波を含む)を必ず持ちます。そこで標本化の直前に、fs/2 未満の成分だけを通しそれ以上を落とすアナログのローパスフィルタ(アンチエイリアスフィルタ)を挿入します。折り返して混ざる前に高域を除去してしまえば、以後の複製は重なりません。
| 対策の段 | 配置 | できること | 限界 |
|---|---|---|---|
| アンチエイリアスフィルタ | 標本化の前(アナログ) | fs/2 以上を除去し折り返しを未然に防ぐ | 理想の急峻なカットは不可、遷移帯域が要る |
| 標本化レートを上げる | ADC の設定 | ナイキスト周波数を上げ余裕(ガードバンド)を作る | データ量と処理コストが増える |
| デジタルフィルタ | 標本化の後 | 帯域内の整形・ノイズ低減 | 既に折り返した成分は分離できない |
現実のフィルタは通過域と阻止域の間に有限の遷移帯域を持つため、fs/2 ちょうどで垂直に切ることはできません。そこで実務では、通過域上端(信号帯域 B)と fs/2 の間にガードバンドと呼ぶ余白を設けます。オーバーサンプリング(必要より高い fs で取る)を使うと fs/2 が高域へ離れ、緩やかななだらかなアナログフィルタで済むうえ、デジタル側で急峻な間引きフィルタをかけてから間引く(デシメーション)ことができます。アナログの負担をデジタルへ移せるのが、現代の ADC がオーバーサンプリングを多用する理由です。
標本化と再構成は表裏一体
標本列から連続信号を戻す再構成も、周波数領域で見れば「複製のうちベースバンドだけを取り出す」操作です。理想は fs/2 で完全に切れるローパスフィルタで、その時間領域での応答(インパルス応答)は sinc 関数です。したがって理想再構成は、各標本を sinc で内挿して足し合わせる補間になります。
x(t) = Σ_k x[k] · sinc( (t − k·Ts) / Ts )
sinc(u) = sin(π·u) / (π·u)
sinc は無限に広がり因果的でないため、そのままでは実装できません。実際の DAC はゼロ次ホールド(各標本値を次の標本まで保持する階段状の出力)で近似します。ゼロ次ホールドの周波数特性は sinc 型の緩やかな低下(アパーチャ効果)を持ち、高域がだれるため、後段のアナログフィルタやデジタルの補正で平坦化します。標本化(前段の帯域制限)と再構成(後段のローパス)は、同じナイキスト周波数を境界に対を成す処理だと捉えると全体が一本の線で繋がります。
- 標本化定理の条件は
fs > 2B(等号ではなく厳密に超える)。境界2Bはナイキストレート、fs/2はナイキスト周波数。 - 標本化はスペクトルを
fs間隔で複製する。複製が重なるとエイリアシング。 - エイリアス周波数は
fs/2を折り返し点としてf_alias = |f − n·fs|で求まる(0〜fs/2 に写る)。 - 折り返した成分は事後のデジタル処理では除去不能。対策は標本化前のアナログフィルタ。
- 理想再構成は sinc 補間、実機はゼロ次ホールド+補正(アパーチャ効果)。
まとめ
標本化定理は、帯域 B の信号なら fs > 2B(ナイキスト周波数 fs/2 が B を上回る)で標本列から完全に復元できることを保証します。要点は、(1) 標本化は周波数軸でスペクトルを fs 間隔に複製する操作であること、(2) fs/2 を超える成分は複製の重なりでベースバンドへ折り返り、偽周波数(エイリアス)として混入すること、(3) 折り返しは標本化後には分離できず、標本化前のアンチエイリアスフィルタで防ぐしかないこと、(4) 再構成は sinc 補間が理論解でゼロ次ホールドが実用近似であること。この原理は音声・画像・通信からセンサ計測、電源 のスイッチング測定、AI の特徴量前処理まで、離散化を伴うあらゆる場面で共通の土台になります。
信号処理・制御 Article
標本化定理とエイリアシングを実務で読む
TL;DRは入口です。実際に選ぶ・使う段階では、何を解決するか、何と比較するか、導入後にどこで詰まるかまで見る必要があります。
解決すること
標本化定理
比較で見る軸
難易度: advanced / カテゴリ: 信号処理・制御 / タグ数: 5
導入後に効く点
fs/2 を超える成分は、標本化でスペクトルが fs 間隔に複製される際に低域へ写り込み、別名(エイリアス)の偽周波数として現れる。一度折り返した成分は事後処理では分離できない。
先に潰すリスク
用語だけ覚えても、設計・実装・運用でどこに効くかを確認しないと判断を誤る。
- 難易度
- advanced
- カテゴリ
- 信号処理・制御
- タグ数
- 5
判断チェックリスト
- 自社の用途が「標本化定理 / ナイキスト周波数」に近いか確認する。
- 強みである「帯域が B までの信号は、標本化周波数 fs が 2B を超えれば標本列から完全に復元できる(標本化定理)。fs/2 をナイキスト周波数と呼び、これが表現できる周波数の上限。」が本当に評価軸になるか確認する。
- 注意点の「用語だけ覚えても、設計・実装・運用でどこに効くかを確認しないと判断を誤る。」を運用で吸収できるか確認する。
- 公開値や仕様値は、対象プラン・対象機種・対象リージョンまで確認する。
- 既存システム、ID、ネットワーク、監視、バックアップとの接続方法を先に洗い出す。
- 小さく試してから、本番移行、権限設計、障害時手順、コスト監視を決める。