Finite Differences
Finite Differences
P. Sam Johnson
January 9, 2015
Finite Differences
January 9, 2015
1 / 24
Overview
Assume that we have a table of values
(xi , yi ),
i = 0, 1, 2, . . . , n
Finite Differences
January 9, 2015
2 / 24
i = 0, 1, 2, . . . , n
i = 0, 1, 2, . . . , n.
Finite Differences
January 9, 2015
3 / 24
x = x0 , x1 , . . . , xn1 .
Finite Differences
January 9, 2015
4 / 24
Exercises
1 Find the successive differences of f (x) = ab c x and sum the first n
differences.
2 Find the function (with suitable h) whose first difference is
ax + b
sin x
ex
e a+bx
Finite Differences
January 9, 2015
5 / 24
The following result can be proved by the method of induction and it helps
in expressing, any powers of operated on y (x), in terms of yi s.
Theorem
Suppose that we have a table of values (xi , yi ), i = 0, 1, 2, . . . , n of any
function y = f (x), the values of x being equally spaced with h as the
interval of differencing. Then for any non-negative integer k,
k
X
ki k
yr =
(1)
yr +i .
i
k
i=0
Finite Differences
January 9, 2015
6 / 24
Shift Operator, E
An operator E defined as
Ef (x) = f (x + h)
is called a shift operator which results in increasing the argument by the
interval of differencing.
E n stands for the operation E being carried n times. Then
E n f (x) = f (x + nh).
The inverse operator E 1 is defined as
E 1 f (x) = f (x h).
Finite Differences
January 9, 2015
7 / 24
Finite Differences
January 9, 2015
8 / 24
With the help of the operators E and , we can prove the above result
which is restated below.
For any non-negative integer k,
k
X
ki k
yr =
(1)
yr +i .
i
k
i=0
Finite Differences
January 9, 2015
9 / 24
x
x0
y = f (x)
y0
x0 + h
y1
f (x)
2 f (x)
3 f (x)
4 f (x)
y0
2 y0
3 y0
y1
x1 + 2h
2 y1
y2
y2
x1 + 3h
y3
x1 + 4h
y4
4 y0
3 y1
2 y2
y3
Finite Differences
January 9, 2015
10 / 24
Exercise
An operator T is said to be linear if
T [af (x) + bg (x)] = aT [f (x)] + bT [g (x)],
where a, b are constants. Prove that the operators E and are linear.
The following theorem says that any value of the function f (x) can be
expressed in terms of leading term and the leading differences of an
ordinary difference table. It can be proved by the method of induction.
Theorem
For all integral values of n, f (a + nh) =
Pn
n
r =0 r
r f (a).
Exercise
Prove the above result with the help of the operators E and .
Finite Differences
January 9, 2015
11 / 24
Finite Differences
January 9, 2015
12 / 24
Exercises
1 Prove that
e =
2
E
ex .
Ee x
,
2 e x
Prove that
n1
X
yx = yn y0 .
x=0
3
Prove that y4 = y3 + y2 + 2 y1 + 3 y1 .
Show that
y4 = y0 + 4y0 + 62 y1 + 103 y1
if fourth and higher differences are zero.
Finite Differences
January 9, 2015
13 / 24
Exercises
1 Prove that
n+1
n+1
n+1
y0 + y1 + + yn =
y0 +
y0 + +
n y0 .
1
2
n+1
2
X
x=0
y2x
1X
1
=
yx +
2
4
x=0
4
1 +
2
4
y0 .
Finite Differences
January 9, 2015
14 / 24
Finite Differences
January 9, 2015
15 / 24
Averaging Operator,
We can observe the following
y0 = y1 = y1/2 ,
3 y2 = 3 y5 = 3 y7/2 .
Finite Differences
January 9, 2015
16 / 24
Differential Operator, D
d
(f (x)).
dx
Exercise
Deduce the following:
1
E e hD
hD
2
2
3
3
Finite Differences
January 9, 2015
17 / 24
Factorial Notation
[n]
x!
x
=
= n!
(x > n).
(x n)!
n
Finite Differences
January 9, 2015
18 / 24
Finite Differences
January 9, 2015
19 / 24
We now discuss the factorial function x [n] when n is any negative integer.
A simple calculation shows that for any positive integer n
x [n] = (x [n 1]h)x [n1] .
(1)
1
x+h .
In general,
x [n] =
1
(x + nh)[n]
Finite Differences
January 9, 2015
20 / 24
Exercises
1 Express the function
x 4 12x 3 + 24x 2 30x + 9
in factorial notation, the interval of differencing being unity.
2
Finite Differences
January 9, 2015
21 / 24
Exercises
1 Find the relation between , and in order that
+ x + x 2
may be expressible in one term in the factorial notation.
Finite Differences
January 9, 2015
22 / 24
Exercises
1 Using the method of separation of symbols, prove the following
identities:
1
2
3
1
25(c b) + 3(a c)
c+
2
256
where a = y0 + y5 , b = y1 + y4 , c = y2 + y5 .
2
x
r nn = nr
, r <n
n xn = 1
x [r ] (x rh)[n] = x [r +n] .
Finite Differences
January 9, 2015
23 / 24
References
Richard L. Burden and J. Douglas Faires, Numerical Analysis Theory and Applications, Cengage Learning, Singapore.
David Kincaid and Ward Cheney, Numerical Analysis Mathematics of Scientific Computing, American Mathematical
Society, Providence, Rhode Island.
Finite Differences
January 9, 2015
24 / 24