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.
The fundamental equation that appears in the Wiener-Hopf method is of the form
A(\alpha)\Xi+(\alpha)+B(\alpha)\Psi-(\alpha)+C(\alpha)=0,
A
B
C
\Xi+(\alpha)
\Psi-(\alpha)
\tau-<ak{Im}(\alpha)<\tau+
\alpha
\Xi+(\alpha)
\Psi-(\alpha)
The key step in many Wiener–Hopf problems is to decompose an arbitrary function
\Phi
\Phi\pm
\Phi+(\alpha)=
1 | |
2\pii |
\int | |
C1 |
\Phi(z)
dz | |
z-\alpha |
and
\Phi-(\alpha)=-
1 | |
2\pii |
\int | |
C2 |
\Phi(z)
dz | |
z-\alpha |
,
where the contours
C1
C2
z=\alpha
Similarly, arbitrary scalar functions may be decomposed into a product of +/− functions, i.e.
K(\alpha)=K+(\alpha)K-(\alpha)
\boldsymbol{L}xyf(x,y)=0,
\boldsymbol{L}xy
f=g(x)forx\leq0, fy=0whenx>0
\boldsymbol{x} → infty
Taking a Fourier transform with respect to results in the following ordinary differential equation
\boldsymbol{L}y\widehat{f}(k,y)-P(k,y)\widehat{f}(k,y)=0,
\boldsymbol{L}y
infty | |
\widehat{f}(k,y)=\int | |
-infty |
f(x,y)e-ikxrm{d}x.
If a particular solution of this ordinary differential equation which satisfies the necessary decay at infinity is denoted, a general solution can be written as
\widehat{f}(k,y)=C(k)F(k,y),
The key idea is to split
\widehat{f}
\widehat{f}+
\widehat{f}-
\widehat{f}+
infty | |
(k,y)=\int | |
0 |
f(x,y)e-ikxrm{d}x,
\widehat{f}-
0 | |
(k,y)=\int | |
-infty |
f(x,y)e-ikxrm{d}x.
The boundary conditions then give
\widehat{g}(k)+\widehat{f}+(k,0)=\widehat{f}-(k,0)+\widehat{f}+(k,0)=\widehat{f}(k,0)=C(k)F(k,0)
y
\widehat{f}'-(k,0)=\widehat{f}'-(k,0)+\widehat{f}'+(k,0)=\widehat{f}'(k,0)=C(k)F'(k,0).
Eliminating
C(k)
\widehat{g}(k)+\widehat{f}+(k,0)=\widehat{f}'-(k,0)/K(k),
K(k)= | F'(k,0) |
F(k,0) |
.
Now
K(k)
K-
K+
To be precise,
K(k)=K+(k)K-(k),
logK-=
1 | |
2\pii |
infty | |
\int | |
-infty |
log(K(z)) | |
z-k |
rm{d}z, \operatorname{Im}k>0,
logK+=-
1 | |
2\pii |
infty | |
\int | |
-infty |
log(K(z)) | |
z-k |
rm{d}z, \operatorname{Im}k<0.
K
1
k → infty
K+\widehat{g}
G+
G-
K+(k)\widehat{g}(k)=G+(k)+G-(k).
This can be done in the same way that we factorised
K(k).
G+(k)+K+(k)\widehat{f}+(k,0)=\widehat{f}'-(k,0)/K-(k)-G-(k).
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 continuation 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, an application of Liouville's theorem shows that this entire function is identically zero, therefore
\widehat{f}+(k,0)=-
G+(k) | |
K+(k) |
,
C(k)=
K+(k)\widehat{g | |
(k)-G |
+(k)}{K+(k)F(k,0)}.