Wiener–Hopf method

Wiener–Hopf method

The Wiener–Hopf method is a mathematical technique widely used in applied mathematics. It was initially developed by Norbert Wiener and Eberhard Hopf as a method to solve systems of integral equations, but has found wider use in solving two-dimensional partial differential equations with mixed boundary conditions on the same boundary. In general, the method works by exploiting the complex-analytical properties of transformed functions. Typically, the standard Fourier transform is used, but examples exist using other transforms, such as the Mellin transform.

In general, the governing equations and boundary conditions are transformed and these transforms are used to define a pair of complex functions (typically denoted with '+' and '-' subscripts) which are respectively analytic in the upper and lower halves of the complex plane, and have growth no faster than polynomials in these regions. These two functions will also coincide on some region of the complex plane, typically, a thin strip containing the real line. Analytic continuation guarantees that these two functions define a single function analytic in the entire complex plane, and Liouville's theorem implies that this function is an unknown polynomial, which is often zero or constant. Analysis of the conditions at the edges and corners of the boundary allows one to determine the degree of this polynomial.

Wiener–Hopf decomposition

The key step in many Wiener–Hopf problems is to decompose an arbitrary function Phi into two functions Phi_{pm} with the desired properties outlined above. In general, this can be done by writing

: Phi_+(alpha) = frac{1}{2pi i} int_{C_1} Phi(z) frac{dz}{z-alpha}

and

: Phi_-(alpha) = - frac{1}{2pi i} int_{C_2} Phi(z) frac{dz}{z-alpha},

where the contours C_1 and C_2 are parallel to the real line, but pass above and below the point z=alpha, respectively.

Similarly, arbitrary scalar functions may be decomposed into a product of +/- functions, i.e. K(alpha) = K_+(alpha)K_-(alpha), by first taking the logarithm, and then performing a sum decomposition. Product decompositions of matrix functions (which occur in coupled multi-modal systems such as elastic waves) are considerably more problematic since the logarithm is not well defined, and any decomposition might be expected to be non-commutative. A small subclass of commutative decompositions were obtained by Khrapkov, and various approximate methods have also been developed.

Example

Let us consider the linear partial differential equation

:oldsymbol{L}_{xy}f(x,y)=0,

where oldsymbol{L}_{xy} is a linear operator which contains derivatives with respect to x and y, subject to the mixed conditions on y=0, for some prescribed function g(x),

:f=g(x) for xleq 0, quad f_{y}=0 when x>0.

and decay at infinity i.e. f ightarrow 0 as oldsymbol{x} ightarrow infty. Taking a Fourier transform with respect to x results in the following ordinary differential equation

: oldsymbol{L}_{y}hat{f}(k,y)-P(k,y)hat{f}(k,y)=0,

where oldsymbol{L}_{y} is a linear operator containing y derivatives only, P(k,y) is a known function of y and k and

: hat{f}(k)=int_{-infty}^{infty} f(x,y)e^{-ikx} extrm{d}x.

If a particular solution of this ordinary differential equation which satisfies the necessary decay at infinity is denoted hat{F}(k,y), a general solution can be written as

: hat{f}=C(k)hat{F}(k,y),

where C(k) is an unknown function to be determined by the boundary conditions on y=0.

The key idea is to split hat{f} into two separate functions, hat{f}_{+} and hat{f}_{-} which are analytic in the lower- and upper-halves of the complex plane, respectively

: hat{f}_{+}(k,y)=int_{0}^{infty} f(x,y)e^{-ikx} extrm{d}x,

: hat{f}_{-}(k,y)=int_{-infty}^{0} f(x,y)e^{-ikx} extrm{d}x.

The boundary conditions then give

: hat{f}_{-}(k,0)+hat{f}_{+}(k,0)=hat{g}(k)+hat{f}_{+}(k,0)=C(k)F(k,0)

and, on taking derivatives with respect to y,

: hat{f}'_{-}(k,0)+hat{f}'_{+}(k,0) =hat{f}'_{-}(k,0) = C(k)F'(k,0).

Eliminating C(k) yields

: hat{g}(k)+hat{f}_{+}(k,0) - hat{f}'_{-}(k,0)/K(k) = 0,

where

: K(k)=frac{F'(k,0)}{F(k,0)}.

Now K(k) can be decomposed into the product of functions K^{-} and K^{+} which analytic in the upper-half plane or lower-half plane, respectively

: K(k)=K^{+}(k)K^{-}(k),

: hbox{log}K^{-}=frac{1}{2pi i}int_{-infty}^{infty} frac{hbox{log}(K(z))}{z-k} extrm{d}z,quadhbox{Im}k>0,

: hbox{log}K^{+}=-frac{1}{2pi i}int_{-infty}^{infty} frac{hbox{log}(K(z))}{z-k} extrm{d}z,quadhbox{Im}k<0.

Consequently,

: K_{+}(k)hat{g}_{+}(k) + K_{+}(k)hat{f}_{+}(k,0) =hat{f}'_{-}(k,0)/K_{-}(k) - K_{+}(k)hat{g}_{-}(k),

where it has been assumed that g can be broken down into functions analytic in the lower-half plane g_{+} and upper-half plane g_{-}, respectively. Now, as the left-hand side of the above equation is analytic in the lower-half plane, whilst the right-hand side is analytic in the upper-half plane, analytic continution guarantees existence of an entire function which coincides with the left- or right-hand sides in their respective half-planes. Furthermore, since it can be shown that the functions on either side of the above equation decay at large k, an application of Liouville's theorem shows that this entire function is identically zero, therefore

:hat{f}_{+}(k,0) = -hat{g}_{+}(k),

and so

: C(k) = frac{hat{g}(k)-hat{g}_{+}(k)}{F(k,0)}.

See also

* Wiener filter

External links

*
* [http://www.wikiwaves.org/index.php/Category:Wiener-Hopf Wiener-Hopf method] at Wikiwaves


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать реферат

Look at other dictionaries:

  • Eberhard Hopf — Eberhard Frederich Ferdinand Hopf (April 4, 1902 Salzburg, Austria – July 24, 1983 Bloomington, Indiana) was a mathematician and astronomer, one of the founding fathers of ergodic theory and a pioneer of bifurcation theory who also made… …   Wikipedia

  • List of mathematics articles (W) — NOTOC Wad Wadge hierarchy Wagstaff prime Wald test Wald Wolfowitz runs test Wald s equation Waldhausen category Wall Sun Sun prime Wallenius noncentral hypergeometric distribution Wallis product Wallman compactification Wallpaper group Walrasian… …   Wikipedia

  • Butcher group — In mathematics, the Butcher group, named after the New Zealand mathematician John C. Butcher by Hairer Wanner (1974), is an infinite dimensional group first introduced in numerical analysis to study solutions of non linear ordinary differential… …   Wikipedia

  • automata theory — Body of physical and logical principles underlying the operation of any electromechanical device (an automaton) that converts information input in one form into another, or into some action, according to an algorithm. Norbert Wiener and Alan M.… …   Universalium

  • List of partial differential equation topics — This is a list of partial differential equation topics, by Wikipedia page. Contents 1 General topics 2 Specific partial differential equations 3 Numerical methods for PDEs 4 …   Wikipedia

  • Гохберг, Израиль Цудикович — Израиль Цудикович Гохберг Дата рождения: 23 августа 1928(1928 08 23) …   Википедия

  • Гохберг И. — Израиль Цудикович Гохберг Дата рождения: 23 августа 1928 Место рождения: Тарутино (Бессарабия) Израиль Цудикович Гохберг (р. 23 август …   Википедия

  • Гохберг И. Ц. — Израиль Цудикович Гохберг Дата рождения: 23 августа 1928 Место рождения: Тарутино (Бессарабия) Израиль Цудикович Гохберг (р. 23 август …   Википедия

  • Гохберг, Израиль — Израиль Цудикович Гохберг Дата рождения: 23 августа 1928 Место рождения: Тарутино (Бессарабия) Израиль Цудикович Гохберг (р. 23 август …   Википедия

  • Гохберг Израиль Цудикович — Израиль Цудикович Гохберг Дата рождения: 23 августа 1928 Место рождения: Тарутино (Бессарабия) Израиль Цудикович Гохберг (р. 23 август …   Википедия

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”