साधारण अवकल समीकरण के हल की संख्यात्मक विधियाँ

मुक्त ज्ञानकोश विकिपीडिया से
नेविगेशन पर जाएँ खोज पर जाएँ
विभिन्न संख्यात्मक विधियों द्वारा अवकल समीकरण <math>y'=y, y(0)=1</math> का हल ;
नीला: आयलर विधि द्वारा,
हरा: मध्यबिन्दु विधि द्वारा,
लाल : यथार्थ हल, <math>y=e^t</math>
यहाँ संख्यात्मक विधियों में स्टेप साइज <math>h=1.0</math> लिया गया है।
<math>h=0.25</math> लेकर उपरोक्त समीकरण का हल। इससे स्पष्ट है कि मध्यबिन्दु विधि, आयलर विधि की अपेक्षा तेजी से कन्वर्ज होती है।

साधारण अवकल समीकरणों के हल की संख्यात्मक विधियाँ उन समीकरणों का भी हल दे सकतीं हैं जो वैश्लेषिक रीति से सम्भव नहीं होते। कभी-कभी इसे 'संख्यात्मक समाकलन' भी कहते हैं (यद्यपि 'संख्यात्मक समाकलन' का मतलब संख्यात्मक रीति से निश्चित समाकल का मान निकालना भी होता है।)।

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

समस्या

प्रथम ऑर्डर के अवकल समीकरण का हल वस्तुतः एक आरम्भिक मान समस्या (Initial value problem) है जिसका स्वरूप निम्नलिखित होता है-

<math>y'(t) = f(t,y(t)), \qquad y(t_0)=y_0, \qquad\qquad (1)</math>

विधियाँ

एक चरण वाली विधियाँ

  • (१) आयलर की स्पष्ट (इक्प्लिसिट) विधि
<math> y_{n+1} = y_n + hf(t_n,y_n) </math>
  • (२) आयलर की अस्पष्ट (इम्प्लिसिट) विधि, या आयलर की बैकवर्ड विधि
<math> y_{n+1} = y_n + hf(t_{n+1},y_{n+1}) </math>
  • (३) ट्रैपीजॉयडल विधि
<math>y_{n+1} = y_n + \frac{h}{2}(f_n + f_{n+1})=y_n + \frac{h}{2}(y' _n + y'_{n+1})</math>
  • (४) ह्यून (Heun) की विधि
सबसे पहले यह गणना करें:
<math> y_{n+1}^{(0)}=y_n+h f_n </math>
फिर यह निकालें:
<math> y_{n+1}=y_n+\frac{h}{2}(f(x_n,y_n)+f(x_n+h,y_{n+1}^{(0)})) </math>

बहुचरणी विधियाँ

  • (१) आदम्स-बाशफोर्थ (Adams-Bashforth) विधि
<math>y_{n+1}=y_n+\frac{h}{12} (23 f_n-16 f_{n-1}+5 f_{n-2}) </math>
  • (२) आदम्स-मोल्टन (Adams-Moulton) विधि
<math>y_{n+1}=y_n+\frac{h}{12} (5 f_{n+1}+ 8 f_{n}- f_{n-1}) </math>

प्राक्कलन-संशोधन विधि

इस विधि को 'प्रेडिक्टर-करेक्टर विधि' भी कहते हैं। इसमें एक स्पष्ट विधि (प्रेडिक्टर) द्वारा पहले एक सन्निकट (approx value) निकालते हैं और फिर दूसरी अस्पष्ट विधि (करेक्टर) का उपयोग करते हुए अधिक शुद्ध मान निकाला जाता है।

पॉवर श्रेणी के रूप में सन्निकटन द्वारा

इन्हें भी देखें