रैखिक मल्टीस्टेप विधि: Difference between revisions

From Vigyanwiki
(Created page with "{{Short description|Class of iterative numerical methods for solving differential equations}} {{redirect|Adams' method|the electoral apportionment method|Method of smallest di...")
 
(text)
Line 1: Line 1:
{{Short description|Class of iterative numerical methods for solving differential equations}}
{{Short description|Class of iterative numerical methods for solving differential equations}}
{{redirect|Adams' method|the electoral apportionment method|Method of smallest divisors}}
{{redirect|एडम्स की विधि|चुनावी विभाजन विधि|सबसे छोटे भाजक की विधि}}
{{no footnotes|date=June 2017}}
 
[[संख्यात्मक साधारण अंतर समीकरण]]ों के लिए रैखिक मल्टीस्टेप विधियों का उपयोग किया जाता है। वैचारिक रूप से, एक संख्यात्मक विधि एक प्रारंभिक बिंदु से शुरू होती है और फिर अगले समाधान बिंदु को खोजने के लिए समय में एक छोटा कदम आगे बढ़ाती है। समाधान निकालने के लिए प्रक्रिया बाद के चरणों के साथ जारी रहती है। एकल-चरण विधियाँ (जैसे यूलर की विधि) वर्तमान मूल्य निर्धारित करने के लिए केवल एक पिछले बिंदु और उसके व्युत्पन्न को संदर्भित करती हैं। रनगे-कुट्टा विधियां|रंज-कुट्टा जैसी विधियां उच्च क्रम विधि प्राप्त करने के लिए कुछ मध्यवर्ती कदम (उदाहरण के लिए, आधा कदम) लेती हैं, लेकिन फिर दूसरा कदम उठाने से पहले सभी पिछली जानकारी को त्याग देती हैं। मल्टीस्टेप विधियाँ पिछले चरणों की जानकारी को त्यागने के बजाय उसे बनाए रखने और उसका उपयोग करके दक्षता हासिल करने का प्रयास करती हैं। नतीजतन, मल्टीस्टेप विधियां कई पिछले बिंदुओं और व्युत्पन्न मूल्यों को संदर्भित करती हैं। ''रैखिक'' मल्टीस्टेप विधियों के मामले में, पिछले बिंदुओं और व्युत्पन्न मूल्यों के एक [[रैखिक संयोजन]] का उपयोग किया जाता है।
[[संख्यात्मक साधारण अंतर समीकरण]] के लिए रैखिक बहुपदीय विधियों का उपयोग किया जाता है। वैचारिक रूप से, एक संख्यात्मक विधि एक प्रारंभिक बिंदु से प्रारम्भ होती है और फिर अगले समाधान बिंदु को खोजने के लिए समय में एक छोटा कदम आगे बढ़ाती है। समाधान निकालने के लिए प्रक्रिया बाद के चरणों के साथ जारी रहती है। एकल-चरण विधियाँ (जैसे यूलर की विधि) वर्तमान मूल्य निर्धारित करने के लिए केवल एक पिछले बिंदु और उसके व्युत्पन्न को संदर्भित करती हैं। रंज-कुट्टा जैसी विधियां उच्च क्रम विधि प्राप्त करने के लिए कुछ मध्यवर्ती कदम (उदाहरण के लिए, आधा कदम) लेती हैं, लेकिन फिर दूसरा कदम उठाने से पहले सभी पिछली जानकारी को त्याग देती हैं। बहुपदीय विधियाँ पिछले चरणों की जानकारी को त्यागने के स्थान पर उसे बनाए रखने और उसका उपयोग करके दक्षता प्राप्त करने का प्रयास करती हैं। नतीजतन, बहुपदीय विधियां कई पिछले बिंदुओं और व्युत्पन्न मूल्यों को संदर्भित करती हैं। रैखिक बहुपदीय विधियों की स्तिथि में, पिछले बिंदुओं और व्युत्पन्न मूल्यों के एक [[रैखिक संयोजन]] का उपयोग किया जाता है।


== परिभाषाएँ ==
== परिभाषाएँ ==
साधारण अंतर समीकरणों के लिए संख्यात्मक विधियाँ फॉर्म की [[प्रारंभिक मूल्य समस्या]]ओं का अनुमानित समाधान करती हैं
साधारण अंतर समीकरणों के लिए संख्यात्मक विधियाँ विधि की [[प्रारंभिक मूल्य समस्या|प्रारंभिक मान समस्या]] का अनुमानित समाधान करती हैं
<math display="block"> y' = f(t,y), \quad y(t_0) = y_0. </math>
<math display="block"> y' = f(t,y), \quad y(t_0) = y_0. </math>
परिणाम के मूल्य के लिए अनुमान है <math> y(t) </math> अलग-अलग समय पर <math> t_i </math>:
परिणाम के मूल्य के लिए अनुमान <math> y(t) </math> अलग-अलग समय पर <math> t_i </math> है:
<math display="block"> y_i \approx y(t_i) \quad\text{where}\quad t_i = t_0 + i h, </math>
<math display="block"> y_i \approx y(t_i) \quad\text{where}\quad t_i = t_0 + i h, </math>
कहाँ <math> h </math> समय चरण है (कभी-कभी इसे कहा जाता है)। <math> \Delta t </math>) और <math>i</math> एक पूर्णांक है.
जहाँ <math> h </math> समय चरण है (कभी-कभी इसे <math> \Delta t </math> कहा जाता है) और <math>i</math> एक पूर्णांक है।


मल्टीस्टेप विधियाँ पिछली जानकारी का उपयोग करती हैं <math> s </math> अगले मान की गणना करने के चरण. विशेष रूप से, एक रैखिक मल्टीस्टेप विधि एक रैखिक संयोजन का उपयोग करती है <math> y_i </math> और <math> f(t_i,y_i) </math> के मूल्य की गणना करने के लिए <math> y </math> वांछित वर्तमान चरण के लिए. इस प्रकार, एक रैखिक मल्टीस्टेप विधि फॉर्म की एक विधि है
बहुपदीय विधियाँ अगले मान की गणना करने के लिए पिछले चरणों <math> s </math> की जानकारी का उपयोग करती हैं। विशेष रूप से, एक रैखिक बहुपदीय विधि वांछित वर्तमान चरण के लिए <math> y </math> के मान की गणना करने के लिए <math> y_i </math> और <math> f(t_i,y_i) </math> के रैखिक संयोजन का उपयोग करती है। इस प्रकार, एक रैखिक बहुपदीय विधि रूप की एक विधि है
<math display="block"> \begin{align}
<math display="block"> \begin{align}
& y_{n+s} + a_{s-1} \cdot y_{n+s-1} + a_{s-2} \cdot y_{n+s-2} + \cdots + a_0 \cdot y_n \\
& y_{n+s} + a_{s-1} \cdot y_{n+s-1} + a_{s-2} \cdot y_{n+s-2} + \cdots + a_0 \cdot y_n \\
Line 17: Line 17:
& \Leftrightarrow \sum_{j=0}^s a_jy_{n+j} = h\sum_{j=0}^sb_jf(t_{n+j},y_{n+j}),
& \Leftrightarrow \sum_{j=0}^s a_jy_{n+j} = h\sum_{j=0}^sb_jf(t_{n+j},y_{n+j}),
\end{align} </math>
\end{align} </math>
साथ <math>a_s=1</math>. गुणांक <math> a_0, \dotsc, a_{s-1} </math> और <math> b_0, \dotsc, b_s </math> विधि निर्धारित करें. विधि का डिजाइनर लागू करने में आसान विधि प्राप्त करने की इच्छा के विरुद्ध सही समाधान के लिए एक अच्छा अनुमान प्राप्त करने की आवश्यकता को संतुलित करते हुए, गुणांक का चयन करता है। विधि को सरल बनाने के लिए अक्सर कई गुणांक शून्य होते हैं।
<math>a_s=1</math> के साथ है। गुणांक <math> a_0, \dotsc, a_{s-1} </math> और <math> b_0, \dotsc, b_s </math> विधि निर्धारित करें। विधि का अभिकल्पक लागू करने में आसान विधि प्राप्त करने की इच्छा के विरुद्ध सही समाधान के लिए एक अच्छा अनुमान प्राप्त करने की आवश्यकता को संतुलित करते हुए, गुणांक का चयन करता है। विधि को सरल बनाने के लिए प्रायः कई गुणांक शून्य होते हैं।


कोई भी स्पष्ट और अंतर्निहित तरीकों के बीच अंतर कर सकता है। अगर <math> b_s = 0 </math>, तो विधि को स्पष्ट कहा जाता है, क्योंकि सूत्र सीधे गणना कर सकता है <math> y_{n+s} </math>. अगर <math> b_s \ne 0 </math> तो विधि को अंतर्निहित कहा जाता है, क्योंकि इसका मान है <math> y_{n+s} </math> के मूल्य पर निर्भर करता है <math> f(t_{n+s}, y_{n+s}) </math>, और समीकरण को हल किया जाना चाहिए <math> y_{n+s} </math>. अंतर्निहित सूत्र को हल करने के लिए अक्सर न्यूटन की विधि जैसी पुनरावृत्तीय विधियों का उपयोग किया जाता है।
कोई भी स्पष्ट और अंतर्निहित तरीकों के बीच अंतर कर सकता है। अगर <math> b_s = 0 </math>, तो विधि को स्पष्ट कहा जाता है, क्योंकि सूत्र <math> y_{n+s} </math> सीधे गणना कर सकता है। अगर <math> b_s \ne 0 </math> तो विधि को अंतर्निहित कहा जाता है, क्योंकि इसका मान <math> y_{n+s} </math> के मूल्य <math> f(t_{n+s}, y_{n+s}) </math> पर निर्भर करता है, और समीकरण को हल <math> y_{n+s} </math> किया जाना चाहिए। अंतर्निहित सूत्र को हल करने के लिए प्रायः न्यूटन की विधि जैसी पुनरावृत्तीय विधियों का उपयोग किया जाता है।


कभी-कभी मूल्य की भविष्यवाणी करने के लिए एक स्पष्ट मल्टीस्टेप विधि का उपयोग किया जाता है <math> y_{n+s} </math>. फिर उस मान को सही करने के लिए एक अंतर्निहित सूत्र में उपयोग किया जाता है। परिणाम एक भविष्यवक्ता-सुधारक विधि है।
कभी-कभी मूल्य की भविष्यवाणी करने के लिए एक स्पष्ट बहुपदीय विधि <math> y_{n+s} </math> का उपयोग किया जाता है। फिर उस मान को सही करने के लिए एक अंतर्निहित सूत्र में उपयोग किया जाता है। परिणाम एक भविष्यवक्ता-सुधारक विधि है।


== उदाहरण ==
== उदाहरण ==
Line 27: Line 27:
उदाहरण के लिए समस्या पर विचार करें
उदाहरण के लिए समस्या पर विचार करें
<math display="block"> y' = f(t,y)=y, \quad y(0) = 1. </math>
<math display="block"> y' = f(t,y)=y, \quad y(0) = 1. </math>
सटीक समाधान है <math> y(t) = e^t </math>.
सटीक समाधान <math> y(t) = e^t </math> है।


=== वन-स्टेप यूलर ===
=== वन-चरण यूलर ===
एक सरल संख्यात्मक विधि यूलर की विधि है:
एक सरल संख्यात्मक विधि यूलर की विधि है:
<math display="block"> y_{n+1} = y_n + hf(t_n, y_n). </math>
<math display="block"> y_{n+1} = y_n + hf(t_n, y_n). </math>
यूलर की विधि को एक चरण के विकृत मामले के लिए एक स्पष्ट मल्टीस्टेप विधि के रूप में देखा जा सकता है।
यूलर की विधि को एक चरण के विकृत स्तिथि के लिए एक स्पष्ट बहुपदीय विधि के रूप में देखा जा सकता है।


यह विधि, चरण आकार के साथ लागू होती है <math> h = \tfrac{1}{2} </math> समस्या पर <math> y' = y </math>, निम्नलिखित परिणाम देता है:
समस्या <math> y' = y </math> पर चरण आकार <math> h = \tfrac{1}{2} </math> के साथ लागू की गई यह विधि निम्नलिखित परिणाम देती है:
<math display="block"> \begin{align}
<math display="block"> \begin{align}
   y_1 &= y_0 + hf(t_0, y_0) = 1 + \tfrac{1}{2} \cdot 1 = 1.5, \\
   y_1 &= y_0 + hf(t_0, y_0) = 1 + \tfrac{1}{2} \cdot 1 = 1.5, \\
Line 46: Line 46:
यूलर की विधि एक चरणीय विधि है। एक सरल बहुचरणीय विधि दो-चरणीय एडम्स-बैशफोर्थ विधि है
यूलर की विधि एक चरणीय विधि है। एक सरल बहुचरणीय विधि दो-चरणीय एडम्स-बैशफोर्थ विधि है
<math display="block"> y_{n+2} = y_{n+1} + \tfrac{3}{2} hf(t_{n+1},y_{n+1}) - \tfrac{1}{2} hf(t_n,y_n). </math>
<math display="block"> y_{n+2} = y_{n+1} + \tfrac{3}{2} hf(t_{n+1},y_{n+1}) - \tfrac{1}{2} hf(t_n,y_n). </math>
इस विधि के लिए दो मानों की आवश्यकता है, <math> y_{n+1} </math> और <math> y_n </math>, अगले मान की गणना करने के लिए, <math> y_{n+2} </math>. हालाँकि, प्रारंभिक मूल्य समस्या केवल एक मान प्रदान करती है, <math> y_0 = 1 </math>. इस समस्या को हल करने की एक संभावना का उपयोग करना है <math> y_1 </math> यूलर की विधि द्वारा दूसरे मान के रूप में गणना की गई। इस विकल्प के साथ, एडम्स-बैशफोर्थ विधि उत्पन्न होती है (चार अंकों तक पूर्णांकित):
इस विधि के लिए दो मानों <math> y_{n+1} </math> और <math> y_n </math> अगले मान <math> y_{n+2} </math> की गणना करने की आवश्यकता है, हालाँकि, प्रारंभिक मूल्य समस्या केवल एक मान <math> y_0 = 1 </math> प्रदान करती है। इस समस्या को हल करने की एक संभावना यूलर की विधि द्वारा गणना किए गए <math> y_1 </math> को दूसरे मान के रूप में उपयोग करना है। इस विकल्प के साथ, एडम्स-बैशफोर्थ विधि उत्पन्न होती है (चार अंकों तक पूर्णांकित):
<math display="block"> \begin{align}
<math display="block"> \begin{align}
   y_2 &= y_1 + \tfrac 3 2 hf(t_1, y_1) - \tfrac 1 2 hf(t_0, y_0) = 1.5 + \tfrac 3 2 \cdot \tfrac 1 2 \cdot 1.5 - \tfrac 1 2 \cdot \tfrac 1 2 \cdot 1 = 2.375, \\
   y_2 &= y_1 + \tfrac 3 2 hf(t_1, y_1) - \tfrac 1 2 hf(t_0, y_0) = 1.5 + \tfrac 3 2 \cdot \tfrac 1 2 \cdot 1.5 - \tfrac 1 2 \cdot \tfrac 1 2 \cdot 1 = 2.375, \\
Line 52: Line 52:
   y_4 &= y_3 + \tfrac 3 2 hf(t_3, y_3) - \tfrac 1 2 hf(t_2, y_2) = 3.7812 + \tfrac 3 2 \cdot \tfrac 1 2 \cdot 3.7812 - \tfrac 1 2 \cdot \tfrac 1 2 \cdot 2.375 = 6.0234.
   y_4 &= y_3 + \tfrac 3 2 hf(t_3, y_3) - \tfrac 1 2 hf(t_2, y_2) = 3.7812 + \tfrac 3 2 \cdot \tfrac 1 2 \cdot 3.7812 - \tfrac 1 2 \cdot \tfrac 1 2 \cdot 2.375 = 6.0234.
\end{align} </math>
\end{align} </math>
पर सटीक समाधान <math> t = t_4 = 2 </math> है <math> e^2 = 7.3891\ldots </math>, इसलिए दो-चरणीय एडम्स-बैशफोर्थ विधि यूलर की विधि से अधिक सटीक है। यदि चरण का आकार काफी छोटा है तो यह हमेशा मामला होता है।
<math> t = t_4 = 2 </math> पर सटीक समाधान <math> e^2 = 7.3891\ldots </math> है, इसलिए दो-चरणीय एडम्स-बैशफोर्थ विधि यूलर की विधि से अधिक सटीक है। यदि चरण का आकार काफी छोटा है तो यह हमेशा स्तिथि होती है।


== मल्टीस्टेप विधियों के परिवार ==
== बहुपदीय विधियों के समूह ==
रैखिक मल्टीस्टेप विधियों के तीन परिवार आमतौर पर उपयोग किए जाते हैं: एडम्स-बैशफोर्थ विधियां, एडम्स-मौल्टन विधियां, और पिछड़े भेदभाव सूत्र (बीडीएफ)।
रैखिक बहुपदीय विधियों के तीन समूह सामान्यतः उपयोग किए जाते हैं: एडम्स-बैशफोर्थ विधियां, एडम्स-मौल्टन विधियां, और पिछड़े भेदभाव सूत्र (बीडीएफ)।


=== एडम्स-बैशफोर्थ विधियाँ ===
=== एडम्स-बैशफोर्थ विधियाँ ===
एडम्स-बैशफोर्थ विधियाँ स्पष्ट विधियाँ हैं। गुणांक हैं <math> a_{s-1}=-1 </math> और <math> a_{s-2} = \cdots = a_0 = 0 </math>, जब <math> b_j </math> ऐसे चुना जाता है कि विधियों का क्रम s हो (यह विधियों को विशिष्ट रूप से निर्धारित करता है)।
एडम्स-बैशफोर्थ विधियाँ स्पष्ट विधियाँ हैं। <math> a_{s-1}=-1 </math> और <math> a_{s-2} = \cdots = a_0 = 0 </math> गुणांक हैं, जब <math> b_j </math> ऐसे चुना जाता है कि विधियों का क्रम s हो (यह विधियों को विशिष्ट रूप से निर्धारित करता है)।


एडम्स-बैशफोर्थ विधियाँ s = 1, 2, 3, 4, 5 के साथ हैं ({{harvnb|Hairer|Nørsett|Wanner|1993|loc=§III.1}}; {{harvnb|Butcher|2003|p=103}}):
एडम्स-बैशफोर्थ विधियाँ s = 1, 2, 3, 4, 5 के साथ हैं ({{harvnb|हेयरर|नॉरसेट|वानर|1993|loc=§III.1}}; {{harvnb|बुचर|2003|p=103}}):
<math display="block"> \begin{align}
<math display="block"> \begin{align}
y_{n+1} &= y_n + hf(t_n, y_n) , \qquad\text{(This is the Euler method)} \\
y_{n+1} &= y_n + hf(t_n, y_n) , \qquad\text{(This is the Euler method)} \\
Line 68: Line 68:
y_{n+5} &= y_{n+4} + h\left( \frac{1901}{720} f(t_{n+4}, y_{n+4}) - \frac{2774}{720} f(t_{n+3}, y_{n+3}) + \frac{2616}{720} f(t_{n+2}, y_{n+2}) - \frac{1274}{720} f(t_{n+1}, y_{n+1}) + \frac{251}{720} f(t_n, y_n) \right) .
y_{n+5} &= y_{n+4} + h\left( \frac{1901}{720} f(t_{n+4}, y_{n+4}) - \frac{2774}{720} f(t_{n+3}, y_{n+3}) + \frac{2616}{720} f(t_{n+2}, y_{n+2}) - \frac{1274}{720} f(t_{n+1}, y_{n+1}) + \frac{251}{720} f(t_n, y_n) \right) .
\end{align} </math>
\end{align} </math>
गुणांक <math> b_j </math> निम्नानुसार निर्धारित किया जा सकता है। घात का बहुपद p ज्ञात करने के लिए [[बहुपद प्रक्षेप]] का उपयोग करें <math> s-1 </math> ऐसा है कि
गुणांक <math> b_j </math> निम्नानुसार निर्धारित किया जा सकता है। <math> s-1 </math> घात का बहुपद p ज्ञात करने के लिए [[बहुपद प्रक्षेप]] का उपयोग करें, यह ऐसा है कि
<math display="block"> p(t_{n+i}) = f(t_{n+i}, y_{n+i}), \qquad \text{for } i=0,\ldots,s-1. </math>
<math display="block"> p(t_{n+i}) = f(t_{n+i}, y_{n+i}), \qquad \text{for } i=0,\ldots,s-1. </math>
बहुपद प्रक्षेप पैदावार के लिए [[लैग्रेंज बहुपद]]
बहुपद प्रक्षेप उपज के लिए [[लैग्रेंज बहुपद]]
<math display="block"> p(t) = \sum_{j=0}^{s-1} \frac{(-1)^{s-j-1}f(t_{n+j}, y_{n+j})}{j!(s-j-1)!h^{s-1}} \prod_{i=0 \atop i\ne j}^{s-1} (t-t_{n+i}). </math>
<math display="block"> p(t) = \sum_{j=0}^{s-1} \frac{(-1)^{s-j-1}f(t_{n+j}, y_{n+j})}{j!(s-j-1)!h^{s-1}} \prod_{i=0 \atop i\ne j}^{s-1} (t-t_{n+i}). </math>
बहुपद p स्थानीय रूप से अवकल समीकरण के दाएँ पक्ष का एक अच्छा सन्निकटन है <math> y' = f(t,y) </math> इसे हल करना है, इसलिए समीकरण पर विचार करें <math> y' = p(t) </math> बजाय। इस समीकरण को बिल्कुल हल किया जा सकता है; समाधान केवल p का अभिन्न अंग है। यह लेने का सुझाव देता है
बहुपद p स्थानीय रूप से अवकल समीकरण <math> y' = f(t,y) </math> के दाएँ पक्ष का एक अच्छा सन्निकटन इसे हल करना है, इसलिए समीकरण <math> y' = p(t) </math> स्थान पर विचार करें। इस समीकरण को बिल्कुल हल किया जा सकता है; समाधान केवल p का अभिन्न अंग है। यह निम्न लेने का सुझाव देता है
<math display="block"> y_{n+s} = y_{n+s-1} + \int_{t_{n+s-1}}^{t_{n+s}} p(t)\,\mathrm dt. </math>
<math display="block"> y_{n+s} = y_{n+s-1} + \int_{t_{n+s-1}}^{t_{n+s}} p(t)\,\mathrm dt. </math>
एडम्स-बैशफोर्थ विधि तब उत्पन्न होती है जब पी के लिए सूत्र प्रतिस्थापित किया जाता है। गुणांक <math> b_j </math> द्वारा दिया जाना निकला
एडम्स-बैशफोर्थ विधि तब उत्पन्न होती है जब p के लिए सूत्र प्रतिस्थापित किया जाता है। गुणांक <math> b_j </math> निम्न द्वारा दिए गए हैं
<math display="block"> b_{s-j-1} = \frac{(-1)^j}{j!(s-j-1)!} \int_0^1 \prod_{i=0 \atop i\ne j}^{s-1} (u+i) \,\mathrm du, \qquad \text{for } j=0,\ldots,s-1. </math>
<math display="block"> b_{s-j-1} = \frac{(-1)^j}{j!(s-j-1)!} \int_0^1 \prod_{i=0 \atop i\ne j}^{s-1} (u+i) \,\mathrm du, \qquad \text{for } j=0,\ldots,s-1. </math>
की जगह <math> f(t, y) </math> इसके इंटरपोलेंट पी द्वारा ऑर्डर एच की त्रुटि उत्पन्न होती है<sup>एस</sup>, और यह इस प्रकार है कि एस-स्टेप एडम्स-बैशफोर्थ विधि में वास्तव में ऑर्डर एस है {{harv|Iserles|1996|loc=§2.1}}
<math> f(t, y) </math> की जगह इसके इंटरपोलेंट पी द्वारा क्रम H<sup>s</sup> की त्रुटि उत्पन्न होती है, और यह इस प्रकार है कि एस-चरण एडम्स-बैशफोर्थ विधि में वास्तव में क्रम s {{harv|इसरल्स|1996|loc=§2.1}} है


एडम्स-बैशफोर्थ विधियों को [[जॉन काउच एडम्स]] द्वारा [[फ्रांसिस बैशफोर्थ]] के कारण केशिका क्रिया मॉडलिंग के अंतर समीकरण को हल करने के लिए डिजाइन किया गया था। {{harvtxt|Bashforth|1883}} ने उनके सिद्धांत और एडम्स की संख्यात्मक पद्धति को प्रकाशित किया {{harv|Goldstine|1977}}.
एडम्स-बैशफोर्थ विधियों को [[जॉन काउच एडम्स]] द्वारा [[फ्रांसिस बैशफोर्थ]] के कारण केशिका क्रिया प्रतिरूपण के अंतर समीकरण को हल करने के लिए अभिकल्पित किया गया था। {{harvtxt|बैशफोर्थ|1883}} ने उनके सिद्धांत और एडम्स की संख्यात्मक पद्धति {{harv|गोल्डस्टाइन|1977}} को प्रकाशित किया।


=== एडम्स-मौलटन विधियाँ ===
=== एडम्स-मौलटन विधियाँ ===
एडम्स-मौलटन विधियाँ एडम्स-बैशफोर्थ विधियों के समान हैं, उनमें भी हैं <math> a_{s-1} = -1 </math> और <math> a_{s-2} = \cdots = a_0 = 0 </math>. उच्चतम संभव क्रम प्राप्त करने के लिए फिर से b गुणांक को चुना जाता है। हालाँकि, एडम्स-मौल्टन विधियाँ अंतर्निहित विधियाँ हैं। उस प्रतिबंध को हटाकर <math> b_s = 0 </math>, एक एस-स्टेप एडम्स-मौलटन विधि ऑर्डर तक पहुंच सकती है <math> s+1 </math>, जबकि एस-स्टेप एडम्स-बैशफोर्थ विधियों में केवल ऑर्डर एस है।
एडम्स-मौलटन विधियाँ एडम्स-बैशफोर्थ विधियों के समान हैं, उनमें <math> a_{s-1} = -1 </math> और <math> a_{s-2} = \cdots = a_0 = 0 </math> भी हैं। उच्चतम संभव क्रम प्राप्त करने के लिए फिर से b गुणांक को चुना जाता है। हालाँकि, एडम्स-मौल्टन विधियाँ अंतर्निहित विधियाँ हैं। उस प्रतिबंध <math> b_s = 0 </math> को हटाकर, एक एस-चरण एडम्स-मौलटन विधि क्रम <math> s+1 </math> तक पहुंच सकती है, जबकि एस-चरण एडम्स-बैशफोर्थ विधियों में केवल क्रम एस है।


s = 0, 1, 2, 3, 4 के साथ एडम्स-मौलटन विधियाँ हैं ({{harvnb|Hairer|Nørsett|Wanner|1993|loc=§III.1}}; {{harvnb|Quarteroni|Sacco|Saleri|2000}}) सूचीबद्ध है, जहां पहले दो तरीके क्रमशः [[बैकवर्ड यूलर विधि]] और ट्रेपेज़ॉइडल नियम (अंतर समीकरण) हैं:
s = 0, 1, 2, 3, 4 के साथ एडम्स-मौलटन विधियाँ ({{harvnb|हेयरर|नॉरसेट|वानर|1993|loc=§III.1}}; {{harvnb|क्वार्टरोनी|सैको|सालेरी|2000}}) सूचीबद्ध हैं, जहां पहले दो तरीके क्रमशः [[बैकवर्ड यूलर विधि]] और ट्रेपेज़ॉइडल नियम (अंतर समीकरण) हैं:
<math display="block"> \begin{align}
<math display="block"> \begin{align}
y_{n} &= y_{n-1} + h f(t_{n},y_{n}), \\
y_{n} &= y_{n-1} + h f(t_{n},y_{n}), \\
Line 91: Line 91:
y_{n+4} &= y_{n+3} + h \left( \frac{251}{720} f(t_{n+4},y_{n+4}) + \frac{646}{720} f(t_{n+3},y_{n+3}) - \frac{264}{720} f(t_{n+2},y_{n+2}) + \frac{106}{720} f(t_{n+1},y_{n+1}) - \frac{19}{720} f(t_n,y_n) \right) .
y_{n+4} &= y_{n+3} + h \left( \frac{251}{720} f(t_{n+4},y_{n+4}) + \frac{646}{720} f(t_{n+3},y_{n+3}) - \frac{264}{720} f(t_{n+2},y_{n+2}) + \frac{106}{720} f(t_{n+1},y_{n+1}) - \frac{19}{720} f(t_n,y_n) \right) .
\end{align} </math>
\end{align} </math>
एडम्स-मौलटन पद्धति की व्युत्पत्ति एडम्स-बैशफोर्थ पद्धति के समान है; हालाँकि, प्रक्षेप बहुपद न केवल बिंदुओं का उपयोग करता है <math>t_{n-1},\dots, t_{n-s} </math>, जैसा कि ऊपर है, लेकिन यह भी <math> t_n </math>. गुणांक द्वारा दिए गए हैं
एडम्स-मौलटन पद्धति की व्युत्पत्ति एडम्स-बैशफोर्थ पद्धति के समान है; हालाँकि, प्रक्षेप बहुपद न केवल ऊपर दिए गए बिंदुओं <math>t_{n-1},\dots, t_{n-s} </math>t का उपयोग करता है, बल्कि <math> t_n </math> का भी उपयोग करता है। गुणांक निम्न द्वारा दिए गए हैं
<math display="block"> b_{s-j} = \frac{(-1)^j}{j!(s-j)!} \int_0^1 \prod_{i=0 \atop i\ne j}^{s} (u+i-1) \,\mathrm du, \qquad \text{for } j=0,\ldots,s. </math>
<math display="block"> b_{s-j} = \frac{(-1)^j}{j!(s-j)!} \int_0^1 \prod_{i=0 \atop i\ne j}^{s} (u+i-1) \,\mathrm du, \qquad \text{for } j=0,\ldots,s. </math>
एडम्स-बैशफोर्थ विधियों की तरह, एडम्स-मौल्टन विधियाँ पूरी तरह से जॉन काउच एडम्स के कारण हैं। [[वन रे मौलटन]] का नाम इन विधियों के साथ जुड़ गया क्योंकि उन्हें एहसास हुआ कि इन्हें एडम्स-बैशफोर्थ विधियों के साथ मिलकर [[भविष्यवक्ता-सुधारक विधि]]|प्रीडिक्टर-करेक्टर जोड़ी के रूप में इस्तेमाल किया जा सकता है। {{harv|Moulton|1926}}; {{harvtxt|Milne|1926}} का भी यही विचार था. एडम्स ने अंतर्निहित समीकरण को हल करने के लिए न्यूटन की विधि का उपयोग किया {{harv|Hairer|Nørsett|Wanner|1993|loc=§III.1}}.
एडम्स-बैशफोर्थ विधियों की तरह, एडम्स-मौल्टन विधियाँ पूरी तरह से जॉन काउच एडम्स के कारण हैं। [[वन रे मौलटन]] का नाम इन विधियों के साथ जुड़ गया क्योंकि उन्हें एहसास हुआ कि इन्हें एडम्स-बैशफोर्थ विधियों के साथ मिलकर [[भविष्यवक्ता-सुधारक विधि]] जोड़ी के रूप में इस्तेमाल किया जा सकता है। {{harv|मौलटन|1926}}; {{harvtxt|मिलन|1926}} का भी यही विचार था। एडम्स ने अंतर्निहित समीकरण को हल करने के लिए न्यूटन की विधि {{harv|हेयरर|नॉरसेट|वानर|1993|loc=§III.1}} का उपयोग किया।


=== पिछड़ा विभेदन सूत्र (पीडीएफ) ===
=== पिछड़ा विभेदन सूत्र (पीडीएफ) ===
{{main|Backward differentiation formula}}
{{main|पिछड़ा विभेदन सूत्र}}
बीडीएफ विधियां अंतर्निहित विधियां हैं <math> b_{s-1} = \cdots = b_0 = 0 </math> और अन्य गुणांक इस प्रकार चुने गए कि विधि क्रम s (अधिकतम संभव) प्राप्त कर ले। इन विधियों का प्रयोग विशेष रूप से कठोर समीकरणों के समाधान के लिए किया जाता है।
 
बीडीएफ विधियां अंतर्निहित विधियां <math> b_{s-1} = \cdots = b_0 = 0 </math> हैं और अन्य गुणांक इस प्रकार चुने गए कि विधि क्रम s (अधिकतम संभव) प्राप्त कर ले। इन विधियों का प्रयोग विशेष रूप से कठोर समीकरणों के समाधान के लिए किया जाता है।


== विश्लेषण ==
== विश्लेषण ==
रैखिक मल्टीस्टेप विधियों के विश्लेषण में केंद्रीय अवधारणाएं, और वास्तव में अंतर समीकरणों के लिए किसी भी संख्यात्मक विधि, संख्यात्मक साधारण अंतर समीकरण # विश्लेषण | अभिसरण, क्रम और स्थिरता हैं।
रैखिक बहुपदीय विधियों के विश्लेषण में केंद्रीय अवधारणाएं, और वास्तव में अंतर समीकरणों के लिए किसी भी संख्यात्मक विधि, संख्यात्मक साधारण अंतर समीकरण अभिसरण, क्रम और स्थिरता हैं।


=== संगति और क्रम ===
=== संगति और क्रम ===
Line 108: Line 109:
& \qquad {} = h \bigl( b_s f(t_{n+s},y_{n+s}) + b_{s-1} f(t_{n+s-1},y_{n+s-1}) + \cdots + b_0 f(t_n,y_n) \bigr),
& \qquad {} = h \bigl( b_s f(t_{n+s},y_{n+s}) + b_{s-1} f(t_{n+s-1},y_{n+s-1}) + \cdots + b_0 f(t_n,y_n) \bigr),
\end{align} </math>
\end{align} </math>
अंतर समीकरण का एक अच्छा सन्निकटन <math> y' = f(t,y) </math>? अधिक सटीक रूप से, एक मल्टीस्टेप विधि सुसंगत होती है यदि स्थानीय ट्रंकेशन त्रुटि चरण आकार h की तुलना में तेजी से शून्य हो जाती है क्योंकि h शून्य पर चला जाता है, जहां स्थानीय ट्रंकेशन त्रुटि को परिणाम के बीच अंतर के रूप में परिभाषित किया जाता है <math>y_{n+s}</math> विधि का, यह मानते हुए कि पिछले सभी मान <math>y_{n+s-1}, \ldots, y_n</math> सटीक हैं, और समय पर समीकरण का सटीक समाधान हैं <math>t_{n+s}</math>. [[टेलर श्रृंखला]] का उपयोग करते हुए एक गणना से पता चलता है कि एक रैखिक मल्टीस्टेप विधि सुसंगत है यदि और केवल यदि
अंतर समीकरण का एक अच्छा सन्निकटन <math> y' = f(t,y) </math> है ? अधिक सटीक रूप से, एक बहुपदीय विधि सुसंगत होती है यदि स्थानीय खंडन त्रुटि चरण आकार h की तुलना में तीव्रता से शून्य हो जाती है क्योंकि h शून्य पर चला जाता है, जहां स्थानीय खंडन त्रुटि को परिणाम <math>y_{n+s}</math> के बीच अंतर के रूप में परिभाषित किया जाता है, यह मानते हुए कि पिछले सभी मान <math>y_{n+s-1}, \ldots, y_n</math> सटीक हैं, और <math>t_{n+s}</math> समय पर समीकरण का सटीक समाधान हैं। [[टेलर श्रृंखला]] का उपयोग करते हुए एक गणना से पता चलता है कि एक रैखिक बहुपदीय विधि सुसंगत है यदि और केवल यदि
<math display="block"> \sum_{k=0}^{s-1} a_k = -1 \quad\text{and}\quad \sum_{k=0}^s b_k = s + \sum_{k=0}^{s-1} k a_k. </math>
<math display="block"> \sum_{k=0}^{s-1} a_k = -1 \quad\text{and}\quad \sum_{k=0}^s b_k = s + \sum_{k=0}^{s-1} k a_k. </math>
ऊपर उल्लिखित सभी विधियाँ सुसंगत हैं {{harv|Hairer|Nørsett|Wanner|1993|loc=§III.2}}.
ऊपर उल्लिखित सभी विधियाँ {{harv|हेयरर|नॉरसेट|वानर|1993|loc=§III.2}} सुसंगत हैं।


यदि विधि सुसंगत है, तो अगला प्रश्न यह है कि संख्यात्मक विधि को परिभाषित करने वाला अंतर समीकरण कितनी अच्छी तरह अंतर समीकरण का अनुमान लगाता है। यदि स्थानीय त्रुटि क्रम की है तो मल्टीस्टेप विधि को ऑर्डर पी कहा जाता है <math>O(h^{p+1})</math> जैसे ही h शून्य पर जाता है। यह विधियों के गुणांकों पर निम्नलिखित शर्त के बराबर है:
यदि विधि सुसंगत है, तो अगला प्रश्न यह है कि संख्यात्मक विधि को परिभाषित करने वाला अंतर समीकरण कितनी अच्छी तरह अंतर समीकरण का अनुमान लगाता है। यदि स्थानीय त्रुटि क्रम की है तो बहुपदीय विधि को क्रम पी कहा जाता है <math>O(h^{p+1})</math> जैसे ही h शून्य पर जाता है। यह विधियों के गुणांकों पर निम्नलिखित परिस्थिति के बराबर है:
<math display="block"> \sum_{k=0}^{s-1} a_k = -1 \quad\text{and}\quad q \sum_{k=0}^s k^{q-1} b_k = s^q + \sum_{k=0}^{s-1} k^q a_k \text{ for } q=1,\ldots,p. </math>
<math display="block"> \sum_{k=0}^{s-1} a_k = -1 \quad\text{and}\quad q \sum_{k=0}^s k^{q-1} b_k = s^q + \sum_{k=0}^{s-1} k^q a_k \text{ for } q=1,\ldots,p. </math>
एस-स्टेप एडम्स-बैशफोर्थ विधि में ऑर्डर एस है, जबकि एस-स्टेप एडम्स-मौल्टन विधि में ऑर्डर है <math>s+1</math> {{harv|Hairer|Nørsett|Wanner|1993|loc=§III.2}}.
एस-चरण एडम्स-बैशफोर्थ विधि में क्रम एस है, जबकि एस-चरण एडम्स-मौल्टन विधि में क्रम <math>s+1</math> {{harv|हेयरर|नॉरसेट|वानर|1993|loc=§III.2}} है।


ये स्थितियां अक्सर विशिष्ट बहुपदों का उपयोग करके तैयार की जाती हैं
ये स्थितियां प्रायः विशिष्ट बहुपदों का उपयोग करके तैयार की जाती हैं
<math display="block"> \rho(z) = z^s + \sum_{k=0}^{s-1} a_k z^k \quad\text{and}\quad \sigma(z) = \sum_{k=0}^s b_k z^k. </math>
<math display="block"> \rho(z) = z^s + \sum_{k=0}^{s-1} a_k z^k \quad\text{and}\quad \sigma(z) = \sum_{k=0}^s b_k z^k. </math>
इन बहुपदों के संदर्भ में, क्रम p रखने की विधि के लिए उपरोक्त शर्त बन जाती है
इन बहुपदों के संदर्भ में, क्रम p रखने की विधि के लिए उपरोक्त परिस्थिति बन जाती है
<math display="block"> \rho(e^h) - h\sigma(e^h) = O(h^{p+1}) \quad \text{as } h\to 0. </math>
<math display="block"> \rho(e^h) - h\sigma(e^h) = O(h^{p+1}) \quad \text{as } h\to 0. </math>
विशेष रूप से, विधि सुसंगत है यदि इसमें कम से कम एक ऑर्डर है, जो कि मामला है <math>\rho(1)=0</math> और <math>\rho'(1)=\sigma(1)</math>.
विशेष रूप से, विधि सुसंगत है यदि इसमें कम से कम एक क्रम है, जो कि स्तिथि  <math>\rho(1)=0</math> और <math>\rho'(1)=\sigma(1)</math> है।


===स्थिरता और अभिसरण ===
===स्थिरता और अभिसरण ===
<!-- [[Zero stability]] redirects here -->
एक-चरणीय विधि का संख्यात्मक समाधान प्रारंभिक स्थिति <math> y_0 </math> पर निर्भर करता है, लेकिन एस-चरण विधि का संख्यात्मक समाधान सभी प्रारम्भिक मानों <math> y_0, y_1, \ldots, y_{s-1} </math> पर निर्भर करता है। इस प्रकार यह रुचि का विषय है कि क्या प्रारंभिक मूल्यों में गड़बड़ी के संबंध में संख्यात्मक समाधान स्थिर है। एक रैखिक बहुपदीय विधि किसी निश्चित समय अंतराल पर एक निश्चित अंतर समीकरण के लिए शून्य-स्थिर है, यदि आकार ε के प्रारम्भिक मूल्यों में गड़बड़ी के कारण उस समय अंतराल पर संख्यात्मक समाधान K के कुछ मूल्य के लिए Kε से अधिक नहीं बदलता है जो चरण आकार h पर निर्भर नहीं करता है। इसे शून्य-स्थिरता कहा जाता है क्योंकि यह अंतर समीकरण <math> y' = 0 </math> {{harv|सुली|मेयर्स|2003|p=332}} की स्थिति की जांच करने के लिए पर्याप्त है।
एक-चरणीय विधि का संख्यात्मक समाधान प्रारंभिक स्थिति पर निर्भर करता है <math> y_0 </math>, लेकिन एस-स्टेप विधि का संख्यात्मक समाधान सभी शुरुआती मानों पर निर्भर करता है, <math> y_0, y_1, \ldots, y_{s-1} </math>. इस प्रकार यह रुचि का विषय है कि क्या प्रारंभिक मूल्यों में गड़बड़ी के संबंध में संख्यात्मक समाधान स्थिर है। एक रैखिक मल्टीस्टेप विधि किसी निश्चित समय अंतराल पर एक निश्चित अंतर समीकरण के लिए शून्य स्थिरता | शून्य-स्थिर है, यदि आकार ε के शुरुआती मूल्यों में गड़बड़ी के कारण उस समय अंतराल पर संख्यात्मक समाधान K के कुछ मूल्य के लिए Kε से अधिक नहीं बदलता है जो चरण आकार h पर निर्भर नहीं करता है। इसे शून्य-स्थिरता कहा जाता है क्योंकि यह अंतर समीकरण की स्थिति की जांच करने के लिए पर्याप्त है <math> y' = 0 </math> {{harv|Süli|Mayers|2003|p=332}}.


यदि विशिष्ट बहुपद ρ के मूलों का मापांक 1 से कम या उसके बराबर है और मापांक 1 के मूल गुणनफल 1 के हैं, तो हम कहते हैं कि मूल स्थिति संतुष्ट है। एक रैखिक मल्टीस्टेप विधि शून्य-स्थिर है यदि और केवल तभी जब मूल स्थिति संतुष्ट हो {{harv|Süli|Mayers|2003|p=335}}.
यदि विशिष्ट बहुपद ρ के मूलों का मापांक 1 से कम या उसके बराबर है और मापांक 1 के मूल गुणनफल 1 के हैं, तो हम कहते हैं कि मूल स्थिति संतुष्ट है। एक रैखिक बहुपदीय विधि शून्य-स्थिर है यदि और केवल तभी जब मूल स्थिति {{harv|सुली|मेयर्स|2003|p=335}} संतुष्ट हो।


अब मान लीजिए कि एक सुसंगत रैखिक मल्टीस्टेप विधि को पर्याप्त रूप से सुचारू अंतर समीकरण और शुरुआती मूल्यों पर लागू किया जाता है <math> y_1, \ldots, y_{s-1}</math> सभी प्रारंभिक मूल्य पर एकत्रित होते हैं <math> y_0 </math> जैसा <math> h \to 0 </math>. फिर, संख्यात्मक समाधान सटीक समाधान में परिवर्तित हो जाता है <math> h \to 0 </math> यदि और केवल यदि विधि शून्य-स्थिर है। इस परिणाम को डाहलक्विस्ट तुल्यता प्रमेय के रूप में जाना जाता है, जिसका नाम [[जर्मुंड डहलक्विस्ट]] के नाम पर रखा गया है; यह प्रमेय आत्मा में [[परिमित अंतर विधि]]यों के लिए [[लैक्स तुल्यता प्रमेय]] के समान है। इसके अलावा, यदि विधि में ऑर्डर पी है, तो [[वैश्विक ट्रंकेशन त्रुटि]] (एक निश्चित समय पर संख्यात्मक समाधान और सटीक समाधान के बीच का अंतर) है <math> O(h^p) </math> {{harv|Süli|Mayers|2003|p=340}}.
अब मान लीजिए कि एक सुसंगत रैखिक मल्टीस्टेप विधि को पर्याप्त रूप से सुचारू अंतर समीकरण पर लागू किया जाता है और प्रारंभिक मान <math> y_1, \ldots, y_{s-1}</math> सभी प्रारंभिक मान <math> y_0 </math> में <math> h \to 0 </math> के रूप में परिवर्तित हो जाते हैं। फिर, संख्यात्मक समाधान सटीक समाधान <math> h \to 0 </math> में परिवर्तित हो जाता है, यदि और केवल यदि विधि शून्य-स्थिर है। इस परिणाम को डाहलक्विस्ट तुल्यता प्रमेय के रूप में जाना जाता है, जिसका नाम [[जर्मुंड डहलक्विस्ट]] के नाम पर रखा गया है; यह प्रमेय तत्परता में [[परिमित अंतर विधि]]यों के लिए [[लैक्स तुल्यता प्रमेय]] के समान है। इसके अतिरिक्त, यदि विधि में क्रम पी है, तो [[वैश्विक ट्रंकेशन त्रुटि|वैश्विक खंडन त्रुटि]] {{harv|सुली|मेयर्स|2003|p=340}}(एक निश्चित समय पर संख्यात्मक समाधान और सटीक समाधान के बीच का अंतर) <math> O(h^p) </math> है।


इसके अलावा, यदि विधि अभिसरण है, तो विधि को दृढ़ता से स्थिर कहा जाता है <math>z=1</math> मापांक 1 का एकमात्र मूल है। यदि यह अभिसरण है और मापांक 1 की सभी जड़ें दोहराई नहीं जाती हैं, लेकिन ऐसे एक से अधिक मूल हैं, तो इसे अपेक्षाकृत स्थिर कहा जाता है। ध्यान दें कि विधि को अभिसरण करने के लिए 1 को मूल होना चाहिए; इस प्रकार अभिसरण विधियाँ हमेशा इन दोनों में से एक होती हैं।
इसके अतिरिक्त, यदि विधि अभिसरण है, तो विधि को दृढ़ता से स्थिर कहा जाता है, <math>z=1</math> मापांक 1 का एकमात्र मूल है। यदि यह अभिसरण है और मापांक 1 की सभी घात दोहराई नहीं जाती हैं, लेकिन ऐसे एक से अधिक मूल हैं, तो इसे अपेक्षाकृत स्थिर कहा जाता है। ध्यान दें कि विधि को अभिसरण करने के लिए 1 को मूल होना चाहिए; इस प्रकार अभिसरण विधियाँ हमेशा इन दोनों में से एक होती हैं।


कठोर समीकरणों पर रैखिक मल्टीस्टेप विधियों के प्रदर्शन का आकलन करने के लिए, रैखिक परीक्षण समीकरण y' = λy पर विचार करें। चरण आकार h के साथ इस अंतर समीकरण पर लागू एक मल्टीस्टेप विधि विशेषता बहुपद के साथ एक रैखिक [[पुनरावृत्ति संबंध]] उत्पन्न करती है
कठोर समीकरणों पर रैखिक बहुपदीय विधियों के प्रदर्शन का आकलन करने के लिए, रैखिक परीक्षण समीकरण y' = λy पर विचार करें। चरण आकार h के साथ इस अंतर समीकरण पर लागू एक बहुपदीय विधि विशेषता बहुपद के साथ एक रैखिक [[पुनरावृत्ति संबंध]] उत्पन्न करती है
<math display="block"> \pi(z; h\lambda) = (1 - h\lambda\beta_s) z^s + \sum_{k=0}^{s-1} (\alpha_k - h\lambda\beta_k) z^k = \rho(z) - h\lambda\sigma(z). </math>
<math display="block"> \pi(z; h\lambda) = (1 - h\lambda\beta_s) z^s + \sum_{k=0}^{s-1} (\alpha_k - h\lambda\beta_k) z^k = \rho(z) - h\lambda\sigma(z). </math>
इस बहुपद को मल्टीस्टेप विधि का स्थिरता बहुपद कहा जाता है। यदि इसकी सभी जड़ों का मापांक एक से कम है तो मल्टीस्टेप विधि का संख्यात्मक समाधान शून्य में परिवर्तित हो जाएगा और मल्टीस्टेप विधि को hλ के उस मान के लिए बिल्कुल स्थिर कहा जाता है। विधि को ए-स्थिर कहा जाता है यदि यह नकारात्मक वास्तविक भाग वाले सभी hλ के लिए बिल्कुल स्थिर है। पूर्ण स्थिरता का क्षेत्र सभी hλ का समुच्चय है जिसके लिए मल्टीस्टेप विधि बिल्कुल स्थिर है {{harv|Süli|Mayers|2003|pp=347 & 348}}. अधिक विवरण के लिए, कठोर समीकरण#मल्टीस्टेप विधियों पर अनुभाग देखें।
इस बहुपद को बहुपदीय विधि का स्थिरता बहुपद कहा जाता है। यदि इसकी सभी घात का मापांक एक से कम है तो बहुपदीय विधि का संख्यात्मक समाधान शून्य में परिवर्तित हो जाएगा और बहुपदीय विधि को hλ के उस मान के लिए बिल्कुल स्थिर कहा जाता है। विधि को ए-स्थिर कहा जाता है यदि यह नकारात्मक वास्तविक भाग वाले सभी hλ के लिए बिल्कुल स्थिर है। पूर्ण स्थिरता का क्षेत्र सभी hλ का समुच्चय है जिसके लिए बहुपदीय विधि बिल्कुल स्थिर है {{harv|सुली|मेयर्स|2003|pp=347 & 348}}अधिक विवरण के लिए, कठोर समीकरण बहुपदीय विधियों पर अनुभाग देखें।


=== उदाहरण ===
=== उदाहरण ===
एडम्स-बैशफोर्थ तीन-चरणीय विधि पर विचार करें
एडम्स-बैशफोर्थ तीन-चरणीय विधि पर विचार करें<math display="block">y_{n+3} = y_{n+2} + h\left( {23\over 12} f(t_{n+2}, y_{n+2}) - {4 \over 3} f(t_{n+1}, y_{n+1}) + {5\over 12}f(t_{n}, y_{n})\right).</math>
<!-- save expression that computes y_n+1 rather than y_n+s
<math display="block">y_{n+1} = y_n + h\left( {23\over 12} f(t_n, y_n) - {16 \over 12} f(t_{n-1}, y_{n-1}) + {5\over 12}f(t_{n-2}, y_{n-2})\right).</math>
-->
<math display="block">y_{n+3} = y_{n+2} + h\left( {23\over 12} f(t_{n+2}, y_{n+2}) - {4 \over 3} f(t_{n+1}, y_{n+1}) + {5\over 12}f(t_{n}, y_{n})\right).</math>
इस प्रकार एक अभिलक्षणिक बहुपद है
इस प्रकार एक अभिलक्षणिक बहुपद है
<math display="block">\rho(z) = z^3-z^2 = z^2(z-1)</math>
<math display="block">\rho(z) = z^3-z^2 = z^2(z-1)</math>
जिसकी जड़ें हैं <math>z=0, 1</math>, और उपरोक्त शर्तें पूरी होती हैं। जैसा <math>z=1</math> मापांक 1 का एकमात्र मूल है, विधि अत्यधिक स्थिर है।
जिसकी घात <math>z=0, 1</math> हैं, और उपरोक्त स्तिथियाँ पूरी होती हैं। जैसे <math>z=1</math> मापांक 1 का एकमात्र मूल है, विधि अत्यधिक स्थिर है।


अन्य विशेषता बहुपद है
अन्य विशेषता बहुपद निम्न है
<math display="block">\sigma(z) = \frac{23}{12} z^2 - \frac{4}{3} z + \frac{5}{12} </math>
<math display="block">\sigma(z) = \frac{23}{12} z^2 - \frac{4}{3} z + \frac{5}{12} </math>




==पहली और दूसरी डहलक्विस्ट बाधाएँ==
==पहली और दूसरी डहलक्विस्ट बाधाएँ==
ये दो परिणाम जर्मुंड डहलक्विस्ट द्वारा सिद्ध किए गए थे और अभिसरण के क्रम के लिए और एक रैखिक मल्टीस्टेप विधि के कठोर समीकरण # ए-स्थिरता | ए-स्थिरता के लिए एक महत्वपूर्ण सीमा का प्रतिनिधित्व करते हैं। पहला डहलक्विस्ट बैरियर सिद्ध हुआ था {{harvtxt|Dahlquist|1956}} और दूसरे में {{harvtxt|Dahlquist|1963}}.
ये दो परिणाम जर्मुंड डहलक्विस्ट द्वारा सिद्ध किए गए थे और अभिसरण के क्रम के लिए और एक रैखिक बहुपदीय विधि के कठोर समीकरण ए-स्थिरता के लिए एक महत्वपूर्ण सीमा का प्रतिनिधित्व करते हैं। पहला डहलक्विस्ट अवरोध {{harvtxt|डहलक्विस्ट|1956}} और दूसरे में {{harvtxt|डहलक्विस्ट|1963}} सिद्ध हुआ था।


===पहला डहलक्विस्ट बैरियर===
===पहला डहलक्विस्ट अवरोध===


पहला डहलक्विस्ट बैरियर बताता है कि एक शून्य-स्थिर और रैखिक q-स्टेप मल्टीस्टेप विधि q + 1 से अधिक अभिसरण का क्रम प्राप्त नहीं कर सकती है यदि q विषम है और यदि q सम है तो q + 2 से अधिक है। यदि विधि भी स्पष्ट है, तो यह q से अधिक ऑर्डर प्राप्त नहीं कर सकती है {{harv|Hairer|Nørsett|Wanner|1993|loc=Thm III.3.5}}.
पहला डहलक्विस्ट अवरोध बताता है कि एक शून्य-स्थिर और रैखिक q-चरण बहुपदीय विधि q + 1 से अधिक अभिसरण का क्रम प्राप्त नहीं कर सकती है यदि q विषम है और यदि q सम है तो q + 2 से अधिक है। यदि विधि भी स्पष्ट है, तो यह q से अधिक क्रम प्राप्त नहीं कर सकती है {{harv|हेयरर|नॉरसेट|वानर|1993|loc=Thm III.3.5}}


===दूसरा डहलक्विस्ट अवरोध===
===दूसरा डहलक्विस्ट अवरोध===


दूसरा डहलक्विस्ट बैरियर बताता है कि कोई भी स्पष्ट रैखिक मल्टीस्टेप विधियां कठोर समीकरण#ए-स्थिरता|ए-स्थिर नहीं हैं। इसके अलावा, एक (अंतर्निहित) ए-स्थिर रैखिक मल्टीस्टेप विधि का अधिकतम क्रम 2 है। क्रम 2 के ए-स्थिर रैखिक मल्टीस्टेप तरीकों में, ट्रैपेज़ॉइडल नियम में सबसे छोटी त्रुटि स्थिरांक है {{harv|Dahlquist|1963|loc=Thm 2.1 and 2.2}}.
दूसरा डहलक्विस्ट अवरोध बताता है कि कोई भी स्पष्ट रैखिक बहुपदीय विधियां कठोर समीकरण ए-स्थिर नहीं हैं। इसके अतिरिक्त, एक (अंतर्निहित) ए-स्थिर रैखिक बहुपदीय विधि का अधिकतम क्रम 2 है। क्रम 2 के ए-स्थिर रैखिक बहुपदीय तरीकों में, समलंबी नियम में सबसे छोटी त्रुटि स्थिरांक है {{harv|डहलक्विस्ट|1963|loc=टीएचएम 2.1 and 2.2}}.


== यह भी देखें ==
== यह भी देखें ==
*[[डिजिटल ऊर्जा लाभ]]
*[[डिजिटल ऊर्जा लाभ|अंकीय ऊर्जा लाभ]]


== संदर्भ ==
== संदर्भ ==
* {{citation | first1 = Francis | last1 = Bashforth | year = 1883 | title = An Attempt to test the Theories of Capillary Action by comparing the theoretical and measured forms of drops of fluid. With an explanation of the method of integration employed in constructing the tables which give the theoretical forms of such drops, by J. C. Adams | location = Cambridge }}.
* {{citation | first1 = Francis | last1 = Bashforth | year = 1883 | title = An Attempt to test the Theories of Capillary Action by comparing the theoretical and measured forms of drops of fluid. With an explanation of the method of integration employed in constructing the tables which give the theoretical forms of such drops, by J. C. Adams | location = Cambridge }}.
* {{Citation | last1=Butcher | first1=John C. | author1-link = John C. Butcher | title=Numerical Methods for Ordinary Differential Equations | publisher=John Wiley | isbn=978-0-471-96758-3 | year=2003}}.
* {{Citation | last1=Butcher | first1=John C. | author1-link = John C. Butcher | title=Numerical Methods for Ordinary Differential Equations | publisher=John Wiley | isbn=978-0-471-96758-3 | year=2003}}.
* {{Citation | last1=Dahlquist | first1=Germund | author1-link=Germund Dahlquist | title=Convergence and stability in the numerical integration of ordinary differential equations | year=1956 | journal=Mathematica Scandinavica | volume=4 | pages=33–53| doi=10.7146/math.scand.a-10454 | doi-access=free }}.
* {{Citation | last1=Dahlquist | first1=Germund | author1-link=Germund Dahlquist | title=साधारण अंतर समीकरणों के संख्यात्मक एकीकरण में अभिसरण और स्थिरता | year=1956 | journal=Mathematica Scandinavica | volume=4 | pages=33–53| doi=10.7146/math.scand.a-10454 | doi-access=free }}.
* {{Citation | last1=Dahlquist | first1=Germund | author1-link=Germund Dahlquist | title=A special stability problem for linear multistep methods | doi=10.1007/BF01963532 | year=1963 | journal=BIT | issn=0006-3835 | volume=3 | pages=27–43| s2cid=120241743 }}.
* {{Citation | last1=Dahlquist | first1=Germund | author1-link=Germund Dahlquist | title=A special stability problem for linear multistep methods | doi=10.1007/BF01963532 | year=1963 | journal=BIT | issn=0006-3835 | volume=3 | pages=27–43| s2cid=120241743 }}.
* {{citation | first1 = Herman H. | last1 = Goldstine | author1-link = Herman Goldstine | year = 1977 | title = A History of Numerical Analysis from the 16th through the 19th Century | publisher = Springer-Verlag | location = New York | isbn = 978-0-387-90277-7 }}.
* {{citation | first1 = Herman H. | last1 = Goldstine | author1-link = Herman Goldstine | year = 1977 | title = A History of Numerical Analysis from the 16th through the 19th Century | publisher = Springer-Verlag | location = New York | isbn = 978-0-387-90277-7 }}.

Revision as of 12:49, 26 July 2023

संख्यात्मक साधारण अंतर समीकरण के लिए रैखिक बहुपदीय विधियों का उपयोग किया जाता है। वैचारिक रूप से, एक संख्यात्मक विधि एक प्रारंभिक बिंदु से प्रारम्भ होती है और फिर अगले समाधान बिंदु को खोजने के लिए समय में एक छोटा कदम आगे बढ़ाती है। समाधान निकालने के लिए प्रक्रिया बाद के चरणों के साथ जारी रहती है। एकल-चरण विधियाँ (जैसे यूलर की विधि) वर्तमान मूल्य निर्धारित करने के लिए केवल एक पिछले बिंदु और उसके व्युत्पन्न को संदर्भित करती हैं। रंज-कुट्टा जैसी विधियां उच्च क्रम विधि प्राप्त करने के लिए कुछ मध्यवर्ती कदम (उदाहरण के लिए, आधा कदम) लेती हैं, लेकिन फिर दूसरा कदम उठाने से पहले सभी पिछली जानकारी को त्याग देती हैं। बहुपदीय विधियाँ पिछले चरणों की जानकारी को त्यागने के स्थान पर उसे बनाए रखने और उसका उपयोग करके दक्षता प्राप्त करने का प्रयास करती हैं। नतीजतन, बहुपदीय विधियां कई पिछले बिंदुओं और व्युत्पन्न मूल्यों को संदर्भित करती हैं। रैखिक बहुपदीय विधियों की स्तिथि में, पिछले बिंदुओं और व्युत्पन्न मूल्यों के एक रैखिक संयोजन का उपयोग किया जाता है।

परिभाषाएँ

साधारण अंतर समीकरणों के लिए संख्यात्मक विधियाँ विधि की प्रारंभिक मान समस्या का अनुमानित समाधान करती हैं

परिणाम के मूल्य के लिए अनुमान अलग-अलग समय पर है:
जहाँ समय चरण है (कभी-कभी इसे कहा जाता है) और एक पूर्णांक है।

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

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

कोई भी स्पष्ट और अंतर्निहित तरीकों के बीच अंतर कर सकता है। अगर , तो विधि को स्पष्ट कहा जाता है, क्योंकि सूत्र सीधे गणना कर सकता है। अगर तो विधि को अंतर्निहित कहा जाता है, क्योंकि इसका मान के मूल्य पर निर्भर करता है, और समीकरण को हल किया जाना चाहिए। अंतर्निहित सूत्र को हल करने के लिए प्रायः न्यूटन की विधि जैसी पुनरावृत्तीय विधियों का उपयोग किया जाता है।

कभी-कभी मूल्य की भविष्यवाणी करने के लिए एक स्पष्ट बहुपदीय विधि का उपयोग किया जाता है। फिर उस मान को सही करने के लिए एक अंतर्निहित सूत्र में उपयोग किया जाता है। परिणाम एक भविष्यवक्ता-सुधारक विधि है।

उदाहरण

उदाहरण के लिए समस्या पर विचार करें

सटीक समाधान है।

वन-चरण यूलर

एक सरल संख्यात्मक विधि यूलर की विधि है:

यूलर की विधि को एक चरण के विकृत स्तिथि के लिए एक स्पष्ट बहुपदीय विधि के रूप में देखा जा सकता है।

समस्या पर चरण आकार के साथ लागू की गई यह विधि निम्नलिखित परिणाम देती है:


दो-चरणीय एडम्स-बैशफोर्थ

यूलर की विधि एक चरणीय विधि है। एक सरल बहुचरणीय विधि दो-चरणीय एडम्स-बैशफोर्थ विधि है

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

बहुपदीय विधियों के समूह

रैखिक बहुपदीय विधियों के तीन समूह सामान्यतः उपयोग किए जाते हैं: एडम्स-बैशफोर्थ विधियां, एडम्स-मौल्टन विधियां, और पिछड़े भेदभाव सूत्र (बीडीएफ)।

एडम्स-बैशफोर्थ विधियाँ

एडम्स-बैशफोर्थ विधियाँ स्पष्ट विधियाँ हैं। और गुणांक हैं, जब ऐसे चुना जाता है कि विधियों का क्रम s हो (यह विधियों को विशिष्ट रूप से निर्धारित करता है)।

एडम्स-बैशफोर्थ विधियाँ s = 1, 2, 3, 4, 5 के साथ हैं (हेयरर, नॉरसेट & वानर 1993, §III.1; बुचर 2003, p. 103):