यूलर पद्धति: Difference between revisions
(Created page with "{{Short description|Approach to finding numerical solutions of ordinary differential equations}} {{for multi|integrating with respect to the Euler characteristic|Euler calculu...") |
|||
| Line 214: | Line 214: | ||
* | |||
* | |||
== बाहरी संबंध == | == बाहरी संबंध == | ||
{{wikibooks|1=Calculus|2=Euler's Method}} | {{wikibooks|1=Calculus|2=Euler's Method}} | ||
Revision as of 11:41, 28 November 2022
| अंतर समीकरण |
|---|
| Navier–Stokes differential equations used to simulate airflow around an obstruction |
| दायरा |
| वर्गीकरण |
| समाधान |
| लोग |
गणित और कम्प्यूटेशनल विज्ञान में, यूलर विधि (जिसे फॉरवर्ड यूलर विधि भी कहा जाता है) एक प्रारंभिक मूल्य समस्या के साथ सामान्य अंतर समीकरणों (ओडीई) को हल करने के लिए प्रथम क्रम संख्यात्मक विश्लेषण प्रक्रिया है। यह संख्यात्मक साधारण अंतर समीकरणों के लिए सबसे बुनियादी स्पष्ट और निहित तरीके हैं और सबसे सरल रनगे-कुट्टा विधि है। यूलर विधि का नाम लियोनहार्ड यूलर के नाम पर रखा गया है, जिन्होंने अपनी पुस्तक इंटीग्रल कैलकुस संस्थान (प्रकाशित 1768-1870) में इसका इलाज किया था।[1]
यूलर विधि एक प्रथम-क्रम विधि है, जिसका अर्थ है कि स्थानीय त्रुटि (प्रति चरण त्रुटि) चरण आकार के वर्ग के समानुपाती होती है, और वैश्विक त्रुटि (किसी निश्चित समय पर त्रुटि) चरण आकार के समानुपाती होती है। यूलर विधि अक्सर अधिक जटिल विधियों के निर्माण के आधार के रूप में कार्य करती है, उदाहरण के लिए, भविष्यवक्ता-सुधारक विधि।
अनौपचारिक ज्यामितीय विवरण
एक अज्ञात वक्र के आकार की गणना करने की समस्या पर विचार करें जो एक दिए गए बिंदु से शुरू होता है और एक दिए गए अवकल समीकरण को संतुष्ट करता है। यहां, एक अंतर समीकरण को एक सूत्र के रूप में माना जा सकता है जिसके द्वारा उस बिंदु की स्थिति की गणना करने के बाद वक्र पर स्पर्शरेखा रेखा की ढलान की गणना वक्र पर किसी भी बिंदु पर की जा सकती है।
विचार यह है कि जबकि वक्र प्रारंभ में अज्ञात है, इसका प्रारंभिक बिंदु, जिसे हम निरूपित करते हैं ज्ञात है (शीर्ष दाईं ओर चित्र देखें)। फिर, अंतर समीकरण से, ढलान से वक्र तक गणना की जा सकती है, और इसलिए, स्पर्श रेखा।
उस स्पर्श रेखा के साथ एक बिंदु तक एक छोटा कदम उठाएं इस छोटे कदम के साथ, ढलान बहुत ज्यादा नहीं बदलता है, इसलिए वक्र के करीब होगा। अगर हम ऐसा दिखावा करते हैं अभी भी वक्र पर है, वही तर्क जो बिंदु के लिए है ऊपर इस्तेमाल किया जा सकता है। कई चरणों के बाद, एक बहुभुज वक्र गणना की जाती है। सामान्य तौर पर, यह वक्र मूल अज्ञात वक्र से बहुत दूर नहीं जाता है, और दो वक्रों के बीच की त्रुटि को छोटा किया जा सकता है यदि चरण का आकार काफी छोटा है और गणना का अंतराल परिमित है:[2]
कोई मान चुनें प्रत्येक चरण और सेट के आकार के लिए . अब, यूलर विधि का एक चरण प्रति है:[3]
का मूल्य समय पर ODE के समाधान का एक अनुमान है : . यूलर विधि स्पष्ट और अंतर्निहित विधियाँ हैं, अर्थात समाधान का एक स्पष्ट कार्य है के लिये .
जबकि यूलर विधि पहले क्रम के ODE, आदेश के किसी भी ODE को एकीकृत करती है प्रथम क्रम ODEs की एक प्रणाली के रूप में प्रतिनिधित्व किया जा सकता है: समीकरण का इलाज करने के लिए
हम सहायक चर पेश करते हैं और प्राप्त करें समतुल्य समीकरण:
यह चर में प्रथम-क्रम प्रणाली है और यूलर की विधि या, वास्तव में, प्रथम-क्रम प्रणालियों के लिए किसी अन्य योजना द्वारा नियंत्रित किया जा सकता है।[4]
उदाहरण
प्रारंभिक मूल्य समस्या को देखते हुए
हम अनुमानित करने के लिए यूलर विधि का उपयोग करना चाहेंगे .[5]
1 के बराबर चरण आकार का उपयोग करना (h = 1)
छवि: संख्यात्मक एकीकरण चित्रण, चरण =1.svg|right|thumb|समीकरण के लिए संख्यात्मक एकीकरण का चित्रण नीला यूलर विधि है; हरा, मध्यबिंदु विधि; लाल, सटीक समाधान, स्टेप साइज है यूलर विधि है
इसलिए पहले हमें गणना करनी चाहिए . इस सरल अंतर समीकरण में, फ़ंक्शन द्वारा परिभाषित किया गया है . हमारे पास है
उपरोक्त चरण को करके, हमने उस रेखा का ढलान पाया है जो बिंदु पर समाधान वक्र की स्पर्शरेखा है . याद रखें कि ढलान को परिवर्तन के रूप में परिभाषित किया गया है में परिवर्तन से विभाजित , या .
अगला चरण उपरोक्त मान को चरण आकार से गुणा करना है , जिसे हम यहाँ एक के बराबर लेते हैं:
चूंकि स्टेप साइज में बदलाव है , जब हम चरण आकार और स्पर्शरेखा के ढलान को गुणा करते हैं, तो हमें एक परिवर्तन मिलता है मूल्य। यह मान तब प्रारंभिक में जोड़ा जाता है गणना के लिए उपयोग किए जाने वाले अगले मान को प्राप्त करने के लिए मूल्य।
उपरोक्त चरणों को खोजने के लिए दोहराया जाना चाहिए , तथा .
इस एल्गोरिथ्म की दोहरावदार प्रकृति के कारण, त्रुटियों से बचने के लिए, जैसा कि नीचे देखा गया है, गणनाओं को चार्ट के रूप में व्यवस्थित करना सहायक हो सकता है।
0 1 0 1 1 1 2 1 2 1 2 1 2 4 2 4 2 4 1 4 8 3 8 3 8 1 8 16
इस गणना का निष्कर्ष यह है कि . अवकल समीकरण का सटीक हल है , इसलिए . हालांकि इस विशिष्ट मामले में यूलर पद्धति का सन्निकटन बहुत सटीक नहीं था, विशेष रूप से एक बड़े मूल्य चरण आकार के कारण , इसका व्यवहार गुणात्मक रूप से सही है जैसा कि चित्र में दिखाया गया है।
अन्य चरण आकारों का उपयोग
छवि: संख्यात्मक एकीकरण उदाहरण चरण =0.25.svg|right|thumb|के लिए एक ही दृष्टांत जैसा कि प्रस्तावना में सुझाया गया है, कदम आकार के मामले में यूलर विधि अधिक सटीक है छोटा है। नीचे दी गई तालिका विभिन्न चरणों के आकार के साथ परिणाम दिखाती है। शीर्ष पंक्ति पिछले अनुभाग में उदाहरण से मेल खाती है, और दूसरी पंक्ति चित्र में सचित्र है।
step size result of Euler's method error 1 16.00 38.60 0.25 35.53 19.07 0.1 45.26 9.34 0.05 49.56 5.04 0.025 51.98 2.62 0.0125 53.26 1.34
तालिका के अंतिम कॉलम में दर्ज की गई त्रुटि सटीक समाधान के बीच का अंतर है और यूलर सन्निकटन। तालिका के निचले भाग में, पिछली पंक्ति में चरण का आकार आधा है, और त्रुटि भी पिछली पंक्ति में त्रुटि का लगभग आधा है। इससे पता चलता है कि त्रुटि चरण आकार के लगभग आनुपातिक है, कम से कम चरण आकार के काफी छोटे मूल्यों के लिए। यह व्यापक रूप से अन्य समीकरणों के लिए भी सत्य है; अधिक विवरण के लिए खंड #वैश्विक ट्रंकेशन त्रुटि देखें।
अन्य विधियाँ, जैसे कि मिडपॉइंट विधि भी आंकड़ों में सचित्र है, अधिक अनुकूल व्यवहार करती है: मिडपॉइंट विधि की वैश्विक त्रुटि लगभग चरण आकार के वर्ग के समानुपाती होती है। इस कारण से, यूलर विधि को प्रथम कोटि की विधि कहा जाता है, जबकि मध्य बिंदु विधि को दूसरी कोटि की विधि कहा जाता है।
हम उपरोक्त तालिका से एक्सट्रपलेशन कर सकते हैं कि तीन दशमलव स्थानों तक सही उत्तर प्राप्त करने के लिए आवश्यक चरण आकार लगभग 0.00001 है, जिसका अर्थ है कि हमें 400,000 चरणों की आवश्यकता है। इतनी बड़ी संख्या में चरणों में उच्च कम्प्यूटेशनल लागत होती है। इस कारण से, उच्च-क्रम विधियों को नियोजित किया जाता है जैसे रनगे-कुट्टा विधियों या रैखिक मल्टीस्टेप विधियों, विशेष रूप से यदि उच्च सटीकता वांछित है।[6]
व्युत्पत्ति
यूलर विधि को कई तरीकों से प्राप्त किया जा सकता है। सबसे पहले, ऊपर ज्यामितीय विवरण है।
फ़ंक्शन के टेलर विस्तार पर विचार करने की एक और संभावना है चारों ओर :
अंतर समीकरण बताता है कि . यदि इसे टेलर विस्तार में प्रतिस्थापित किया जाता है और द्विघात और उच्च-क्रम की शर्तों को अनदेखा किया जाता है, तो यूलर विधि उत्पन्न होती है।[7] यूलर विधि द्वारा की गई त्रुटि का विश्लेषण करने के लिए टेलर विस्तार का उपयोग नीचे किया गया है, और इसे रनगे-कुट्टा विधियों का उत्पादन करने के लिए बढ़ाया जा सकता है।
व्युत्पन्न के लिए आगे परिमित अंतर सूत्र को प्रतिस्थापित करने के लिए एक निकटता से संबंधित व्युत्पत्ति है,
अंतर समीकरण में . दोबारा, यह यूलर विधि उत्पन्न करता है।[8] इसी तरह की गणना मिडपॉइंट विधि और बैकवर्ड यूलर विधि की ओर ले जाती है।
अंत में, कोई अंतर समीकरण को एकीकृत कर सकता है प्रति और कलन की मौलिक प्रमेय को प्राप्त करने के लिए लागू करें:
अब बाएँ हाथ की आयत विधि (केवल एक आयत के साथ) द्वारा अभिन्न का अनुमान लगाएं:
दोनों समीकरणों को मिलाकर, फिर से यूलर विधि प्राप्त होती है।[9] विभिन्न रेखीय मल्टीस्टेप विधियों पर पहुंचने के लिए विचार की इस पंक्ति को जारी रखा जा सकता है।
स्थानीय ट्रंकेशन त्रुटि
यूलर विधि की स्थानीय ट्रंकेशन त्रुटि एक चरण में की गई त्रुटि है। यह एक कदम के बाद संख्यात्मक समाधान के बीच का अंतर है, , और समय पर सटीक समाधान . द्वारा संख्यात्मक समाधान दिया गया है
सटीक समाधान के लिए, हम ऊपर #Derivation सेक्शन में उल्लिखित टेलर विस्तार का उपयोग करते हैं:
यूलर विधि द्वारा शुरू की गई स्थानीय ट्रंकेशन त्रुटि (एलटीई) इन समीकरणों के बीच के अंतर से दी गई है:
यह परिणाम मान्य है यदि एक सीमित तीसरा व्युत्पन्न है।[10] इससे पता चलता है कि छोटे के लिए , स्थानीय ट्रंकेशन त्रुटि लगभग आनुपातिक है . यह यूलर विधि को कम सटीक बनाता है (छोटे के लिए ) रनगे-कुट्टा विधियों और रैखिक मल्टीस्टेप विधियों जैसी अन्य उच्च-क्रम तकनीकों की तुलना में, जिसके लिए स्थानीय ट्रंकेशन त्रुटि चरण आकार की उच्च शक्ति के समानुपाती होती है।
टेलर के प्रमेय में शेष अवधि के लिए लैग्रेंज फॉर्म का उपयोग करके स्थानीय ट्रंकेशन त्रुटि के लिए थोड़ा अलग फॉर्मूलेशन प्राप्त किया जा सकता है। यदि लगातार दूसरा व्युत्पन्न है, तो वहां मौजूद है ऐसा है कि
त्रुटि के लिए उपरोक्त भावों में, अज्ञात सटीक समाधान का दूसरा व्युत्पन्न अंतर समीकरण के दाईं ओर शामिल एक अभिव्यक्ति द्वारा प्रतिस्थापित किया जा सकता है। दरअसल, यह समीकरण से चलता है वह[12]
वैश्विक ट्रंकेशन त्रुटि
वैश्विक ट्रंकेशन त्रुटि एक निश्चित समय पर त्रुटि है आरंभिक समय से उस समय तक पहुंचने के लिए विधि को जितने भी कदम उठाने होंगे, उसके बाद। वैश्विक ट्रंकेशन त्रुटि प्रत्येक चरण में की गई स्थानीय ट्रंकेशन त्रुटियों का संचयी प्रभाव है।[13] चरणों की संख्या आसानी से निर्धारित की जाती है , जो आनुपातिक है , और प्रत्येक चरण में की गई त्रुटि आनुपातिक है (पिछला भाग देखें)। इस प्रकार, यह उम्मीद की जानी चाहिए कि वैश्विक ट्रंकेशन त्रुटि आनुपातिक होगी .[14] इस सहज तर्क को सटीक बनाया जा सकता है। यदि समाधान एक बंधा हुआ दूसरा व्युत्पन्न है और लिप्सचिट्ज़ निरंतरता अपने दूसरे तर्क में है, तो वैश्विक ट्रंकेशन त्रुटि (जीटीई) से घिरा है
कहाँ पे के दूसरे व्युत्पन्न पर एक ऊपरी सीमा है दिए गए अंतराल पर और का लिपशिट्ज स्थिरांक है .[15] इस बाउंड का सटीक रूप थोड़ा व्यावहारिक महत्व का है, क्योंकि ज्यादातर मामलों में बाउंड यूलर विधि द्वारा की गई वास्तविक त्रुटि को बहुत अधिक बढ़ा देता है।[16] जो महत्वपूर्ण है वह यह दर्शाता है कि वैश्विक ट्रंकेशन त्रुटि (लगभग) आनुपातिक है . इस कारण से, यूलर विधि को प्रथम कोटि की विधि कहा जाता है।[17]
संख्यात्मक स्थिरता
यूलर विधि संख्यात्मक रूप से संख्यात्मक स्थिरता भी हो सकती है, विशेष रूप से कठोर समीकरणों के लिए, जिसका अर्थ है कि संख्यात्मक समाधान उन समीकरणों के लिए बहुत बड़ा हो जाता है जहां सटीक समाधान नहीं होता है। इसे रैखिक समीकरण का उपयोग करके चित्रित किया जा सकता है
अचूक उपाय है , जो शून्य के रूप में घटता है . हालाँकि, यदि इस समीकरण पर चरण आकार के साथ यूलर विधि लागू की जाती है , तो संख्यात्मक समाधान गुणात्मक रूप से गलत है: यह दोलन करता है और बढ़ता है (चित्र देखें)। अस्थिर होने का यही अर्थ है। उदाहरण के लिए, यदि छोटे चरण आकार का उपयोग किया जाता है , तो संख्यात्मक समाधान क्षय से शून्य हो जाता है।
यदि यूलर विधि को रैखिक समीकरण पर लागू किया जाता है , तो संख्यात्मक समाधान अस्थिर है अगर उत्पाद क्षेत्र के बाहर है
दाईं ओर चित्रित। इस क्षेत्र को (रैखिक) स्थिरता क्षेत्र कहा जाता है।[18] उदाहरण में, , तो अगर फिर जो स्थिरता क्षेत्र के बाहर है, और इस प्रकार संख्यात्मक समाधान अस्थिर है।
यह सीमा - त्रुटि के धीमे अभिसरण के साथ — का मतलब है कि संख्यात्मक एकीकरण के एक साधारण उदाहरण को छोड़कर, यूलर विधि का अक्सर उपयोग नहीं किया जाता है[citation needed].
राउंडिंग एरर
चरण में यूलर विधि में, गोलाई त्रुटि मोटे तौर पर परिमाण की है कहाँ पे मशीन एप्सिलॉन है। यह मानते हुए कि राउंडिंग त्रुटियां स्वतंत्र यादृच्छिक चर हैं, अपेक्षित कुल राउंडिंग त्रुटि आनुपातिक है .[19] इस प्रकार, चरण आकार के अत्यंत छोटे मानों के लिए ट्रंकेशन त्रुटि छोटी होगी लेकिन राउंडिंग त्रुटि का प्रभाव बड़ा हो सकता है। राउंडिंग एरर के अधिकांश प्रभाव को आसानी से टाला जा सकता है यदि यूलर विधि के सूत्र में मुआवजा योग का उपयोग किया जाता है।[20]
संशोधन और एक्सटेंशन
यूलर विधि का एक सरल संशोधन जो नोट की गई स्थिरता की समस्याओं को समाप्त करता है # संख्यात्मक स्थिरता पश्च यूलर विधि है:
यह उस कार्य में (मानक, या आगे) यूलर विधि से भिन्न है प्रारंभिक बिंदु के बजाय चरण के अंत बिंदु पर मूल्यांकन किया जाता है। बैकवर्ड यूलर विधि एक स्पष्ट और निहित विधि है, जिसका अर्थ है कि बैकवर्ड यूलर विधि का सूत्र है दोनों तरफ, इसलिए पश्चगामी यूलर विधि को लागू करते समय हमें एक समीकरण को हल करना होगा। यह कार्यान्वयन को और अधिक महंगा बनाता है।
यूलर विधि के अन्य संशोधन जो स्थिरता में मदद करते हैंघातीय यूलर विधि विधि या अर्ध-अंतर्निहित यूलर विधि उत्पन्न करते हैं।
अधिक जटिल विधियां उच्च क्रम (और अधिक सटीकता) प्राप्त कर सकती हैं। अधिक फ़ंक्शन मूल्यांकनों का उपयोग करने की एक संभावना है। यह मिडपॉइंट विधि द्वारा सचित्र है जिसका उल्लेख इस लेख में पहले ही किया जा चुका है:
- .
यह रनगे-कुट्टा विधियों के परिवार की ओर जाता है।
दूसरी संभावना अधिक पिछले मूल्यों का उपयोग करने की है, जैसा कि दो-चरणीय एडम्स-बैशफोर्थ विधि द्वारा दिखाया गया है:
यह रैखिक मल्टीस्टेप विधियों के परिवार की ओर जाता है। ऐसे अन्य संशोधन हैं जो मेमोरी उपयोग को कम करने के लिए कंप्रेसिव सेंसिंग की तकनीकों का उपयोग करते हैं[21]
लोकप्रिय संस्कृति में
फिल्म छिपे हुए आंकड़े में, कैथरीन गोबल ने पृथ्वी की कक्षा से अंतरिक्ष यात्री जॉन ग्लेन के पुन: प्रवेश की गणना करने के लिए यूलर विधि का सहारा लिया।[22]
यह भी देखें
- क्रैंक-निकोलसन विधि
- ढतला हुआ वंश इसी तरह परिमित चरणों का उपयोग करता है, यहाँ कार्यों की न्यूनतमता खोजने के लिए
- रनगे-कुट्टा विधियों की सूची
- लीनियर मल्टीस्टेप विधि
- संख्यात्मक एकीकरण (निश्चित समाकलों की गणना के लिए)
- साधारण अंतर समीकरणों के लिए संख्यात्मक तरीके
टिप्पणियाँ
- ↑ Butcher 2003, p. 45; Hairer, Nørsett & Wanner 1993, p. 35
- ↑ Atkinson 1989, p. 342; Butcher 2003, p. 60
- ↑ Butcher 2003, p. 45; Hairer, Nørsett & Wanner 1993, p. 36
- ↑ Butcher 2003, p. 3; Hairer, Nørsett & Wanner 1993, p. 2
- ↑ See also Atkinson 1989, p. 344
- ↑ Hairer, Nørsett & Wanner 1993, p. 40
- ↑ Atkinson 1989, p. 342; Hairer, Nørsett & Wanner 1993, p. 36
- ↑ Atkinson 1989, p. 342
- ↑ Atkinson 1989, p. 343
- ↑ Butcher 2003, p. 60
- ↑ Atkinson 1989, p. 342
- ↑ Stoer & Bulirsch 2002, p. 474
- ↑ Atkinson 1989, p. 344
- ↑ Butcher 2003, p. 49
- ↑ Atkinson 1989, p. 346; Lakoba 2012, equation (1.16)
- ↑ Iserles 1996, p. 7
- ↑ Butcher 2003, p. 63
- ↑ Butcher 2003, p. 70; Iserles 1996, p. 57
- ↑ Butcher 2003, pp. 74–75
- ↑ Butcher 2003, pp. 75–78
- ↑ Unni, M. P.; Chandra, M. G.; Kumar, A. A. (March 2017). "कंप्रेसिव सेंसिंग का उपयोग करके डिफरेंशियल इक्वेशन के न्यूमेरिकल सॉल्यूशन के लिए मेमोरी रिडक्शन". 2017 IEEE 13th International Colloquium on Signal Processing Its Applications (CSPA): 79–84. doi:10.1109/CSPA.2017.8064928. ISBN 978-1-5090-1184-1. S2CID 13082456.
- ↑ Khan, Amina (9 January 2017). "मिलिए 'हिडन फिगर्स' गणितज्ञ से जिन्होंने अमेरिकियों को अंतरिक्ष में भेजने में मदद की". Los Angeles Times. Retrieved 12 February 2017.
संदर्भ
- Atkinson, Kendall A. (1989). An Introduction to Numerical Analysis (2nd ed.). New York: John Wiley & Sons. ISBN 978-0-471-50023-0.
- Ascher, Uri M.; Petzold, Linda R. (1998). Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations. Philadelphia: Society for Industrial and Applied Mathematics. ISBN 978-0-89871-412-8.
- Butcher, John C. (2003). Numerical Methods for Ordinary Differential Equations. New York: John Wiley & Sons. ISBN 978-0-471-96758-3.
- Hairer, Ernst; Nørsett, Syvert Paul; Wanner, Gerhard (1993). Solving ordinary differential equations I: Nonstiff problems. Berlin, New York: Springer-Verlag. ISBN 978-3-540-56670-0.
- Iserles, Arieh (1996). A First Course in the Numerical Analysis of Differential Equations. Cambridge University Press. ISBN 978-0-521-55655-2.
- Stoer, Josef; Bulirsch, Roland (2002). Introduction to Numerical Analysis (3rd ed.). Berlin, New York: Springer-Verlag. ISBN 978-0-387-95452-3.
- Lakoba, Taras I. (2012), Simple Euler method and its modifications (PDF) (Lecture notes for MATH334), University of Vermont, retrieved 29 February 2012
- Unni, M P. (2017). "Memory reduction for numerical solution of differential equations using compressive sensing". 2017 IEEE 13th International Colloquium on Signal Processing & its Applications (CSPA). IEEE CSPA. pp. 79–84. doi:10.1109/CSPA.2017.8064928. ISBN 978-1-5090-1184-1. S2CID 13082456.
बाहरी संबंध
- File:Commons-logo.svg Media related to यूलर पद्धति at Wikimedia Commons
- Euler method implementations in different languages by Rosetta Code
- "Euler method", Encyclopedia of Mathematics, EMS Press, 2001 [1994]