Fast and stable computation of highly oscillatory and/or exponentially decaying integrals using a Clenshaw-Curtis product-integration rule

तकनीकी सारांश: अत्यधिक दोलनशील और/या तेजी से क्षय होने वाले समाकलों की तीव्र और स्थिर गणना

लेखक: विक्टर डोमिंगुएज़ (Víctor Domínguez)
दिनांक: 10 फरवरी, 2026
विषय क्षेत्र: संख्यात्मक विश्लेषण (Numerical Analysis) (क्वाड्रचर नियम, दोलनशील समाकल)

1. समस्या

यह शोध पत्र निम्नलिखित रूप के समाकलों के संख्यात्मक मूल्यांकन को संबोधित करता है:
$$I_z(f) = \int_{0}^{2} f(s) \exp(zs) \, ds$$
जहाँ $z \in \mathbb{C}$ एक जटिल पैरामीटर है। मुख्य चुनौती तब उत्पन्न होती है जब $z$ का काल्पनिक भाग (imaginary part) बड़ा होता है (उच्च दोलन का कारण), या वास्तविक भाग (real part) का ऋणात्मक मान बहुत अधिक होता है (तेजी से घातीय क्षय का कारण), या दोनों स्थितियाँ एक साथ होती हैं।

मानक क्वाड्रचर नियम (जैसे Gauss-Legendre) इन परिदृश्यों में विफल हो जाते हैं क्योंकि आवश्यक नोड्स (nodes) की संख्या $|z|$ के साथ बढ़ती है, जिससे गणना की लागत अत्यधिक बढ़ जाती है। हालांकि मौजूदा विधियाँ जैसे Filon-type या Steepest Descent मौजूद हैं, वे अक्सर तब संघर्ष करती हैं जब $z$ व्यापक रूप से बदलता है या जब समाकल (integrand) में क्षय और दोलन दोनों का संयोजन होता है।

2. कार्यप्रणाली: प्रोडक्ट क्लेंसवॉ-कर्टिस (Product Clenshaw-Curtis - PCC) नियम

लेखक एक प्रोडक्ट क्लेंसवॉ-कर्टिस (PCC) नियम प्रस्तावित करते हैं। कार्यप्रणाली निम्नलिखित चरणों का पालन करती है:

• अंतर्वेशन (Interpolation): समाकल के सुचारू भाग, $f(s)$, को शिफ्ट किए गए चेबिशेव नोड्स (Chebyshev nodes) पर एक बहुपद $p_L(s)$ का उपयोग करके अंतर्वेशित किया जाता है।
• भारों का सटीक एकीकरण (Exact Integration of Weights): पूरे समाकल का अनुमान लगाने के बजाय, यह नियम अंतर्वेशित बहुपद और घातीय पद के गुणनफल का सटीक एकीकरण करता है:
  $$I_{L,z}(f) = \sum_{\ell=0}^{L} \!'' \alpha_{\ell,L}(f) \omega_\ell(z)$$
  जहाँ $\omega_\ell(z) = \int_{0}^{2} T_\ell(s-1) \exp(zs) \, ds$ प्रथम श्रेणी के चेबिशेव बहुपदों से जुड़े "भार" (weights) हैं।
• तीन-चरणीय एल्गोरिदम कार्यान्वयन: भार $\omega_\ell(z)$ की गणना स्थिर और कुशल तरीके से सुनिश्चित करने के लिए, लेखक एक सुदृढ़ तीन-चरणीय एल्गोरिदम पेश करते हैं:
  1. चरण 1 (पुनरावृत्ति - Recurrence): छोटे सूचकांकों $n \lesssim |z|^{1/2}$ के लिए प्रत्यक्ष पुनरावृत्ति संबंध का उपयोग करता है।
  2. चरण 2 (त्रिकोणीय प्रणाली - Tridiagonal System): बड़े सूचकांकों के लिए जहाँ पुनरावृत्ति अस्थिर हो जाती है, समस्या को एक त्रिकोणीय रैखिक प्रणाली के रूप में पुनर्गठित किया जाता है, जिसे थॉमस एल्गोरिदम का उपयोग करके $O(L)$ समय में हल किया जाता है।
  3. चरण 3 (चरम गुणांक - Extreme Coefficient): त्रिकोणीय प्रणाली को बंद करने के लिए आवश्यक सीमा गुणांक $\rho_{n_1+1}(z)$ की गणना करने के लिए एक विशेष विधि का उपयोग करता है।

3. मुख्य योगदान

• एकीकृत ढांचा (Unified Framework): पिछले कार्यों के विपरीत जो शुद्ध रूप से दोलनशील ($\Re z = 0$) और शुद्ध रूप से क्षयकारी ($\Im z = 0$) मामलों को अलग-अलग मानते हैं, यह विधि सामान्य जटिल मामले के लिए एक एकल, स्थिर ढांचा प्रदान करती है।
• स्थिरता विश्लेषण (Stability Analysis): लेखक प्रमाण प्रदान करते हैं कि यह एल्गोरिदम मानक पुनरावृत्ति संबंधों के विशिष्ट घातीय त्रुटि विकास (exponential error growth) से बचता है। त्रिकोणीय प्रणाली की कंडीशन नंबर (condition number) को सुव्यवस्थित दिखाया गया है।
• जटिलता (Complexity): इसकी गणना लागत अनिवार्य रूप से $|z|$ से स्वतंत्र है और नोड्स की संख्या $L$ के साथ रैखिक रूप से स्केल करती है, जो इसे "स्टिफ" (stiff) समाकलों के लिए अत्यधिक कुशल बनाती है।
• विशेष फलनों का परिहार (Avoidance of Special Functions): इसका कार्यान्वयन जटिल विशेष फलनों के बजाय बुनियादी अंकगणित और रैखिक बीजगणित पर निर्भर करता है, जो संख्यात्मक सुदृढ़ता को बढ़ाता है।

4. परिणाम और सत्यापन

शोध पत्र कई कठोर परीक्षणों के माध्यम से विधि को सत्यापित करता है:

• अभिसरण दरें (Convergence Rates): सैद्धांतिक विश्लेषण दिखाता है कि सुचारू $f$ के लिए, त्रुटि $L$ के संबंध में सुपर-अल्जेब्रिक (superalgebraic) रूप से घटती है और निश्चित $L$ के लिए $O(|z|^{-2})$ (Filon-type convergence) का अनुसरण करती है।
• संख्यात्मक प्रयोग:
  • सुचारू फलन: मशीन परिशुद्धता (machine precision) तक सुपर-अल्जेब्रिक अभिसरण प्रदर्शित किया गया।
  • विलक्षण समाकल (Singular Integrands): एंडपॉइंट विलक्षणताओं (जैसे $\alpha = 1/2, 3/2$) वाले फलनों का परीक्षण किया गया। हालांकि विलक्षणताएं $L$ के अभिसरण दर को कम करती हैं, फिर भी यह विधि अत्यधिक प्रभावी बनी रहती है।
  • स्थिरता परीक्षण: पुष्टि की गई कि तीन-चरणीय एल्गोरिदम राउंड-ऑफ त्रुटियों के कारण विफल होने वाली प्रत्यक्ष पुनरावृत्ति संबंध के बावजूद परिशुद्धता बनाए रखता है।
• वास्तविक दुनिया का अनुप्रयोग: लेखक लैपलेस ट्रांसफॉर्म व्युत्क्रमण (Laplace transform inversion) के माध्यम से फ्रैक्शनल-ऑर्डर इवोल्यूशन समीकरणों (Riemann-Liouville फ्रैक्शनल डेरिवेटिव्स) को हल करने के लिए इस नियम की उपयोगिता प्रदर्शित करते हैं। यह क्वाड्रचर नियम ब्रॉमविच कंटूर इंटीग्रेशन (Bromwich contour integration) के लिए आवश्यक $z$ के विस्तृत रेंज को सफलतापूर्वक संभालता है।

5. महत्व

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