Differentiation, numerical

Finding the derivative of a function by numerical methods. Such differentiation is resorted to when the methods of differential calculus are inapplicable (the function is obtained from tables) or involves considerable difficulties (the analytic expression of the function is complicated).

Let a function be defined on an interval and let the nodal points , , be given. The totality of points , , is known as a table. The result of numerical differentiation of the table is the function which approximates, in some sense, the -th derivative of the function on some sets of points . The use of numerical differentiation will be expedient if only an insignificant amount of computational effort is required to obtain the function for each . Linear methods of numerical differentiation are commonly employed; the result thus obtained is written in the form

(1)

where are functions defined on . The most popular method for obtaining formulas (1) is as follows: One constructs the function

interpolating , and assumes that

The accuracy of the algorithms based on the interpolation formulas of Lagrange, Newton and others strongly depends on the selection of the manner of interpolation, and may sometimes be very low even if the function is sufficiently smooth and if the number of nodal points is large [1]. Algorithms of numerical differentiation involving spline-interpolation [2] are often free from this disadvantage. If all that is needed is the computation of the approximate values of the derivative at the nodal points only, formula (1) assumes the form

(2)

and is fully defined by specifying a coefficient matrix for a given . Formulas such as (2) are known as difference formulas for numerical differentiation. The coefficients of such formulas are determined from the condition that the difference

has highest order of smallness with respect to . As a rule, formulas (2) are very simple and easy to handle. Thus, if , they assume the form

Numerical differentiation algorithms are often employed with tables in which the values of are given (or obtained) inaccurately. In such cases there is need for a preliminary smoothing, since a direct application of the formula may result in large errors in the results [3].

References

[1]	I.S. Berezin, N.P. Zhidkov, "Computing methods" , Pergamon (1973) (Translated from Russian)
[2]	J.H. Ahlberg, E.N. Nilson, J.F. Walsh, "Theory of splines and their applications" , Acad. Press (1967)
[3]	V.A. Morozov, "The differentiation problem and certain algorithms of approximation by experimental information" , Computing methods and programming , 14 , Moscow (1970) pp. 46–62 (In Russian)

Comments

Explicit differentiation formulas are given in, e.g., [a1], [a2].

References

[a1]	A. Segun, M. Abramowitz, "Handbook of mathematical functions" , Appl. Math. Ser. , 55 , Nat. Bur. Standards (1970)
[a2]	F.B. Hildebrand, "Introduction to numerical analysis" , McGraw-Hill (1974)