Struve function

Graph of ${\mathrm{H}}_n(x)$ for $n\in [0,1,2,3,4,5]$

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

$x^2\frac{d^{2}y}{d{x}^2}+x\frac{dy}{dx}+(x^2-{\alpha }^2)y=\frac{4{\left(\frac{x}{2}\right)}^{\alpha +1}}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}$

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 ${\mathbf{K}}_{\alpha }(x)$ as ${\mathbf{K}}_{\alpha }(x)={\mathbf{H}}_{\alpha }(x)-Y_{\alpha }(x)$.

The modified Struve functions L_α(x) are equal to −ie^{−iαπ / 2}H_α(ix), are solutions y(x) of the non-homogeneous Bessel's differential equation:

$x^2\frac{d^{2}y}{d{x}^2}+x\frac{dy}{dx}-(x^2+{\alpha }^2)y=\frac{4{\left(\frac{x}{2}\right)}^{\alpha +1}}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}$

And further defined its second-kind version ${\mathbf{M}}_{\alpha }(x)$ as ${\mathbf{M}}_{\alpha }(x)={\mathbf{L}}_{\alpha }(x)-I_{\alpha }(x)$.

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

${\mathbf{H}}_{\alpha }(z)=\sum _{m=0}^{\mathrm{\infty }}\frac{(-1{)}^m}{\mathrm{\Gamma }(m+\frac{3}{2})\mathrm{\Gamma }(m+\alpha +\frac{3}{2})}{\left(\frac{z}{2}\right)}^{2m+\alpha +1},$

where Γ(z) is the gamma function.

The modified Struve functions, denoted L_ν(z), have the following power series form

${\mathbf{L}}_{\nu }(z)={\left(\frac{z}{2}\right)}^{\nu +1}\sum _{k=0}^{\mathrm{\infty }}\frac{1}{\mathrm{\Gamma }(\frac{3}{2}+k)\mathrm{\Gamma }(\frac{3}{2}+k+\nu )}{\left(\frac{z}{2}\right)}^{2k}.$

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:

${\mathbf{H}}_{\alpha }(x)=\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{1}{\int }}(1-t^2{)}^{\alpha -\frac{1}{2}}\sin xtdt=\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\frac{\pi }{2}}{\int }}\sin (x\cos \tau ){\sin }^{}\tau d\tau =\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\frac{\pi }{2}}{\int }}\sin (x\sin \tau ){\cos }^{}\tau d\tau$

${\mathbf{K}}_{\alpha }(x)=\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\mathrm{\infty }}{\int }}(1+t^2{)}^{\alpha -\frac{1}{2}}{e}^{-xt}dt=\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\mathrm{\infty }}{\int }}{e}^{-x\sinh \tau }{\cosh }^{2\alpha }\tau d\tau$

${\mathbf{L}}_{\alpha }(x)=\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{1}{\int }}(1-t^2{)}^{\alpha -\frac{1}{2}}\sinh xtdt=\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\frac{\pi }{2}}{\int }}\sinh (x\cos \tau ){\sin }^{}\tau d\tau =\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\frac{\pi }{2}}{\int }}\sinh (x\sin \tau ){\cos }^{}\tau d\tau$

${\mathbf{M}}_{\alpha }(x)=-\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{1}{\int }}(1-t^2{)}^{\alpha -\frac{1}{2}}{e}^{-xt}dt=-\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\frac{\pi }{2}}{\int }}{e}^{-x\cos \tau }{\sin }^{2\alpha }\tau d\tau =-\frac{2{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}\underset{0}{\overset{\frac{\pi }{2}}{\int }}{e}^{-x\sin \tau }{\cos }^{2\alpha }\tau d\tau$

Asymptotic forms

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

For large x, one obtains:

${\mathbf{H}}_{\alpha }(x)-Y_{\alpha }(x)=\frac{{\left(\frac{x}{2}\right)}^{\alpha -1}}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{1}{2})}+O\left({\left(\frac{x}{2}\right)}^{\alpha -3}\right),$

where Y_α(x) is the Neumann function.

Properties

The Struve functions satisfy the following recurrence relations:

$\begin{array}{rl}{\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 }\mathrm{\Gamma }(\alpha +\frac{3}{2})},\\ {\mathbf{H}}_{\alpha -1}(x)-{\mathbf{H}}_{\alpha +1}(x)& =2\frac{d}{dx}({\mathbf{H}}_{\alpha }(x))-\frac{{\left(\frac{x}{2}\right)}^{\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{3}{2})}.\end{array}$

Relation to other functions

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

$\begin{array}{rl}{\mathbf{E}}_n(z)& =\frac{1}{\pi }\sum _{k=0}^{\lfloor \frac{n-1}{2}\rfloor }\frac{\mathrm{\Gamma }(k+\frac{1}{2}){\left(\frac{z}{2}\right)}^{n-2k-1}}{\mathrm{\Gamma }(n-k+\frac{1}{2})}-{\mathbf{H}}_n(z),\\ {\mathbf{E}}_{-n}(z)& =\frac{(-1{)}^{n+1}}{\pi }\sum _{k=0}^{\lfloor \frac{n-1}{2}\rfloor }\frac{\mathrm{\Gamma }(n-k-\frac{1}{2}){\left(\frac{z}{2}\right)}^{-n+2k+1}}{\mathrm{\Gamma }(k+\frac{3}{2})}-{\mathbf{H}}_{-n}(z).\end{array}$

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

${\mathbf{H}}_{-n-\frac{1}{2}}(z)=(-1{)}^n{J}_{n+\frac{1}{2}}(z),$

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 ₁F₂ (which is not the Gauss hypergeometric function ₂F₁):

${\mathbf{H}}_{\alpha }(z)=\frac{z^{\alpha +1}}{2^{\alpha }\sqrt{\pi }\mathrm{\Gamma }(\alpha +\frac{3}{2})}{}_1{F}_2(1,\frac{3}{2},\alpha +\frac{3}{2},-\frac{z^2}{4}).$

References

• R. M. Aarts and Augustus J. E. M. Janssen (2003). "Approximation of the Struve function H₁ occurring in impedance calculations". J. Acoust. Soc. Am. 113 (5): 2635–2637. Bibcode:2003ASAJ..113.2635A. doi:10.1121/1.1564019. PMID 12765381.
• R. M. Aarts and Augustus J. E. M. Janssen (2016). "Efficient approximation of the Struve functions H_n occurring in the calculation of sound radiation quantities". J. Acoust. Soc. Am. 140 (6): 4154–4160. Bibcode:2016ASAJ..140.4154A. doi:10.1121/1.4968792. PMID 28040027.
• Abramowitz, Milton; Stegun, Irene Ann, eds. (1983) [June 1964]. "Chapter 12". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. p. 496. ISBN 978-0-486-61272-0. LCCN 64-60036. MR 0167642. LCCN 65-12253.
• Ivanov, A. B. (2001) [1994], "Struve function", Encyclopedia of Mathematics, EMS Press
• Paris, R. B. (2010), "Struve function", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248
• Struve, H. (1882). "Beitrag zur Theorie der Diffraction an Fernröhren". Annalen der Physik und Chemie. 17 (13): 1008–1016. Bibcode:1882AnP...253.1008S. doi:10.1002/andp.18822531319.

External links

• Struve functions at the Wolfram functions site.