Product of convolution

In Mathematical, the produces convolution of two real functions or complex F and G generally notes “ \ ast” and is written:

(F \ ast G) (X) = \ int_ {- \ infty} ^ {+ \ infty} F (x-t) \ cdot G (T) \ cdot dt = \ int_ {- \ infty} ^ {+ \ infty} F (T) \ cdot G (x-t) \ cdot dt

One can regard this formula as a generalization of the Moving average idea of .

So that this definition has a direction, it is necessary to carry out certain assumptions on F and G , for example if these two functions are summable S their product of convolution is defined for almost all X and is itself summable.

Properties of the product of convolution

  • the product of convolution is commutative.
(F \ ast G) (X) \ stackrel {\ mathrm {def.}} {=} \ int_ {- \ infty} ^ {+ \ infty} F (x-t) \ cdot G (T) \, dt

\ int_ {- \ infty} ^ {+ \ infty} F (T) \ cdot G (x-T) \, D (x-T)

\ int_ {- \ infty} ^ {+ \ infty} F (T) \ cdot G (x-T) \, dT \;

\ stackrel {\ mathrm {def.}} {=} (G \ ast F) (X)

Where T=x-t, is t=x-T.

  • the product of convolution is distributive

(F \ ast (g+h)) (X) \ stackrel {\ mathrm {def.}} {=} \ int_ {- \ infty} ^ {+ \ infty} F (x-t) \ cdot (G (T) +h (T)) \, dt

\ int_ {- \ infty} ^ {+ \ infty} \ cdot G (T) + F (x-t) \ cdot H (T) \, dt

= \ int_ {- \ infty} ^ {+ \ infty} F (x-t) \ cdot G (T) \, dt + \ int_ {- \ infty} ^ {+ \ infty} F (x-t) \ cdot H (T) \, dt
\ stackrel {\ mathrm {def.}} {=} (F \ ast G) (X) + (F \ ast H) (X)
  • the product of convolution corresponds to the multiplication of the transformed of Fourier of the functions

F \ ast G = \ mathcal {F} ^* \ left (\ mathcal {F} (F) \ cdot \ mathcal {F} (G) \ right)

where \ mathcal {F} indicates the transformation of Fourier and \ mathcal {F} ^* the opposite transformation of Fourier. The principal interest of the calculation of the product of convolution by transforms of Fourier is that these operations are less expensive in time for a computer than the direct calculation of the integral.

Use of the product of convolution

The product of convolution is used in the treatment of the signal, when filters are used (low-pass, high-pass, band pass). If there is an entering signal S (T) and a filter element having a transfer transfer function H (T) then the output signal S_s (T) will be the convolution of these two functions:

S_s (T) =S \ ast H

Another use of the products of convolution is in the field of quantum mechanics, where one carries out products of convolution starting from the functions of wave will bra and ket.

In a general way, one can write the linear differential equations corresponding to many problems physique in the shape of the product of convolution of an operator by a function describing the system. One can then solve in a generic way the problem by determining the reverse of convolution of the operator (called Fonction of Green). Joseph Fourier was at the origin of this method when he sought to solve the equation of heat. Its modern formulation has to await the arrival of the theory of the distributions introduced by Laurent Schwartz.

The product of convolution spreads with many algebras of a group, for example with the algebras of a group finished. So moreover the group is abelian, then the theory of the harmonic Analyze on an abelian group finished allows to establish all the traditional results of the product of convolution.

Approaches popularized

The manner simplest to represent the product of convolution consists in considering the Fonction δ of Dirac δ has ( X ); this function is worth 0 if X ≠ has , and its integral is worth 1. This can seem odd at first sight, one can imagine it like the limit of a succession of functions, bell-shaped curves or rectangles having all same surface 1, but increasingly fine (thus increasingly high); when the width of the curves tends towards 0, its height tends towards +∞, but surface remains equal to 1. For practical reasons, one often represents Dirac as a stick positioned in is and height 1.


Dirac: limit of a succession of functions

Because of his form, one calls also sometimes Dirac “function impulse”. The product of convolution by Dirac δ has corresponds to a translation of the initial function of a value of has

f \ ast \ delta_a (X) = F (x-a)


Produces convolution of a function by Dirac

One sees that δ0 leaves invariant a function, it is the neutral element of the product of convolution

f \ ast \ delta_0 (X) = F (x-0)

If one now considers the product of convolution by a balanced Somme of two diracs (α.δ has + β.δ B ), one obtains the superposition of two dilated curves.


Produces convolution of a function by a balanced sum of two diracs

Let us consider a function now carries Pa, b ; it is a function which is worth 1 ( Ba ) between has and B , and 0 elsewhere (its integral is worth 1). This function can be seen like a succession of diracs. The convolution of F by Pa, b thus will be obtained while making slip F on the interval. One obtains a “widening” of F .


Produit convolution of a function by a function carries

If one considers an unspecified function now G , one can see G like a succession of diracs balanced by the value of G at the point considered. The product of convolution of F by G is thus obtained while making slip the function F and by dilating it according to the value of G .


Produces convolution of a function by an unspecified function

The product of convolution and filtering

The product of convolution is related to the concept of filtering under two conditions, namely the linearity and independence of the filter with respect to time (Système invariant). From these two conditions, the operator of convolution can be built. The convolution corresponds to the response of the filter to a given entry (noted e (T) ). The filter is entirely characterized by its impulse response h (T) . Put in equation, the answer of the filter is s (T) = H (T) \ ast E (T) .

The construction of the operator of convolution is worked out in the following way. First of all, one is interested in the two conditions imposed on the filter. One notes f (E) the filtering carried out by the filter on the entry e. The linearity of the filter implies that:

f (\ lambda E) = \ lambda F (E)

f (e_1+e_2) =f (e_1) +f (e_2)

One can note that the response of the filter to a null signal is null. The independence of time is summarized by:

f (e^d) = (F (E))^d

Where e^d is the signal e delayed d quantity.

From there, one can build the response of the filter linear and independent of time to the entry e (T) . Indeed, as the filter is linear, one can break up the signal e (T) into independent parts, using a whole of signals e_i with disjoined supports compact so that e (T) = \ Sigma_i e_i (T) . One injects each part of the signal in the filter then one summons the various answers. Thus filtering will give: f (E) = \ Sigma_i F (e_i) . This temporal decomposition of e (T) can be carried out in a recursive way on the signals e_i (T) . At the end, one obtains a succession of signals whose support is summarized at a point. These signals, elementary because nondecomposable temporally, correspond each one of them to the distribution of Dirac \ delta (T \ tau) centered in \ tau with an amplitude e (\ tau) , the impulse is written \ delta (T \ tau) .e (\ tau) . It is enough to summon all the impulses according to the variable \ tau to obtain the signal e (T) :

e (T) = \ int \ delta (T \ tau) .e (\ tau) D \ tau

One applies the operation of filtering to e (T) . As the filter is linear and independent of time, we have:

f (E) = \ int F (\ delta (T \ tau) .e (\ tau)) D \ tau

= \ int E (\ tau) .f (\ delta (T \ tau)) D \ tau \ \ (linearity)
= \ int E (\ tau). (F \ circ \ delta) (T \ tau) D \ tau \ \ (independence of time)

The response of the filter f to the impulse \ delta (T) is named the impulse response of the filter h (T) . Finally one a:

f (E) = \ int E (\ tau) .h (T \ tau) D \ tau

who is only other than the product of convolution.

In conclusion: if the filter is linear and independent of time, then it is entirely characterized by its answer h (T) and the answer of the filter to the entry e (T) is given by the operator of convolution.

Another fundamental conclusion of the filters linear and independent of time: if one enters a signal e (T) =e^ {2 \ pi \ jmath F T} , the output signal will be:

s (T) = \ int e^ {2 \ pi \ jmath F \ tau} H (T \ tau) D \ tau = \ int e^ {2 \ pi \ jmath F (T \ tau)} H (T) D \ tau

s (T) = e^ {2 \ pi \ jmath F T} \ int e^ {- 2 \ pi \ jmath F \ tau} H (\ tau) D \ tau

s (T) = e^ {2 \ pi \ jmath F T} H (F)

s (T) will be also a signal of the form e^ {2 \ pi \ jmath F T} with the factor H (F) near. This factor is, neither more, nor less, that the Transformée of Fourier of h (T) .

Random links:Ultima IV | Brenne | Conjecture Monstrous Moonshine | Deliver bedside | Alex Gerry | Wikify