畳み込み

2つの正方形による畳み込み。解として得る波形は三角波となる。黄色の領域で示されている面積が2つの方形波の合成積である。
正方形がRC回路に入力された場合の出力信号波形を得るために、RC回路のインパルス応答と方形波の畳み込みを行っている。 黄色の領域で示されている面積が合成積である。
畳み込み（たたみこみ、convolution）とは、関数 g を平行移動しながら関数 f に重ね足し合わせる二項演算である。あるいはコンボリューションとも呼ばれる。

**定義**

**一次元**

**連続**
連続関数 f, g の畳み込み f ∗ g は以下のように定義される：

積分を用いて2つの関数を合わせることから畳み込み積分、合成積、重畳積分とも呼ばれる。
積分範囲は関数の定義域に依存する。通常は区間 (−∞, +∞) で定義される関数を扱うことが多いので、積分範囲は −∞ から +∞ で計算されることが多い。一方 f, g が有限区間でしか定義されない場合には、g(t − τ) が定義域内に入るように f, g を周期関数と見なして計算される。この周期関数と見なして畳み込みをすることを循環畳み込み（じゅんかんたたみこみ、cyclic convolution）と呼ぶ。

**離散**
離散信号 f, g の畳み込み f ∗ g は以下のように定義される：

すなわち積分のかわりに総和を使って同様に定義される。そのため畳み込み和・重畳和とも呼ばれる。
総和の範囲も関数の定義域に依存し、関数が有限区間でしか定義されていない場合は周期関数とみなして畳み込み演算が行われる。また定義域外の値を 0 と定義し直した関数での畳み込みがよく行われる。これを線形畳み込み（せんけいたたみこみ、linear convolution）あるいは直線畳み込み（ちょくせんたたみこみ）と呼ぶ。

**高次元**
上の複素数値函数 fと g の畳み込みは、それ自身が Rd 上の複素数値函数として

で定義されるものであるが、右辺の積分が存在してこれが定義可能となるには、fと g が無限遠において十分Vanish at infinity必要がある。とはいえ、たとえば g が無限遠において爆発するとしても、その影響は f が十分に急減少であれば容易に打ち消すことができるから、この積分の存在条件は込み入ったものも考え得る。この問題をクリアする函数の条件としてよく用いられる場合を以下に挙げる。

**コンパクト台付き函数**
函数 f と g がともにコンパクト台連続函数ならば、それらの畳み込みは存在して、やはりコンパクト台連続函数となるloc=Chapter 1。より一般に、一方がコンパクト台、他方が局所可積分函数ならば、畳み込み f ∗ g が定義されて連続である。
上では両者が局所自乗可積分の場合、あるいは両者がともに半無限区間  (あるいはともに ) に台を持つ場合でも畳み込みが定まる。

**可積分函数**
函数 f と g がともにL1(Rd)に属するルベーグ可積分函数ならば、それらの畳み込み f ∗ g が存在してやはり可積分であるloc=Theorem 1.3。これはトネリの定理の帰結である。このことは  に属する数列の離散畳み込みや、より一般の群上の  の畳み込みでも成立する。
同様にして、 f ∈ L1(Rd) と g ∈ Lp(Rd) が 1 ≤ p ≤ ∞ のとき、 f ∗ g ∈ Lp(Rd) かつ
 
を満たす。特に 1=p = 1 のとき、これにより L1 は畳み込みを積としてバナッハ代数を成す(また、等号成立は f と g がともに殆ど至る所非負のときである)。

より一般に、畳み込みに対するヤングの不等式により、畳み込み積は適当な -空間上の連続双線型演算となることが従う。具体的に書けば、 1 ≤ p,q,r ≤ ∞ が
 
なる関係を満足するとして、
 
となるから、畳み込み積は  なる連続双線型写像を定めている。
畳み込みに対するヤングの不等式、循環畳み込みや離散畳み込みなどほかの文脈でも成立する。また、 R 上では先に掲げた不等式はより厳しく評価できる: 先と同様の関係を持つ 1 1(T) を固定して、ヒルベルト空間 L2(T) に作用するよく知られた作用素:
 
がとれる。作用素 T はコンパクト作用素である。直接計算により、その随伴作用素 T* は  による畳み込みであることが示せる。上で掲げた可換性により、 T は正規作用素 (1=T*T = TT* である。また T は平行移動作用素とも可換である。そのような畳み込み作用素と平行移動作用素全体の成す作用素族を S とすれば、 S は正規作用素からなる可換族である。ヒルベルト空間上のスペクトル論に従えば、 S を同時対角化する正規直交基底  が存在して、これが円周上の畳み込みを特徴付ける。具体的には
 
がちょうど T の指標の全体の成す集合に一致する。この基底に属する各畳み込み作用素がコンパクト乗算作用素であることが、上で述べた循環畳み込みに対する畳み込み定理としてみることができる。
離散畳み込みの例は位数 n の有限巡回群をとる。この場合の畳み込み作用素は巡回行列によって表現され、離散フーリエ変換によって対角化することができる。
同様の結果が (アーベルとは限らない) コンパクト群 G に対しても知られている: 有限次元ユニタリ表現の行列要素の全体が  の正規直交基底を成し (Peter–Weyl theorem)、適当な意味での畳み込み定理が (フーリエ変換に基づく調和解析の他の多くの側面とともに) 引き続き満足される。

**群上の測度の畳み込み**
位相群 G 上の有限ボレル測度 μ と ν に対し、それらの畳み込み μ ∗ ν は G の各可測部分集合 E に対して
 
で定義され、やはり有限測度となる。実際、全変動に関するヤングの不等式

が満足される。 G が局所コンパクト群で左ハール測度 λ を持ち、 μ と ν が λ に関して絶対連続で各々密度函数を持つ場合には、畳み込み μ ∗ ν もまた絶対連続で、その密度函数は各測度の密度函数の (通常の函数としての) 畳み込みに一致する。
考える位相群が実数の加法群 (R, +) のとき、その上の確率測度 μ と ν をとれば、測度の畳み込み μ ∗ ν は、分布 μ および ν に従う独立確率変数 X および Y の和 X + Y の確率分布に対応する。

**性質**
積分演算に由来する性質として以下の性質がある。- 交換律: - 結合律: - 分配律: - スカラー倍: , a は任意の複素数。- 微分: , D は微分演算子（離散系の場合は ）- redirect=1: ,  はフーリエ変換

**畳み込み定理**
は次の式で示される。  

ここで  はフーリエ変換である。この定理によりフーリエ変換を使って畳み込み演算を単純な掛け算に変換することが出来る。この定理はラプラス変換・Z変換やメリン変換といった変換に対しても適用できる。

**計算**
畳み込み演算を実際に計算する際には様々な技法が利用される。

**離散時間信号**
離散時間信号の畳み込み和は定義通りの畳み込み計算でなく、redirect=1を利用して計算される場合が多い。
離散時間信号では高速フーリエ変換 (FFT) を用いることで少ない計算量でフーリエ変換を実行できる。ゆえに関数 f, g をFFTにより周波数表現へ変換し、その積を取って、最後にIFFTで時間表現へ戻すことで高速に畳み込み和を計算できる。

**有限直線畳み込み**
有限長信号の直線畳み込みは「範囲外が0埋めされた関数」と見なすことで通常の直線畳み込み演算を適用できる。しかし値が0と事前に確定している領域を実際に計算するのは無駄であり、実践的には有限区間に限って計算がおこなわれる。畳み込みは2つの関数を重ね合わせる演算であるため「0埋めした領域を畳み込んだ場合に出力扱いするか」にバリエーションが存在する。以下はその一例である。- valid: 両関数が非ゼロの領域のみ出力 (L=N-M+1)- full: 両関数が非ゼロでない限り出力 (L=N+M-1)- same: 中央部分を入力と同じ長さだけ出力 (L=N)
fullは「両側 M-1 ゼロパディング+valid直線畳み込み」と同義であり、sameは「両側 (M-1)/2 ゼロパディング+valid直線畳み込み」と同義である。

**応用**

**確率測度における畳み込み**
集合関数の一種である確率測度の畳み込みは次のように表現される。確率測度 μ1, μ2 において任意のボレル集合 B に対し、

と表現される。ここで1BはBの定義関数である。これは μ1, μ2 を集合関数として捉えて、変数変換することで求まる。これにより、μ1, μ2 を分布に持つ確率変数 X, Y においてその和 X + Y の分布が畳み込みにあたることが分かる。

**多項式の掛け算**
多項式の掛け算の結果の係数列は、元の多項式の係数列の線形畳み込みになる。実際

であり、掛け算の結果の係数が a*b となる。

**線形システム**
電気回路といった古典的な時不変（シフト不変）線形システムは、任意の入力 x(t) に対する出力 y(t) が x(t) とインパルス応答 h(t) の畳み込みで記述できる：

ここで特に、入力 x(t) がデルタ関数 δ(t) のとき出力は h(t) そのものになる。
ここで上式の両辺をフーリエ変換もしくはラプラス変換（離散系の場合はZ変換）すると畳み込み定理より下式のようになる。

ここで、
 
を伝達関数といい、この式は古典制御論の基礎となっている。

**音響学**
エコーは元の音波と、音を反射するさまざまな物体に因る特性（インパルス応答）との畳み込みで記述される。カラオケやシンセサイザーに搭載されているエコー機能は、この畳み込みの効果を電気回路もしくはコンピュータでシミュレートすることで実現している。

**光学および画像処理**
撮像時のブレなどの多くのぶれ (blur) は畳み込みで記述できる。例えば、ピントがぼけた写真は、ピントがあった仮想的な画像と、絞りの特性を示す円との畳み込みである。また被写体等の動きによるブレも、静止した仮想的な画像と動きの特性との畳み込みであり、グラフィックソフトウェアのモーションブラーはこの畳み込み演算を計算によりシミュレートすることで実現している。
画像認識においても、異なるスケールの画像を認識するにあたり、畳み込みでぶれをつくってから、画像処理することがある。

**統計学**
がそれぞれ独立な連続型確率変数とすると、和の  の確率密度関数は畳み込みによって与えられる。X, Y の確率密度関数をそれぞれ  と表記すると、S の密度関数は以下の式で与えられる。

**歴史**
畳み込み積分が用いられた最初期の例の一つは d'Alembert (1754) Recherches sur différents points importants du système du monde におけるテイラーの定理の導出にあるp=2。また
 
の形の式は label=Lacroix Treatise on differences and series百科辞典シリーズ Traité du calcul différentiel et du calcul intégral, Chez Courcier, Paris, 1797-1800. の最後の三巻 の505頁で用いられp=4、そのすぐ後にLaplace、Fourier、Poissonらの研究に畳み込み演算が現れている。名称自体が広く用いられるようになるには1950年代あるいは1960年代を待たなければならない。それに先立ってはfaltung（「畳み込み」）、composition product（「合成積」）、 superposition integral（「重ね合わせ積分」）などとも呼ばれ、あるいはカールソンの積分 (Carson's integral)
とも言った。現代的な定義がより古い用例に馴染むわけでもないが、それでも早くは1903年ごろには出現している。
合成積の特別の場合としての演算
 
は Volterra (1913)
"Leçons sur les fonctions de lignes" にある

**脚注**

**注釈**

**出典**

**参考文献**
 Dominguez-Torres, Alejandro (Nov 2, 2010). "Origin and history of convolution". 41 pgs. http://www.slideshare.net/Alexdfar/origin-adn-history-of-convolution. Cranfield, Bedford MK43 OAL, UK. Retrieved Mar 13, 2013.- isbn=3-540-12104-8 - isbn=0-691-08078-X
*mr=0493420-  isbn=0-471-52364-X.

**関連項目**
逆畳み込み
*ブラインド・デコンボリューション
*インパルス応答 - 伝達関数
*フーリエ変換 - ラプラス変換
*軟化子
*基本解
*自己回帰移動平均モデル
*巡回行列
*アナログ信号処理
*コーシー積
*積和演算

**外部リンク**
 合成積（畳み込み）の意味と応用３つ- title=Convolution- title=convolution- title=Convolution of functions- title=Convolution transform-  The Joy of Convolution Java Applet を使った視覚的な畳み込みの説明-  Examples of sampled impulse responses to be used in convolution reverbs (Fokke Van Saane)-  Examples of impulse responses synthesized from oscillator spectra, to be used in convolution reverbs (Emmanuel Deruty)-  BruteFIR; A software for applying long FIR filters to multi-channel digital audio, either offline or in realtime.-  Freeverb3 Reverb Impulse Response Processor: DSP library with convolution engines.