Struve-functie

In de wiskunde is de Struve-functie een speciale functie die in 1882 werd geïntroduceerd door de astronoom Hermann Struve tijdens zijn theoretisch onderzoek van diffractieverschijnselen[1] in de optica. De functie heeft inmiddels toepassingen gevonden in de wiskunde, de optica, de hydrodynamica en de akoestiek. De functie wordt meestal voorgesteld door H ν ( z ) {\displaystyle H_{\nu }(z)} waarin ν {\displaystyle \nu } de orde aangeeft. De Struve-functie beschrijft oplossingen van de Besselse differentiaalvergelijking.

Definitie

Verloop van de functie H ν ( z ) {\displaystyle H_{\nu }(z)} voor ν = 0 , 1 , 2 , 3 , 4 , 5 {\displaystyle \nu =0,1,2,3,4,5}

De Struve-functie H ν ( z ) {\displaystyle H_{\nu }(z)} van de eerste soort is een particuliere oplossing van de volgende inhomogene besselse differentiaalvergelijking met speciaal tweede lid

z 2 d 2 y d z 2 + z d y d z + ( z 2 ν 2 ) y = 4 ( z 2 ) ν + 1 π Γ ( ν + 1 2 ) {\displaystyle z^{2}{\frac {{\rm {d}}^{2}y}{{\rm {d}}z^{2}}}+z{\frac {{\rm {d}}y}{{\rm {d}}z}}+\left(z^{2}-\nu ^{2}\right)y={\frac {4\left({\frac {z}{2}}\right)^{\nu +1}}{{\sqrt {\pi }}\,\Gamma \left(\nu +{\frac {1}{2}}\right)}}}

waarin Γ {\displaystyle \Gamma } de Gammafunctie voorstelt. Het complexe getal ν {\displaystyle \nu } geeft de orde aan van de Struve-functie en is meestal een geheel getal.

Reeksontwikkeling

Struve-functies, aangeduid als H ν ( z ) {\displaystyle H_{\nu }(z)} , kunnen weergegeven worden door de volgende machtreeksen.

H ν ( z ) = k = 0 ( 1 ) k Γ ( k + 3 2 ) Γ ( k + ν + 3 2 ) ( z 2 ) 2 k + ν + 1 {\displaystyle H_{\nu }(z)=\sum _{k=0}^{\infty }{\frac {(-1)^{k}}{\Gamma \left(k+{\frac {3}{2}}\right)\Gamma \left(k+\nu +{\frac {3}{2}}\right)}}\left({\frac {z}{2}}\right)^{2k+\nu +1}}

waarin Γ ( x ) {\displaystyle \Gamma (x)} de Gammafunctie is.

H 0 ( z ) = 2 π ( z z 3 1 2 3 2 + z 5 1 2 3 2 5 2 ) {\displaystyle H_{0}(z)={\frac {2}{\pi }}\left(z-{\frac {z^{3}}{1^{2}\cdot 3^{2}}}+{\frac {z^{5}}{1^{2}\cdot 3^{2}\cdot 5^{2}}}-\ldots \right)}
H 1 ( z ) = 2 π ( z 2 1 2 3 z 4 1 2 3 2 5 + z 6 1 2 3 2 5 2 7 ) {\displaystyle H_{1}(z)={\frac {2}{\pi }}\left({\frac {z^{2}}{1^{2}\cdot 3}}-{\frac {z^{4}}{1^{2}\cdot 3^{2}\cdot 5}}+{\frac {z^{6}}{1^{2}\cdot 3^{2}\cdot 5^{2}\cdot 7}}-\ldots \right)}

Integraalvoorstelling

Struve-functie kunnen ook door een integraal voorgesteld worden voor ( ν ) > 1 / 2 {\displaystyle \Re (\nu )>-1/2}

H ν ( z ) = 2 ( z 2 ) ν π Γ ( ν + 1 2 ) 0 π 2 sin ( z cos θ ) sin 2 ν ( θ ) d θ . {\displaystyle H_{\nu }(z)={\frac {2\left({\frac {z}{2}}\right)^{\nu }}{{\sqrt {\pi }}\,\Gamma \left(\nu +{\frac {1}{2}}\right)}}\int _{0}^{\frac {\pi }{2}}\sin(z\cos \theta )\sin ^{2\nu }(\theta ){\rm {d}}\theta .}

Een andere voorstelling krijgt men door de substitutie t = cos ( θ ) {\displaystyle t=\cos(\theta )} .

H ν ( z ) = 2 ( z 2 ) ν π Γ ( ν + 1 2 ) 0 1 ( 1 t 2 ) ν 1 2 sin ( z t ) d t {\displaystyle H_{\nu }(z)={\frac {2\left({\frac {z}{2}}\right)^{\nu }}{{\sqrt {\pi }}\,\Gamma \left(\nu +{\frac {1}{2}}\right)}}\int _{0}^{1}(1-t^{2})^{\nu -{\frac {1}{2}}}\sin \left(zt\right){\rm {d}}t}

Recursieformules

De Struve functie voldoet aan de volgende recursierelaties:

H ν 1 ( z ) + H ν + 1 ( z ) = 2 ν z H ν ( z ) + ( z 2 ) ν π Γ ( ν + 3 2 ) , H ν 1 ( z ) H ν + 1 ( z ) = 2 d d z H ν ( z ) ( z 2 ) ν π Γ ( ν + 3 2 ) , {\displaystyle {\begin{aligned}H_{\nu -1}(z)+H_{\nu +1}(z)&={\frac {2\nu }{z}}H_{\nu }(z)+{\frac {\left({\frac {z}{2}}\right)^{\nu }}{{\sqrt {\pi }}\,\Gamma \left(\nu +{\frac {3}{2}}\right)}},\\H_{\nu -1}(z)-H_{\nu +1}(z)&=2{\frac {\rm {d}}{{\rm {d}}z}}H_{\nu }(z)-{\frac {\left({\frac {z}{2}}\right)^{\nu }}{{\sqrt {\pi }}\,\Gamma \left(\nu +{\frac {3}{2}}\right)}},\end{aligned}}}
d H 0 d z = 1 2 π H 1 {\displaystyle {\frac {{\rm {d}}H_{0}}{{\rm {d}}z}}={\tfrac {1}{2}}\pi -H_{1}}
d d z ( z ν H ν ) = z ν H ν 1 {\displaystyle {\frac {\rm {d}}{{\rm {d}}z}}\left(z^{\nu }H_{\nu }\right)=z^{\nu }H_{\nu -1}}
d d z ( z ν H ν ) = 1 π 2 ν Γ ( ν + 3 2 ) z ν H ν + 1 ( z ) {\displaystyle {\frac {\rm {d}}{{\rm {d}}z}}\left(z^{-\nu }H_{\nu }\right)={\frac {1}{{\sqrt {\pi }}\,2^{\nu }\Gamma \left(\nu +{\frac {3}{2}}\right)}}-z^{-\nu }H_{\nu +1}(z)}

Integralen

Uit de vierde recursieformule volgt onmiddellijk de integraal:

0 a z H 0 ( z ) d z = a H 1 ( a ) {\displaystyle \int _{0}^{a}z\,H_{0}(z)\,{\rm {d}}z=aH_{1}(a)} .
0 t 1 H 0 ( t ) d t = 1 2 π {\displaystyle \int _{0}^{\infty }t^{-1}H_{0}(t)\,{\rm {d}}t={\tfrac {1}{2}}\pi }

De volgende integraal wordt ook wel de Struve-integraal genoemd:

4 π z t 2 H 1 ( t ) d t = 2 π z H 1 ( z ) + 2 π z t 1 H 0 ( t ) d t {\displaystyle {\frac {4}{\pi }}\int _{z}^{\infty }t^{-2}\,H_{1}(t)\,{\rm {d}}t={\frac {2}{\pi z}}H_{1}(z)+{\frac {2}{\pi }}\int _{z}^{\infty }t^{-1}H_{0}(t)\,{\rm {d}}t}
2 π z t 1 H 0 ( t ) d t = 1 4 π 2 ( z z 3 1 2 3 2 3 + z 5 1 2 3 2 5 2 5 ) {\displaystyle {\frac {2}{\pi }}\int _{z}^{\infty }t^{-1}H_{0}(t)\,{\rm {d}}t=1-{\frac {4}{\pi ^{2}}}\left(z-{\frac {z^{3}}{1^{2}\cdot 3^{2}\cdot 3}}+{\frac {z^{5}}{1^{2}\cdot 3^{2}\cdot 5^{2}\cdot 5}}-\ldots \right)}

Asymptotische ontwikkeling

Voor kleine z {\displaystyle z} gebruikt men de gegeven bovenstaande machtreeksontwikkeling.

Voor grote z {\displaystyle z} verkrijgt men:

H ν ( z ) Y ν ( z ) ( z 2 ) ν 1 π Γ ( ν + 1 2 ) + O ( ( z 2 ) ν 3 ) , {\displaystyle H_{\nu }(z)-Y_{\nu }(z)\to {\frac {\left({\frac {z}{2}}\right)^{\nu -1}}{{\sqrt {\pi }}\,\Gamma \left(\nu +{\frac {1}{2}}\right)}}+O\left(\left({\tfrac {z}{2}}\right)^{\nu -3}\right),}

waarin Y ν ( z ) {\displaystyle Y_{\nu }(z)} de besselfunctie van de tweede soort van de orde ν {\displaystyle \nu } is. Alleen de eerste term van deze expansie is weergegeven.

Benaderende formules

In de wetenschappelijke literatuur vindt men vele benaderingsformules voor de Bessel- en de Struve-functies. De meeste daarvan splitsen het gebied van z {\displaystyle z} op in een gebied waar z {\displaystyle z} klein is en een gebied waar z {\displaystyle z} groot is. Zo publiceerde J. Newman[2] nauwkeurige veeltermbenaderingen voor z < 3 {\displaystyle z<3} en voor z > 3 {\displaystyle z>3} . Een effectieve benadering[3] voor H 1 {\displaystyle H_{1}} , die geldig is voor alle waarden van z 0 {\displaystyle z\geq 0} met een maximale absolute fout van 0,0049, wordt gegeven door de volgende uitdrukking:

H 1 ( z ) 2 π J 0 ( z ) + ( 16 π 5 ) sin z z + ( 12 36 π ) 1 cos z z 2 {\displaystyle H_{1}(z)\approx {\frac {2}{\pi }}-J_{0}(z)+\left({\frac {16}{\pi }}-5\right){\frac {\sin z}{z}}+\left(12-{\frac {36}{\pi }}\right){\frac {1-\cos z}{z^{2}}}}

waarin J 0 ( z ) {\displaystyle J_{0}(z)} de besselfunctie van de eerste soort en orde 0 is. Met behulp van de vierde recursieformule kan men dan een benadering voor H 0 {\displaystyle H_{0}} krijgen met een maximale absolute fout van 0,0063. Door gebruik te maken van een soortgelijke maar verbeterde methode kon de nauwkeurigheid voor H 0 {\displaystyle H_{0}} en H 1 {\displaystyle H_{1}} opgevoerd worden tot een maximale absolute fout van respectievelijk 0,00125 en 0,00185.[4] De verbeterde nauwkeurigheid opent de weg naar nauwkeuriger benaderingen voor de hogere orde Struve-functies H ν ( z ) {\displaystyle H_{\nu }(z)} met ν = 2 , 3 , {\displaystyle \nu =2,3,\ldots } en dit met behulp van de eerste recursieformule.

Literatuur

Milton Abramowitz & Irene Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover Publications, New York; isbn 978-0-486-61272-0 1972.

  • Struve functies op the Wolfram functies website.
  • Struve functions and related functions
Bronnen, noten en/of referenties
  1. Struve, H., Beitrag zur Theorie der Diffraction an Fernröhren Ann. Physik Chemie 17 (13): 1008-1016, 1882.
  2. J. Newman, Approximations for the Bessel and Struve functions, Mathematics of Computing, vol.43, nr. 168, pp. 551-556, oktober 1984
  3. *Ronald M. Aarts en Augustus J.E.M. Janssen, Approximation of the Struve function H 1 {\displaystyle H_{1}} occurring in impedance calculations, J. Acoust. Soc. Am.,vol. 113 , nr.5 ,pp. 2635–2637, mei 2003
  4. *Ronald M. Aarts en Augustus J.E.M. Janssen, Efficient approximation of the Struve functions H n {\displaystyle H_{n}} occurring in the calculation of sound radiation quantities, J. Acoust. Soc. Am.,vol. 140 , nr.6 ,pp. 4154-4160, december 2016