畳み込み積分の理解


<div>
    <p>
        本記事では畳み込み積分について説明します. <br>
        基本的には連続な関数に対しての畳み込みを説明しますが, 離散な関数に対しての畳み込みも同様の考え方で大丈夫でしょう. <br>
        読み物感覚で見てもらって間違いとかあったらぜひ指摘ください. <br>
        うまく見えてなかったら<a href="https://study1ng.github.io/rinkou/blog.html"> github pages </a>を見てください.
    </p>
    <h2>前提知識</h2>
    <p>
        基本的に本記事では理論的な話をふんわりとした感じで書くので
        積分を区分求積法のイメージで理解できればだいたい理解できると思います。<br>
    </p>
    <h2>
        対象読者
    </h2>
    <p>
        畳み込み積分をはじめてきいた人<br>
        授業で習ったけどよくわからない人
    </p>
    <h2>フィルタ処理によるノイズ除去</h2>
    <p>

![soundwave.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2859784/50b9c4ff-dc04-4e87-a315-7eaad2015aab.png)
<br>
        音声信号を例にとってみます. 一般に音声をグラフに示すときには横軸に時間$t$を, 縦軸に振幅$x(t)$をとります.
        上の図をご覧いただきたい. 左側には私の声を示すグラフがあります. きれいな正弦波ですね. <br>
        しかし, 一般に録音を行うとどうしても環境音などの細かなノイズが入ってきてしまいます. そうすると, 右のグラフのようにノイズが乗っかかってきて
        ギザギザなグラフになってしまいます. <br>
        どうにかしてこのノイズを取り除いて, 左側の真の私の声を取り出す方法はないだろうか. <br>
        賢いsさんは考えました. <br>
        <blockquote>
            ノイズの平均はゼロとなるはずだから, 振幅を知りたい音声の時間を$t'$としたときに, $t'$の前後一秒間の振幅の平均をとれば
            ノイズは打ち消しあって, 音声の振幅が得られるのではないか.
        </blockquote>
        連続関数$x(t)$に対して, $t'$の前後一秒間の振幅の平均をとるというのは, 式に直せば
        $$
            \frac{1}{2} \int_{-1}^{+1} x(t' + t) dt
        $$
        となります. <br>
        つまり, $x(t)$の$[t'-1, t'+1]$の区間の合計$\int_{-1}^{+1} x(t' + t) dt$を区間の長さ$2$で割ったものです. <br>
        実際にこの処理をかけてみると, 
        
![mean_filter.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2859784/6051e74e-2673-4c42-982e-c2299824b007.png)
        <br>
        のようになります. <br>
        確かにグラフのギザギザはなくなりました. しかし, オリジナルの音声のグラフの振幅の絶対値の最大値は1.0であったのに対して, こちらは0.75程度となっています. <br>
        これではオリジナルの音声を復元できたとは言えないでしょう. <br>
        sさんはこの誤差の原因についてこのように考えました. <br>
        <blockquote>
            ノイズの関数をn(t)とし, オリジナルの音声の関数をx(t)としたときに, 録音された音声の関数はx(t) + n(t)となる. <br>
            単純に平均をとっただけではn(t)が平均されて0になるだけでなく, x(t)も平均されてしまうではないか. <br>
            しかし, 録音された音声のn(t)だけを取り出すことはできない. <br>
            ならばせめて影響を小さくするために, 前後一秒間ではなく, 前後0.1秒間の平均をとるのはどうだろうか.
        </blockquote>
        このとき, ノイズ除去された音声の関数は
        $$
            \frac{1}{0.2} \int_{-0.1}^{+0.1} x(t' + t) dt
        $$
        となります. <br>
        これの処理をかけた結果は
        ![smaller_mean_filter.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2859784/fd2a5d7e-99ea-44c2-a96b-1b30cff83ccd.png)
        <br>となります. <br>
        最大値はたしかに1.0程度となりましたが, グラフががくがくしていて正弦波には見えません. <br>
        sさんの苦戦を見ていたgさんはこのように言いました. <br>
        <blockquote>
            フィルタをかける時間区間を変えてもうまくいくかもしれないけど, <br>
            いい感じの時間区間を見つけるためには何度も実験して確認しないといけないよね. <br>
            データが違えばそのデータにあった時間区間を実験で見つけないといけないなんてめんどくさすぎる. <br>
            私なら時間は前後一秒間に固定して, フィルタを変えてみるかな. <br>
            例えば, フィルタの中央の値の影響を大きくして, その周りの値の影響を小さくするフィルタなんてどうかな. <br>
        </blockquote>
        そんなことを言いながら, gさんは
        $$
        \frac{3}{\sqrt{2\pi}} e^{-\frac{9t^2}{2}}
        $$
        を持ち出してみました. <br>
        この関数でフィルタを行うと, $$
        \int_{-1}^{+1} x(t' + t) \frac{3}{\sqrt{2\pi}} e^{-\frac{9(t)^2}{2}} dt
        $$というようになります. <br>
        このグラフの概形は下図のようになります. 
        ![gauss.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2859784/7e16a611-d9d1-4330-ab20-d9677fdad4ee.png)
        <br>グラフの$[-1, 1]$の区間の外ではほとんど関数の値は0に近い値になっています. 実際, このグラフの$[-1, 1]$の範囲を積分してみると
        0.9973程度になります. <br>
        この関数で処理してみると,
        
![gauss_filter.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2859784/d0d19405-8fce-45c8-8811-84372e2bed90.png)
        
のようになります. <br>
        単に平均値をとる処理よりかはきれいですね.
    </p>
    <h2>フィルタの一般化</h2>
    <p>
        このように, 信号$f(t)$に対しては, 用途や必要な精度に応じて, フィルタの関数や適用区間を変えることがあります. <br>
        どのようなフィルタや時間区間にも適用できるように, フィルタ処理を一般化すると
        $$
        filtered(t') = \int\limits_{S} f(t' + t) h(t) dt
        $$となります. <br>
        ただし, $S$というのは実数上の区間で, 特に$[-1, 1]$である必要はなく, 例えばフィルタの特性や, リアルタイム処理で未来のデータが得られない場合には$[-1, 0]$のようなフィルタを用いることもあるでしょう. <br>
ここで, 
$$
\tilde{f}_ S(t) = 
\begin{cases}
    f(t) & (t \in S) \\
    0 & (t \notin S)
\end{cases}
$$
$$
\underline{h}(x) = h(-x)
$$
という新たな記法を導入します. 一般的な数学の記号ではありませんので他のところでは使わないほうがいいです. <br>
        この一般化されたフィルタをさらに変形してみましょう. <br>
        $$
        \begin{align*}
            filtered(t') &= \int\limits_{S} f(t' + t) h(t) dt \\
            &= -\int\limits_{-S} f(t' - t) h(-t) dt \\ 
            &= -\int\limits_{-S} f(t' - t) \underline{h}(t)dt \\
            &= -\int_{-\infty}^{\infty} \tilde{f} _ S(t' - t) \underline{h}(t) dt \\
            &= -(\tilde{f} _ S * \underline{h})(t')
        \end{align*}
        $$
        なじみのある畳み込みの形になりました. <br>
        ここで, $\tilde{f} _ {S} $は単に信号を0でパディングした関数, $\underline{h}$はフィルタの関数を反転させたものなので, 
        畳み込みの正体というのはフィルタ処理といえるでしょう.
    </p>
    <h2>周波数領域での話</h2>
    <p>
        ノイズとは何でしょう. 再び最初の図の右のグラフを見てみましょう. とてもギザギザですよね. <br>
        ギザギザというのはすなわち, 値の変動が激しい高周波数成分が含まれているということです. <br>
        このノイズが含まれたギザギザのグラフにフィルタ処理を行うと, なめらかなグラフになります. <br>
        なめらかというのは値の変動が少ない低周波数成分の影響が強いということです. <br>
        つまり, 私たちが今までやってきたフィルタ処理というのは高周波数成分を除去して, 低周波数成分を強調する処理を行っているということです. <br>
        このような効果をもたらすフィルタをローパスフィルタと呼びます. <br>
        ただし, フィルタの中には高周波数成分を強調するフィルタもあります. <br>
        例えば, 画像処理でエッジを強調するフィルタなどは高周波数成分を強調するフィルタです. <br>
        このようなフィルタをハイパスフィルタと呼びます. <br>
        周波数の話をするならフーリエ変換$$
        \mathcal{F}(f(t)) = F(\omega) = \int_{-\infty}^{\infty} f(t) e^{-i\omega t} dt
        $$は通れないでしょう. <br>
        $F(\omega)$は一般に複素数となり, 信号に含まれる$a sin(\omega t + \phi)$のような成分の様々な情報を表します. <br>
        $F(\omega)$の虚部は成分の位相を表し, 実部は成分の振幅を表します. <br>
        フーリエ変換を畳み込み積分に適用すると, $$
        \mathcal{F}(f(t) * g(t)) = \mathcal{F}(f(t)) \cdot \mathcal{F}(g(t))
        $$となります. 
        つまり, 畳み込みによってフィルタをかけているとき, 周波数領域でもフィルタをかけているのです. <br>
        証明は教科書に書いてあると思います. 
    </p>
    <h2>高次元での話</h2>
    <p>
        ここまで説明してきた話は高次元でも成り立ちます. <br>
        フーリエ変換と畳み込み積分の関係の証明は一次元のときと同様に積分の順序交換を行えばいいと思います. <br>
        まあ現実的にはせいぜい三次元動画像の属する四次元くらいまでじゃないんですかね使うことがあるとしても.
    </p>
</div>