Poisson wavelet explained

In mathematics, in functional analysis, several different wavelets are known by the name Poisson wavelet. In one context, the term "Poisson wavelet" is used to denote a family of wavelets labeled by the set of positive integers, the members of which are associated with the Poisson probability distribution. These wavelets were first defined and studied by Karlene A. Kosanovich, Allan R. Moser and Michael J. Piovoso in 1995–96.^{[1] [2]} In another context, the term refers to a certain wavelet which involves a form of the Poisson integral kernel.^[3] In still another context, the terminology is used to describe a family of complex wavelets indexed by positive integers which are connected with the derivatives of the Poisson integral kernel.^[4]

Wavelets associated with Poisson probability distribution

Definition

For each positive integer n the Poisson wavelet

\psi _{n(t)} is defined by\psi _{n(t)} = \begin{cases} \left( (t - n)/(n!) \right) t ^{n-1} e ^{-t} & for t\ge 0 \ 0 & for t< 0 .\end{cases}

To see the relation between the Poisson wavelet and the Poisson distribution let X be a discrete random variable having the Poisson distribution with parameter (mean) t and, for each non-negative integer n, let Prob(X = n) = p_n(t). Then we have

p _{n(t)} = (t ^{n})/(n!) e ^{-t} .

The Poisson wavelet

\psi _{n(t)} is now given by\psi _{n(t)} =- (d)/(dt) p _{n(t).}

Basic properties

\psi _{n(t)} is the backward difference of the values of the Poisson distribution: \psi _{n(t)=p} _{n(t)-p} _{n-1} (t).

• The "waviness" of the members of this wavelet family follows from

\int ^{inf ty} _{-inf ty} \psi _{n(t)dt} =0 .

• The Fourier transform of

\psi _{n(t)} is given\Psi(\omega) = (- i\omega)/(( 1+ i\omega) ^{n+1}) .

• The admissibility constant associated with

\psi _{n(t)} isC _{\psi _{n}} = \int ^{inf ty} _{-inf ty} (\left|\Psi ^{2} _{n(\omega)\right|})/(|\omega|) d\omega = (1)/(n) .

• Poisson wavelet is not an orthogonal family of wavelets.

Poisson wavelet transform

The Poisson wavelet family can be used to construct the family of Poisson wavelet transforms of functions defined the time domain. Since the Poisson wavelets satisfy the admissibility condition also, functions in the time domain can be reconstructed from their Poisson wavelet transforms using the formula for inverse continuous-time wavelet transforms.

If f(t) is a function in the time domain its n-th Poisson wavelet transform is given by

(W _{nf)(a,b)} = (1)/(\sqrt{|a|) }\int_^ f(t)\psi_n\left(\frac\right)\, dt

In the reverse direction, given the n-th Poisson wavelet transform

(W _{nf)(a,b)} of a function f(t) in the time domain, the function f(t) can be reconstructed as follows:f(t) = (1)/(C _{\psi _{n}}) \int ^{inf ty} _{-inf ty} \left[ \int ^{inf ty} _{-inf ty} \left\{(W _{nf)(a,b) (1)/(\sqrt{|a|)} } \psi_n\left(\frac\right)\,\right\} db\right] \frac

Applications

Poisson wavelet transforms have been applied in multi-resolution analysis, system identification, and parameter estimation. They are particularly useful in studying problems in which the functions in the time domain consist of linear combinations of decaying exponentials with time delay.

Wavelet associated with Poisson kernel

Definition

The Poisson wavelet is defined by the function^[3]

\psi(t) = (1)/(\pi) (1- t ^{2})/(( 1+ t ^{2)} ^{2})

This can be expressed in the form

\psi(t) = P(t) + t (d)/(dt) P(t) where P(t) = (1)/(\pi) (1)/(1+ t ^{2}) .

Relation with Poisson kernel

The function

P(t) appears as an integral kernel in the solution of a certain initial value problem of the Laplace operator.

This is the initial value problem: Given any

s(x) in L ^{p(R)} , find a harmonic function \phi(x,y) defined in the upper half-plane satisfying the following conditions:\int ^{inf ty} _{-inf ty} |\phi(x,y)| ^{p} dx \le c < inf ty , and\phi(x,y) → s(x) as y → 0 in L ^{p(R)} .

The problem has the following solution: There is exactly one function

\phi(x,y) satisfying the two conditions and it is given by\phi(t,y) = P _{y(t)\star} s(t)

where

P _{y(t) = (1)/(y)} P\left( (t)/(\right) = y) (1)/(\pi) (y)/(t ^{2+y} ^{2}) and where "\star " denotes the convolution operation. The function P _{y(t)} is the integral kernel for the function \phi(x,y) . The function \phi(x,y) is the harmonic continuation of s(x) into the upper half plane.

Properties

• The "waviness" of the function follows from

\int ^{inf ty} _{-inf ty} \psi(t) dt =0 .

• The Fourier transform of

\psi(t) is given by\Psi(\omega) = |\omega|e ^{-|\omega|} .

• The admissibility constant is

C _{\psi} = \int ^{inf ty} _{-inf ty} (\left|\Psi(\omega)\right| ^{2})/(|\omega|) d\omega =2 .

A class of complex wavelets associated with the Poisson kernel

Definition

The Poisson wavelet is a family of complex valued functions indexed by the set of positive integers and defined by^{[4] [5]}

\psi _{n(t) = (1)/(2 \pi)} ( 1- it) ^{-(n+1)} where n =1 , 2 , 3 , \ldots

Relation with Poisson kernel

The function

\psi _{n(t)} can be expressed as an n-th derivative as follows:\psi _{n(t)=} (1)/(2 \pi) (1)/(n! i ^{n}) (d ^{n})/(dt ^{n}) \left(( 1- it) ^{-1} \right) Writing the function ( 1- it) ^{-1} in terms of the Poisson integral kernel P(t) = (1)/(1+ t ^{2}) as ( 1- it) ^{-1} = P(t) + itP(t) we have\psi _{n(t) = (1)/(2 \pi)} (1)/(n!i ^{n}) (d ^{n})/(dt ^{n}) P(t) + i \left( (1)/(2 \pi) (1)/(n! i ^{n}) (d ^{n})/(dt ^{n}) \left(tP(t)\right)\right) Thus \psi _{n(t)} can be interpreted as a function proportional to the derivatives of the Poisson integral kernel.

Properties

The Fourier transform of

\psi _{n(t)} is given by\Psi _{n(\omega)} = (1)/(\Gamma(n +1 )) \omega ^{n} e ^{-\omega} u(\omega) where u(\omega) is the unit step function.

Notes and References

1. Karlene A. Kosanovich, Allan R. Moser and Michael J. Piovoso. The Poisson wavelet transform. Chemical Engineering Communications. 1996. 146. 1. 131–138. 10.1080/00986449608936485 .
2. Karlene A. Kosanovich, Allan R. Moser and Michael J. Piovoso. A new family of wavelets: the Poisson wavelet transform. Computers in Chemical Engineering. 1997. 21. 6. 601–620. 10.1016/S0098-1354(96)00294-3 .
3. Book: Roland. Klees. Roger. Haagmans . Wavelets in the Geosciences. 2000. Springer. Berlin. 18–20.
4. Book: Abdul J. Jerri. The Gibbs Phenomenon in Fourier Analysis, Splines and Wavelet Approximations. limited. 1998. Springer Science+Business Media. Dordrech. 978-1-4419-4800-7. 222–224.
5. Book: Wojbor A. Woyczynski. Distributions in the Physical and Engineering Sciences: Distributional and Fractal Calculus, Integral Transforms and Wavelets, Volume 1. 1997. Springer Science & Business Media. 9780817639242. 223.