рдмреЗрд╕реЗрд▓ Y рдлрд▓рди рдЯреЗрдмрд▓ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рд╣реИ?
рдпрд╣ рдЯреВрд▓ рджреНрд╡рд┐рддреАрдп рдкреНрд░рдХрд╛рд░ рдХреЗ рдмреЗрд╕реЗрд▓ рдлрд▓рди рдХреА рдЯреЗрдмрд▓ рддреИрдпрд╛рд░ рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕реЗ рд╡реЗрдмрд░ рдлрд▓рди рдпрд╛ рдиреНрдпреВрдореИрди рдлрд▓рди рднреА рдХрд╣рд╛ рдЬрд╛рддрд╛ рд╣реИ рдФрд░ рдЬрд┐рд╕реЗ \(Y_{\nu}(x)\) рдХреЗ рд░реВрдк рдореЗрдВ рд▓рд┐рдЦрд╛ рдЬрд╛рддрд╛ рд╣реИред рдпрд╣ рдмреЗрд╕реЗрд▓ рдХреЗ рдЕрд╡рдХрд▓ рд╕рдореАрдХрд░рдг (Bessel's differential equation) рдХрд╛ рджреВрд╕рд░рд╛ рд░реИрдЦрд┐рдХ рд╕реНрд╡рддрдВрддреНрд░ рд╣рд▓ рд╣реИред рдХрд┐рд╕реА рдирд┐рд╢реНрдЪрд┐рдд рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдХреЛрдЯрд┐ \(\nu\) рдХреЗ рд▓рд┐рдП, рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ \(x\) рдХреЗ рдЙрди рдорд╛рдиреЛрдВ рдХреА рд╢реНрд░реГрдВрдЦрд▓рд╛ рдкрд░ \(Y_{\nu}(x)\) рдХреА рдЧрдгрдирд╛ рдХрд░рддрд╛ рд╣реИ рдЬреЛ рдПрдХ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди, рдПрдХ рд╡реГрджреНрдзрд┐ (increment) рдФрд░ рдмрд┐рдВрджреБрдУрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ рд╕реЗ рддрдп рд╣реЛрддреЗ рд╣реИрдВ, рдФрд░ рдЗрд╕ рдкреНрд░рдХрд╛рд░ рдПрдХ рдкреВрд░реА рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рдЯреЗрдмрд▓ рдмрдирд╛ рджреЗрддрд╛ рд╣реИред
рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ
рдХреЛрдЯрд┐ \(\nu\) рджрд░реНрдЬ рдХрд░реЗрдВ (рдЬреЛ рдЕрдкреВрд░реНрдгрд╛рдВрдХ рдпрд╛ рдЛрдгрд╛рддреНрдордХ рднреА рд╣реЛ рд╕рдХрддреА рд╣реИ), \(x\) рдХрд╛ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди, рдмрд┐рдВрджреБрдУрдВ рдХреЗ рдмреАрдЪ рдХреА рд╡реГрджреНрдзрд┐ (рд╕реНрдЯреЗрдк), рдФрд░ рджреЛрд╣рд░рд╛рд╡реЛрдВ рдпрд╛рдиреА рдкрдВрдХреНрддрд┐рдпреЛрдВ (rows) рдХреА рд╕рдВрдЦреНрдпрд╛ рднрд░реЗрдВред рдХреИрд▓рдХреБрд▓реЗрдЯрд░ \(x_i = \text{startX} + i \cdot \text{stepX}\) (рдЬрд╣рд╛рдБ \(i = 0\) рд╕реЗ \(\text{pointCount} - 1\) рддрдХ) рдХреЗ рдЕрдиреБрд╕рд╛рд░ рдорд╛рди рдмрдирд╛рддрд╛ рд╣реИ рдФрд░ рд╣рд░ рдорд╛рди рдХреЗ рд▓рд┐рдП \(Y_{\nu}(x)\) рд╕реВрдЪреАрдмрджреНрдз рдХрд░рддрд╛ рд╣реИред рдзреНрдпрд╛рди рджреЗрдВ рдХрд┐ \(x = 0\) рдкрд░ \(Y_{\nu}(x)\) рдЛрдгрд╛рддреНрдордХ рдЕрдирдВрдд (тИТтИЮ) рдХреА рдУрд░ рдЕрдкрд╕рд░рд┐рдд рд╣реЛрддрд╛ рд╣реИ рдФрд░ рдпрд╣ рдХреЗрд╡рд▓ \(x > 0\) рдХреЗ рд▓рд┐рдП рд╣реА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╣реЛрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП \(x \le 0\) рд╡рд╛рд▓реА рдХрд┐рд╕реА рднреА рдкрдВрдХреНрддрд┐ рдХреЛ рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд (undefined) рдХреЗ рд░реВрдк рдореЗрдВ рдЪрд┐рд╣реНрдирд┐рдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред
рд╕реВрддреНрд░
рдЕрдкреВрд░реНрдгрд╛рдВрдХ рдХреЛрдЯрд┐ рдХреЗ рд▓рд┐рдП:
$$Y_{\nu}(x) = \frac{J_{\nu}(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}$$рдкреВрд░реНрдгрд╛рдВрдХ рдХреЛрдЯрд┐ \(n\) рдХреЗ рд▓рд┐рдП, рд╕реАрдорд╛ (limit) рд▓реЗрдиреЗ рдкрд░ рдПрдХ рдмрдВрдж рд░реВрдк (closed form) рдкреНрд░рд╛рдкреНрдд рд╣реЛрддрд╛ рд╣реИ рдЬрд┐рд╕рдореЗрдВ \(J_n(x)\cdot\ln(x/2)\) рд╡рд╛рд▓рд╛ рдПрдХ рд▓рдШреБрдЧрдгрдХреАрдп (logarithmic) рдкрдж, рдПрдХ рдкрд░рд┐рдорд┐рдд рдШрд╛рдд-рд╢реНрд░реЗрдгреА рд╕реБрдзрд╛рд░, рдФрд░ рдПрдХ рдбрд╛рдЗрдЧрд╛рдорд╛ (digamma) рд╢реНрд░реЗрдгреА рд╢рд╛рдорд┐рд▓ рд╣реЛрддреА рд╣реИред рдкреНрд░рдердо рдкреНрд░рдХрд╛рд░ рдХреЗ рдлрд▓рди \(J_{\nu}(x)\) рдХреА рдЧрдгрдирд╛ рдЙрд╕рдХреА рдШрд╛рдд-рд╢реНрд░реЗрдгреА рд╕реЗ рдпреЛрдЧ рдХрд░рдХреЗ рдХреА рдЬрд╛рддреА рд╣реИ, рдФрд░ рдЧрд╛рдорд╛ рдлрд▓рди рдХрд╛ рдорд╛рди рд▓рд╛рдВрдХрдЬрд╝реЛрд╕ (Lanczos) рд╕рдиреНрдирд┐рдХрдЯрди рд╕реЗ рдирд┐рдХрд╛рд▓рд╛ рдЬрд╛рддрд╛ рд╣реИред
рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг
рдорд╛рди рд▓реЗрдВ \(\nu = 0\), \(\text{startX} = 0\), \(\text{stepX} = 0.2\), \(\text{pointCount} = 51\) тАФ рддрдм рдкрдВрдХреНрддрд┐рдпрд╛рдБ \(x = 0.0\) рд╕реЗ \(10.0\) рддрдХ рдЪрд▓реЗрдВрдЧреАред \(Y_0(0)\) рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд рд╣реИ (тИТтИЮ), \(Y_0(0.2) \approx -1.0811\), \(Y_0(1.0) \approx 0.0883\), \(Y_0(2.0) \approx 0.5104\), рдФрд░ \(Y_0(10.0) \approx 0.0557\)ред рдореБрдЦреНрдп рдкрд░рд┐рдгрд╛рдо рдореЗрдВ "рдкрд╣рд▓рд╛ рдкрд░рд┐рдорд┐рдд рдорд╛рди" рдХреЗ рд░реВрдк рдореЗрдВ \(-1.0811\) рджрд┐рдЦрд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред
рдкрд░рд┐рднрд╛рд╖рд╛рдПрдБ рдФрд░ рд╢рдмреНрджрдХреЛрд╢
- рдХреНрд░рдо \(\nu\)
-
рдкреИрд░рд╛рдореАрдЯрд░ (the
orderрдлреАрд▓реНрдб) рдЬреЛ рдмреЗрд╕рд▓ рдлрдВрдХреНрд╢рди рдХреЗ рдкрд░рд┐рд╡рд╛рд░ рдХреЛ рдЕрдиреБрдХреНрд░рдорд┐рдд рдХрд░рддрд╛ рд╣реИред рдпрд╣ рдХреЛрдИ рднреА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╣реЛ рд╕рдХрддрд╛ рд╣реИред рдкреВрд░реНрдгрд╛рдВрдХ рдХреНрд░рдо (0, 1, 2, тАж) рдмреЗрд▓рдирд╛рдХрд╛рд░ рд╕рдорд░реВрдкрддрд╛ рд╡рд╛рд▓реА рднреМрддрд┐рдХ рд╕рдорд╕реНрдпрд╛рдУрдВ рдореЗрдВ рд╕рдмрд╕реЗ рд╕рд╛рдорд╛рдиреНрдп рд╣реИрдВ; рдЕрд░реНрдз-рдкреВрд░реНрдгрд╛рдВрдХ рдХреНрд░рдо рдЧреЛрд▓рд╛рдХрд╛рд░ рдмреЗрд╕рд▓ рдлрдВрдХреНрд╢рди рджреЗрддреЗ рд╣реИрдВред - рджреВрд╕рд░реЗ рдкреНрд░рдХрд╛рд░ рдХрд╛ рдмреЗрд╕рд▓ рдлрдВрдХреНрд╢рди \(Y_\nu(x)\)
- рд╡реЗрдмрд░ рдпрд╛ рдиреНрдпреВрдорд╛рди рдлрдВрдХреНрд╢рди (рдХрднреА-рдХрднреА \(N_\nu\) рдХреЗ рд░реВрдк рдореЗрдВ рд▓рд┐рдЦрд╛ рдЬрд╛рддрд╛ рд╣реИ) рднреА рдХрд╣рд╛ рдЬрд╛рддрд╛ рд╣реИред рдпрд╣ рдмреЗрд╕рд▓ рдХреЗ рд╕рдореАрдХрд░рдг рдХрд╛ рдПрдХ рд╕рдорд╛рдзрд╛рди рд╣реИ рдЬреЛ рдореВрд▓ рдкрд░ рдЕрд╕реАрдорд┐рдд (рд╡рд┐рд▓рдХреНрд╖рдг) рд╣реИред рдЧреИрд░-рдкреВрд░реНрдгрд╛рдВрдХ \(\nu\) рдХреЗ рд▓рд┐рдП рдкрд░рд┐рднрд╛рд╖рд┐рдд рд╣реИ: \(Y_\nu(x) = \dfrac{J_\nu(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}\), рдкреВрд░реНрдгрд╛рдВрдХ рд╕реНрдерд┐рддрд┐ рд╕реАрдорд╛ рдХреЗ рд░реВрдк рдореЗрдВ рдкреНрд░рд╛рдкреНрдд рд╣реЛрддреА рд╣реИред
- \(J_\nu\) рдмрдирд╛рдо \(Y_\nu\)
- \(J_\nu(x)\) (рдкреНрд░рдердо рдкреНрд░рдХрд╛рд░) \(x=0\) рдкрд░ рдкрд░рд┐рдорд┐рдд рд╣реИ; \(Y_\nu(x)\) (рджреНрд╡рд┐рддреАрдп рдкреНрд░рдХрд╛рд░) \(x\to 0^+\) рдХреЗ рд░реВрдк рдореЗрдВ \(-\infty\) рддрдХ рд╡рд┐рдЪрд▓рд┐рдд рд╣реЛрддрд╛ рд╣реИред рд╕рд╛рде рдореЗрдВ рд╡реЗ рдмреЗрд╕рд▓ рдХреЗ рд╕рдореАрдХрд░рдг рдХреЗ рд╕реНрд╡рддрдВрддреНрд░ рд╕рдорд╛рдзрд╛рдиреЛрдВ рдХреА рдПрдХ рд╕рдВрдкреВрд░реНрдг рдЬреЛрдбрд╝реА рдмрдирд╛рддреЗ рд╣реИрдВред
- рдмреЗрд╕рд▓ рдХрд╛ рдЕрд╡рдХрд▓ рд╕рдореАрдХрд░рдг
- рд░реИрдЦрд┐рдХ ODE \(x^2 y'' + x y' + (x^2 - \nu^2) y = 0\)ред рдЗрд╕рдХрд╛ рд╕рд╛рдорд╛рдиреНрдп рд╕рдорд╛рдзрд╛рди рд╣реИ \(y = c_1 J_\nu(x) + c_2 Y_\nu(x)\)ред
- рдЧрд╛рдорд╛ рдлрдВрдХреНрд╢рди \(\Gamma(z)\)
- рднрд╛рдЬреНрдп рдХрд╛ рд╕рддрдд рд╡рд┐рд╕реНрддрд╛рд░, \(\Gamma(n+1) = n!\), \(J_\nu\) рдФрд░ \(Y_\nu\) рдХреА рд╢реНрд░реЗрдгреА рдЧреБрдгрд╛рдВрдХреЛрдВ рдореЗрдВ рдкреНрд░рдХрдЯ рд╣реЛрддрд╛ рд╣реИред
- рдбрд┐рдЧрд╛рдорд╛ рдлрдВрдХреНрд╢рди \(\psi(z)\)
- рд▓рдШреБрдЧрдгрдХреАрдп рд╡реНрдпреБрддреНрдкрдиреНрди \(\psi(z) = \Gamma'(z)/\Gamma(z)\)ред рдпрд╣ рдкреВрд░реНрдгрд╛рдВрдХ-рдХреНрд░рдо \(Y_n(x)\) рдХреА рд╢реНрд░реЗрдгреА рдореЗрдВ рд╕реНрдкрд╖реНрдЯ рд░реВрдк рд╕реЗ рдкреНрд░рдХрдЯ рд╣реЛрддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рдПрдХ рд▓рдШреБрдЧрдгрдХреАрдп рдкрдж \(\tfrac{2}{\pi}\ln(x/2)J_n(x)\) рдФрд░ рдбрд┐рдЧрд╛рдорд╛-рднрд╛рд░рд┐рдд рдЧреБрдгрд╛рдВрдХ рд╣реЛрддреЗ рд╣реИрдВред
- рд▓реИрдирдЬрд╝реЛрд╕ рд╕рдиреНрдирд┐рдХрдЯрди
- рдЬрдЯрд┐рд▓ рдпрд╛ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рддрд░реНрдХ рдХреЗ рд▓рд┐рдП рдЧрд╛рдорд╛ рдлрдВрдХреНрд╢рди \(\Gamma(z)\) рдХреА рдЧрдгрдирд╛ рдХрд░рдиреЗ рдХреА рдПрдХ рдЕрддреНрдпрдзрд┐рдХ рд╕рдЯреАрдХ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╡рд┐рдзрд┐, рдЬрд┐рд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рд╕рд╛рдорд╛рдиреНрдпрддрдГ рдмреЗрд╕рд▓-рдлрдВрдХреНрд╢рди рджрд┐рдирдЪрд░реНрдпрд╛ рдХреЗ рдЕрдВрджрд░ рд╢реНрд░реЗрдгреА рдЧреБрдгрд╛рдВрдХреЛрдВ рдХреА рдЧрдгрдирд╛ рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдП рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред
- рд░реИрдЦрд┐рдХрддрдГ рд╕реНрд╡рддрдВрддреНрд░ рд╕рдорд╛рдзрд╛рди
- рдПрдХ рджреВрд╕рд░рд╛ рд╕рдорд╛рдзрд╛рди рдЬреЛ рдкреНрд░рдердо рдХреЗ рдПрдХ рдЕрдЪрд░ рдЧреБрдгрдЬ рдХреЗ рд░реВрдк рдореЗрдВ рд╡реНрдпрдХреНрдд рдирд╣реАрдВ рд╣реИред рдХреНрдпреЛрдВрдХрд┐ \(J_\nu\) рдЕрдХреЗрд▓реЗ рдореВрд▓ рдкрд░ рд╡рд┐рд▓рдХреНрд╖рдг рд╕рдорд╛рдзрд╛рдиреЛрдВ рдХрд╛ рдкреНрд░рддрд┐рдирд┐рдзрд┐рддреНрд╡ рдирд╣реАрдВ рдХрд░ рд╕рдХрддрд╛, \(Y_\nu\) рд╕рд╛рдорд╛рдиреНрдп рд╕рдорд╛рдзрд╛рди рдХреЗ рд▓рд┐рдП рдЖрд╡рд╢реНрдпрдХ рд╕реНрд╡рддрдВрддреНрд░ рд╕рд╛рдереА рдХреА рдЖрдкреВрд░реНрддрд┐ рдХрд░рддрд╛ рд╣реИред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рд╕рд╡рд╛рд▓
рдкрд╣рд▓реА рдкрдВрдХреНрддрд┐ рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд рдХреНрдпреЛрдВ рд╣реЛрддреА рд╣реИ? \(Y_{\nu}(x)\) рдХрд╛ \(x = 0\) рдкрд░ рдПрдХ рд╡рд┐рдЪрд┐рддреНрд░рддрд╛ (singularity) рд╣реИ рдФрд░ рдпрд╣ \(-\infty\) рдХреА рдУрд░ рдЕрдкрд╕рд░рд┐рдд рд╣реЛрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рд╡рд╣рд╛рдБ рдЗрд╕рдХрд╛ рдХреЛрдИ рдкрд░рд┐рдорд┐рдд рдорд╛рди рдирд╣реАрдВ рд╣реЛрддрд╛ред
рдХреНрдпрд╛ рдХреЛрдЯрд┐ рдЛрдгрд╛рддреНрдордХ рд╣реЛ рд╕рдХрддреА рд╣реИ? рд╣рд╛рдБред рдЛрдгрд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХ рдХреЛрдЯрд┐ рдХреЗ рд▓рд┐рдП рд╕рдордорд┐рддрд┐ \(Y_{-n}(x) = (-1)^{n}Y_{n}(x)\) рд▓рд╛рдЧреВ рд╣реЛрддреА рд╣реИ; рдЛрдгрд╛рддреНрдордХ рдЕрдкреВрд░реНрдгрд╛рдВрдХ рдХреЛрдЯрд┐ рдХреЗ рд▓рд┐рдП рд╕рд╛рдорд╛рдиреНрдп рд╕реВрддреНрд░ рд╕реАрдзреЗ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред
рдпрд╣ рдХрд┐рддрдирд╛ рд╕рдЯреАрдХ рд╣реИ? рд╢реНрд░реЗрдгрд┐рдпреЛрдВ рдХрд╛ рдпреЛрдЧ рддрдм рддрдХ рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ рдЬрдм рддрдХ рдкрдж рдорд╢реАрди рд╕рд╣рдирд╢реАрд▓рддрд╛ (machine tolerance) рд╕реЗ рдиреАрдЪреЗ рди рдЖ рдЬрд╛рдПрдВ, рдЬрд┐рд╕рд╕реЗ рдордзреНрдпрдо \(x\) рдорд╛рдиреЛрдВ рдХреЗ рд▓рд┐рдП рд▓рдЧрднрдЧ 6-7 рд╕рд╛рд░реНрдердХ рдЕрдВрдХреЛрдВ рддрдХ рдХреА рд╕рдЯреАрдХрддрд╛ рдорд┐рд▓рддреА рд╣реИред