Dublin Core The Dublin Core metadata element set is common to all Omeka records, including items, files, and collections. For more information see, http://dublincore.org/documents/dces/. Title A name given to the resource note Description An account of the resource 日常のノート note 日常のノートの記録 note テキストの自由記述 <div> <div><span># 平均値フィルタによる雑音抑制</span></div> <br /> <div><span>平均値フィルタ（移動平均フィルタ）は、音声信号などに含まれる高周波数の雑音成分を抑制し、低周波数の主要成分を残すために用いられる基本的なディジタルフィルタである。本資料では、フィルタの概要から計算式、雑音抑制の理論、線形システムとしての位置づけを整理する。</span></div> <br /> <div><span>---</span></div> <br /> <div><span>## 1. ディジタルフィルタの概要</span></div> <br /> <div><span>### ディジタルフィルタのしくみ</span></div> <div><span>-</span><span> ディジタル入力信号 x(n) → </span><span>**ディジタルフィルタ**</span><span> → ディジタル出力信号 y(n)</span></div> <div><span>-</span><span> 足し算・掛け算を主体とした演算の組み合わせ</span></div> <div><span>-</span><span> この演算の組み合わせを適切に与えることによって、雑音除去が実現できる</span></div> <br /> <div><span>#### 図1: 平均値フィルタの信号の流れ</span></div> <div><span>```mermaid</span></div> <div><span>flowchart LR</span></div> <div><span> X["入力信号 x(n)"] --> A["加算・平均化"]</span></div> <div><span> A --> Y["出力信号 y(n)"]</span></div> <div><span> subgraph フィルタ内部</span></div> <div><span> A</span></div> <div><span> end</span></div> <div><span> %% M 個のサンプルを足し合わせて 1/M を掛ける</span></div> <div><span>```</span></div> <br /><br /> <div><span>### 今回扱うフィルタ</span></div> <div><span>**平均値フィルタ（移動平均フィルタ）**</span></div> <br /><br /> <div><span>## 2. 平均値フィルタ</span></div> <br /> <div><span>### M点の平均値フィルタの計算式</span></div> <br /> <div><span>$$</span></div> <div><span>y(n) = </span><span>\</span><span>frac</span><span>{</span><span>1</span><span>}{M} \{x(n) + x(n-</span><span>1</span><span>) + x(n-</span><span>2</span><span>) + </span><span>\</span><span>cdots + x(n-M+</span><span>1</span><span>)\}</span></div> <div><span> = </span><span>\</span><span>frac</span><span>{</span><span>1</span><span>}{M} </span><span>\</span><span>sum_</span><span>{k=</span><span>0</span><span>}^{M-</span><span>1</span><span>} x(n-k)</span></div> <div><span>$$</span></div> <br /> <div><span>### 説明</span></div> <div><span>出力信号のn番目のサンプルy(n)の値を、「入力信号のn番目からn-(M-1)番目までにわたる合計M個のサンプル（すなわちx(n), x(n-1), ⋯, x(n-M+1)）の平均値」として与える。</span></div> <br /> <div><span>### フィルタ係数</span></div> <div><span>平均値フィルタでは、すべてのフィルタ係数の値が </span><span>**1/M**</span><span> となる。</span></div> <br /> <div><span>$$</span></div> <div><span>y(n) = </span><span>\</span><span>frac</span><span>{</span><span>1</span><span>}{M}x(n) + </span><span>\</span><span>frac</span><span>{</span><span>1</span><span>}{M}x(n-</span><span>1</span><span>) + </span><span>\</span><span>frac</span><span>{</span><span>1</span><span>}{M}x(n-</span><span>2</span><span>) + </span><span>\</span><span>cdots + </span><span>\</span><span>frac</span><span>{</span><span>1</span><span>}{M}x(n-(M-</span><span>1</span><span>))</span></div> <div><span>$$</span></div> <br /> <div><span>（フィルタ係数は合計M個）</span></div> <br /> <div><span>---</span></div> <br /> <div><span>## 3. 平均値フィルタによる雑音除去の原理</span></div> <br /> <div><span>### 3.1 時間領域に基づく説明</span></div> <br /> <div><span>#### 雑音の特性</span></div> <div><span>-</span><span> 雑音は不規則な形であるため、激しく時間変動する</span></div> <br /> <div><span>#### 移動平均の効果</span></div> <div><span>-</span><span> 移動平均の演算によって、激しい時間変動の成分が抑制され、滑らかになった（平滑化された）出力信号が得られる</span></div> <div><span>-</span><span> 移動平均のサンプル数を増やすことによって、平滑化の度合いが強くなる</span></div> <div><span>-</span><span> </span><span>**注意**</span><span>：平均のサンプル数を増やし過ぎると、雑音とともに所望の信号の成分も過剰に失われる</span></div> <br /> <div><span>※参考：新型コロナ感染者の7日間移動平均も、これと同じ原理である</span></div> <br /> <div><span>### 3.2 周波数領域（スペクトル）に基づく説明</span></div> <br /> <div><span>#### 基本概念</span></div> <div><span>-</span><span> </span><span>**周波数領域**</span><span>：時間領域をフーリエ変換したもの</span></div> <div><span>-</span><span> </span><span>**スペクトル**</span><span>：時間領域の信号波形をフーリエ変換したもの</span></div> <br /> <div><span>#### 注目すべき3つのスペクトル</span></div> <div><span>1.</span><span> </span><span>**入力信号（音声+雑音）の振幅スペクトル |X(f)|**</span></div> <div><span>2.</span><span> </span><span>**平均値フィルタの振幅スペクトル |H(f)|**</span></div> <div><span>3.</span><span> </span><span>**出力信号（平均値フィルタの出力）の振幅スペクトル |Y(f)|**</span></div> <br /> <div><span>#### 平均値フィルタの振幅スペクトル |H(f)| とは</span></div> <div><span>平均の演算で使っているM個のフィルタ係数（1/M, 1/M, ⋯, 1/M）を「信号」と見なして、その「信号」のフーリエ変換の絶対値をとったもの。</span></div> <br /> <div><span>#### 図2: 平均値フィルタの典型的な振幅特性（概念図）</span></div> <div><span>```</span></div> <div><span>|H(f)|</span></div> <div><span>1.0 ┤───────╮</span></div> <div><span> │ │╲</span></div> <div><span>0.5 ┤ │ ╲</span></div> <div><span> │ │ ╲</span></div> <div><span>0.0 ┼───────┴───╲── f</span></div> <div><span> 低周波高周波</span></div> <div><span>```</span></div> <div><span>低周波域では係数和がほぼ 1 のため透過し、高周波域ではサンプル同士が打ち消し合うため減衰する。</span></div> <br /> <div><span>#### 重要な関係式</span></div> <div><span>$$</span></div> <div><span>|Y(f)| = |H(f)| </span><span>\</span><span>times |X(f)|</span></div> <div><span>$$</span></div> <br /> <div><span>※この理論は、平均値フィルタを含めた任意の線形システムに対して成り立つ</span></div> <br /> <div><span>### 3.3 雑音抑制できる理由（周波数領域の観点から）</span></div> <br /> <div><span>#### フィルタ入力信号 |X(f)| について</span></div> <div><span>-</span><span> </span><span>**声の主要成分**</span><span>：低い周波数領域に集中している</span></div> <div><span>-</span><span> </span><span>**雑音の成分**</span><span>：高い周波数領域に集中している</span></div> <br /> <div><span>#### 平均値フィルタ |H(f)| について</span></div> <div><span>-</span><span> </span><span>**低い周波数領域**</span><span>：|H(f)|の値は </span><span>**1に近い**</span><span>（通過させる）</span></div> <div><span>-</span><span> </span><span>**高い周波数領域**</span><span>：|H(f)|の値は </span><span>**0に近い**</span><span>（抑制する）</span></div> <br /> <div><span>#### フィルタ出力信号 |Y(f)| について</span></div> <div><span>上記の特徴と、|Y(f)| = |H(f)| × |X(f)| より：</span></div> <div><span>-</span><span> </span><span>**低い周波数領域の成分（声）**</span><span>：ほぼ残る</span></div> <div><span>-</span><span> </span><span>**高い周波数領域の成分（雑音）**</span><span>：抑制される</span></div> <br /> <div><span>### まとめ：平均値フィルタで雑音抑制できる理由</span></div> <div><span>1.</span><span> 人の声は低周波領域に集中している</span></div> <div><span>2.</span><span> 雑音は高周波領域に集中しやすい</span></div> <div><span>3.</span><span> 平均値フィルタは、低周波領域の成分を通過させ、高周波領域の成分を抑制させる特性をもつ（</span><span>**ローパスフィルタ**</span><span>）</span></div> <br /> <div><span>---</span></div> <br /> <div><span>## 4. 補足：線形システムについて</span></div> <br /> <div><span>### 4.1 線形システムとは</span></div> <div><span>以下に述べる線形性が成立するシステムのこと。</span></div> <br /> <div><span>#### 線形性</span></div> <div><span>関数f(n), f₁(n), f₂(n)がそれぞれ与えられているとする。このとき、任意の実数aとbに対して以下の2つの関係式を成り立たせる作用素Sは、線形性を満足する作用素である。</span></div> <br /> <div><span>$$</span></div> <div><span>S[af(n)] = aS[f(n)]</span></div> <div><span>$$</span></div> <br /> <div><span>$$</span></div> <div><span>S[af_</span><span>1</span><span>(n) + bf_</span><span>2</span><span>(n)] = aS[f_</span><span>1</span><span>(n)] + bS[f_</span><span>2</span><span>(n)]</span></div> <div><span>$$</span></div> <br /> <div><span>**平均値フィルタは、線形性を満足するシステムである。**</span></div> <br /> <div><span>### 4.2 線形性を満たさないシステムの例</span></div> <div><span>以下の作用素Sによって与えられるシステム</span></div> <br /> <div><span>$$</span></div> <div><span>y(n) = S[x(n)] = x^</span><span>2</span><span>(n)</span></div> <div><span>$$</span></div> <br /> <div><span>は、線形性を満たさない（すなわち</span><span>**非線形システム**</span><span>である）。</span></div> <br /> <div><span>#### 理由</span></div> <div><span>入力信号をax(n) と与えたとき、出力は</span></div> <br /> <div><span>$$</span></div> <div><span>y(n) = S[ax(n)] = a^</span><span>2</span><span> x^</span><span>2</span><span>(n)</span></div> <div><span>$$</span></div> <br /> <div><span>となり、y(n) ≠ aS[</span><span>x(n)</span><span>] と結論づけられるので、このシステムは線形性を満たさない。</span></div> <br /> <div><span>---</span></div> <br /> <div><span>## 5. まとめ</span></div> <br /> <div><span>### 平均値フィルタにおいて成立する理論（振幅スペクトルに基づいた理論）</span></div> <div><span>$$</span></div> <div><span>|Y(f)| = |H(f)| </span><span>\</span><span>times |X(f)|</span></div> <div><span>$$</span></div> <br /> <div><span>※この理論は、平均値フィルタを含めた任意の線形システムに対して成り立つ</span></div> <br /> <div><span>### 平均値フィルタを使って雑音抑制できる理由（振幅スペクトルの観点から）</span></div> <div><span>1.</span><span> 人の声は低周波領域に集中している</span></div> <div><span>2.</span><span> 雑音は高周波領域に集中しやすい</span></div> <div><span>3.</span><span> 平均値フィルタは、低周波領域の成分を通過させ、高周波領域の成分を抑制させる特性をもつ（</span><span>**ローパスフィルタ**</span><span>）</span></div> </div> Dublin Core The Dublin Core metadata element set is common to all Omeka records, including items, files, and collections. For more information see, http://dublincore.org/documents/dces/. Title A name given to the resource デジタル信号処理１０ Creator An entity primarily responsible for making the resource unjuno Date Created Date of creation of the resource. 2025-11-26 Rights Information about rights held in and over the resource CC BY-NC 4.0 - 非営利目的のみ許可