Finite Differences
Finite Differences
SANDLI SRIVASTAVA
Department of Mathematics. Class - B.sc III
Roll No –
1
CERTIFICATE FROM THE GUIDE
This is to certify that the work embodied in this project
report entitled “FINITE DIFFERENCES” submitted to
CSJM University, Kanpur, for a award of the degree of
B.Sc has been carried out by SANDLI SRIVASTAVA
under my supervision at Department of Mathematics ,
PPN COLLEGE.
Date -……………………..
2
ACKNOWLEDGMENT
The success and final outcome of this project report required a lot
of guidance and assistance from many people and I am extremely
fortunate to have got this all along the completion of my project
work. I would like to extend my sincere thanks to all of them.
SUBMITTED BY-
SANDLI SRIVASTAVA
BSc III year
3
INDEX
1. Introduction 5
2. Shift Operator 7
3. Forward Differences 9
4. Backward Differences 13
5. Average Operator 15
9. Factorial Polynomials 19
4
INTRODUCTION
Numerical Analysis is a branch of Mathematics which leads to approximate
solution by repeated applications of four basic operations of Algebra. The
knowledge of finite differences is essential for the study of Numerical
Analysis.
PRELIMINARIES –
• PARTITION: - Let [a, b] be a closed interval. Then a set points:
P = {a=x0, x1, x2, ……...,xn-1 ,xn}
a=x0<x1< x2< ……...<xn-1 <xn = b
is called of partition of [a, b]
• Grid point: - The elements of a partition are called its grid points.
• Step size: - If the grid points of the partition are equally space (h) then h is
called the step size of the partition. Xi – Xi-1 = h, where i = 1, 2, 3...n. The last
one of these relations can be also defined as
h = (b-a)/n .
5
FINITE DIFFERENCES
Consider the arguments x0, x1, x2…xn and the entries y0, y1, y2…yn. Let y=f(x) be
a function of x. Let us assume that the values of x are in increasing order and
equally spaced with a space length h. Then the values of x may be taken to be
x0, x0+h, x0+2h....x0+nh and the function assumes the values f (x0 ),
f (x0 + h ), f (x 0 + 2h),…f (x0 + nh) . Here we study some of the finite
differences of the function y = f (x)
Arguments x₀ x₁ x₂ …. …. xₑ
x =x₀+h =x₀+2h = x₀+eh
Entry y₀ y₁ y₂ …. …. yₑ
y = f(x) = f(x) =f(x₀+h) =f(x₀+2h) =f(x₀+eh)
Suppose that we are required to recover / estimate the value of f(x) or f'(x) for
some intermediate value of x in the interval x ≤ x ≤ x+nh. The methods for
solution of these and many more related problems are based on the concept of
finite differences of the function.
The following three such types of differences are commonly used in numerical
analysis:
1. Forward differences.
2. Backward differences.
3. Central differences.
6
SHIFT OPERATOR
Let y = f (x) be a given function of x and x0 , x0+h , x0+2h , x0+3h ,….,
x0+nh be the consecutive values of x. Then the operator E is defined as-
E[f(x0)] = f(x0+h)
Enf(x) =f(x+nh)
7
Properties of the shift operator -
1. Linearity property: The shift operator E is linear.
E[C1f1(x) +C2f2(x)] = C1 E[f1(x)] + C2 E[f2(x)]
2. Law of indices: If m and n are positive integers, then
EmEn = Em+n
3. Inverse shift operator(E-1): It is defined as
E-1 f(x) = f(x-h)
(E-1)-1 = E.
Similarly, the general formula for the inverse shift operator will
be -
E-nf(x) = f(x-nh)
8
FORWORD DIFFERENCE OPERATOR
(∆)
Let y = f(x) be a given function of x. Let y0, y1 , y2 ,….. yn be the value
of y at = x0 , x1 , x2 ,…. xn respectively.Then y1 − y0, y2 − y1, y 3 − y 2, …. yn − yn−1 are
called the first (forward) differences of the function y. They are denoted
by Δy0 , Δy1 , Δy2 ,..., Δyn−1 respectively.
(i.e) Δy 0 = y1 − y0 , Δy1 = y2 − y1, Δy2 = y3 − y2 ,..., Δyn −1 = yn − yn−1
In general, Δyn = yn+1 − yn , n = 0,1,2,3,...
The symbol Δ is called the forward difference operator and pronounced as delta.
The forward difference operator ∆ can also be defined
as Df ( x) = f ( x + h ) − f ( x), h is the equal interval of spacing.
9
= Δ f (x) + Δ g (x)
Similarly we can show
that Δ [f ( x) − g ( x)] = f (x) − Δ [g(x)]general, Δ[ f1 (x) + f 2 (x)...... + f n (x)] = Δf
1 (x)+ Δ f 2 (x) + .... + Δ fn (x).
The differences of the first differences denoted by Δ2y0, Δ2y1, ….. Δ2yn, are
10
In general kth differences of yn is
11
The forward difference table for f(x) is given below.
12
BACKWORD DIFFERENCE OPERATOR (∇)
13
Second (backward) differences:
∇ 2 y n = ∇ y n − ∇yn+1 , n = 1,2,3,…
Third (backward) differences:
∇ 3 y n = ∇ 2 yn − ∇2 yn−1 n = 1,2,3,…
In general, kth (backward) differences:
14
First differences: ∇ f (x + h) = f (x + h) − f (x)
∇ f (x + 2h) = f (x + 2h) − f (x + h),..., h is the interval of spacing.
Second differences:
∇ 2 f (x + h) = ∇(∇f (x + h) = ∇( f (x + h) − f (x)) = ∇ f (x + h) − ∇f (x)
∇ 2 f (x + 2h) = ∇ f (x + 2h ) − ∇f (x + h)
Third differences:
∇ 3 f (x + h) = ∇ 2 f (x + h ) − ∇2 f (x)
∇ 3 f (x + 2h) = ∇ 2 f (x + 2h ) − ∇2 f (x + h)
Here we note that,
∇ f (x + h) = f (x + h ) − f (x) = Δf (x)
∇ f (x + 2h) = f (x + 2h ) − f (x + h) = Δf (x + h)
∇ 2 f (x + 2h) = ∇ f (x + 2h ) − ∇f (x + h) = Δf (x + h ) − Δf (x)
= Δ2 f (x)
1 1
𝜇𝑓(𝑥 ) = 1/2 {𝑓 (𝑥 + ℎ) + 𝑓(𝑥 − ℎ)}
2 2
1 1
𝛿𝑓(𝑥 ) = 𝑓 (𝑥 + ℎ) − 𝑓(𝑥 − ℎ)
2 2
16
Δ f (x) = f (x + h ) − f (x)
Δ f (x) = Ef (x ) − f (x)
⇒ Δ f (x) = (E −1) f (x)
Δ = E −1
∴ E=1+Δ
∆=E-1
∇ f (x) = f(x) - f (x - h )
Δ f (x) = f (x) – E -1f (x )
⇒ Δ f (x) = (E-1 −1) f (x)
∇ =∇ 1and
(III) Relation between E, – E∆-1
Proof:
(E∇) f(x) = E [∇f(x) ]= E [ f(x)- f(x-h) ]
= f(x+h ) – f(x)
= Δf(x)
(∇E ) f(x) = ∇[E f(x)]
= f(x+h) – f(x) = Δf(x)
17
We conclude that
E∇ = Δ = ∇ E
∇ = E −1 / E
∇f (x) = f (x) - f (x-h)
= f (x) – E-1f (x)
= (1- E-1) f (x)
∇ = (1- E-1)
∇ = 1 – 1/E
Hence ∇ = [E – 1]/E
4. Δ∇ = ∇Δ= Δ - ∇
5. (1+ Δ) (1-∇) = 1
Using the difference operators and shifting operator we can able to find the
missing terms.
Example :
From the following table find the missing value
18
Solution:
Since only four values of f(x) are given, the polynomial which fits the data is
of degree three. Hence fourth differences are zeros.
19
FACTORIAL POLYNOMIALS
X(n) = x(x - h)(x – 2h)…(x – n – 1h), x ≠ 0
Factorial Polynomials of n is a positive integer and h > 0, the function x (n)
defined by the product is called a factorial polynomial (or factorial function) of
degree n.
By convention, we define x (0) =1. Clearly
21
DIFFERENCES OF ZERO
Let y0 be the true values of a function y f(x) at the equally spaced points x0, xt,
.t9. Suppose the value be affected with an error e such that the erroneous value
of is e. Then the effect of this error in the successive differences is as shown in
the following table :
23
(i) The error increases with the order of the differences.
(ii) The error is maximum (in magnitude) along the horizontal line
through the erroneous tabulated value.
(iii) In the kth difference column, the coefficients of errors are the
binomial coefficients
in the expansion of (1 − x)k... In particular, the errors in the second
difference
column are ε, −2ε, ε, in the third difference column these are ε, −3ε, 3ε,
−ε, and
so on.
(iv) The algebraic sum of errors in any complete column is zero.
If there is any error in a single entry of the table, then we can detect and
correct
it from the difference table. The position of the error in an entry can be
identified by
performing the following steps.
(i) If at any stage, the differences do not follow a smooth pattern,
then there is an error
(ii) If the differences of some order (it is generally happens in
higher order) becomes alternating in sign then the middle entry
contains an error.
24
REFERENCE
The below followings are used to taken help and
references for the following Project Report.
WEBSITES
www.google.com
www.youtube.com
www.Wikipedia.org
BOOKS
25