Approximation by short exponential sums with geometric error decay based on Gauss quadrature

तकनीकी सारांश: गॉस क्वाड्रचर (Gauss Quadrature) पर आधारित ज्यामितीय त्रुटि क्षय के साथ लघु घातांकीय योगों द्वारा सन्निकटन

समस्या विवरण
यह शोध पत्र एक निश्चित, उपयोगकर्ता-निर्धारित लंबाई $N$ वाले लघु घातांकीय योगों (short exponential sums) का उपयोग करके दो विशिष्ट विश्लेषणात्मक फलनों को सन्निकट करने की चुनौती को संबोधित करता है:

1. डोमेन $[0, \infty)$ पर फलन $f_1(x) = \frac{1}{a+x}$ जहाँ $a > 0$।
2. डोमेन $\mathbb{R}$ पर फलन $f_2(x) = e^{-x^2/2\sigma}$ जहाँ $\sigma > 0$।

प्राथमिक उद्देश्य एक ज्यामितीय त्रुटि क्षय (geometric error decay) $\rho^{-2N}$ प्राप्त करना है (जहाँ $\rho > 1$ एक उपयोगकर्ता-निर्धारित पैरामीटर है) और साथ ही संख्यात्मक स्थिरता (numerical stability) बनाए रखना है। लेखक नोट करते हैं कि सभी $x \in [0, \infty)$ के लिए एक ही निश्चित घातांकीय योग के साथ इस क्षय दर को प्राप्त करना असंभव है। फलस्वरूप, समस्या को $\rho$ के वांछित क्षय पैरामीटर पर निर्भर अंतराल $[b_j, b_{j+1})$ पर सन्निकटन बनाने के लिए पुनर्गठित किया गया है।

कार्यप्रणाली
प्रस्तावित विधि लक्षित फलनों को पैरामीट्रिक समाकल निरूपणों (parametric integral representations) के माध्यम से प्रस्तुत करने पर निर्भर करती है जो सीधे गॉसियन क्वाड्रचर सूत्रों के प्रत्यक्ष अनुप्रयोग की अनुमति देते हैं।

• $f_1(x)$ के लिए: लेखक समाकल $\frac{1}{x+a} = \frac{1}{ab} \int_0^\infty e^{-s(\frac{x+a}{ab}-1)} e^{-s} ds$ का उपयोग करते हैं। वे इस समाकल पर गॉस-लैगुएर (Gauss-Laguerre) क्वाड्रचर लागू करते हैं। एक प्रमुख सैद्धांतिक योगदान इंटीग्रैंड की विशिष्ट संरचना के अनुकूल एक नया, सटीक त्रुटि अनुमान (प्रमेय 2.1) प्राप्त करना है, जो मानक पीनो कर्नेल (Peano kernel) अनुमानों में सुधार करता है।
• $f_2(x)$ के लिए: गॉसियन को फूरियर रूपांतरण और पहचान $e^{-x^2/2\sigma} = \frac{1}{\sqrt{\pi}} \int_{-\infty}^\infty e^{-u^2} \cos(\frac{x}{\sqrt{\sigma}}\sqrt{2}u) du$ का उपयोग करके दर्शाया गया है। लेखक इस समाकल पर गॉस-हर्मिट (Gauss-Hermite) क्वाड्रचर लागू करते हैं। आवर्ती फलन के मामले के समान, वे हर्मिट बहुपदों (Hermite polynomials) में इंटीग्रैंड के विस्तार के आधार पर विशिष्ट त्रुटि सीमाएँ व्युत्पन्न करते हैं।

प्रमुख योगदान और निर्माण

1. अंतराल रणनीति (Interval Strategy): विधि डोमेन को क्रमिक अंतरालों $[b_j, b_{j+1})$ में विभाजित करती है।
   • $f_1$ के लिए, अंतराल घातांकीय रूप से बढ़ते हैं: $b_{j+1} - b_j = \frac{2a}{\rho-1} (\frac{\rho+1}{\rho-1})^j$।
   • $f_2$ के लिए, अंतराल समदूरस्थ (equidistant) होते हैं: $b_{j+1} - b_j = \frac{2}{\rho} \sqrt{\frac{N\sigma}{e}}$।
2. पैरामीटर पुनरावृत्ति (Parameter Recurrence): इस निर्माण की एक महत्वपूर्ण विशेषता यह है कि किसी भी अंतराल $j > 0$ पर घातांकीय योग के पैरामीटर (भार और घातांक) प्रारंभिक अंतराल $[b_0, b_1)$ के पैरामीटर्स से सरल समन्वय रूपांतरणों (coordinate transformations) के माध्यम से आसानी से संगणना योग्य हैं। यह सुनिश्चित करता है कि विधि प्रत्येक अंतराल के लिए क्वाड्रचर नोड्स और भार को पुन: संगणित किए बिना, गणनात्मक रूप से कुशल और संख्यात्मक रूप से स्थिर है।
3. त्रुटि नियंत्रण: यह तकनीक सभी अंतरालों में एक सुसंगत सापेक्ष त्रुटि संरचना की गारंटी देती है। पूर्ण त्रुटि (absolute error) $c^{(j)} \rho^{-2N}$ द्वारा सीमित होती है, जहाँ अनुक्रम $c^{(j)}$ घातांकीय रूप से घटता है, जिससे यह सुनिश्चित होता है कि सन्निकटन $x$ के बढ़ने पर भी अत्यधिक सटीक बना रहे।
4. अनुप्रयोग:
   • लॉगारिदम (Logarithm): $\frac{1}{t+1}$ के सन्निकटन को एकीकृत करके, लेखक $[1, \infty)$ पर $\log(x)$ के लिए एक लघु घातांकीय योग सन्निकटन प्राप्त करते हैं, जिसमें पूर्वानुमानित ज्यामितीय त्रुटि क्षय होता है।
   • त्रुटि फलन (Error Function): गॉसियन के सन्निकटन को एकीकृत करके, वे त्रुटि फलन $\text{erf}(x)$ के लिए एक अत्यधिक सटीक सन्निकटन का निर्माण करते हैं।

परिणाम और संख्यात्मक प्रयोग
शोध पत्र डबल-प्रिसिजन अंकगणंड (double-precision arithmetic) का उपयोग करते हुए योग लंबाई $N=8$ और $N=10$ के लिए संख्यात्मक उदाहरण प्रदान करता है।

• व्युत्क्रम फलन (Reciprocal Function): $N=8$ और $\rho=4$ के लिए, विधि पहले अंतराल पर $10^{-13}$ से कम की पूर्ण त्रुटि प्राप्त करती है, जिसमें पूर्व-संगणित सुधार पदों (रिमार्क 2.2) का उपयोग करने पर $10^{-15}$ तक के बेहतर सीमाएँ प्राप्त होती हैं।
• गॉसियन फलन: $N=8$ और $\rho=4$ के लिए, $e^{-x^2/2}$ के सन्निकटन से कई अंतरालों में $10^{-15}$ की सीमा में त्रुटियाँ प्राप्त होती हैं।
• लॉगारिदम और त्रुटि फलन: $\log(x)$ और $\text{erf}(x)$ के लिए प्राप्त सन्निकटन $10^{-15}$ की सीमा तक पहुँचने वाली त्रुटियों का प्रदर्शन करते हैं, जो प्रभावी रूप से डबल-प्रिसिजन अंकगणंड की पूर्ण क्षमता का उपयोग करते हैं।
• लेखक इस बात पर जोर देते हैं कि $N \le 10$ के लिए, विधि मशीन-प्रिसिजन सटीकता प्राप्त करने के लिए पर्याप्त है, और सभी आवश्यक पैरामीटर मानक गॉस क्वाड्रचर रूटीन का उपयोग करके उच्च सटीकता के साथ संगणित किए जा सकते हैं।

महत्व और दावे
लेखक दावा करते हैं कि उनका दृष्टिकोण मौजूदा विधियों (जैसे प्रोंनी की विधि, रेमेज़ एल्गोरिदम या विभेदक समीकरण दृष्टिकोण) से भिन्न है क्योंकि यह उपयोगकर्ता को वांछित त्रुटि सीमा और योग लंबाई $N$ को अग्रिम रूप से निर्धारित करने की अनुमति देता है। कॉम्पैक्ट अंतरालों पर विशिष्ट नॉर्म्स (norms) के लिए अनुकूलित करने वाली विधियों के विपरीत, यह तकनीक गारंटीकृत ज्यामितीय त्रुटि क्षय के साथ अर्ध-अनंत या अनंत डोमेन तक सन्निकटन को विस्तारित करने का एक व्यवस्थित तरीका प्रदान करती है।

यह शोध पत्र दावा करता है कि विधि:

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

यह कार्य निष्कर्ष निकालता है कि यह क्वाड्रचर-आधारित रणनीति उन अनुप्रयोगों के लिए एक मजबूत विकल्प प्रदान करती है जिनमें उच्च-परिशुद्धता वाले लघु घातांकीय योगों के रूप में कर्नेल और ग्रीन फलन (Green's functions) का प्रतिनिधित्व करने की आवश्यकता होती है।