2025-11-19T12:16:13.706625

Uniform Asymptotic approximation and numerical evaluation of the Reverse Generalized Bessel Polynomial zeros

Dunster, Gil, Ruiz-Antolin et al.
Uniform asymptotic expansions are derived for the zeros of the reverse generalized Bessel polynomials of large degree $n$ and real parameter $a$. It is assumed that $-Δ_{1} n+\frac{3}{2} \leq a \leq Δ_{2} n$ for fixed arbitrary $Δ_{1} \in (0,1)$ and bounded positive $Δ_{2}$. For this parameter range at most one of the zeros is real, with the rest being complex conjugates. The new expansions are uniformly valid for all the zeros, and are shown to be highly accurate for moderate or large values of $n$. They are consequently used as initial values in a very efficient numerical algorithm designed to obtain the remaining complex zeros using Taylor series.
academic

रिवर्स जनरलाइज्ड बेसेल पॉलीनोमियल शून्यों का एकसमान अनंतस्पर्शी सन्निकटन और संख्यात्मक मूल्यांकन

मूल जानकारी

  • पेपर ID: 2510.14295
  • शीर्षक: रिवर्स जनरलाइज्ड बेसेल पॉलीनोमियल शून्यों का एकसमान अनंतस्पर्शी सन्निकटन और संख्यात्मक मूल्यांकन
  • लेखक: T. M. Dunster, A. Gil, D. Ruiz-Antolin, J. Segura
  • वर्गीकरण: math.CA (शास्त्रीय विश्लेषण और ODEs)
  • प्रकाशन तिथि: 16 अक्टूबर 2025
  • पेपर लिंक: https://arxiv.org/abs/2510.14295

सारांश

यह पेपर बड़ी घात nn और वास्तविक पैरामीटर aa के लिए रिवर्स जनरलाइज्ड बेसेल पॉलीनोमियल शून्यों के एकसमान अनंतस्पर्शी विस्तार को प्राप्त करता है। पैरामीटर श्रेणी Δ1n+32aΔ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n मानते हुए, जहाँ Δ1(0,1)\Delta_1 \in (0,1) एक निश्चित मनमाना स्थिरांक है और Δ2\Delta_2 एक परिबद्ध धनात्मक स्थिरांक है। इस पैरामीटर श्रेणी में, अधिकतम एक शून्य वास्तविक है, शेष सभी जटिल संयुग्म युग्म हैं। नए विस्तार सभी शून्यों के लिए एकसमान रूप से प्रभावी हैं, nn के मध्यम या बड़े मानों पर उच्च सटीकता प्रदर्शित करते हैं। इन विस्तारों का उपयोग एक कुशल संख्यात्मक एल्गोरिथ्म के लिए प्रारंभिक मानों के रूप में किया जाता है, जो शेष जटिल शून्यों की गणना के लिए टेलर श्रृंखला का उपयोग करता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या विवरण

जनरलाइज्ड बेसेल पॉलीनोमियल को इस प्रकार परिभाषित किया जाता है: yn(z;a)=k=0n(nk)(n+a1)k(12z)ky_n(z; a) = \sum_{k=0}^{n} \binom{n}{k} (n+a-1)_k \left(\frac{1}{2z}\right)^k

जहाँ (α)k=Γ(α+k)/Γ(α)(α)_k = Γ(α + k)/Γ(α) पोचहैमर प्रतीक है। इन पॉलीनोमियल के शून्य आमतौर पर जटिल-मूल्यवान होते हैं, और अनुप्रयुक्त गणित और इंजीनियरिंग में विविध अनुप्रयोग होते हैं।

अनुसंधान प्रेरणा

  1. व्यावहारिक अनुप्रयोग की आवश्यकता: जनरलाइज्ड बेसेल पॉलीनोमियल के शून्य अनुप्रयुक्त गणित और इंजीनियरिंग क्षेत्र में महत्वपूर्ण अनुप्रयोग हैं
  2. कम्प्यूटेशनल चुनौती: इन शून्यों की गणना के लिए मौजूदा तरीके दक्षता और सटीकता में सीमित हैं
  3. सैद्धांतिक पूर्णता: सभी शून्यों के व्यवहार का वर्णन करने के लिए एकसमान रूप से प्रभावी अनंतस्पर्शी सिद्धांत स्थापित करने की आवश्यकता है
  4. संख्यात्मक एल्गोरिथ्म अनुकूलन: जटिल शून्यों की गणना के लिए कुशल संख्यात्मक तरीकों की आवश्यकता है

मौजूदा तरीकों की सीमाएँ

  • पूर्ववर्ती कार्य 2, 13 में कम्प्यूटेशनल सटीकता और दक्षता में सीमाएँ हैं
  • सभी शून्यों के लिए एकसमान रूप से प्रभावी अनंतस्पर्शी विस्तार का अभाव है
  • संख्यात्मक एल्गोरिथ्म के प्रारंभिक मानों का चयन पर्याप्त रूप से अनुकूलित नहीं है

मुख्य योगदान

  1. एकसमान अनंतस्पर्शी विस्तार प्राप्त किया: रिवर्स जनरलाइज्ड बेसेल पॉलीनोमियल के शून्यों के लिए सभी शून्यों के लिए एकसमान रूप से प्रभावी अनंतस्पर्शी विस्तार स्थापित किया
  2. टर्निंग पॉइंट सिद्धांत स्थापित किया: WKB सिद्धांत और Liouville-Green विधि के आधार पर, टर्निंग पॉइंट वाले अवकल समीकरणों को संभाला
  3. कुशल संख्यात्मक एल्गोरिथ्म का निर्माण किया: अनंतस्पर्शी विस्तार और टेलर श्रृंखला को जोड़कर, जटिल शून्यों की गणना के लिए एक कुशल एल्गोरिथ्म विकसित किया
  4. उच्च सटीकता सन्निकटन प्रदान किया: विस्तार nn के मध्यम या बड़े मानों पर 101510^{-15} परिमाण स्तर की सापेक्ष त्रुटि प्राप्त करता है

विधि विवरण

कार्य परिभाषा

रिवर्स जनरलाइज्ड बेसेल पॉलीनोमियल θn(z;a)=znyn(z1;a)θ_n(z; a) = z^n y_n(z^{-1}; a) के शून्य tm(u,a)t_m(u, a) खोजना, जहाँ u=n+12u = n + \frac{1}{2}, पैरामीटर श्रेणी Δ1n+32aΔ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n के तहत।

मुख्य सैद्धांतिक ढाँचा

1. अवकल समीकरण परिवर्तन

स्केलिंग फ़ंक्शन के माध्यम से: wn(0)(z;a)=2na+1z1na/2ezθn(z;a)w_n^{(0)}(z; a) = 2^{-n-a+1} z^{1-n-a/2} e^{-z} θ_n(z; a)

समस्या को अवकल समीकरण को हल करने में परिवर्तित किया जाता है: d2wdz2={1+a2z+(2n+a)(2n+a2)4z2}w\frac{d^2w}{dz^2} = \left\{1 + \frac{a-2}{z} + \frac{(2n+a)(2n+a-2)}{4z^2}\right\} w

2. Liouville-Green सिद्धांत

पैरामीटर परिभाषित करें: u=n+12,α=a2uu = n + \frac{1}{2}, \quad α = \frac{a-2}{u}

अवकल समीकरण को इस प्रकार पुनः लिखा जाता है: d2wdz2={u2f(α,z)+g(z)}w\frac{d^2w}{dz^2} = \{u^2 f(α, z) + g(z)\} w

जहाँ: f(α,z)=(z+12α)2+1+αz2,g(z)=14z2f(α, z) = \frac{(z + \frac{1}{2}α)^2 + 1 + α}{z^2}, \quad g(z) = -\frac{1}{4z^2}

3. टर्निंग पॉइंट विश्लेषण

टर्निंग पॉइंट स्थित हैं: z1,2(α)=±iσ12α,σ=1+αz_{1,2}(α) = ±iσ - \frac{1}{2}α, \quad σ = \sqrt{1 + α}

Liouville चर ζζ और LG चर ξξ के माध्यम से: 23ζ3/2=ξ=z1(α)zf1/2(α,t)dt\frac{2}{3}ζ^{3/2} = ξ = \int_{z_1(α)}^z f^{1/2}(α, t) dt

शून्यों का अनंतस्पर्शी विस्तार

शून्य निहित समीकरण को संतुष्ट करते हैं: Z(u,a,u1tm(u,a))=u2/3amZ(u, a, u^{-1}t_m(u, a)) = u^{-2/3}a_m

जहाँ ama_m Airy फ़ंक्शन Ai(x)Ai(x) का mm-वाँ ऋणात्मक शून्य है।

विस्तार रूप है: tm(u,a)us=0τm,s(α)u2s(u)t_m(u, a) \sim u \sum_{s=0}^{\infty} \frac{τ_{m,s}(α)}{u^{2s}} \quad (u → ∞)

मुख्य गुणांक गणना

  • अग्रणी पद τm,0τ_{m,0}: जटिल निहित समीकरण (3.16) द्वारा निर्धारित
  • बाद के पद: पुनरावर्ती संबंध (3.20)-(3.23) के माध्यम से गणना की जाती है:

τm,1=Υm,1ζm,0τ_{m,1} = -\frac{Υ_{m,1}}{ζ'_{m,0}}

τm,2=12ζm,0{τm,12ζm,0+2τm,1Υm,1+2Υm,2}τ_{m,2} = -\frac{1}{2ζ'_{m,0}}\{τ_{m,1}^2 ζ''_{m,0} + 2τ_{m,1} Υ'_{m,1} + 2Υ_{m,2}\}

संख्यात्मक एल्गोरिथ्म डिज़ाइन

1. पुनरावृत्ति फ़ंक्शन

Tn(a,z)=z1Ωn(a,z)arctan(Ωn(a,z)wn(0)(z;a)wn(0)(z;a)/z)T_n(a, z) = z - \frac{1}{\sqrt{Ω_n(a, z)}} \arctan\left(\frac{\sqrt{Ω_n(a, z)}w_n^{(0)}(z; a)}{\partial w_n^{(0)}(z; a)/\partial z}\right)

2. स्टेप साइज़ फ़ंक्शन

Hn(a,z)=z+πΩn(a,z)H_n(a, z) = z + \frac{π}{\sqrt{Ω_n(a, z)}}

3. टेलर श्रृंखला विस्तार

पुनरावर्ती संबंध (4.7) का उपयोग करके उच्च-क्रम व्युत्पन्न की गणना करें, फ़ंक्शन मानों की सीधी गणना से बचें।

प्रायोगिक सेटअप

परीक्षण पैरामीटर

  • घात श्रेणी: n=15,30,50,200,500,1000,2000n = 15, 30, 50, 200, 500, 1000, 2000
  • पैरामीटर मान: a=1.01,20.2a = 1.01, 20.2 आदि
  • शून्य सूचकांक: m=1,3,10,15,25m = 1, 3, 10, 15, 25 आदि

मूल्यांकन संकेतक

  • सापेक्ष त्रुटि: Maple उच्च-सटीकता कार्यान्वयन के साथ तुलना
  • कम्प्यूटेशन समय: CPU निष्पादन समय
  • अभिसरण: पुनरावृत्ति एल्गोरिथ्म की अभिसरण सटीकता ε=1012ε = 10^{-12} पर सेट है

कार्यान्वयन विवरण

  • अनंतस्पर्शी विस्तार पहले 5 पद गुणांकों का उपयोग करता है
  • Maple कार्यान्वयन 60-बिट सटीकता का उपयोग करता है
  • Matlab एल्गोरिथ्म Dell Latitude 7410 पर परीक्षण किया गया

प्रायोगिक परिणाम

मुख्य परिणाम

1. उच्च सटीकता सन्निकटन

तालिका 1 और 2 द्वारा दर्शित सापेक्ष त्रुटि:

  • a=1.01a = 1.01 के लिए: सापेक्ष त्रुटि 101810^{-18} से 101510^{-15} परिमाण में है
  • a=20.2a = 20.2 के लिए: सापेक्ष त्रुटि 101810^{-18} से 101510^{-15} परिमाण में है
  • सभी परीक्षण मामलों में सापेक्ष त्रुटि 101510^{-15} के निकट या बेहतर है

2. पैरामीटर निर्भरता विश्लेषण

चित्र 1 दर्शाता है कि पैरामीटर aa बढ़ने के साथ सापेक्ष त्रुटि घटती है, अधिकतम सापेक्ष त्रुटि 4×10154×10^{-15} से कम है।

3. एल्गोरिथ्म दक्षता

तालिका 3 के कम्प्यूटेशन समय परिणाम:

  • n=30n = 30 (15 शून्य): 2.8×1032.8×10^{-3} सेकंड
  • n=2000n = 2000 (1000 शून्य): 3.7×1023.7×10^{-2} सेकंड
  • शून्य संख्या में 67 गुना वृद्धि, कम्प्यूटेशन समय में केवल 13 गुना वृद्धि, एल्गोरिथ्म की उच्च दक्षता प्रदर्शित करता है

अभिसरण विश्लेषण

चित्र 4 दर्शाता है कि बड़ी संख्या में शून्यों की गणना करते समय भी, सापेक्ष त्रुटि अच्छी तरह से नियंत्रित रहती है, एल्गोरिथ्म उत्कृष्ट संख्यात्मक स्थिरता प्रदर्शित करता है।

केस सत्यापन

चित्र 2 F(w)=0F(w) = 0 समीकरण समाधान का दृश्य प्रदर्शन करता है, Newton विधि शून्यों को सटीकता से स्थानीयकृत कर सकती है।

संबंधित कार्य

ऐतिहासिक विकास

  1. Carpenter (1992) 2: जनरलाइज्ड बेसेल पॉलीनोमियल शून्यों का प्रारंभिक अनंतस्पर्शी विश्लेषण
  2. Pasquini (2000) 13: जनरलाइज्ड बेसेल पॉलीनोमियल शून्यों की सटीक गणना के लिए संख्यात्मक विधि
  3. Segura (2013) 14: विशेष फ़ंक्शन जटिल शून्यों की गणना के लिए सामान्य विधि
  4. de Bruin आदि (1981) 3,4: जनरलाइज्ड बेसेल पॉलीनोमियल शून्यों के वितरण डोमेन का अध्ययन

इस पेपर के लाभ

  • एकसमानता: विस्तार सभी शून्यों के लिए एकसमान रूप से प्रभावी है
  • सटीकता: मशीन सटीकता स्तर तक पहुँचता है
  • दक्षता: एल्गोरिथ्म जटिलता मौजूदा तरीकों से काफी बेहतर है
  • प्रयोज्यता: व्यापक पैरामीटर श्रेणी को कवर करता है

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. सैद्धांतिक योगदान: रिवर्स जनरलाइज्ड बेसेल पॉलीनोमियल शून्यों का पूर्ण अनंतस्पर्शी सिद्धांत स्थापित किया
  2. संख्यात्मक उपलब्धि: उच्च सटीकता, उच्च दक्षता वाले शून्य गणना एल्गोरिथ्म विकसित किया
  3. व्यावहारिक मूल्य: संबंधित अनुप्रयोगों के लिए विश्वसनीय कम्प्यूटेशनल उपकरण प्रदान किया

सीमाएँ

  1. पैरामीटर प्रतिबंध: विधि Δ1n+32aΔ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n पैरामीटर श्रेणी के लिए लागू है
  2. जटिलता: अनंतस्पर्शी विस्तार की व्युत्पत्ति गणितीय रूप से काफी जटिल है
  3. कार्यान्वयन आवश्यकताएँ: उच्च-सटीकता अंकगणितीय संचालन समर्थन की आवश्यकता है

भविष्य की दिशाएँ

  1. अधिक सामान्य पैरामीटर श्रेणी तक विस्तार करें
  2. अन्य विशेष फ़ंक्शन शून्यों के लिए समान विधियों का अध्ययन करें
  3. एल्गोरिथ्म की कम्प्यूटेशन दक्षता को अनुकूलित करें
  4. वास्तविक इंजीनियरिंग समस्याओं में अनुप्रयोग की खोज करें

गहन मूल्यांकन

शक्तियाँ

  1. सैद्धांतिक कठोरता: टर्निंग पॉइंट सिद्धांत और WKB विधि पर आधारित कठोर गणितीय व्युत्पत्ति
  2. विधि नवाचार: अनंतस्पर्शी विस्तार को संख्यात्मक एल्गोरिथ्म के साथ कुशलतापूर्वक संयोजित करना
  3. असाधारण सटीकता: 101510^{-15} परिमाण स्तर की सापेक्ष त्रुटि प्राप्त करना
  4. उत्कृष्ट दक्षता: पारंपरिक तरीकों की तुलना में एल्गोरिथ्म जटिलता में सुधार
  5. मजबूत व्यावहारिकता: पूर्ण कार्यान्वयन और ओपन-सोर्स कोड प्रदान करना

कमियाँ

  1. सैद्धांतिक जटिलता: विधि का गणितीय आधार काफी गहन है, समझने की सीमा अधिक है
  2. पैरामीटर निर्भरता: विधि की प्रभावशीलता विशिष्ट पैरामीटर श्रेणी मान्यताओं पर निर्भर है
  3. कार्यान्वयन जटिलता: कई जटिल गुणांकों और व्युत्पन्नों की गणना की आवश्यकता है

प्रभाव

  1. शैक्षणिक मूल्य: विशेष फ़ंक्शन सिद्धांत के लिए महत्वपूर्ण सैद्धांतिक उपकरण प्रदान करता है
  2. अनुप्रयोग संभावना: सिग्नल प्रोसेसिंग, नियंत्रण सिद्धांत आदि क्षेत्रों में संभावित अनुप्रयोग
  3. पद्धति महत्व: समान समस्याओं के अनुसंधान के लिए एक उदाहरण प्रदान करता है

प्रयोज्य परिदृश्य

  1. वैज्ञानिक कम्प्यूटिंग: बेसेल पॉलीनोमियल शून्यों की उच्च-सटीकता गणना की आवश्यकता वाली स्थितियाँ
  2. इंजीनियरिंग अनुप्रयोग: फ़िल्टर डिज़ाइन, नियंत्रण प्रणाली विश्लेषण आदि
  3. सैद्धांतिक अनुसंधान: विशेष फ़ंक्शन और अनंतस्पर्शी विश्लेषण का आगे का अनुसंधान

संदर्भ

यह पेपर 14 महत्वपूर्ण संदर्भों का हवाला देता है, जो अनंतस्पर्शी विश्लेषण, विशेष फ़ंक्शन, संख्यात्मक एल्गोरिथ्म आदि संबंधित क्षेत्रों के शास्त्रीय और अग्रणी कार्यों को शामिल करते हैं, जो अनुसंधान के लिए एक ठोस सैद्धांतिक आधार प्रदान करते हैं।


कुल मूल्यांकन: यह एक उच्च-गुणवत्ता वाला गणितीय विश्लेषण पेपर है, जो सैद्धांतिक व्युत्पत्ति और संख्यात्मक कार्यान्वयन दोनों में बहुत उच्च मानदंड तक पहुँचता है। विधि की नवाचारिता, परिणामों की सटीकता और एल्गोरिथ्म की दक्षता सभी प्रशंसनीय हैं, जो संबंधित क्षेत्रों के अनुसंधान और अनुप्रयोग के लिए मूल्यवान उपकरण प्रदान करते हैं।