example of under-determined polynomial interpolation


Consider the following interpolationMathworldPlanetmath problem:

Given x1,y1,x2,y2∈R with x1≠x2 to determine all cubic polynomials

p⁢(x)=a⁢x3+b⁢x2+c⁢x+d,x,a,b,c,d∈ℝ

such that

p⁢(x1)=y1,p⁢(x2)=y2.

This is a linear problem. Let 𝒫3 denote the vector spaceMathworldPlanetmath of cubic polynomials. The underlying linear mapping is the multi-evaluation mapping

E:𝒫3→ℝ2,

given by

p↦(p⁢(x1)p⁢(x2)),p∈𝒫3

The interpolation problem in question is represented by the equation

E⁢(p)=(y1y2)

where p∈𝒫3 is the unknown. One can recast the problem into the traditional form by taking standard bases of 𝒫3 and ℝ2 and then seeking all possible a,b,c,d∈ℝ such that

((x1)3(x1)2x11(x2)3(x2)2x21)⁢(abcd)=(y1y2)

However, it is best to treat this problem at an abstract level, rather than mucking about with row reduction. The Lagrange interpolation formula gives us a particular solution, namely the linear polynomial

p0⁢(x)=x-x1x2-x1⁢y1+x-x2x1-x2⁢y2,x∈ℝ

The general solution of our interpolation problem is therefore given as p0+q, where q∈𝒫3 is a solution of the homogeneousPlanetmathPlanetmath problem

E⁢(q)=0.

A basis of solutions for the latter is, evidently,

q1⁢(x)=(x-x1)⁢(x-x2),q2⁢(x)=x⁢q1⁢(x),x∈ℝ

The general solution to our interpolation problem is therefore

p⁢(x)=x-x1x2-x1⁢y1+x-x2x1-x2⁢y2+(a⁢x+b)⁢(x-x1)⁢(x-x2),x∈ℝ,

with a,b∈ℝ arbitrary. The general under-determined interpolation problem is treated in an entirely analogous manner.

Title example of under-determined polynomial interpolation
Canonical name ExampleOfUnderdeterminedPolynomialInterpolation
Date of creation 2013-03-22 12:35:22
Last modified on 2013-03-22 12:35:22
Owner rmilson (146)
Last modified by rmilson (146)
Numerical id 5
Author rmilson (146)
Entry type Example
Classification msc 15A06