संख्यात्मक विभेदन

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

परिमित भिन्नताएँ

Further information: परिमित भिन्नताएँ

इस प्रकार से अधिक सरल विधि परिमित अंतर सन्निकटन का उपयोग करना है।

एक सरल दो-बिंदु अनुमान बिंदुओं (x, f(x)) और (x + h, f (x + h)) के माध्यम से समीप की व्युत्पन्न रेखा के प्रवणता की गणना करना है। छोटी संख्या h चुनना, h x में छोटे परिवर्तन को दर्शाता है, और यह या तो धनात्मक या ऋणात्मक हो सकता है। इस रेखा को रूप देना है

${\frac {f(x+h)-f(x)}{h}}.$

यह अभिव्यक्ति आइजैक न्यूटन का अंतर भागफल है (जिसे प्रथम-क्रम विभाजित अंतर के रूप में भी जाना जाता है)।

इस व्युत्पन्न रेखा का प्रवणता स्पर्शरेखा रेखा के प्रवणता से उस मात्रा में भिन्न होता है जो लगभग h के समानुपाती होता है। जैसे-जैसे h शून्य के समीप पहुंचता है, व्युत्पन्न रेखा का प्रवणता स्पर्शरेखा रेखा के प्रवणता के समीप पहुंचता है। इसलिए, 'f' 'at' 'x' का वास्तविक 'व्युत्पन्न अंतर भागफल के मान की सीमा है क्योंकि व्युत्पन्न रेखाएं स्पर्श रेखा होने के समीप और समीप आती हैं:

$f'(x)=\lim _{h\to 0}{\frac {f(x+h)-f(x)}{h}}.$

चूँकि h के स्थान पर तुरंत 0 प्रतिस्थापन (तर्क) पर ${\frac {0}{0}}$ अनिश्चित रूप, प्राप्त होता है, सीधे व्युत्पन्न की गणना करना सहज ज्ञान युक्त नहीं हो सकता है।

समान रूप से, प्रवणता का अनुमान पदों (x - h) और x का उपयोग करके लगाया जा सकता है।

अतः अन्य दो-बिंदु सूत्र बिंदुओं (x - h, f(x - h)) और (x + h, f (x + h)) के माध्यम से समीप की व्युत्पन्न रेखा के प्रवणता की गणना करना है। इस रेखा का प्रवणता है

${\frac {f(x+h)-f(x-h)}{2h}}.$

इस सूत्र को सममित अंतर भागफल के रूप में जाना जाता है। इस मामले में प्रथम-क्रम त्रुटियाँ रद्द हो जाती हैं, इसलिए इन व्युत्पन्न रेखाओं का प्रवणता स्पर्शरेखा रेखा के प्रवणता से उस मात्रा में भिन्न होता है जो लगभग $h^{2}$ आनुपातिक होता है . इसलिए h के छोटे मानों के लिए यह एकतरफा अनुमान की तुलना में स्पर्शरेखा रेखा का अधिक स्पष्ट अनुमान है। चूंकि , प्रवणता की गणना x पर की जा रही है, जहाँ x पर फलन का मान सम्मिलित नहीं है।

इस प्रकार से अनुमान त्रुटि द्वारा दिया गया है

$R={\frac {-f^{(3)}(c)}{6}}h^{2}$,

जहाँ $c$ के मध्य कुछ बिंदु $x-h$ और $x+h$ है.

अतः संख्याओं का प्रतिनिधित्व करने और सीमित परिशुद्धता में गणना किए जाने के कारण इस त्रुटि में पूर्णांकन त्रुटि सम्मिलित नहीं है।

सममित अंतर भागफल को TI-82, TI-83, TI-84, TI-85 सहित अनेक कैलकुलेटरों में व्युत्पन्न का अनुमान लगाने की विधि के रूप में नियोजित किया जाता है, जो सभी h = 0.001 के साथ इस विधि का उपयोग करते हैं।

चरण आकार

See also: अनुकूली चरण का आकार

पूर्णांकन त्रुटि और सूत्र त्रुटि दोनों के कारण $h$ चुनने में कठिनाई दर्शाने वाला उदाहरण है

वास्तविक में महत्वपूर्ण विचार जब फलन की गणना परिमित परिशुद्धता के फ़्लोटिंग-पॉइंट अंकगणित का उपयोग करके की जाती है, तब चरण आकार का विकल्प h होता है, यदि अधिक छोटा चुना जाता है, तो घटाव में उच्च वृत्ताकार त्रुटि उत्पन्न होती है। वास्तव में, सभी परिमित-अंतर सूत्र असंबद्ध हैं और यदि h अधिक छोटा है तो रद्दीकरण के कारण शून्य का मान उत्पन्न होता है। यदि बहुत बड़ा है, तब व्युत्पन्न रेखा के प्रवणता की गणना अधिक स्पष्ट रूप से की जाती है, किन्तु व्युत्पन्न का उपयोग करके स्पर्शरेखा के प्रवणता का अनुमान व्यर्थ हो सकता है।

इस प्रकार से मूलभूत केंद्रीय अंतरों के लिए, इष्टतम चरण मशीन ईपीएसलॉन का क्यूब-रूट है।

जहाँ x और x + h पर मूल्यांकित संख्यात्मक व्युत्पन्न सूत्र के लिए, h के लिए विकल्प जो उच्च वृत्ताकार त्रुटि उत्पन्न किए बिना छोटा है ${\sqrt {\varepsilon }}x$ (चूंकि तब नहीं जब x = 0), जहां मशीन ईपीएसलॉन ε समान्यतः 2.2×10⁻¹⁶ के क्रम का होता है दोहरी परिशुद्धता के लिए। h के लिए सूत्र जो इष्टतम स्पष्ट ता के लिए सेकेंट त्रुटि के विरुद्ध गोलाई त्रुटि को संतुलित करता है

$h=2{\sqrt {\varepsilon \left|{\frac {f(x)}{f''(x)}}\right|}}$

(चूंकि जब नहीं $f''(x)=0$), और इसे नियोजित करने के लिए फलन के ज्ञान की आवश्यकता होती है।

कंप्यूटर गणनाओं के लिए समस्याएँ और भी बढ़ जाती हैं क्योंकि, जहाँ x आवश्यक रूप से कुछ परिशुद्धता (32 या 64-बिट, आदि) में प्रतिनिधित्व योग्य फ़्लोटिंग-पॉइंट संख्या रखता है, इस प्रकार से जहाँ x + h लगभग निश्चित रूप से उस परिशुद्धता में बिल्कुल प्रतिनिधित्व करने योग्य नहीं है। इसका अर्थ यह है कि x + h को समीप की मशीन-प्रतिनिधित्व योग्य संख्या में (व्रत या निष्कोणन करके) परिवर्तन किया जाता है, जिसके परिणामस्वरूप (x + h) − x, h के समान नहीं होगा; दो फलन मूल्यांकन बिल्कुल अलग नहीं होंगे। इस संबंध में, चूंकि अधिकांश दशमलव अंश बाइनरी में आवर्ती अनुक्रम होते हैं (जैसे 1/3 दशमलव में होता है) प्रतीत होता है कि गोल चरण जैसे कि h = 0.1 बाइनरी में गोल संख्या नहीं होगी; यह 0.000110011001100...₂ है संभावित दृष्टिकोण इस प्रकार है:

h := sqrt(eps) * x;
 xph := x + h;
 dx := xph - x;
 slope := (F(xph) - F(x)) / dx;

चूंकि , कंप्यूटर के साथ, कंपाइलर अनुकूलन सुविधाएं वास्तविक कंप्यूटर अंकगणित के विवरण में भाग लेने में विफल हो सकती हैं और इसके अतिरिक्त यह अनुमान लगाने के लिए गणित के सिद्धांतों को प्रस्तुत करती हैं कि dx और h समान हैं। सी (प्रोग्रामिंग लैंग्वेज) और समान लैंग्वेज के साथ, निर्देश कि xph अस्थिर वेरिएबल है, इसे रोक देता है।

अन्य विधियाँ

उच्च-क्रम विधियाँ

Further information: परिमित अंतर गुणांक

व्युत्पन्न का अनुमान लगाने के लिए उच्च-क्रम विधियाँ, साथ ही उच्च व्युत्पन्न के लिए विधियाँ उपस्तिथ हैं।

प्रथम व्युत्पन्न (एक आयाम में पांच-बिंदु स्टैंसिल) के लिए पांच-बिंदु विधि नीचे दी गई है:

$f'(x)={\frac {-f(x+2h)+8f(x+h)-8f(x-h)+f(x-2h)}{12h}}+{\frac {h^{4}}{30}}f^{(5)}(c),$

जहाँ $c\in [x-2h,x+2h]$.

अन्य स्टैंसिल कॉन्फ़िगरेशन और व्युत्पन्न आदेशों के लिए, परिमित अंतर गुणांक कैलकुलेटर उपकरण है जिसका उपयोग किसी भी व्युत्पन्न क्रम के साथ किसी भी स्टैंसिल के लिए व्युत्पन्न सन्निकटन विधियों को उत्पन्न करने के लिए (निःशंदेह कोई समाधान उपस्तिथ हो) किया जा सकता है।

उच्चतर डेरिवेटिव

इस प्रकार से न्यूटन के अंतर भागफल का उपयोग करते हुए,

$f'(x)=\lim _{h\to 0}{\frac {f(x+h)-f(x)}{h}}$

निम्नलिखित दिखाया जा सकता है (n>0 के लिए):

$f^{(n)}(x)=\lim _{h\to 0}{\frac {1}{h^{n}}}\sum _{k=0}^{n}(-1)^{k+n}{\binom {n}{k}}f(x+kh)$

समष्टि -परिवर्तनीय विधियाँ

इस प्रकार से संख्यात्मक विभेदन के लिए शास्त्रीय परिमित-अंतर सन्निकटन व्यर्थ स्थिति वाले हैं। चूंकि , यदि $f$ होलोमोर्फिक फलन है, जो की वास्तविक रेखा पर वास्तविक-मूल्यवान है, जिसका मूल्यांकन समीप के समष्टि विमान में बिंदुओं $x$ पर किया जा सकता है , फिर संख्यात्मक स्थिरता विधियाँ हैं। अतः उदाहरण के लिए, प्रथम व्युत्पन्न की गणना समष्टि -चरण व्युत्पन्न सूत्र द्वारा की जा सकती है:

$f^{\prime }(x)={\frac {\Im (f(x+\mathrm {i} h))}{h}}+O(h^{2}),\quad \mathrm {i^{2}} :=-1.$

विभिन्न स्थितियों के लिए स्पष्ट डेरिवेटिव प्राप्त करने के लिए अनुशंसित चरण आकार $h=10^{-200}$ है . यह सूत्र टेलर श्रृंखला विस्तार द्वारा प्राप्त किया जा सकता है:

$f(x+\mathrm {i} h)=f(x)+\mathrm {i} hf^{\prime }(x)-h^{2}f''(x)/2!-\mathrm {i} h^{3}f^{(3)}(x)/3!+\cdots .$

समष्टि -चरण व्युत्पन्न सूत्र केवल प्रथम-क्रम व्युत्पन्न की गणना के लिए मान्य है। किसी भी क्रम के डेरिवेटिव की गणना के लिए उपरोक्त का सामान्यीकरण मल्टीकॉम्प्लेक्स संख्याएँ को नियोजित करता है, जिसके परिणामस्वरूप मल्टीकॉम्प्लेक्स डेरिवेटिव होते हैं।

$f^{(n)}(x)\approx {\frac {{\mathcal {C}}_{n^{2}-1}^{(n)}(f(x+\mathrm {i} ^{(1)}h+\ldots +\mathrm {i} ^{(n)}h))}{h^{n}}}$

जहां $\mathrm {i} ^{(k)}$ मल्टीकॉम्प्लेक्स काल्पनिक इकाइयों को दर्शाता है; $\mathrm {i} ^{(1)}\equiv \mathrm {i}$, ${\mathcal {C}}_{k}^{(n)}$ ऑपरेटर स्तर ${\mathcal {C}}_{0}^{(n)}$ की मल्टीकॉम्प्लेक्स संख्या के $k$th घटक को निकालता है उदाहरण के लिए $n$ वास्तविक घटक को निकालता है और ${\mathcal {C}}_{n^{2}-1}^{(n)}$अंतिम, "अधिक काल्पनिक" घटक को निकालता है। इस विधि को मिश्रित डेरिवेटिव पर प्रस्तुत किया जा सकता है, जैसे दूसरे क्रम के व्युत्पन्न के लिए

${\frac {\partial ^{2}f(x,y)}{\partial x\,\partial y}}\approx {\frac {{\mathcal {C}}_{3}^{(2)}(f(x+\mathrm {i} ^{(1)}h,y+\mathrm {i} ^{(2)}h))}{h^{2}}}$

अतः मल्टीकॉम्प्लेक्स अंकगणित का C++ कार्यान्वयन उपलब्ध है।

सामान्य रूप से, किसी भी क्रम के व्युत्पन्न की गणना कॉची के अभिन्न सूत्र का उपयोग करके की जा सकती है:

$f^{(n)}(a)={\frac {n!}{2\pi i}}\oint _{\gamma }{\frac {f(z)}{(z-a)^{n+1}}}\,\mathrm {d} z,$

जहां एकीकरण संख्यात्मक एकीकरण किया जाता है।

इस प्रकार से संख्यात्मक विभेदन के लिए समष्टि वेरिएबल का उपयोग 1967 में लिनेस और मोलर द्वारा प्रारंभ किया गया था। उनका एल्गोरिदम उच्च-क्रम डेरिवेटिव पर प्रस्तुत होता है।

समष्टि लाप्लास परिवर्तन के संख्यात्मक व्युत्क्रम पर आधारित विधि एबेट और डबनेर द्वारा विकसित की गई थी। एल्गोरिदम जिसका उपयोग विधि या फलन के चरित्र के पश्चात में ज्ञान की आवश्यकता के बिना किया जा सकता है, फोर्नबर्ग द्वारा विकसित किया गया था।

विभेदक चतुर्भुज

इस प्रकार से विभेदक चतुर्भुज फलन मानों के भारित योग का उपयोग करके डेरिवेटिव का अनुमान है। विभेदक चतुर्भुज व्यावहारिक रुचि का है क्योंकि यह किसी को ध्वनि वाले डेटा से डेरिवेटिव की गणना करने की अनुमति देता है। अतः नाम चतुर्भुज के अनुरूप है, जिसका अर्थ है संख्यात्मक एकीकरण है, जहां भारित रकम का उपयोग सिम्पसन की विधि या ट्रेपेज़ॉइडल नियम जैसी विधियों में किया जाता है। वज़न गुणांक निर्धारित करने के लिए विभिन्न विधियाँ हैं, अतः उदाहरण के लिए, सविट्ज़की-गोले फ़िल्टर है। चूंकि आंशिक अवकल समीकरणों को हल करने के लिए विभेदक चतुर्भुज का उपयोग किया जाता है। ध्वनि वाले डेटा से डेरिवेटिव की गणना के लिए अनेक विधि हैं।

यह भी देखें

• स्वचालित विभेदीकरण
• पांच-बिंदु स्टेंसिल
• सवित्ज़की-व्रत फ़िल्टर
• संख्यात्मक एकीकरण – Family of algorithms for finding the definite integral of a function
• संख्यात्मक साधारण अवकल समीकरण
• संख्यात्मक चौरसाई और विभेदन
• संख्यात्मक-विश्लेषण सॉफ्टवेयर की सूची

बाहरी संबंध

• Numerical Differentiation from wolfram.com
• Numerical Differentiation Resources: Textbook notes, PPT, Worksheets, Audiovisual YouTube Lectures at Numerical Methods for STEM Undergraduate
• Fortran code for the numerical differentiation of a function using Neville's process to extrapolate from a sequence of simple polynomial approximations.
• NAG Library numerical differentiation routines
• http://graphulator.com Online numerical graphing calculator with calculus function.
• Boost. Math numerical differentiation, including finite differencing and the complex step derivative
• Complex Step Differentiation
• Differentiation With(out) a Difference by Nicholas Higham, SIAM News.
• findiff Python project