स्पर्श रेखा सन्निकटन
ज्यामितीय विचार
एक साधारण अंतर समीकरण dy/dx = f(x,y) (x,y) समतल में प्रत्येक बिंदु को एक ढलान निर्दिष्ट करता है — एक दिशा क्षेत्र। सच्चा समाधान y(x) एक वक्र है जो हर जगह उन निर्दिष्ट ढलानों का अनुसरण करता है।
यूलर की विधि सतत दिशा क्षेत्र को एक असतत चलने में परिवर्तित करती है:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
बिंदु (xₙ, yₙ) से, स्पर्श रेखा के साथ दूरी h की ओर जाएं। एक अनुमानित अगला बिंदु पर पहुंचें। दोहराएं।
ज्यामितीय त्रुटि: (xₙ, yₙ) पर स्पर्श रेखा की ढलान f(xₙ, yₙ) है, लेकिन सच्चे वक्र की अंतराल [xₙ, xₙ + h] के साथ हर बिंदु पर एक अलग ढलान है। यूलर चरण पूरे समय बाएं अंतबिंदु पर ढलान का उपयोग करता है — 'जो ढलान थी।' प्रति चरण त्रुटि h² के रूप में बढ़ती है।
संचित त्रुटि
एक निश्चित अंतबिंदु x = a तक पहुंचने के लिए N चरणों पर, h = a/N के साथ:
- स्थानीय ट्रंकेशन त्रुटि प्रति चरण: O(h²)
- चरणों की संख्या: N = a/h
- वैश्विक त्रुटि: O(h²) × (a/h) = O(h) — प्रथम-क्रम सटीकता
यूलर की विधि प्रथम-क्रम है: h को आधा करने से वैश्विक त्रुटि आधी हो जाती है।
यूलर की विधि चलाना
dy/dx = y पर विचार करें, प्रारंभिक स्थिति y(0) = 1 के साथ। सच्चा समाधान: y(x) = eˣ, तो y(1) = e ≈ 2.71828।
h = 0.5 के साथ यूलर की विधि लागू करें, x = 0 से x = 1 तक (2 चरण):
चरण 1: y₁ = y₀ + h·f(x₀, y₀) = 1 + 0.5·(1) = 1.5। नया बिंदु: (0.5, 1.5)।
चरण 2: y₂ = y₁ + h·f(x₁, y₁) = 1.5 + 0.5·(1.5) = 2.25। नया बिंदु: (1.0, 2.25)।
यूलर 2.25 देता है बनाम सच्चा मान 2.71828। त्रुटि: 0.468। सापेक्ष त्रुटि: ~17%।
यूलर की स्थिरता क्षेत्र व्युत्पन्न करना
परीक्षण समीकरण dy/dx = λy के लिए (जहां λ एक जटिल संख्या है), यूलर की विधि देती है:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
प्रति चरण प्रवर्धन कारक: z = 1 + hλ।
स्थिरता स्थिति: गणना किया गया समाधान सीमित रहता है यदि और केवल यदि |z| ≤ 1, अर्थात्, |1 + hλ| ≤ 1।
यह जटिल hλ-समतल में एक ज्यामितीय स्थिति है: बिंदु hλ (-1, 0) पर केंद्रित त्रिज्या 1 के वृत्त के अंदर होना चाहिए।
यूलर की स्थिरता क्षेत्र: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
वास्तविक, नकारात्मक λ के लिए (एक क्षय ODE जैसे dy/dx = -2y): hλ वास्तविक अक्ष पर अंतराल (-2, 0] में होना चाहिए। λ = -2 और h = 0.5 के साथ: hλ = -1। यह ठीक स्थिरता सीमा पर है — विधि सीमांत रूप से स्थिर है, जो पहले के उदाहरण में गुणात्मक विफलता की व्याख्या करता है।
h = 1 और λ = -2 के साथ: hλ = -2, हमें स्थिरता क्षेत्र के बाहर रखता है। समाधान बढ़ते हुए आयाम के साथ दोलन करता है।
स्थिरता सीमा ढूंढना
Runge-Kutta 4 (RK4) के पास यूलर की तुलना में एक बड़ा स्थिरता क्षेत्र है, जो इसे अधिकांश समस्याओं के लिए पसंद किए जाने का एक कारण है।
वास्तविक नकारात्मक λ के लिए, RK4 hλ को वास्तविक अक्ष पर लगभग -2.785 तक की अनुमति देता है (यूलर की -2 सीमा बनाम)।
बहुत अलग परिमाण पर eigenvalues के साथ कठोर समीकरणों के लिए — कहो λ₁ = -1 और λ₂ = -1000 — स्थिरता के लिए hλ₂ क्षेत्र के अंदर रहना आवश्यक है। वास्तविक अक्ष पर RK4 के लिए: h·(-1000) ≥ -2.785, तो h ≤ 0.002785।
यह छोटा चरण आकार, कठोर eigenvalue λ₂ द्वारा निर्दिष्ट, सिमुलेशन को महंगा बनाता है भले ही धीमा घटक λ₁ h = 2 का उपयोग कर सकता है।
निश्चित बिंदु & आकर्षण बेसिन
dy/dx = f(y) पर लागू यूलर की विधि एक असतत मानचित्र को परिभाषित करती है: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ)।
इस मानचित्र का एक निश्चित बिंदु: y जैसे कि g(y) = y। dy/dx = f(y) पर यूलर के लिए, निश्चित बिंदु f(y) = 0 को संतुष्ट करते हैं — ODE के संतुलन।
निश्चित बिंदु की स्थिरता: यदि |g'(y)| < 1, पास के पुनरावृत्ति y में अभिसरण करते हैं। यदि |g'(y*)| > 1, वे विचलित होते हैं।
g'(y) = 1 + h·f'(y)। एक निश्चित बिंदु y पर: स्थिरता के लिए |1 + h·f'(y)| < 1।
यह बिल्कुल λ = f'(y*) के साथ यूलर स्थिरता स्थिति है — संतुलन पर ODE का रैखिकीकरण।
आकर्षण बेसिन: प्रारंभिक स्थितियों का सेट जो यूलर मानचित्र के तहत y* में अभिसरण करते हैं। अरैखिक प्रणालियों के लिए, बेसिन सीमा परिभाषित करता है कि सिमुलेशन ODE संतुलन को विश्वसनीयतापूर्वक ट्रैक करेगा या दूसरे आकर्षक की ओर विचलित होगा।
सिमुलेशन लूप एक असतत गतिशील प्रणाली है। इसका गुणात्मक व्यवहार — अभिसरण, दोलन, विचलन — ODE की दिशा क्षेत्र की ज्यामिति के सापेक्ष चरण आकार h पर निर्भर करता है।
ज्यामिति को सिमुलेशन डिजाइन से जोड़ना
संख्यात्मक सिमुलेशन की ज्यामिति तीन सवालों तक आती है:
1. स्थिरता क्षेत्र कहां है? यूलर के लिए: डिस्क |1 + hλ| ≤ 1। RK4 के लिए बड़ा। निहित तरीकों के लिए अनबाउंड (पूरा बाएं आधा-समतल)।
2. ODE की eigenvalues कहां हैं? प्रत्येक बिंदु पर f के Jacobian की eigenvalues λ निर्धारित करते हैं कि किस स्थिरता क्षेत्र में hλ होना चाहिए।
3. कौन सा h सभी hλ को क्षेत्र के अंदर रखता है? अधिकतम अनुमत h = (स्थिरता क्षेत्र त्रिज्या) / max|λ|।
कठोर प्रणालियों के लिए: max|λ| विशाल है, स्पष्ट तरीकों के लिए छोटा h बजबरदस्ती करता है। निहित तरीके प्रति चरण महंगे हैं लेकिन बड़े h की अनुमति देते हैं — वे प्रति-चरण लागत को स्थिरता के लिए व्यापार करते हैं।
हैमिंग का अंतर्दृष्टि अनुवाद करता है: संख्यात्मक विधि की पसंद ODE की eigenvalue स्पेक्ट्रम की ज्यामिति के बारे में एक दांव को एन्कोड करता है। उस दांव को स्पष्ट बनाना किसी भी सिमुलेशन में पहला डिजाइन निर्णय है।