ट्रंकेशन त्रुटि (संख्यात्मक एकीकरण): Difference between revisions
(Created page with "{{Short description|Errors arising in numerical integration}} संख्यात्मक एकीकरण में ट्रंकेशन त्रुटिया...") |
No edit summary |
||
Line 2: | Line 2: | ||
[[संख्यात्मक एकीकरण]] में ट्रंकेशन त्रुटियाँ दो प्रकार की होती हैं: | [[संख्यात्मक एकीकरण]] में ट्रंकेशन त्रुटियाँ दो प्रकार की होती हैं: | ||
* ''स्थानीय | * ''स्थानीय खंडन त्रुटियाँ'' - एक पुनरावृत्ति के कारण होने वाली त्रुटि, और | ||
* '' | * ''ग्लोबल ट्रंकेशन त्रुटियां'' - कई पुनरावृत्तियों के कारण होने वाली संचयी त्रुटि। | ||
== परिभाषाएँ == | == परिभाषाएँ == | ||
Line 9: | Line 9: | ||
: <math> y' = f(t,y), \qquad y(t_0) = y_0, \qquad t \geq t_0 </math> | : <math> y' = f(t,y), \qquad y(t_0) = y_0, \qquad t \geq t_0 </math> | ||
और हम | |||
और हम अलग-अलग समय चरणों <math> t_1,t_2,\ldots,t_N </math> पर वास्तविक समाधान <math> y(t_n) </math> के एक अनुमान <math> y_n </math> की गणना करना चाहते हैं। सरलता के लिए, मान लें कि समय चरण समान दूरी पर हैं: | |||
: <math> h = t_n - t_{n-1}, \qquad n=1,2,\ldots,N. </math> | : <math> h = t_n - t_{n-1}, \qquad n=1,2,\ldots,N. </math> | ||
मान लीजिए हम अनुक्रम | मान लीजिए कि हम प्रपत्र की एक-चरणीय विधि से अनुक्रम <math> y_n </math> की गणना करते हैं | ||
: <math> y_n = y_{n-1} + h A(t_{n-1}, y_{n-1}, h, f). </math> | : <math> y_n = y_{n-1} + h A(t_{n-1}, y_{n-1}, h, f). </math> | ||
फलन <math> A </math> को इंक्रीमेंट फलन कहा जाता है, और इसकी व्याख्या स्लोप <math> \frac{y(t_n)-y(t_{n-1})}{h} </math> के अनुमान के रूप में की जा सकती है। | |||
=== स्थानीय ट्रंकेशन त्रुटि === | === स्थानीय ट्रंकेशन त्रुटि === | ||
स्थानीय | स्थानीय ट्रंकेशन त्रुटि <math> \tau_n</math> यह त्रुटि है कि हमारा वेतन वृद्धि फ़ंक्शन, <math> A </math>, एक ही पुनरावृत्ति के दौरान कारण, पिछले पुनरावृत्ति में सही समाधान का सही ज्ञान मानता है। | ||
अधिक औपचारिक रूप से, | अधिक औपचारिक रूप से, चरण <math> n </math> पर स्थानीय ट्रंकेशन त्रुटि ,<math> \tau_n </math> की गणना वेतन वृद्धि <math> y_n \approx y_{n-1} + h A(t_{n-1}, y_{n-1}, h, f) </math> के लिए समीकरण के बाएँ और दाएँ पक्ष के बीच के अंतर से की जाती है। | ||
:<math> \tau_n = y(t_n) - y(t_{n-1}) - h A(t_{n-1}, y(t_{n-1}), h, f). </math><ref>{{cite journal|last=Gupta|first=G. K.|last2=Sacks-Davis |first2=R. |last3=Tischer |first3=P. E. |title=ODE को हल करने में हाल के विकास की समीक्षा|journal=Computing Surveys|date=March 1985|volume=17|issue=1|pages=5–47|citeseerx = 10.1.1.85.783|doi=10.1145/4078.4079}}</ref><ref>{{harvnb|Süli|Mayers|2003|p=317}}, calls <math> \tau_n/h </math> the truncation error.</ref> | :<math> \tau_n = y(t_n) - y(t_{n-1}) - h A(t_{n-1}, y(t_{n-1}), h, f). </math><ref>{{cite journal|last=Gupta|first=G. K.|last2=Sacks-Davis |first2=R. |last3=Tischer |first3=P. E. |title=ODE को हल करने में हाल के विकास की समीक्षा|journal=Computing Surveys|date=March 1985|volume=17|issue=1|pages=5–47|citeseerx = 10.1.1.85.783|doi=10.1145/4078.4079}}</ref><ref>{{harvnb|Süli|Mayers|2003|p=317}}, calls <math> \tau_n/h </math> the truncation error.</ref> | ||
यदि स्थानीय ट्रंकेशन त्रुटि | यदि स्थानीय ट्रंकेशन त्रुटि <math> o(h) </math> है तो संख्यात्मक विधि सुसंगत है (इसका अर्थ है कि प्रत्येक <math> \varepsilon > 0 </math> के लिए एक <math> H </math> उपस्थित है जैसे कि <math> |\tau_n| < \varepsilon h </math> सभी <math> h < H </math> के लिए; देखें लिटिल-ओ संकेतन)। यदि वृद्धि फ़ंक्शन <math> A </math> निरंतर है, तो विधि सुसंगत है यदि, और केवल यदि <math> A(t,y,0,f) = f(t,y) </math> है।<ref>{{harvnb|Süli|Mayers|2003|pp=321 & 322}}</ref> | ||
इसके अतिरिक्त , हम कहते हैं कि संख्यात्मक विधि में ऑर्डर <math> p </math> है यदि प्रारंभिक मूल्य समस्या के किसी भी पर्याप्त रूप से सुचारू समाधान के लिए, स्थानीय ट्रंकेशन त्रुटि <math> O(h^{p+1}) </math> है (जिसका अर्थ है कि स्थिरांक <math> C </math> और <math> H </math> उपस्थित हैं जैसे वह <math> |\tau_n| < Ch^{p+1} </math> सभी <math> h < H </math> के लिए है ।<ref>{{harvnb|Iserles|1996|p=8}}; {{harvnb|Süli|Mayers|2003|p=323}}</ref> | |||
=== ग्लोबल ट्रंकेशन त्रुटि === | === ग्लोबल ट्रंकेशन त्रुटि === | ||
ग्लोबल ट्रंकेशन त्रुटि सभी पुनरावृत्तियों पर ''स्थानीय ट्रंकेशन त्रुटि'' का संचय है, जो प्रारंभिक समय चरण में सही समाधान का सही ज्ञान मानती है।{{Citation needed|date=April 2013}} | |||
अधिक औपचारिक रूप से, | अधिक औपचारिक रूप से, ग्लोबल ट्रंकेशन त्रुटि,<math> e_n </math> समय <math> t_n </math> पर परिभाषित की गई है: | ||
: <math> | : <math> | ||
Line 38: | Line 39: | ||
\end{align} | \end{align} | ||
</math><ref>{{harvnb|Süli|Mayers|2003|p=317}}</ref> | </math><ref>{{harvnb|Süli|Mayers|2003|p=317}}</ref> | ||
यदि चरण आकार शून्य हो जाता है तो | यदि चरण आकार शून्य हो जाता है तो ग्लोबल ट्रंकेशन त्रुटि शून्य हो जाती है तो संख्यात्मक विधि अभिसरण होती है; दूसरे शब्दों में, संख्यात्मक समाधान स्पष्ट समाधान में परिवर्तित हो जाता है: <math> \lim_{h\to0} \max_n |e_n| = 0 </math>.<ref>{{harvnb|Iserles|1996|p=5}}</ref> | ||
== स्थानीय और | == स्थानीय और ग्लोबल ट्रंकेशन त्रुटियों के बीच संबंध == | ||
कभी-कभी | कभी-कभी ग्लोबल ट्रंकेशन त्रुटि पर ऊपरी सीमा की गणना करना संभव है, अगर हम पहले से ही स्थानीय ट्रंकेशन त्रुटि जानते हैं। इसके लिए आवश्यक है कि हमारा वेतन वृद्धि कार्य पर्याप्त रूप से अच्छा हो। | ||
ग्लोबल ट्रंकेशन त्रुटि पुनरावृत्ति संबंध को संतुष्ट करती है: | |||
:<math> e_{n+1} = e_n + h \Big( A(t_n, y(t_n), h, f) - A(t_n, y_n, h, f) \Big) + \tau_{n+1}. </math> | :<math> e_{n+1} = e_n + h \Big( A(t_n, y(t_n), h, f) - A(t_n, y_n, h, f) \Big) + \tau_{n+1}. </math> | ||
यह परिभाषाओं से तुरंत अनुसरण करता है। अब मान लें कि वृद्धि फ़ंक्शन दूसरे तर्क में | यह परिभाषाओं से तुरंत अनुसरण करता है। अब मान लें कि वेतन वृद्धि फ़ंक्शन दूसरे तर्क में लिप्सचिट्ज़ निरंतर है, अर्थात, एक स्थिरांक <math>L</math> उपस्थित है जैसे कि सभी <math>t</math> और <math>y_1</math> और <math>y_2</math> के लिए, हमारे पास है: | ||
:<math> | A(t,y_1,h,f) - A(t,y_2,h,f) | \le L |y_1-y_2|. </math> | :<math> | A(t,y_1,h,f) - A(t,y_2,h,f) | \le L |y_1-y_2|. </math> | ||
तब | तब ग्लोबल त्रुटि बाध्यता को संतुष्ट करती है | ||
:<math> | e_n | \le \frac{\max_j \tau_j}{hL} \left( \mathrm{e}^{L(t_n-t_0)} - 1 \right). </math><ref>{{harvnb|Süli|Mayers|2003|p=318}}</ref> | :<math> | e_n | \le \frac{\max_j \tau_j}{hL} \left( \mathrm{e}^{L(t_n-t_0)} - 1 \right). </math><ref>{{harvnb|Süli|Mayers|2003|p=318}}</ref> | ||
ग्लोबल त्रुटि के लिए उपरोक्त सीमा से यह पता चलता है कि यदि अंतर समीकरण में फ़ंक्शन <math> f </math> पहले तर्क में निरंतर है और लिप्सचिट्ज़ दूसरे तर्क में निरंतर है (पिकार्ड-लिंडेलोफ प्रमेय से स्थिति), और वृद्धि फ़ंक्शन <math> A </math> निरंतर है सभी तर्कों में और दूसरे तर्क में लिप्सचिट्ज़ निरंतर, तो ग्लोबल त्रुटि शून्य हो जाती है क्योंकि चरण आकार <math> h </math> शून्य के समीप पहुंचता है (दूसरे शब्दों में, संख्यात्मक विधि स्पष्ट समाधान में परिवर्तित हो जाती है)।<ref>{{harvnb|Süli|Mayers|2003|p=322}}</ref> | |||
== [[रैखिक मल्टीस्टेप विधि]]यों का विस्तार == | == [[रैखिक मल्टीस्टेप विधि]]यों का विस्तार == | ||
अब सूत्र द्वारा दी गई एक रैखिक मल्टीस्टेप विधि पर विचार करें | अब सूत्र द्वारा दी गई एक रैखिक मल्टीस्टेप विधि पर विचार करें | ||
Line 61: | Line 61: | ||
इस प्रकार, संख्यात्मक समाधान के लिए अगले मान की गणना इसके अनुसार की जाती है | इस प्रकार, संख्यात्मक समाधान के लिए अगले मान की गणना इसके अनुसार की जाती है | ||
: <math> y_{n+s} = - \sum_{k=0}^{s-1} a_{k} y_{n+k} + h \sum_{k=0}^s b_k f(t_{n+k}, y_{n+k}). </math> | : <math> y_{n+s} = - \sum_{k=0}^{s-1} a_{k} y_{n+k} + h \sum_{k=0}^s b_k f(t_{n+k}, y_{n+k}). </math> | ||
एक रैखिक मल्टीस्टेप विधि का अगला पुनरावृत्त पिछले चरण पर निर्भर करता है। इस प्रकार, स्थानीय ट्रंकेशन त्रुटि की परिभाषा में, अब यह माना जाता है कि पिछले सभी पुनरावृत्त | एक रैखिक मल्टीस्टेप विधि का अगला पुनरावृत्त पिछले चरण पर निर्भर करता है। इस प्रकार, स्थानीय ट्रंकेशन त्रुटि की परिभाषा में, अब यह माना जाता है कि पिछले सभी पुनरावृत्त स्पष्ट समाधान के अनुरूप हैं: | ||
: <math> \tau_n = y(t_{n+s}) + \sum_{k=0}^{s-1} a_{k} y(t_{n+k}) - h \sum_{k=0}^s b_k f(t_{n+k}, y(t_{n+k})). </math><ref>{{harvnb|Süli|Mayers|2003|p=337}}, uses a different definition, dividing this by essentially by ''h''</ref> | : <math> \tau_n = y(t_{n+s}) + \sum_{k=0}^{s-1} a_{k} y(t_{n+k}) - h \sum_{k=0}^s b_k f(t_{n+k}, y(t_{n+k})). </math><ref>{{harvnb|Süli|Mayers|2003|p=337}}, uses a different definition, dividing this by essentially by ''h''</ref> | ||
पुनः, | पुनः, यदि <math> \tau_n = o(h) </math> तो विधि सुसंगत है और यदि <math> \tau_n = O(h^{p+1}) </math> है तो इसका क्रम p है। ग्लोबल ट्रंकेशन त्रुटि की परिभाषा भी अपरिवर्तित है। | ||
स्थानीय और ग्लोबल ट्रंकेशन त्रुटियों के बीच का संबंध एक-चरणीय विधियों की सरल सेटिंग से थोड़ा अलग है। रैखिक मल्टीस्टेप विधियों के लिए, स्थानीय और ग्लोबल ट्रंकेशन त्रुटियों के बीच संबंध को समझाने के लिए शून्य-स्थिरता नामक एक अतिरिक्त अवधारणा की आवश्यकता होती है। शून्य-स्थिरता की स्थिति को पूरा करने वाली रैखिक मल्टीस्टेप विधियाँ स्थानीय और ग्लोबल त्रुटियों के बीच एक-चरणीय विधियों के समान संबंध रखती हैं। दूसरे शब्दों में, यदि एक रैखिक मल्टीस्टेप विधि शून्य-स्थिर और सुसंगत है, तो यह अभिसरण करती है। और यदि एक रैखिक मल्टीस्टेप विधि शून्य-स्थिर है और इसमें स्थानीय त्रुटि <math> \tau_n = O(h^{p+1}) </math> है, तो इसकी ग्लोबल त्रुटि <math> e_n = O(h^p) </math> को संतुष्ट करती है।<ref>{{harvnb|Süli|Mayers|2003|p=340}}</ref> | |||
== यह भी देखें == | == यह भी देखें == | ||
* [[सटीकता का क्रम]] | * [[सटीकता का क्रम|स्पष्टता का क्रम]] | ||
* संख्यात्मक एकीकरण | * संख्यात्मक एकीकरण | ||
*[[संख्यात्मक साधारण अंतर समीकरण]] | *[[संख्यात्मक साधारण अंतर समीकरण]] | ||
* [[काट-छाँट त्रुटि]] | * [[काट-छाँट त्रुटि|ट्रंकेशन त्रुटि]] | ||
== टिप्पणियाँ == | == टिप्पणियाँ == |
Revision as of 09:13, 30 July 2023
संख्यात्मक एकीकरण में ट्रंकेशन त्रुटियाँ दो प्रकार की होती हैं:
- स्थानीय खंडन त्रुटियाँ - एक पुनरावृत्ति के कारण होने वाली त्रुटि, और
- ग्लोबल ट्रंकेशन त्रुटियां - कई पुनरावृत्तियों के कारण होने वाली संचयी त्रुटि।
परिभाषाएँ
मान लीजिए हमारे पास एक सतत अवकल समीकरण है
और हम अलग-अलग समय चरणों पर वास्तविक समाधान के एक अनुमान की गणना करना चाहते हैं। सरलता के लिए, मान लें कि समय चरण समान दूरी पर हैं:
मान लीजिए कि हम प्रपत्र की एक-चरणीय विधि से अनुक्रम की गणना करते हैं
फलन को इंक्रीमेंट फलन कहा जाता है, और इसकी व्याख्या स्लोप के अनुमान के रूप में की जा सकती है।
स्थानीय ट्रंकेशन त्रुटि
स्थानीय ट्रंकेशन त्रुटि यह त्रुटि है कि हमारा वेतन वृद्धि फ़ंक्शन, , एक ही पुनरावृत्ति के दौरान कारण, पिछले पुनरावृत्ति में सही समाधान का सही ज्ञान मानता है।
अधिक औपचारिक रूप से, चरण पर स्थानीय ट्रंकेशन त्रुटि , की गणना वेतन वृद्धि के लिए समीकरण के बाएँ और दाएँ पक्ष के बीच के अंतर से की जाती है।
यदि स्थानीय ट्रंकेशन त्रुटि है तो संख्यात्मक विधि सुसंगत है (इसका अर्थ है कि प्रत्येक के लिए एक उपस्थित है जैसे कि सभी के लिए; देखें लिटिल-ओ संकेतन)। यदि वृद्धि फ़ंक्शन निरंतर है, तो विधि सुसंगत है यदि, और केवल यदि है।[3]
इसके अतिरिक्त , हम कहते हैं कि संख्यात्मक विधि में ऑर्डर है यदि प्रारंभिक मूल्य समस्या के किसी भी पर्याप्त रूप से सुचारू समाधान के लिए, स्थानीय ट्रंकेशन त्रुटि है (जिसका अर्थ है कि स्थिरांक और उपस्थित हैं जैसे वह सभी के लिए है ।[4]
ग्लोबल ट्रंकेशन त्रुटि
ग्लोबल ट्रंकेशन त्रुटि सभी पुनरावृत्तियों पर स्थानीय ट्रंकेशन त्रुटि का संचय है, जो प्रारंभिक समय चरण में सही समाधान का सही ज्ञान मानती है।[citation needed]
अधिक औपचारिक रूप से, ग्लोबल ट्रंकेशन त्रुटि, समय पर परिभाषित की गई है:
यदि चरण आकार शून्य हो जाता है तो ग्लोबल ट्रंकेशन त्रुटि शून्य हो जाती है तो संख्यात्मक विधि अभिसरण होती है; दूसरे शब्दों में, संख्यात्मक समाधान स्पष्ट समाधान में परिवर्तित हो जाता है: .[6]
स्थानीय और ग्लोबल ट्रंकेशन त्रुटियों के बीच संबंध
कभी-कभी ग्लोबल ट्रंकेशन त्रुटि पर ऊपरी सीमा की गणना करना संभव है, अगर हम पहले से ही स्थानीय ट्रंकेशन त्रुटि जानते हैं। इसके लिए आवश्यक है कि हमारा वेतन वृद्धि कार्य पर्याप्त रूप से अच्छा हो।
ग्लोबल ट्रंकेशन त्रुटि पुनरावृत्ति संबंध को संतुष्ट करती है:
यह परिभाषाओं से तुरंत अनुसरण करता है। अब मान लें कि वेतन वृद्धि फ़ंक्शन दूसरे तर्क में लिप्सचिट्ज़ निरंतर है, अर्थात, एक स्थिरांक उपस्थित है जैसे कि सभी और और के लिए, हमारे पास है:
तब ग्लोबल त्रुटि बाध्यता को संतुष्ट करती है
ग्लोबल त्रुटि के लिए उपरोक्त सीमा से यह पता चलता है कि यदि अंतर समीकरण में फ़ंक्शन पहले तर्क में निरंतर है और लिप्सचिट्ज़ दूसरे तर्क में निरंतर है (पिकार्ड-लिंडेलोफ प्रमेय से स्थिति), और वृद्धि फ़ंक्शन निरंतर है सभी तर्कों में और दूसरे तर्क में लिप्सचिट्ज़ निरंतर, तो ग्लोबल त्रुटि शून्य हो जाती है क्योंकि चरण आकार शून्य के समीप पहुंचता है (दूसरे शब्दों में, संख्यात्मक विधि स्पष्ट समाधान में परिवर्तित हो जाती है)।[8]
रैखिक मल्टीस्टेप विधियों का विस्तार
अब सूत्र द्वारा दी गई एक रैखिक मल्टीस्टेप विधि पर विचार करें
इस प्रकार, संख्यात्मक समाधान के लिए अगले मान की गणना इसके अनुसार की जाती है
एक रैखिक मल्टीस्टेप विधि का अगला पुनरावृत्त पिछले चरण पर निर्भर करता है। इस प्रकार, स्थानीय ट्रंकेशन त्रुटि की परिभाषा में, अब यह माना जाता है कि पिछले सभी पुनरावृत्त स्पष्ट समाधान के अनुरूप हैं:
पुनः, यदि तो विधि सुसंगत है और यदि है तो इसका क्रम p है। ग्लोबल ट्रंकेशन त्रुटि की परिभाषा भी अपरिवर्तित है।
स्थानीय और ग्लोबल ट्रंकेशन त्रुटियों के बीच का संबंध एक-चरणीय विधियों की सरल सेटिंग से थोड़ा अलग है। रैखिक मल्टीस्टेप विधियों के लिए, स्थानीय और ग्लोबल ट्रंकेशन त्रुटियों के बीच संबंध को समझाने के लिए शून्य-स्थिरता नामक एक अतिरिक्त अवधारणा की आवश्यकता होती है। शून्य-स्थिरता की स्थिति को पूरा करने वाली रैखिक मल्टीस्टेप विधियाँ स्थानीय और ग्लोबल त्रुटियों के बीच एक-चरणीय विधियों के समान संबंध रखती हैं। दूसरे शब्दों में, यदि एक रैखिक मल्टीस्टेप विधि शून्य-स्थिर और सुसंगत है, तो यह अभिसरण करती है। और यदि एक रैखिक मल्टीस्टेप विधि शून्य-स्थिर है और इसमें स्थानीय त्रुटि है, तो इसकी ग्लोबल त्रुटि को संतुष्ट करती है।[10]
यह भी देखें
- स्पष्टता का क्रम
- संख्यात्मक एकीकरण
- संख्यात्मक साधारण अंतर समीकरण
- ट्रंकेशन त्रुटि
टिप्पणियाँ
- ↑ Gupta, G. K.; Sacks-Davis, R.; Tischer, P. E. (March 1985). "ODE को हल करने में हाल के विकास की समीक्षा". Computing Surveys. 17 (1): 5–47. CiteSeerX 10.1.1.85.783. doi:10.1145/4078.4079.
- ↑ Süli & Mayers 2003, p. 317, calls the truncation error.
- ↑ Süli & Mayers 2003, pp. 321 & 322
- ↑ Iserles 1996, p. 8; Süli & Mayers 2003, p. 323
- ↑ Süli & Mayers 2003, p. 317
- ↑ Iserles 1996, p. 5
- ↑ Süli & Mayers 2003, p. 318
- ↑ Süli & Mayers 2003, p. 322
- ↑ Süli & Mayers 2003, p. 337, uses a different definition, dividing this by essentially by h
- ↑ Süli & Mayers 2003, p. 340
संदर्भ
- Iserles, Arieh (1996), A First Course in the Numerical Analysis of Differential Equations, Cambridge University Press, ISBN 978-0-521-55655-2.
- Süli, Endre; Mayers, David (2003), An Introduction to Numerical Analysis, Cambridge University Press, ISBN 0521007941.