Numerical Interpolation
Numerical Interpolation
1
Direct Method of
Interpolation
2
What is Interpolation ?
Given (x0,y0), (x1,y1), …… (xn,yn), find the value of ‘y’ at a
value of ‘x’ that is not given.
4
Direct Method
Given ‘n+1’ data points (x0,y0), (x1,y1),………….. (xn,yn),
pass a polynomial of order ‘n’ through the data as given
below:
y a0 a1 x .................... an x . n
5
Example 1
The upward velocity of a rocket is given as a
function of time in Table 1.
Find the velocity at t=16 seconds using the
direct method for linear interpolation.
Table 1 Velocity as a function
of time.
t , s vt , m/s
0 0
10 227.04
15 362.78
20 517.35
22.5 602.97
30 901.67 Figure 2 Velocity vs. time data for the
rocket example
6
Linear Interpolation
vt a0 a1t y
Hence
vt 100.93 30.914t , 15 t 20.
v16 100.93 30.91416 393.7 m/s
7
Example 2
The upward velocity of a rocket is given as a
function of time in Table 2.
Find the velocity at t=16 seconds using the
direct method for quadratic interpolation.
Table 2 Velocity as a function
of time.
t , s vt , m/s
0 0
10 227.04
15 362.78
20 517.35
22.5 602.97
30 901.67 Figure 5 Velocity vs. time data for the
rocket example
8
Quadratic Interpolation
y
f 2 x
v20 a0 a1 20 a2 20 517.35
2
x0 , y 0
x
ys
f x desired 350
250
227.04 200
10 12 14 16 18 20
10 x s range x desired 20
t , s vt , m/s
0 0
10 227.04
15 362.78
20 517.35
22.5 602.97
30 901.67 Figure 6 Velocity vs. time data for the
rocket example
11
Cubic Interpolation
y
x3 , y3
vt a0 a1t a2t a3t
2 3
x1 , y1
v10 227.04 a0 a1 10 a2 10 a3 10
2 3
f 3 x
x2 , y 2
v15 362.78 a0 a1 15 a2 15 a3 15
2 3
x0 , y0
12
Cubic Interpolation (contd)
vt 4.2540 21.266t 0.13204t 2 0.0054347t 3 , 10 t 22.5
v 16 4 .2540 21 .266 16 0 .13204 16 0.0054347 16
2 3
f x desired
a 100
300
392 .06
227.04 200
10
10
12 14 16 18
x s range x desired
20 22 24
22.5
0.033269 %
13
Comparison Table
Table 4 Comparison of different orders of the polynomial.
Order of
1 2 3
Polynomial
vt 16 m/s 393.7 392.19 392.06
Absolute Relative
---------- 0.38410 % 0.033269 %
Approximate Error
14
Distance from Velocity Profile
Find the distance covered by the rocket from t=11s to t=16s ?
vt 4.3810 21.289t 0.13064t 2 0.0054606t 3 , 10 t 22.5
16
s 16 s 11 vt dt
11
16
4.2540 21.266t 0.13204t 2 0.0054347 t 3 dt
11
16
t2 t3 t4
4.2540t 21.266 0.13204 0.0054347
2 3 4 11
1605 m
15
Acceleration from Velocity Profile
Find the acceleration of the rocket at t=16s given that
t 4.2540 21.266t 0.132042 0.0054347t 3 ,10 t 22.5
d
a t vt
dt
4.2540 21.266t 0.13204t 2 0.0054347t 3
d
dt
21.289 0.26130t 0.016382t 2 , 10 t 22.5
29.665 m/s 2
16
Lagrange Method of
Interpolation
Lagrangian Interpolation
Lagrangian interpolating polynomial is given by
n
f n ( x) Li ( x ) f ( xi )
i 0
where ‘ n ’ in f n (x) stands for the n th order polynomial that approximates the function y f (x)
18
Example
The upward velocity of a rocket is given as a function of
time in Table 1. Find the velocity at t=16 seconds using
the Lagrangian method for linear interpolation.
Table Velocity as a
function of time
t (s) v(t ) (m/s)
0 0
10 227.04
15 362.78
20 517.35
22.5 602.97
30 901.67
1
v(t ) Li (t )v(ti ) 500
i 0
ys
L0 (t )v(t 0 ) L1 (t )v(t1 ) f ( range)
450
f x desired
20
Linear Interpolation (contd)
1 t tj t t1
L0 (t )
j 0 t0 t j t 0 t1
j 0
1 t tj t t0
L1 (t )
j 0 t1 t j t1 t 0
j 1
t t1 t t0 t 20 t 15
v(t ) v(t 0 ) v(t1 ) (362.78) (517.35)
t 0 t1 t1 t 0 15 20 20 15
16 20 16 15
v(16) (362.78) (517.35)
15 20 20 15
0.8(362.78) 0.2(517.35)
393.7 m/s.
21
Quadratic Interpolation
For the second order polynomial interpolation (also called quadratic interpolation), we
choose the velocity given by
2
v (t ) Li ( t ) v(t i )
i 0
L0 (t )v (t 0 ) L1 (t ) v( t1 ) L2 (t ) v( t 2 )
22
Example
The upward velocity of a rocket is given as a function of
time in Table 1. Find the velocity at t=16 seconds using
the Lagrangian method for quadratic interpolation.
Table Velocity as a
function of time
t (s) v(t ) (m/s)
0 0
10 227.04
15 362.78
20 517.35
22.5 602.97
30 901.67
ys
400
f ( range)
2 t tj t t1 t t 2
L0 (t )
f x desired 350
j 0 t0 t j t t
0 1 0 2 t t
j 0
300
2 t t j t t0 t t 2
L1 (t )
j 0 t1 t j t1 t0 t1 t 2 250
j 1
2 t tj t t0 t t1 227.04 200
L2 (t ) 10 12 14 16 18 20
j 0 t2 t j t t
2 0 2 1 t t 10 x s range x desired 20
j 2
24
Quadratic Interpolation (contd)
t t1 t t2 t t0 t t2 t t0 t t1
vt vt0 vt1 vt2
t0 t1 t0 t2 t1 t0 t1 t2 t2 t0 t2 t1
1615 16 20 1610 16 20 1610 1615
v16 227 .04 362 .78 517.35
1015 10 20 1510 15 20 20 10 2015
0.08227.04 0.96362.78 0.12527.35
392.19 m/s
The absolute relative approximate error a obtained between the
results from the first and second order polynomial is
392.19 393.70
a 100
392.19
0.38410%
25
Cubic Interpolation
For the third order polynomial (also called cubic interpolation), we choose the velocity given by
3
v (t ) Li ( t ) v(t i )
i 0
600
ys 500
f ( range)
f x desired
400
300
227.04 200
10 12 14 16 18 20 22 24
10 x s range x desired 22.5
26
Example
The upward velocity of a rocket is given as a function of
time in Table 1. Find the velocity at t=16 seconds using
the Lagrangian method for cubic interpolation.
Table Velocity as a
function of time
t (s) v(t ) (m/s)
0 0
10 227.04
15 362.78
20 517.35
22.5 602.97
30 901.67
700
3 t tj t t 1 t t 2 t t 3 602.97
L0 (t ) ;
j 0 t0 t j t t t t
0 1 0 2 0 3 t t 600
j 0
3 t t j t t0 t t 2 t t 3
L1 (t )
ys
500
t1 t j t1 t 0 1 2 t1 t 3
t t
f ( range)
j 0
j 1
f x desired
400
3 t tj t t 0 t t1 t t 3
L2 (t ) ;
t2 t j t 2 t 0 t 2 t 1 t 2 t 3
300
j 0
j 2
t tj t t 0 t t1 t t 2
227.04 200
3
L3 ( t )
10 12 14 16 18 20 22 24
10 x s range x desired 22.5
j 0 t3 t j t3 t 0 t 3 t1 t3 t 2
j 3
28
Cubic Interpolation (contd)
t t1 t t 2 t t3 t t0 t t2 t t3
vt vt1 vt 2
t t t
0 1 0 2 0 3 t t t t t t
1 0 1 2 1 3 t t t
t t0 t t1 t t3 t t1 t t1 t t 2
vt 2 vt3
t t t
2 0 2 1 2 3 t t t t t t
3 1 3 1 3 2 t t t
16 15 16 20 16 22.5 16 10 16 20 16 22.5
v16 227.04 362.78
10 15 10 20 10 22.5 15 10 15 20 15 22.5
16 10 16 15 16 22.5 16 10 16 15 16 20
517.35 602.97
20 10 20 15 20 22 .5 22 .5 10 22 .5 15 22. 5 20
0.0416227.04 0.832362.78 0.312517.35 0.1024602.97
392.06 m/s
The absolute relative approximate error a obtained between the
results from the first and second order polynomial is
392.06 392.19
a 100
392.06
0.033269%
29
Comparison Table
Order of
1 2 3
Polynomial
v(t=16) m/s 393.69 392.19 392.06
Absolute Relative
-------- 0.38410% 0.033269%
Approximate Error
30
Distance from Velocity Profile
Find the distance covered by the rocket from t=11s to
t=16s ?
v(t ) (t 3 57.5t 2 1087.5t 6750)(0.36326) (t 3 52.5t 2 875t 4500)(1.9348)
(t 3 47.5t 2 712.5t 3375)(4.1388) (t 3 45t 2 650t 3000)(2.5727)
v (t ) 4.245 21.265t 0.13195t 2 0.00544t 3 , 10 t 22.5
16
s(16) s (11) v( t ) dt
11
16
( 4.245 21.265t 0.13195t 2 0.00544t 3 ) dt
11
t2 t3 t 4 16
[ 4.245t 21.265 0.13195 0.00544 ]11
2 3 4
31
1605 m
Acceleration from Velocity Profile
Find the acceleration of the rocket at t=16s given that
dt dt
29.665 m / s 2
32
Newton’s Divided
Difference Method of
Interpolation
Newton’s Divided Difference
Method
Linear interpolation: Given ( x0 , y0 ), ( x1 , y1 ), pass a
linear interpolant through the data
f1 ( x) b0 b1 ( x x0 )
where
b0 f ( x 0 )
f ( x1 ) f ( x 0 )
b1
x1 x 0
34
Example
The upward velocity of a rocket is given as a function of
time in Table 1. Find the velocity at t=16 seconds using
the Newton Divided Difference method for linear
interpolation.
Table. Velocity as a
function of time
v(t ) b0 b1 (t t 0 ) 500
ys
v(t1 ) v(t 0 )
b1 30.914 362.78 350
t1 t 0 10
x s 10
12 14 16 18
x s range x desired
20 22 24
x s 10
0 1
36
Linear Interpolation (contd)
550
517.35
500
ys
f ( range)
450
f x desired
400
362.78 350
10 12 14 16 18 20 22 24
x s 10 x s range x desired x s 10
v (t ) b0 b1 (t t 0 )
0 1
b0 f ( x0 )
f ( x1 ) f ( x0 )
b1
x1 x0
f ( x 2 ) f ( x1 ) f ( x1 ) f ( x0 )
x 2 x1 x1 x0
b2
x2 x0
38
Example
The upward velocity of a rocket is given as a function of
time in Table 1. Find the velocity at t=16 seconds using
the Newton Divided Difference method for quadratic
interpolation.
Table. Velocity as a
function of time
500
450
ys
400
f ( range)
f x desired 350
300
250
227.04 200
10 12 14 16 18 20
10 x s range x desired 20
The absolute relative approximate error a obtained between the results from the first
order and second order polynomial is
392.19 393.69
a x100
392.19
= 0.38502 %
42
General Form
f 2 ( x) b0 b1 ( x x0 ) b2 ( x x0 )( x x1 )
where
b0 f [ x0 ] f ( x0 )
f ( x1 ) f ( x 0 )
b1 f [ x1 , x0 ]
x1 x0
f ( x 2 ) f ( x1 ) f ( x1 ) f ( x0 )
f [ x 2 , x1 ] f [ x1 , x 0 ] x 2 x1 x1 x0
b2 f [ x 2 , x1 , x 0 ]
x 2 x0 x 2 x0
Rewriting
f 2 ( x) f [ x0 ] f [ x1 , x0 ]( x x0 ) f [ x 2 , x1 , x0 ]( x x0 )( x x1 )
43
General Form
Given (n 1) data points, x0 , y 0 , x1 , y1 ,......, x n 1 , y n 1 , x n , y n as
f n ( x) b0 b1 ( x x0 ) .... bn ( x x0 )( x x1 )...( x x n 1 )
where
b0 f [[xx0 ]
b1 f [ x1 , x0 ]
b2 f [ x 2 , x1 , x0 ]
bn 1 f [ x n 1 , x n 2 ,...., x0 ]
bn f [ x n , x n 1 ,...., x0 ]
44
General form
The third order polynomial, given ( x0 , y 0 ), ( x1 , y1 ), ( x 2 , y 2 ), and ( x3 , y 3 ), is
f 3 ( x) f [ x0 ] f [ x1 , x0 ]( x x0 ) f [ x 2 , x1 , x0 ]( x x0 )( x x1 )
f [ x3 , x 2 , x1 , x0 ]( x x0 )( x x1 )( x x 2 )
b0
x0 f ( x0 ) b1
f [ x1 , x 0 ] b2
x1 f ( x1 ) f [ x 2 , x1 , x0 ] b3
f [ x 2 , x1 ] f [ x3 , x 2 , x1 , x 0 ]
x2 f ( x2 ) f [ x3 , x 2 , x1 ]
f [ x3 , x 2 ]
x3 f ( x3 )
45
Example
The upward velocity of a rocket is given as a function of
time in Table 1. Find the velocity at t=16 seconds using
the Newton Divided Difference method for cubic
interpolation.
Table. Velocity as a
function of time
47
Example
b0
t0 10 227.04 b1
27.148 b2
t1 15, 362.78 0.37660 b3
30.914 5.4347 10 3
t 2 20, 517.35 0.44453
34.248
t3 22.5, 602.97
48
Example
Hence
v (t ) b0 b1 (t t 0 ) b2 (t t 0 )( t t1 ) b3 (t t 0 )( t t1 )(t t 2 )
227.04 27.148( t 10) 0.37660(t 10)(t 15)
5.4347 * 10 3 (t 10)( t 15)( t 20)
At t 16,
v (16) 227.04 27.148(16 10) 0.37660(16 10)(16 15)
5.4347 * 10 3 (16 10)(16 15)(16 20)
392.06 m/s
The absolute relative approximate error a obtained is
392.06 392.19
a x100
392.06
= 0.033427 %
49
Comparison Table
Order of 1 2 3
Polynomial
v(t=16) 393.69 392.19 392.06
m/s
Absolute Relative ---------- 0.38502 % 0.033427 %
Approximate Error
50
Distance from Velocity Profile
Find the distance covered by the rocket from t=11s to
t=16s ?
v (t ) 227.04 27.148(t 10) 0.37660( t 10)( t 15)
10 t 22.5
5.4347 * 10 (t 10)( t 15)( t 20)
3
16
16
t2 t3 t4
4.2541t 21.265 0.13204 0.0054347
2 3 4 11
51 1605 m
Acceleration from Velocity Profile
Find the acceleration of the rocket at t=16s given that
29.664 m / s 2
52