Struve function

From HandWiki
Graph of [math]\displaystyle{ \mathrm{H}_n(x) }[/math] for [math]\displaystyle{ n\in [0,1,2,3,4,5] }[/math]

In mathematics, the Struve functions Hα(x), are solutions y(x) of the non-homogeneous Bessel's differential equation:

[math]\displaystyle{ x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} + \left (x^2 - \alpha^2 \right )y = \frac{4\left (\frac{x}{2}\right)^{\alpha+1}}{\sqrt{\pi}\Gamma \left (\alpha+\frac{1}{2} \right )} }[/math]

introduced by Hermann Struve (1882). The complex number α is the order of the Struve function, and is often an integer.

And further defined its second-kind version [math]\displaystyle{ \mathbf{K}_\alpha(x) }[/math] as [math]\displaystyle{ \mathbf{K}_\alpha(x)=\mathbf{H}_\alpha(x)-Y_\alpha(x) }[/math].

The modified Struve functions Lα(x) are equal to ieiαπ / 2Hα(ix), are solutions y(x) of the non-homogeneous Bessel's differential equation:

Plot of the Struve function H n(z) with n=2 in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
Plot of the Struve function H n(z) with n=2 in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
[math]\displaystyle{ x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} - \left (x^2 + \alpha^2 \right )y = \frac{4\left (\frac{x}{2}\right)^{\alpha+1}}{\sqrt{\pi}\Gamma \left (\alpha+\frac{1}{2} \right )} }[/math]

And further defined its second-kind version [math]\displaystyle{ \mathbf{M}_\alpha(x) }[/math] as [math]\displaystyle{ \mathbf{M}_\alpha(x)=\mathbf{L}_\alpha(x)-I_\alpha(x) }[/math].

Definitions

Since this is a non-homogeneous equation, solutions can be constructed from a single particular solution by adding the solutions of the homogeneous problem. In this case, the homogeneous solutions are the Bessel functions, and the particular solution may be chosen as the corresponding Struve function.

Power series expansion

Struve functions, denoted as Hα(z) have the power series form

[math]\displaystyle{ \mathbf{H}_\alpha(z) = \sum_{m=0}^\infty \frac{(-1)^m}{\Gamma \left (m+\frac{3}{2} \right ) \Gamma \left (m+\alpha+\frac{3}{2} \right )} \left({\frac{z}{2}}\right)^{2m+\alpha+1}, }[/math]

where Γ(z) is the gamma function.

The modified Struve functions, denoted Lα(z), have the following power series form

[math]\displaystyle{ \mathbf{L}_\alpha(z) = \sum_{m=0}^\infty \frac{1}{\Gamma \left (m+\frac{3}{2} \right ) \Gamma \left (m+\alpha+\frac{3}{2} \right )} \left(\frac{z}{2}\right)^{2m+\alpha+1}. }[/math]
Plot of the modified Struve function L n(z) with n=2 in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
Plot of the modified Struve function L n(z) with n=2 in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D

Integral form

Another definition of the Struve function, for values of α satisfying Re(α) > − 1/2, is possible expressing in term of the Poisson's integral representation:

[math]\displaystyle{ \mathbf{H}_\alpha(x)=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^1(1-t^2)^{\alpha-\frac{1}{2}}\sin xt~dt=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\frac{\pi}{2}\sin(x\cos\tau)\sin^{2\alpha}\tau~d\tau=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\frac{\pi}{2}\sin(x\sin\tau)\cos^{2\alpha}\tau~d\tau }[/math]

[math]\displaystyle{ \mathbf{K}_\alpha(x)=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\infty(1+t^2)^{\alpha-\frac{1}{2}}e^{-xt}~dt=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\infty e^{-x\sinh\tau}\cosh^{2\alpha}\tau~d\tau }[/math]

[math]\displaystyle{ \mathbf{L}_\alpha(x)=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^1(1-t^2)^{\alpha-\frac{1}{2}}\sinh xt~dt=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\frac{\pi}{2}\sinh(x\cos\tau)\sin^{2\alpha}\tau~d\tau=\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\frac{\pi}{2}\sinh(x\sin\tau)\cos^{2\alpha}\tau~d\tau }[/math]

[math]\displaystyle{ \mathbf{M}_\alpha(x)=-\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^1(1-t^2)^{\alpha-\frac{1}{2}}e^{-xt}~dt=-\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\frac{\pi}{2}e^{-x\cos\tau}\sin^{2\alpha}\tau~d\tau=-\frac{2\left(\frac{x}{2}\right)^\alpha}{\sqrt\pi\Gamma\left(\alpha+\frac{1}{2}\right)}\int_0^\frac{\pi}{2}e^{-x\sin\tau}\cos^{2\alpha}\tau~d\tau }[/math]

Asymptotic forms

For small x, the power series expansion is given above.

For large x, one obtains:

[math]\displaystyle{ \mathbf{H}_\alpha(x) - Y_\alpha(x) = \frac{\left(\frac{x}{2}\right)^{\alpha-1}}{\sqrt{\pi} \Gamma \left (\alpha+\frac{1}{2} \right )} + O\left(\left (\tfrac{x}{2}\right)^{\alpha-3}\right), }[/math]

where Yα(x) is the Neumann function.

Properties

The Struve functions satisfy the following recurrence relations:

[math]\displaystyle{ \begin{align} \mathbf{H}_{\alpha -1}(x) + \mathbf{H}_{\alpha+1}(x) &= \frac{2\alpha}{x} \mathbf{H}_\alpha (x) + \frac{\left (\frac{x}{2}\right)^{\alpha}}{\sqrt{\pi}\Gamma \left (\alpha + \frac{3}{2} \right )}, \\ \mathbf{H}_{\alpha -1}(x) - \mathbf{H}_{\alpha+1}(x) &= 2 \frac{d}{dx} \left (\mathbf{H}_\alpha(x) \right) - \frac{ \left( \frac{x}{2} \right)^\alpha}{\sqrt{\pi}\Gamma \left (\alpha + \frac{3}{2} \right )}. \end{align} }[/math]

Relation to other functions

Struve functions of integer order can be expressed in terms of Weber functions En and vice versa: if n is a non-negative integer then

[math]\displaystyle{ \begin{align} \mathbf{E}_n(z) &= \frac{1}{\pi} \sum_{k=0}^{\left \lfloor \frac{n-1}{2} \right \rfloor} \frac{\Gamma \left (k+ \frac{1}{2} \right) \left (\frac{z}{2} \right )^{n-2k-1}}{\Gamma \left (n- k + \frac{1}{2}\right )} -\mathbf{H}_n(z),\\ \mathbf{E}_{-n}(z) &= \frac{(-1)^{n+1}}{\pi}\sum_{k=0}^{\left \lfloor \frac{n-1}{2} \right \rfloor} \frac{\Gamma(n-k-\frac{1}{2}) \left (\frac{z}{2} \right )^{-n+2k+1}}{\Gamma \left (k+ \frac{3}{2} \right)}-\mathbf{H}_{-n}(z). \end{align} }[/math]

Struve functions of order n + 1/2 where n is an integer can be expressed in terms of elementary functions. In particular if n is a non-negative integer then

[math]\displaystyle{ \mathbf{H}_{-n-\frac{1}{2}} (z) = (-1)^n J_{n+\frac{1}{2}}(z), }[/math]

where the right hand side is a spherical Bessel function.

Struve functions (of any order) can be expressed in terms of the generalized hypergeometric function 1F2:

[math]\displaystyle{ \mathbf{H}_{\alpha}(z) = \frac{z^{\alpha+1}}{2^{\alpha}\sqrt{\pi} \Gamma \left (\alpha+\tfrac{3}{2} \right )} {}_1F_2 \left (1;\tfrac{3}{2}, \alpha+\tfrac{3}{2};-\tfrac{z^2}{4} \right ). }[/math]

Applications

The Struve and Weber functions were shown to have an application to beamforming in.[1], and in describing the effect of confining interface on Brownian motion of colloidal particles at low Reynolds numbers.[2]

References

  1. K. Buchanan, C. Flores, S. Wheeland, J. Jensen, D. Grayson and G. Huff, "Transmit beamforming for radar applications using circularly tapered random arrays," 2017 IEEE Radar Conference (RadarConf), 2017, pp. 0112-0117, doi: 10.1109/RADAR.2017.7944181
  2. B. U. Felderhof, "Effect of the wall on the velocity autocorrelation function and long-time tail of Brownian motion." The Journal of Physical Chemistry B 109.45, 2005, pp. 21406-21412

External links