रिचर्डसन एक्सट्रपलेशन



संख्यात्मक विश्लेषण में, रिचर्डसन एक्सट्रपलेशन एक अनुक्रम त्वरण विधि है जिसका उपयोग कुछ मूल्य $$A^\ast = \lim_{h\to 0} A(h)$$ के अनुमानों के अनुक्रम के अभिसरण की दर में सुधार करने के लिए किया जाता है। संक्षेप में, $$h$$ के कई मानों के लिए $$A(h)$$ का मान दिया गया है, हम अनुमानों को $$h=0$$ पर एक्सट्रपलेशन करके $$A^\ast$$का अनुमान लगा सकते हैं। इसका नाम लुईस फ्राई रिचर्डसन के नाम पर रखा गया है, जिन्होंने 20वीं सदी की प्रारंभ में इस तकनीक की प्रारंभ की थी चूँकि यह विचार क्रिस्टियान ह्यूजेंस को π की गणना में पहले से ही पता था। बिरखॉफ और रोटा के शब्दों में, "व्यावहारिक गणनाओं के लिए इसकी उपयोगिता को संभवतः ही कम करके आंका जा सकता है।"

रिचर्डसन एक्सट्रपलेशन के व्यावहारिक अनुप्रयोगों में रोमबर्ग एकीकरण सम्मिलित है, जो ट्रेपेज़ॉइड नियम पर रिचर्डसन एक्सट्रपलेशन को प्रयुक्त करता है, और सामान्य अंतर समीकरणों को हल करने के लिए बुलिर्श-स्टोअर एल्गोरिदम है।

संकेतन
मान लीजिए कि$$A_0(h)$$ $$A^*$$(स्पष्ट मान) का एक अनुमान है जो प्रपत्र के त्रुटि सूत्र के साथ धनात्मक चरण आकार h पर निर्भर करता है
 * $$ A^* = A_0(h)+a_0h^{k_0} + a_1h^{k_1} + a_2h^{k_2} + \cdots $$

जहां $$a_i$$ अज्ञात स्थिरांक हैं और$$k_i$$ ज्ञात स्थिरांक हैं जैसे कि $$h^{k_i}>h^{k_{i+1}}$$ इसके अतिरिक्त $$O(h^{k_i})$$,$$A_i(h)$$ सन्निकटन की काट-छाँट त्रुटि को इस प्रकार दर्शाता है कि $$A^*=A_i(h)+O(h^{k_i}).$$ इसी प्रकार,$$A^*=A_i(h)+O(h^{k_i}),$$ में सन्निकटन $$A_i(h)$$} को $$O(h^{k_i})$$ सन्निकटन कहा जाता है।

ध्यान दें कि बिग ओ अंकन  के साथ सरलीकरण करके, निम्नलिखित सूत्र समतुल्य हैं:

$$ \begin{align} A^* &= A_0(h) + a_0h^{k_0} + a_1h^{k_1} + a_2h^{k_2} + \cdots \\ A^* &= A_0(h)+ a_0h^{k_0} + O(h^{k_1}) \\ A^* &= A_0(h)+O(h^{k_0}) \end{align} $$

उद्देश्य
रिचर्डसन एक्सट्रपलेशन एक ऐसी प्रक्रिया है जो त्रुटि सूत्र को $$A^*=A_0(h)+O(h^{k_0})$$ को $$A^*=A_1(h)+O(h^{k_1}).$$ में बदलकर $$A^*$$का उत्तम अनुमान लगाती है, इसलिए, $$A_0(h)$$ को $$A_1(h)$$ से बदलने से ट्रंकेशन त्रुटि$$O(h^{k_0}) $$ से कम हो गई है $$O(h^{k_1}) $$ समान चरण आकार $$O(h^{k_1}) $$ के लिए। सामान्य पैटर्न तब होता है जिसमें $$A_i(h)$$, $$i>j$$ होने पर $$A_i(h)$$ की तुलना में अधिक स्पष्ट  अनुमान होता है। इस प्रक्रिया द्वारा, हमने त्रुटि में सबसे बड़े पद जो कि $$O(h^{k_0}) $$ था, को घटाकर $$A^*$$ का उत्तम अनुमान प्राप्त किया है। उत्तम अनुमान प्राप्त करने के लिए अधिक त्रुटि शब्दों को हटाने के लिए इस प्रक्रिया को दोहराया जा सकता है।

प्रक्रिया
कुछ स्थिरांक $$t$$ के लिए चरण आकार $$h$$ और $$h / t$$ का उपयोग करते हुए,$$A^*$$ के लिए दो सूत्र हैं:

$$\begin{align} A^* &= A_0(h)+ a_0h^{k_0} + a_1h^{k_1} + a_2h^{k_2} + O(h^{k_3}) & (1) \\\\ A^* &= A_0\!\left(\frac{h}{t}\right) + a_0\left(\frac{h}{t}\right)^{k_0} + a_1\left(\frac{h}{t}\right)^{k_1} + a_2\left(\frac{h}{t}\right)^{k_2} + O(h^{k_3}) & (2) \end{align}$$ पहले त्रुटि पद को हटाकर अपने सन्निकटन को $$O(h^{k_0})$$ से $$O(h^{k_1})$$ तक सुधारने के लिए, हम दूसरे समीकरण (2) को $$t^{k_0}$$ से गुणा करते हैं और पहले समीकरण (1) को घटाकर हमें प्राप्त होता है$$ (t^{k_0}-1)A^* = \bigg[t^{k_0}A_0\left(\frac{h}{t}\right) - A_0(h)\bigg] + \bigg(t^{k_0}a_1\bigg(\frac{h}{t}\bigg)^{k_1}-a_1h^{k_1}\bigg)+ \bigg(t^{k_0}a_2\bigg(\frac{h}{t}\bigg)^{k_2}-a_2h^{k_2}\bigg) + O(h^{k_3}). $$यह गुणा और घटाव इसलिए किया गया क्योंकि $\big[t^{k_0}A_0\left(\frac{h}{t}\right) - A_0(h)\big]$, $$O(h^{k_1})$$ का $$(t^{k_0}-1)A^*$$ सन्निकटन है। हम $$A^*$$ देने के लिए अपने वर्तमान सूत्र को हल कर सकते हैं$$A^* = \frac{\bigg[t^{k_0}A_0\left(\frac{h}{t}\right) - A_0(h)\bigg]}{t^{k_0}-1}

+ \frac{\bigg(t^{k_0}a_1\bigg(\frac{h}{t}\bigg)^{k_1}-a_1h^{k_1}\bigg)}{t^{k_0}-1}

+\frac{\bigg(t^{k_0}a_2\bigg(\frac{h}{t}\bigg)^{k_2}-a_2h^{k_2}\bigg)}{t^{k_0}-1}

+O(h^{k_3}) $$जिसे सेटिंग करके$$A^* = A_1(h)+O(h^{k_1})$$लिखा जा सकता है


 * $$A_1(h) = \frac{t^{k_0}A_0\left(\frac{h}{t}\right) - A_0(h)}{t^{k_0}-1} .$$

पुनरावृत्ति संबंध
एक सामान्य पुनरावृत्ति संबंध को सन्निकटन के लिए परिभाषित किया जा सकता है
 * $$ A_{i+1}(h) = \frac{t^{k_i}A_i\left(\frac{h}{t}\right) - A_i(h)}{t^{k_i}-1} $$

जहाँ $$k_{i+1}$$ संतुष्ट
 * $$ A^* = A_{i+1}(h) + O(h^{k_{i+1}}) $$.

गुण
रिचर्डसन एक्सट्रपलेशन को एक रैखिक अनुक्रम परिवर्तन के रूप में माना जा सकता है।

इसके अतिरिक्त, सामान्य सूत्र का उपयोग $$k_0$$ (ट्रंकेशन त्रुटि का अग्रणी क्रम चरण आकार व्यवहार) का अनुमान लगाने के लिए किया जा सकता है, जब न तो इसका मान और न ही $$A^*$$ को प्राथमिकता से जाना जाता है। ऐसी तकनीक अभिसरण की अज्ञात दर को मापने के लिए उपयोगी हो सकती है। तीन अलग-अलग चरण आकारों $$h$$,$$h / t$$, से $$A^*$$ का अनुमान दिया गया है, और स्पष्ट संबंध $$h / s$$ दिया गया है$$A^*=\frac{t^{k_0}A_i\left(\frac{h}{t}\right) - A_i(h)}{t^{k_0}-1} + O(h^{k_1}) = \frac{s^{k_0}A_i\left(\frac{h}{s}\right) - A_i(h)}{s^{k_0}-1} + O(h^{k_1})$$एक अनुमानित संबंध उत्पन्न करता है (कृपया ध्यान दें कि यहां संकेतन थोड़ा अस्पष्ट  उत्पन्न कर सकता है, उपरोक्त समीकरण में दिखाई देने वाले दो ओ केवल अग्रणी क्रम चरण आकार के व्यवहार को निरुपित करते हैं किंतु  उनके स्पष्ट रूप अलग-अलग हैं और इसलिए दो ओ शब्दों को समाप्त करना लगभग मान्य है)$$A_i\left(\frac{h}{t}\right) + \frac{A_i\left(\frac{h}{t}\right) - A_i(h)}{t^{k_0}-1} \approx A_i\left(\frac{h}{s}\right) +\frac{A_i\left(\frac{h}{s}\right) - A_i(h)}{s^{k_0}-1}$$जिसे $$h$$,$$s$$, और $$t$$ के कुछ इच्छित विकल्पों के लिए $$k_0$$ का अनुमान लगाने के लिए संख्यात्मक रूप से हल किया जा सकता है।

रिचर्डसन एक्सट्रपलेशन का उदाहरण
मान लीजिए कि हम $$A^*$$ का अनुमान लगाना चाहते हैं, और हमारे पास एक विधि $$A(h)$$) है जो एक छोटे पैरामीटर $$h$$ पर इस तरह निर्भर करती है कि $$A(h) = A^\ast + C h^n + O(h^{n+1}).$$ आइए एक नए फलन को परिभाषित करें$$ R(h,t) := \frac{ t^n A(h/t) - A(h)}{t^n-1} $$जहाँ $$h$$ और $$\frac{h}{t}$$ दो अलग-अलग चरण आकार हैं।

तब$$ R(h, t) = \frac{ t^n ( A^* + C \left(\frac{h}{t}\right)^n + O(h^{n+1}) ) - ( A^* + C h^n + O(h^{n+1}) ) }{ t^n - 1} = A^* + O(h^{n+1}). $$$$ R(h,t) $$ को $$ A(h) $$ का रिचर्डसन एक्सट्रपलेशन कहा जाता है, और इसमें $$ A(h) $$ की तुलना में उच्च-क्रम त्रुटि अनुमान $$ O(h^{n+1}) $$ है।

बहुत बार, बहुत छोटे h' के साथ A(h') के अतिरिक्त R(h) का उपयोग करके दी गई स्पष्टता प्राप्त करना बहुत आसान होता है। जहां A(h') सीमित परिशुद्धता (गोल त्रुटियों) और/या आवश्यक गणनाओं की बढ़ती संख्या के कारण समस्याएं उत्पत्ति कर सकता है (नीचे उदाहरण देखें)।

रिचर्डसन एक्सट्रपलेशन के लिए उदाहरण छद्मकोड कोड
मैटलैब शैली में निम्नलिखित छद्मकोड ट्रैपेज़ॉइडल विधि के साथ ओडीई $$y'(t) = -y^2$$, $$y(0) = 1$$को हल करने में सहायता करने के लिए रिचर्डसन एक्सट्रपलेशन को प्रदर्शित करता है। इस उदाहरण में हम प्रत्येक पुनरावृत्ति में चरण आकार को आधा कर देते हैं और इसलिए ऊपर की चर्चा में हमारे पास वह $$t = 2$$ होगा। ट्रैपेज़ॉइडल विधि की त्रुटि को विषम शक्तियों के संदर्भ में व्यक्त किया जा सकता है जिससे कई चरणों में त्रुटि को सम शक्तियों में व्यक्त किया जा सके; यह हमें $$t$$ को दूसरी घात तक बढ़ाने और छद्मकोड में $$4 = 2^2 = t^2$$ की घातें लेने की ओर ले जाता है। हम $$y(5)$$ का मान ज्ञात करना चाहते हैं, जिसका स्पष्ट समाधान$$\frac{1}{5 + 1} = \frac{1}{6} = 0.1666...$$ है... क्योंकि ODE का स्पष्ट समाधान $$y(t) = \frac{1}{1 + t}$$ है। यह छद्म कोड मानता है कि नामक एक फलन  उपस्थित है जो प्रारंभिक  बिंदु और   और चरण आकार एच के साथ फलन   पर ट्रैपेज़ॉयडल विधि निष्पादित करके   की गणना करने का प्रयास करता है।

ध्यान दें कि प्रारंभिक चरण के आकार को बहुत छोटे से प्रारंभ करने से संभावित रूप से अंतिम समाधान में त्रुटि आ सकती है। चूँकि सर्वोत्तम प्रारंभिक चरण आकार चुनने में सहायता करने के लिए डिज़ाइन की गई विधियाँ हैं, एक विकल्प बड़े चरण आकार के साथ प्रारंभ करना है और फिर रिचर्डसन एक्सट्रपलेशन को प्रत्येक पुनरावृत्ति चरण आकार को कम करने की अनुमति देना है जब तक कि त्रुटि वांछित सहनशीलता तक नहीं पहुंच जाती।

यह भी देखें

 * ऐटकेन की डेल्टा-वर्ग प्रक्रिया
 * केंको ताकेबे
 * रिचर्डसन इटेरशन

संदर्भ

 * Extrapolation Methods. Theory and Practice by C. Brezinski and M. Redivo Zaglia, North-Holland, 1991.
 * Ivan Dimov, Zahari Zlatev, Istvan Farago, Agnes Havasi: Richardson Extrapolation: Practical Aspects and Applications'', Walter de Gruyter GmbH & Co KG, ISBN 9783110533002 (2017).

बाहरी संबंध

 * Fundamental Methods of Numerical Extrapolation With Applications, mit.edu
 * Richardson-Extrapolation
 * Richardson extrapolation on a website of Robert Israel (University of British Columbia)
 * Matlab code for generic Richardson extrapolation.
 * Julia code for generic Richardson extrapolation.