Scribd
Upload
25 views3 pages

Lagrange Interpolation

This document describes Lagrange interpolation, which finds a polynomial of degree n-1 that passes through n data points. It presents the Lagrange interpolation formula, which writes the interpolating polynomial as a linear combination of n basis polynomials, each of which is 0 at n-1 of the points. The formula guarantees a unique interpolating polynomial. An example fits a quadratic polynomial to 3 data points and evaluates it.

Uploaded by

Macloud Kamula
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
25 views3 pages

Lagrange Interpolation

This document describes Lagrange interpolation, which finds a polynomial of degree n-1 that passes through n data points. It presents the Lagrange interpolation formula, which writes the interpolating polynomial as a linear combination of n basis polynomials, each of which is 0 at n-1 of the points. The formula guarantees a unique interpolating polynomial. An example fits a quadratic polynomial to 3 data points and evaluates it.

Uploaded by

Macloud Kamula
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

LAGRANGE INTERPOLATION

DARYL DEFORD

1. Introduction
Polynomial interpolation is a method for solving the following problem:
Given a set of n of data points with distinct x–coordinates {(xi , yi )}ni=1 find a poly-
nomial of degree at most n − 1 that passes through each point.
Example graphs of these polynomials for different data sets are shown below:

(a) (b) (c)

(d) (e) (f)

These interpolating polynomials provide a smooth approximation to the data points allowing
for efficient extrapolation of the data. These methods are often used for constructing numerical
approximations for computing complicated function values and evaluating differential equations.
Methods for constructing interpolating polynomials go back hundreds of years. Any polynomial
interpolation problem can be solved by realizing the coefficients as the solution to a linear system
of the form:

1 x1 x21 · · · xn−1 c0 y1
    
1
1 x2 x22 · · · xn−1
2  c1   y2 
1 x3 x23 · · · x3n−1   c2  =  y3 
    
. .
 .. .. .. .. ..   ..  .
  .. 
. . .  .
2
1 xn xn · · · xnn−1 cn−1 yn

Date: September 25, 2015.


1
The theory of Vandermonde determinants shows that this system has a unique solutions and
hence there exists a unique interpolating polynomial: cn−1 xn−1 + cn−2 xn−2 + · · · + c1 x + c0 . However,
solving this system requires inverting a n × n matrix which can be computationally infeasible and
can introduce distortions due to poor matrix conditioning. The method presented in class, which
is described the section below, is due to a combination of Euler, Waring, and Lagrange, and was
published in the late 1700’s.
2. Formula
The Lagrange interpolation formula writes the interpolating polynomial for {(xi , yi )}ni=1 as a
linear combination of n degree n − 1 polynomials each of which is zero at exactly n − 1 of the points.
The polynomials are given by the general formula:
Y x − xj
`i =
x − xj
0≤j≤n: i
j6=i

and the interpolation polynomial is calculated as:


n
X
L(x) = yi `i (x)
i=1
In Math 1 we looked specifically at the version where n = 3 so we obtain the following quadratic
polynomials for {(x1 , y1 ), (x2 , y2 ), (x3 , y3 )}:
  
x − x2 x − x3
`1 =
x1 − x2 x1 − x 3
  
x − x1 x − x3
`2 =
x2 − x1 x2 − x 3
  
x − x1 x − x2
`3 =
x3 − x1 x3 − x 2

L(x) = y1 `1 + y2 `2 + y3 `3
        
x − x2 x − x3 x − x1 x − x3 x − x1 x − x2
= y1 + y2 + y3
x1 − x2 x1 − x3 x2 − x1 x2 − x3 x3 − x1 x3 − x2

2
3. Example
Fit a quadratic polynomial to the following data points {(−2, 9), (5, −12), (10, 33)}.
  
x−5 x − 10
`1 =
−2 − 5 −2 − 10
  
x−5 x − 10
−7 −12
 2 
x − 15x + 50
84
  
x+2 x − 10
`2 =
5+2 5 − 10
  
x+2 x − 10
7 −5
 2 
x − 8x − 20
−35
  
x+2 x−5
`3 =
10 + 2 10 − 5
  
x+2 x−5
12 5
 2 
x − 3x − 10
60
 2   2   2 
x − 15x + 50 x − 8x − 20 x − 3x − 10
L(x) = y1 + y2 + y3
84 −35 60
 2   2   2 
x − 15x + 50 x − 8x − 20 x − 3x − 10
= 9 + −12 + 33
84 −35 60

= x2 − 6x − 7

(a) (b)

Department of Mathematics, Dartmouth College


E-mail address: [email protected]

You might also like