संख्यात्मक विभेदन: Difference between revisions

From Vigyanwiki
(Created page with "{{Short description|Use of numerical analysis to estimate derivatives of functions}} संख्यात्मक विश्लेषण में, संख्या...")
 
No edit summary
Line 9: Line 9:
सबसे सरल विधि परिमित अंतर सन्निकटन का उपयोग करना है।
सबसे सरल विधि परिमित अंतर सन्निकटन का उपयोग करना है।


एक सरल दो-बिंदु अनुमान बिंदुओं (x, f(x)) और (x + h, f (x + h)) के माध्यम से पास की [[छेदक रेखा]] के ढलान की गणना करना है।<ref>Richard L. Burden, J. Douglas Faires (2000), ''Numerical Analysis'', (7th Ed),  Brooks/Cole. {{isbn|0-534-38216-9}}.</ref> एक छोटी संख्या h चुनना, h x में एक छोटे परिवर्तन को दर्शाता है, और यह या तो सकारात्मक या नकारात्मक हो सकता है। इस रेखा का ढलान है
एक सरल दो-बिंदु अनुमान बिंदुओं (x, f(x)) और (x + h, f (x + h)) के माध्यम से पास की [[छेदक रेखा]] के ढलान की गणना करना है।<ref>Richard L. Burden, J. Douglas Faires (2000), ''Numerical Analysis'', (7th Ed),  Brooks/Cole. {{isbn|0-534-38216-9}}.</ref> छोटी संख्या h चुनना, h x में छोटे परिवर्तन को दर्शाता है, और यह या तो सकारात्मक या नकारात्मक हो सकता है। इस रेखा का ढलान है
: <math>\frac{f(x + h) - f(x)}{h}.</math>
: <math>\frac{f(x + h) - f(x)}{h}.</math>
यह अभिव्यक्ति [[आइजैक न्यूटन]] का [[अंतर भागफल]] है (जिसे प्रथम-क्रम [[विभाजित अंतर]] के रूप में भी जाना जाता है)।
यह अभिव्यक्ति [[आइजैक न्यूटन]] का [[अंतर भागफल]] है (जिसे प्रथम-क्रम [[विभाजित अंतर]] के रूप में भी जाना जाता है)।


इस छेदक रेखा का ढलान स्पर्शरेखा रेखा के ढलान से उस मात्रा में भिन्न होता है जो लगभग h के समानुपाती होता है। जैसे-जैसे h शून्य के करीब पहुंचता है, छेदक रेखा का ढलान स्पर्शरेखा रेखा के ढलान के करीब पहुंचता है। इसलिए, 'f' 'at' 'x' का वास्तविक 'व्युत्पन्न अंतर भागफल के मान की सीमा है क्योंकि छेदक रेखाएं एक स्पर्श रेखा होने के करीब और करीब आती हैं:
इस छेदक रेखा का ढलान स्पर्शरेखा रेखा के ढलान से उस मात्रा में भिन्न होता है जो लगभग h के समानुपाती होता है। जैसे-जैसे h शून्य के करीब पहुंचता है, छेदक रेखा का ढलान स्पर्शरेखा रेखा के ढलान के करीब पहुंचता है। इसलिए, 'f' 'at' 'x' का वास्तविक 'व्युत्पन्न अंतर भागफल के मान की सीमा है क्योंकि छेदक रेखाएं स्पर्श रेखा होने के करीब और करीब आती हैं:
: <math>f'(x) = \lim_{h \to 0} \frac{f(x + h) - f(x)}{h}.</math>
: <math>f'(x) = \lim_{h \to 0} \frac{f(x + h) - f(x)}{h}.</math>
चूँकि h के स्थान पर तुरंत [[प्रतिस्थापन (तर्क)]] 0 का परिणाम आता है <math>\frac{0}{0}</math> [[अनिश्चित रूप]], सीधे व्युत्पन्न की गणना करना सहज ज्ञान युक्त नहीं हो सकता है।
चूँकि h के स्थान पर तुरंत [[प्रतिस्थापन (तर्क)]] 0 का परिणाम आता है <math>\frac{0}{0}</math> [[अनिश्चित रूप]], सीधे व्युत्पन्न की गणना करना सहज ज्ञान युक्त नहीं हो सकता है।
Line 35: Line 35:
==चरण आकार==
==चरण आकार==
{{see also|Adaptive step size}}
{{see also|Adaptive step size}}
[[Image:AbsoluteErrorNumericalDifferentiationExample.png|thumb|300px|चुनने की कठिनाई दर्शाने वाला उदाहरण <math>h</math> पूर्णांकन त्रुटि और सूत्र त्रुटि दोनों के कारण]]व्यवहार में एक महत्वपूर्ण विचार जब फ़ंक्शन की गणना परिमित परिशुद्धता के [[फ़्लोटिंग-पॉइंट अंकगणित]] का उपयोग करके की जाती है, तो चरण आकार का विकल्प होता है, एच। यदि बहुत छोटा चुना जाता है, तो घटाव में बड़ी गोलाई त्रुटि उत्पन्न होगी। वास्तव में, सभी परिमित-अंतर सूत्र असंबद्ध हैं<ref name=Fornberg1>Numerical Differentiation of Analytic Functions, B Fornberg – ACM Transactions on Mathematical Software (TOMS), 1981.</ref> और यदि h काफी छोटा है तो रद्दीकरण के कारण शून्य का मान उत्पन्न होगा।<ref name=SquireTrapp1>Using Complex Variables to Estimate Derivatives of Real Functions, W. Squire, G. Trapp – SIAM REVIEW, 1998.</ref> यदि बहुत बड़ा है, तो छेदक रेखा के ढलान की गणना अधिक सटीक रूप से की जाएगी, लेकिन छेदक का उपयोग करके स्पर्शरेखा के ढलान का अनुमान बदतर हो सकता है।<ref name="edo2021">{{Cite book|url= http://flowlab.groups.et.byu.net/mdobook.pdf|title=इंजीनियरिंग डिज़ाइन अनुकूलन|last1=Martins|first1=Joaquim R. R. A.|last2=Ning|first2=Andrew|date=2021-10-01|publisher=Cambridge University Press|isbn=978-1108833417|language=en}}</ref>
[[Image:AbsoluteErrorNumericalDifferentiationExample.png|thumb|300px|चुनने की कठिनाई दर्शाने वाला उदाहरण <math>h</math> पूर्णांकन त्रुटि और सूत्र त्रुटि दोनों के कारण]]व्यवहार में महत्वपूर्ण विचार जब फ़ंक्शन की गणना परिमित परिशुद्धता के [[फ़्लोटिंग-पॉइंट अंकगणित]] का उपयोग करके की जाती है, तो चरण आकार का विकल्प होता है, एच। यदि बहुत छोटा चुना जाता है, तो घटाव में बड़ी गोलाई त्रुटि उत्पन्न होगी। वास्तव में, सभी परिमित-अंतर सूत्र असंबद्ध हैं<ref name=Fornberg1>Numerical Differentiation of Analytic Functions, B Fornberg – ACM Transactions on Mathematical Software (TOMS), 1981.</ref> और यदि h काफी छोटा है तो रद्दीकरण के कारण शून्य का मान उत्पन्न होगा।<ref name=SquireTrapp1>Using Complex Variables to Estimate Derivatives of Real Functions, W. Squire, G. Trapp – SIAM REVIEW, 1998.</ref> यदि बहुत बड़ा है, तो छेदक रेखा के ढलान की गणना अधिक सटीक रूप से की जाएगी, लेकिन छेदक का उपयोग करके स्पर्शरेखा के ढलान का अनुमान बदतर हो सकता है।<ref name="edo2021">{{Cite book|url= http://flowlab.groups.et.byu.net/mdobook.pdf|title=इंजीनियरिंग डिज़ाइन अनुकूलन|last1=Martins|first1=Joaquim R. R. A.|last2=Ning|first2=Andrew|date=2021-10-01|publisher=Cambridge University Press|isbn=978-1108833417|language=en}}</ref>
बुनियादी केंद्रीय अंतरों के लिए, इष्टतम चरण [[मशीन ईपीएसलॉन]] का क्यूब-रूट है।<ref>Sauer, Timothy (2012). ''Numerical Analysis''. Pearson. p.248.</ref>
बुनियादी केंद्रीय अंतरों के लिए, इष्टतम चरण [[मशीन ईपीएसलॉन]] का क्यूब-रूट है।<ref>Sauer, Timothy (2012). ''Numerical Analysis''. Pearson. p.248.</ref>
x और x + h पर मूल्यांकित संख्यात्मक व्युत्पन्न सूत्र के लिए, h के लिए एक विकल्प जो बड़ी गोलाई त्रुटि उत्पन्न किए बिना छोटा है <math>\sqrt{\varepsilon} x</math> (हालांकि तब नहीं जब x = 0), जहां मशीन ईपीएसलॉन ε आम तौर पर 2.2 के क्रम का होता है{{e|−16}} दोहरी परिशुद्धता के लिए।<ref>Following ''[[Numerical Recipes]] in C'', [http://www.nrbook.com/a/bookcpdf/c5-7.pdf Chapter 5.7].</ref> एच के लिए एक सूत्र जो इष्टतम सटीकता के लिए सेकेंट त्रुटि के विरुद्ध गोलाई त्रुटि को संतुलित करता है<ref>[http://www.uio.no/studier/emner/matnat/math/MAT-INF1100/h10/kompendiet/kap11.pdf p. 263].</ref>
 
x और x + h पर मूल्यांकित संख्यात्मक व्युत्पन्न सूत्र के लिए, h के लिए विकल्प जो बड़ी गोलाई त्रुटि उत्पन्न किए बिना छोटा है <math>\sqrt{\varepsilon} x</math> (हालांकि तब नहीं जब x = 0), जहां मशीन ईपीएसलॉन ε आम तौर पर 2.2 के क्रम का होता है{{e|−16}} दोहरी परिशुद्धता के लिए।<ref>Following ''[[Numerical Recipes]] in C'', [http://www.nrbook.com/a/bookcpdf/c5-7.pdf Chapter 5.7].</ref> एच के लिए सूत्र जो इष्टतम सटीकता के लिए सेकेंट त्रुटि के विरुद्ध गोलाई त्रुटि को संतुलित करता है<ref>[http://www.uio.no/studier/emner/matnat/math/MAT-INF1100/h10/kompendiet/kap11.pdf p. 263].</ref>
: <math>h = 2\sqrt{\varepsilon\left|\frac{f(x)}{f''(x)}\right|}</math>
: <math>h = 2\sqrt{\varepsilon\left|\frac{f(x)}{f''(x)}\right|}</math>
(हालांकि कब नहीं <math>f''(x) = 0</math>), और इसे नियोजित करने के लिए फ़ंक्शन के ज्ञान की आवश्यकता होगी।
(हालांकि कब नहीं <math>f''(x) = 0</math>), और इसे नियोजित करने के लिए फ़ंक्शन के ज्ञान की आवश्यकता होगी।


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


==अन्य विधियाँ==
==अन्य विधियाँ==
Line 58: Line 59:
कहाँ <math>c \in [x - 2h, x + 2h]</math>.
कहाँ <math>c \in [x - 2h, x + 2h]</math>.


अन्य स्टैंसिल कॉन्फ़िगरेशन और व्युत्पन्न आदेशों के लिए, [http://web.media.mit.edu/~crtaylor/calculator.html परिमित अंतर गुणांक कैलकुलेटर] एक उपकरण है जिसका उपयोग किसी भी व्युत्पन्न क्रम के साथ किसी भी स्टैंसिल के लिए व्युत्पन्न सन्निकटन विधियों को उत्पन्न करने के लिए किया जा सकता है (बशर्ते कोई समाधान मौजूद हो)।
अन्य स्टैंसिल कॉन्फ़िगरेशन और व्युत्पन्न आदेशों के लिए, [http://web.media.mit.edu/~crtaylor/calculator.html परिमित अंतर गुणांक कैलकुलेटर] उपकरण है जिसका उपयोग किसी भी व्युत्पन्न क्रम के साथ किसी भी स्टैंसिल के लिए व्युत्पन्न सन्निकटन विधियों को उत्पन्न करने के लिए किया जा सकता है (बशर्ते कोई समाधान मौजूद हो)।


=== उच्चतर डेरिवेटिव ===
=== उच्चतर डेरिवेटिव ===
Line 67: Line 68:


:<math>f^{(n)}(x) = \lim_{h\to 0} \frac{1}{h^n} \sum_{k=0}^n (-1)^{k+n} \binom{n}{k} f(x + kh)</math>
:<math>f^{(n)}(x) = \lim_{h\to 0} \frac{1}{h^n} \sum_{k=0}^n (-1)^{k+n} \binom{n}{k} f(x + kh)</math>
<!-- [the following clearly has a wrong sign for n = 1]
However, a factor seems to be missing:
: <math>f^{(n)}(x)= (-1)^n \lim_{h\to 0} \frac{1}{h^n} \sum_{k=0}^n (-1)^{k+1} \binom{n}{k} f(x + kh)</math>
-->
==जटिल-परिवर्तनीय विधियाँ==
==जटिल-परिवर्तनीय विधियाँ==


संख्यात्मक विभेदन के लिए शास्त्रीय परिमित-अंतर सन्निकटन ख़राब स्थिति वाले हैं। हालांकि, यदि <math>f</math> एक [[होलोमोर्फिक फ़ंक्शन]] है, जो वास्तविक रेखा पर वास्तविक-मूल्यवान है, जिसका मूल्यांकन पास के जटिल विमान में बिंदुओं पर किया जा सकता है <math>x</math>, फिर [[संख्यात्मक स्थिरता]] विधियाँ हैं। उदाहरण के लिए,<ref name=SquireTrapp1/>पहले व्युत्पन्न की गणना जटिल-चरण व्युत्पन्न सूत्र द्वारा की जा सकती है:<ref>{{cite journal | last1 = Martins | first1 = J. R. R. A. | first2 = P. | last2 = Sturdza | first3 = J. J. | last3 = Alonso | year = 2003 | citeseerx=10.1.1.141.8002 | title = जटिल-चरण व्युत्पन्न सन्निकटन| journal = ACM Transactions on Mathematical Software | volume = 29 | issue = 3 | pages = 245–262 | doi=10.1145/838250.838251| s2cid = 7022422 }}</ref><ref>[https://sinews.siam.org/Details-Page/differentiation-without-a-difference Differentiation With(out) a Difference] by [[Nicholas Higham]] </ref><ref>[https://blogs.mathworks.com/cleve/2013/10/14/complex-step-differentiation/ article] from [[MathWorks]] blog, posted by [[Cleve Moler]]</ref>
संख्यात्मक विभेदन के लिए शास्त्रीय परिमित-अंतर सन्निकटन ख़राब स्थिति वाले हैं। हालांकि, यदि <math>f</math> [[होलोमोर्फिक फ़ंक्शन]] है, जो वास्तविक रेखा पर वास्तविक-मूल्यवान है, जिसका मूल्यांकन पास के जटिल विमान में बिंदुओं पर किया जा सकता है <math>x</math>, फिर [[संख्यात्मक स्थिरता]] विधियाँ हैं। उदाहरण के लिए,<ref name=SquireTrapp1/> पहले व्युत्पन्न की गणना जटिल-चरण व्युत्पन्न सूत्र द्वारा की जा सकती है:<ref>{{cite journal | last1 = Martins | first1 = J. R. R. A. | first2 = P. | last2 = Sturdza | first3 = J. J. | last3 = Alonso | year = 2003 | citeseerx=10.1.1.141.8002 | title = जटिल-चरण व्युत्पन्न सन्निकटन| journal = ACM Transactions on Mathematical Software | volume = 29 | issue = 3 | pages = 245–262 | doi=10.1145/838250.838251| s2cid = 7022422 }}</ref><ref>[https://sinews.siam.org/Details-Page/differentiation-without-a-difference Differentiation With(out) a Difference] by [[Nicholas Higham]] </ref><ref>[https://blogs.mathworks.com/cleve/2013/10/14/complex-step-differentiation/ article] from [[MathWorks]] blog, posted by [[Cleve Moler]]</ref>
:<math>f^\prime(x) = \frac{\Im(f(x + \mathrm{i}h))}{h} + O(h^2), \quad \mathrm{i^2}:=-1.</math>
:<math>f^\prime(x) = \frac{\Im(f(x + \mathrm{i}h))}{h} + O(h^2), \quad \mathrm{i^2}:=-1.</math>
विभिन्न स्थितियों के लिए सटीक डेरिवेटिव प्राप्त करने के लिए अनुशंसित चरण आकार है <math>h = 10^{-200}</math>.<ref name="edo2021"/>यह सूत्र [[टेलर श्रृंखला]] विस्तार द्वारा प्राप्त किया जा सकता है:
विभिन्न स्थितियों के लिए सटीक डेरिवेटिव प्राप्त करने के लिए अनुशंसित चरण आकार है <math>h = 10^{-200}</math>.<ref name="edo2021"/> यह सूत्र [[टेलर श्रृंखला]] विस्तार द्वारा प्राप्त किया जा सकता है:


:<math>f(x+\mathrm{i}h) = f(x)+\mathrm{i}hf^\prime(x)-h^2f''(x)/2!-\mathrm{i}h^3f^{(3)}(x)/3!+\cdots.</math>
:<math>f(x+\mathrm{i}h) = f(x)+\mathrm{i}hf^\prime(x)-h^2f''(x)/2!-\mathrm{i}h^3f^{(3)}(x)/3!+\cdots.</math>
जटिल-चरण व्युत्पन्न सूत्र केवल प्रथम-क्रम व्युत्पन्न की गणना के लिए मान्य है। किसी भी क्रम के डेरिवेटिव की गणना के लिए उपरोक्त का सामान्यीकरण [[मल्टीकॉम्प्लेक्स संख्याएँ]] को नियोजित करता है, जिसके परिणामस्वरूप मल्टीकॉम्प्लेक्स डेरिवेटिव होते हैं।<ref>{{Cite web |url=http://russell.ae.utexas.edu/FinalPublications/ConferencePapers/2010Feb_SanDiego_AAS-10-218_mulicomplex.pdf |title=संग्रहीत प्रति|access-date=2012-11-24 |archive-url=https://web.archive.org/web/20140109145840/http://russell.ae.utexas.edu/FinalPublications/ConferencePapers/2010Feb_SanDiego_AAS-10-218_mulicomplex.pdf |archive-date=2014-01-09 |url-status=dead }}</ref>
जटिल-चरण व्युत्पन्न सूत्र केवल प्रथम-क्रम व्युत्पन्न की गणना के लिए मान्य है। किसी भी क्रम के डेरिवेटिव की गणना के लिए उपरोक्त का सामान्यीकरण [[मल्टीकॉम्प्लेक्स संख्याएँ]] को नियोजित करता है, जिसके परिणामस्वरूप मल्टीकॉम्प्लेक्स डेरिवेटिव होते हैं।<ref>{{Cite web |url=http://russell.ae.utexas.edu/FinalPublications/ConferencePapers/2010Feb_SanDiego_AAS-10-218_mulicomplex.pdf |title=संग्रहीत प्रति|access-date=2012-11-24 |archive-url=https://web.archive.org/web/20140109145840/http://russell.ae.utexas.edu/FinalPublications/ConferencePapers/2010Feb_SanDiego_AAS-10-218_mulicomplex.pdf |archive-date=2014-01-09 |url-status=dead }}</ref><ref>{{cite journal | last1 = Lantoine | first1 = G. | last2 = Russell | first2 = R. P. | last3 = Dargent | first3 = Th. | title = उच्च-क्रम डेरिवेटिव की स्वचालित गणना के लिए मल्टीकॉम्प्लेक्स चर का उपयोग करना| journal = ACM Trans. Math. Softw. | volume = 38 | year = 2012 | issue = 3 | pages = 1–21 | doi = 10.1145/2168773.2168774| s2cid = 16253562 }}</ref><ref>{{cite web | last1 = Verheyleweghen | first = A. | title = मल्टी-कॉम्प्लेक्स चरण विधि का उपयोग करके उच्च-क्रम डेरिवेटिव की गणना| year = 2014 | url = http://folk.ntnu.no/preisig/HAP_Specials/AdvancedSimulation_files/2014/AdvSim-2014__Verheule_Adrian_Complex_differenetiation.pdf}}</ref>
<ref>{{cite journal | last1 = Lantoine | first1 = G. | last2 = Russell | first2 = R. P. | last3 = Dargent | first3 = Th. | title = उच्च-क्रम डेरिवेटिव की स्वचालित गणना के लिए मल्टीकॉम्प्लेक्स चर का उपयोग करना| journal = ACM Trans. Math. Softw. | volume = 38 | year = 2012 | issue = 3 | pages = 1–21 | doi = 10.1145/2168773.2168774| s2cid = 16253562 }}</ref>
 
<ref>{{cite web | last1 = Verheyleweghen | first = A. | title = मल्टी-कॉम्प्लेक्स चरण विधि का उपयोग करके उच्च-क्रम डेरिवेटिव की गणना| year = 2014 | url = http://folk.ntnu.no/preisig/HAP_Specials/AdvancedSimulation_files/2014/AdvSim-2014__Verheule_Adrian_Complex_differenetiation.pdf}}</ref>
:<math>f^{(n)}(x) \approx \frac{\mathcal{C}^{(n)}_{n^2-1}(f(x + \mathrm{i}^{(1)} h + \ldots + \mathrm{i}^{(n)} h))}{h^n}</math>
:<math>f^{(n)}(x) \approx \frac{\mathcal{C}^{(n)}_{n^2-1}(f(x + \mathrm{i}^{(1)} h + \ldots + \mathrm{i}^{(n)} h))}{h^n}</math>
जहां <math>\mathrm{i}^{(k)}</math> मल्टीकॉम्प्लेक्स काल्पनिक इकाइयों को निरूपित करें; <math>\mathrm{i}^{(1)} \equiv \mathrm{i}</math>. <math>\mathcal{C}^{(n)}_k</math> h> ऑपरेटर निकालता है <math>k</math>स्तर के एक बहुसंकेतन संख्या का वां घटक <math>n</math>, जैसे, <math>\mathcal{C}^{(n)}_0</math> वास्तविक घटक निकालता है और <math>\mathcal{C}^{(n)}_{n^2-1}</math> अंतिम, "सबसे काल्पनिक" घटक निकालता है। इस विधि को मिश्रित डेरिवेटिव पर लागू किया जा सकता है, जैसे दूसरे क्रम के व्युत्पन्न के लिए
जहां <math>\mathrm{i}^{(k)}</math> मल्टीकॉम्प्लेक्स काल्पनिक इकाइयों को निरूपित करें; <math>\mathrm{i}^{(1)} \equiv \mathrm{i}</math>. <math>\mathcal{C}^{(n)}_k</math> h> ऑपरेटर निकालता है <math>k</math>स्तर के बहुसंकेतन संख्या का वां घटक <math>n</math>, जैसे, <math>\mathcal{C}^{(n)}_0</math> वास्तविक घटक निकालता है और <math>\mathcal{C}^{(n)}_{n^2-1}</math> अंतिम, "सबसे काल्पनिक" घटक निकालता है। इस विधि को मिश्रित डेरिवेटिव पर लागू किया जा सकता है, जैसे दूसरे क्रम के व्युत्पन्न के लिए


:<math>\frac{\partial^2 f(x, y)}{\partial x \,\partial y} \approx \frac{\mathcal{C}^{(2)}_3(f(x + \mathrm{i}^{(1)} h, y + \mathrm{i}^{(2)} h))}{h^2}</math>
:<math>\frac{\partial^2 f(x, y)}{\partial x \,\partial y} \approx \frac{\mathcal{C}^{(2)}_3(f(x + \mathrm{i}^{(1)} h, y + \mathrm{i}^{(2)} h))}{h^2}</math>
मल्टीकॉम्प्लेक्स अंकगणित का C++ कार्यान्वयन उपलब्ध है।<ref>{{cite web | last1 = Bell | first1 =  I. H.| title = एमसीएक्स (मल्टीकॉम्प्लेक्स बीजगणित लाइब्रेरी)| website = [[GitHub]]| year = 2019 | url = https://github.com/ianhbell/mcx}}</ref>
मल्टीकॉम्प्लेक्स अंकगणित का C++ कार्यान्वयन उपलब्ध है।<ref>{{cite web | last1 = Bell | first1 =  I. H.| title = एमसीएक्स (मल्टीकॉम्प्लेक्स बीजगणित लाइब्रेरी)| website = [[GitHub]]| year = 2019 | url = https://github.com/ianhbell/mcx}}</ref>
सामान्य तौर पर, किसी भी क्रम के व्युत्पन्न की गणना कॉची के अभिन्न सूत्र का उपयोग करके की जा सकती है:<ref>Ablowitz, M. J., Fokas, A. S.,(2003). Complex variables: introduction and applications. [[Cambridge University Press]]. Check theorem 2.6.2</ref>
सामान्य तौर पर, किसी भी क्रम के व्युत्पन्न की गणना कॉची के अभिन्न सूत्र का उपयोग करके की जा सकती है:<ref>Ablowitz, M. J., Fokas, A. S.,(2003). Complex variables: introduction and applications. [[Cambridge University Press]]. Check theorem 2.6.2</ref>
: <math>f^{(n)}(a) = \frac{n!}{2\pi i} \oint_\gamma \frac{f(z)}{(z - a)^{n+1}} \,\mathrm{d}z,</math>
: <math>f^{(n)}(a) = \frac{n!}{2\pi i} \oint_\gamma \frac{f(z)}{(z - a)^{n+1}} \,\mathrm{d}z,</math>
Line 95: Line 89:
संख्यात्मक विभेदन के लिए जटिल चर का उपयोग 1967 में लिनेस और मोलर द्वारा शुरू किया गया था।<ref name=LynessMoler1>{{cite journal | first1 = J. N. | last1 = Lyness | first2 = C. B. | last2 = Moler | title = विश्लेषणात्मक कार्यों का संख्यात्मक विभेदन| journal = SIAM J. Numer. Anal. | volume = 4 | year = 1967 | issue = 2 | pages = 202–210 | doi=10.1137/0704019| bibcode = 1967SJNA....4..202L }}</ref> उनका एल्गोरिदम उच्च-क्रम डेरिवेटिव पर लागू होता है।
संख्यात्मक विभेदन के लिए जटिल चर का उपयोग 1967 में लिनेस और मोलर द्वारा शुरू किया गया था।<ref name=LynessMoler1>{{cite journal | first1 = J. N. | last1 = Lyness | first2 = C. B. | last2 = Moler | title = विश्लेषणात्मक कार्यों का संख्यात्मक विभेदन| journal = SIAM J. Numer. Anal. | volume = 4 | year = 1967 | issue = 2 | pages = 202–210 | doi=10.1137/0704019| bibcode = 1967SJNA....4..202L }}</ref> उनका एल्गोरिदम उच्च-क्रम डेरिवेटिव पर लागू होता है।


जटिल [[लाप्लास परिवर्तन]] के संख्यात्मक व्युत्क्रम पर आधारित एक विधि एबेट और डबनेर द्वारा विकसित की गई थी।<ref>{{cite journal | title = कार्यों की शक्ति श्रृंखला विस्तार उत्पन्न करने की एक नई विधि| first1 = J | last1 = Abate | first2 = H | last2 = Dubner | journal = SIAM J. Numer. Anal. | volume =5 | issue = 1 | pages = 102–112 |date=March 1968 | doi=10.1137/0705008| bibcode = 1968SJNA....5..102A }}</ref> एक एल्गोरिदम जिसका उपयोग विधि या फ़ंक्शन के चरित्र के बारे में ज्ञान की आवश्यकता के बिना किया जा सकता है, फोर्नबर्ग द्वारा विकसित किया गया था।<ref name=Fornberg1/>
जटिल [[लाप्लास परिवर्तन]] के संख्यात्मक व्युत्क्रम पर आधारित विधि एबेट और डबनेर द्वारा विकसित की गई थी।<ref>{{cite journal | title = कार्यों की शक्ति श्रृंखला विस्तार उत्पन्न करने की एक नई विधि| first1 = J | last1 = Abate | first2 = H | last2 = Dubner | journal = SIAM J. Numer. Anal. | volume =5 | issue = 1 | pages = 102–112 |date=March 1968 | doi=10.1137/0705008| bibcode = 1968SJNA....5..102A }}</ref> एल्गोरिदम जिसका उपयोग विधि या फ़ंक्शन के चरित्र के बारे में ज्ञान की आवश्यकता के बिना किया जा सकता है, फोर्नबर्ग द्वारा विकसित किया गया था।<ref name=Fornberg1/>
 
 
==विभेदक चतुर्भुज==
==विभेदक चतुर्भुज==
विभेदक चतुर्भुज फ़ंक्शन मानों के भारित योग का उपयोग करके डेरिवेटिव का अनुमान है।<ref>Differential Quadrature and Its Application in Engineering: Engineering Applications, Chang Shu, Springer, 2000, {{isbn|978-1-85233-209-9}}.</ref><ref>Advanced Differential Quadrature Methods, Yingyan Zhang, CRC Press, 2009, {{isbn|978-1-4200-8248-7}}.</ref> विभेदक चतुर्भुज व्यावहारिक रुचि का है क्योंकि यह किसी को शोर वाले डेटा से डेरिवेटिव की गणना करने की अनुमति देता है। नाम चतुर्भुज के अनुरूप है, जिसका अर्थ है संख्यात्मक एकीकरण, जहां भारित रकम का उपयोग सिम्पसन की विधि या ट्रेपेज़ॉइडल नियम जैसी विधियों में किया जाता है। वज़न गुणांक निर्धारित करने के लिए विभिन्न विधियाँ हैं, उदाहरण के लिए, सविट्ज़की-गोले फ़िल्टर। आंशिक अवकल समीकरणों को हल करने के लिए विभेदक चतुर्भुज का उपयोग किया जाता है।
विभेदक चतुर्भुज फ़ंक्शन मानों के भारित योग का उपयोग करके डेरिवेटिव का अनुमान है।<ref>Differential Quadrature and Its Application in Engineering: Engineering Applications, Chang Shu, Springer, 2000, {{isbn|978-1-85233-209-9}}.</ref><ref>Advanced Differential Quadrature Methods, Yingyan Zhang, CRC Press, 2009, {{isbn|978-1-4200-8248-7}}.</ref> विभेदक चतुर्भुज व्यावहारिक रुचि का है क्योंकि यह किसी को शोर वाले डेटा से डेरिवेटिव की गणना करने की अनुमति देता है। नाम चतुर्भुज के अनुरूप है, जिसका अर्थ है संख्यात्मक एकीकरण, जहां भारित रकम का उपयोग सिम्पसन की विधि या ट्रेपेज़ॉइडल नियम जैसी विधियों में किया जाता है। वज़न गुणांक निर्धारित करने के लिए विभिन्न विधियाँ हैं, उदाहरण के लिए, सविट्ज़की-गोले फ़िल्टर। आंशिक अवकल समीकरणों को हल करने के लिए विभेदक चतुर्भुज का उपयोग किया जाता है।
शोर वाले डेटा से डेरिवेटिव की गणना के लिए और भी तरीके हैं।<ref>{{Cite journal|last1=Ahnert|first1=Karsten|last2=Abel|first2=Markus|title=Numerical differentiation of experimental data: local versus global methods|journal=Computer Physics Communications|year=2007 |volume=177|issue=10|pages=764–774|doi=10.1016/j.cpc.2007.03.009|bibcode=2007CoPhC.177..764A |issn=0010-4655}}</ref>
शोर वाले डेटा से डेरिवेटिव की गणना के लिए और भी तरीके हैं।<ref>{{Cite journal|last1=Ahnert|first1=Karsten|last2=Abel|first2=Markus|title=Numerical differentiation of experimental data: local versus global methods|journal=Computer Physics Communications|year=2007 |volume=177|issue=10|pages=764–774|doi=10.1016/j.cpc.2007.03.009|bibcode=2007CoPhC.177..764A |issn=0010-4655}}</ref>
==यह भी देखें==
==यह भी देखें==
*{{annotated link|Automatic differentiation}}
*{{annotated link|Automatic differentiation}}
Line 114: Line 104:
==संदर्भ==
==संदर्भ==
{{reflist}}
{{reflist}}


== बाहरी संबंध ==
== बाहरी संबंध ==

Revision as of 12:06, 24 July 2023

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

Derivative.svg

परिमित अंतर

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

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

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

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

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

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

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

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

अनुमान त्रुटि द्वारा दिया गया है

,

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

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


चरण आकार

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

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

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

x और x + h पर मूल्यांकित संख्यात्मक व्युत्पन्न सूत्र के लिए, h के लिए विकल्प जो बड़ी गोलाई त्रुटि उत्पन्न किए बिना छोटा है (हालांकि तब नहीं जब x = 0), जहां मशीन ईपीएसलॉन ε आम तौर पर 2.2 के क्रम का होता है×10−16 दोहरी परिशुद्धता के लिए।[8] एच के लिए सूत्र जो इष्टतम सटीकता के लिए सेकेंट त्रुटि के विरुद्ध गोलाई त्रुटि को संतुलित करता है[9]

(हालांकि कब नहीं ), और इसे नियोजित करने के लिए फ़ंक्शन के ज्ञान की आवश्यकता होगी।

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

 h := sqrt(eps) * x;
 एक्सपीएच := एक्स + एच;
 dx := xph - x;
 ढलान := (F(xph) - F(x)) / dx;

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

अन्य विधियाँ

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

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

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

कहाँ .

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

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

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

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

जटिल-परिवर्तनीय विधियाँ

संख्यात्मक विभेदन के लिए शास्त्रीय परिमित-अंतर सन्निकटन ख़राब स्थिति वाले हैं। हालांकि, यदि होलोमोर्फिक फ़ंक्शन है, जो वास्तविक रेखा पर वास्तविक-मूल्यवान है, जिसका मूल्यांकन पास के जटिल विमान में बिंदुओं पर किया जा सकता है , फिर संख्यात्मक स्थिरता विधियाँ हैं। उदाहरण के लिए,[5] पहले व्युत्पन्न की गणना जटिल-चरण व्युत्पन्न सूत्र द्वारा की जा सकती है:[12][13][14]

विभिन्न स्थितियों के लिए सटीक डेरिवेटिव प्राप्त करने के लिए अनुशंसित चरण आकार है .[6] यह सूत्र टेलर श्रृंखला विस्तार द्वारा प्राप्त किया जा सकता है:

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

जहां मल्टीकॉम्प्लेक्स काल्पनिक इकाइयों को निरूपित करें; . h> ऑपरेटर निकालता है स्तर के बहुसंकेतन संख्या का वां घटक , जैसे, वास्तविक घटक निकालता है और अंतिम, "सबसे काल्पनिक" घटक निकालता है। इस विधि को मिश्रित डेरिवेटिव पर लागू किया जा सकता है, जैसे दूसरे क्रम के व्युत्पन्न के लिए

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

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

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

संख्यात्मक विभेदन के लिए जटिल चर का उपयोग 1967 में लिनेस और मोलर द्वारा शुरू किया गया था।[20] उनका एल्गोरिदम उच्च-क्रम डेरिवेटिव पर लागू होता है।

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

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

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

यह भी देखें

संदर्भ

  1. Richard L. Burden, J. Douglas Faires (2000), Numerical Analysis, (7th Ed), Brooks/Cole. ISBN 0-534-38216-9.
  2. Katherine Klippert Merseth (2003). Windows on Teaching Math: Cases of Middle and Secondary Classrooms. Teachers College Press. p. 34. ISBN 978-0-8077-4279-2.
  3. Tamara Lefcourt Ruby; James Sellers; Lisa Korf; Jeremy Van Horn; Mike Munn (2014). Kaplan AP Calculus AB & BC 2015. Kaplan Publishing. p. 299. ISBN 978-1-61865-686-5.
  4. 4.0 4.1 Numerical Differentiation of Analytic Functions, B Fornberg – ACM Transactions on Mathematical Software (TOMS), 1981.
  5. 5.0 5.1 Using Complex Variables to Estimate Derivatives of Real Functions, W. Squire, G. Trapp – SIAM REVIEW, 1998.
  6. 6.0 6.1 Martins, Joaquim R. R. A.; Ning, Andrew (2021-10-01). इंजीनियरिंग डिज़ाइन अनुकूलन (PDF) (in English). Cambridge University Press. ISBN 978-1108833417.
  7. Sauer, Timothy (2012). Numerical Analysis. Pearson. p.248.
  8. Following Numerical Recipes in C, Chapter 5.7.
  9. p. 263.
  10. Abramowitz & Stegun, Table 25.2.
  11. Shilov, George. Elementary Real and Complex Analysis.
  12. Martins, J. R. R. A.; Sturdza, P.; Alonso, J. J. (2003). "जटिल-चरण व्युत्पन्न सन्निकटन". ACM Transactions on Mathematical Software. 29 (3): 245–262. CiteSeerX 10.1.1.141.8002. doi:10.1145/838250.838251. S2CID 7022422.
  13. Differentiation With(out) a Difference by Nicholas Higham
  14. article from MathWorks blog, posted by Cleve Moler
  15. "संग्रहीत प्रति" (PDF). Archived from the original (PDF) on 2014-01-09. Retrieved 2012-11-24.
  16. Lantoine, G.; Russell, R. P.; Dargent, Th. (2012). "उच्च-क्रम डेरिवेटिव की स्वचालित गणना के लिए मल्टीकॉम्प्लेक्स चर का उपयोग करना". ACM Trans. Math. Softw. 38 (3): 1–21. doi:10.1145/2168773.2168774. S2CID 16253562.
  17. Verheyleweghen, A. (2014). "मल्टी-कॉम्प्लेक्स चरण विधि का उपयोग करके उच्च-क्रम डेरिवेटिव की गणना" (PDF).
  18. Bell, I. H. (2019). "एमसीएक्स (मल्टीकॉम्प्लेक्स बीजगणित लाइब्रेरी)". GitHub.
  19. Ablowitz, M. J., Fokas, A. S.,(2003). Complex variables: introduction and applications. Cambridge University Press. Check theorem 2.6.2
  20. Lyness, J. N.; Moler, C. B. (1967). "विश्लेषणात्मक कार्यों का संख्यात्मक विभेदन". SIAM J. Numer. Anal. 4 (2): 202–210. Bibcode:1967SJNA....4..202L. doi:10.1137/0704019.
  21. Abate, J; Dubner, H (March 1968). "कार्यों की शक्ति श्रृंखला विस्तार उत्पन्न करने की एक नई विधि". SIAM J. Numer. Anal. 5 (1): 102–112. Bibcode:1968SJNA....5..102A. doi:10.1137/0705008.
  22. Differential Quadrature and Its Application in Engineering: Engineering Applications, Chang Shu, Springer, 2000, ISBN 978-1-85233-209-9.
  23. Advanced Differential Quadrature Methods, Yingyan Zhang, CRC Press, 2009, ISBN 978-1-4200-8248-7.
  24. Ahnert, Karsten; Abel, Markus (2007). "Numerical differentiation of experimental data: local versus global methods". Computer Physics Communications. 177 (10): 764–774. Bibcode:2007CoPhC.177..764A. doi:10.1016/j.cpc.2007.03.009. ISSN 0010-4655.

बाहरी संबंध