Daily Chess Puzzle – Train your tactical vision with fresh puzzles. Click any card, think, and then reveal the solution in the post body.

Gregory-Newton’s backward Interpolation

Newton's Backward Interpolation: An Interactive Masterclass
A person looking backward over their shoulder, representing the backward method

Newton's Backward Interpolation: A Masterclass on Predicting from the Present

Discover the essential counterpart to the forward method and complete your toolkit for interpolating data with maximum accuracy and efficiency.

The Limitation of Looking Forward

In our last masterclass, we built a powerful tool for predicting the unknown: the Gregory-Newton Forward Interpolation formula. It's an efficient, "upgradable" way to create a polynomialA polynomial is a friendly type of math expression with variables raised to positive whole-number powers, like 3x² + 5x - 7. They are very easy for computers to calculate. that connects our data points. But it has a specific, inherent bias: it is most accurate near the beginning of our data set.

Imagine you are analyzing a dataset of a rocket's altitude over time, with data points recorded every second from $t=0$ to $t=10$. The forward method is the perfect tool to estimate the altitude at $t=0.5$. But what if you need to know the altitude at $t=9.5$? Using the forward formula, which is anchored at $t=0$, would be like trying to aim a target by looking at where you started your journey, not where you are now. It can be inaccurate and inefficient.

To solve this, we need a different perspective. We need a formula that is anchored at the end of our data and works backwards. This is the Gregory-Newton Backward Interpolation Formula, and it's the essential other half of our interpolation toolkit.


The Backward Difference: A New Perspective on the Same Data

The engine behind the backward formula is the backward difference operator, denoted by the symbol $\nabla$ (nabla). While the forward operator $\Delta$ looked ahead, the backward operator $\nabla$ looks behind.

The Backward Difference Operator ($\nabla$)

The definition is simple: the backward difference is the current value minus the previous value.

$$ \nabla y_i = y_i - y_{i-1} $$

Just like with forward differences, we can apply this operator repeatedly to find higher-order differences:

  • Second Backward Difference ($\nabla^2 y_i$): $\nabla^2 y_i = \nabla y_i - \nabla y_{i-1}$.
  • Third Backward Difference ($\nabla^3 y_i$): $\nabla^3 y_i = \nabla^2 y_i - \nabla^2 y_{i-1}$.

The Difference Table Revisited

Here is the most important insight: a backward difference table contains the exact same numbers as a forward difference table. The only thing that changes is the notation and which diagonal of coefficients we use for our formula.

For example, the first forward difference $\Delta y_0 = y_1 - y_0$ is identical to the first backward difference at the next point, $\nabla y_1 = y_1 - y_0$. The table is the same; only our perspective changes. The interactive lab below makes this clear.

The Difference Table Generator

Enter your equally-spaced data points. Use the toggle to switch between "Forward" and "Backward" modes to see which diagonal coefficients are used for each formula.


The Gregory-Newton Backward Difference Formula

The backward formula is constructed similarly to the forward one, but it starts from the last data point, $y_n$, and uses the coefficients from the backward-pointing diagonal (highlighted in red in our lab).

We define a normalized variable $s$, but this time, it's relative to the last point, $x_n$:

$$ s = \frac{x - x_n}{h} $$

The Gregory-Newton Backward Difference Formula is then:

$$ P_n(x) = y_n + s\nabla y_n + \frac{s(s+1)}{2!}\nabla^2 y_n + \frac{s(s+1)(s+2)}{3!}\nabla^3 y_n + \dots $$

Notice the subtle but crucial difference in the binomial coefficientA term used in combinatorics, often written as 'n choose k'. Here, it's a way of creating polynomial terms of increasing degree. terms. Instead of $s(s-1)$, we now have $s(s+1)$, and so on. This change in sign is what allows the formula to properly build the polynomial by working backwards from the end of the data set.


A Strategic Choice: Forward, Backward, or Central?

With multiple formulas at our disposal, how do we choose which one to use? The decision is based entirely on where the point we want to estimate, $x$, lies within our dataset.

  • Newton's Forward Formula: Uses the top diagonal of the table ($y_0, \Delta y_0, \dots$). It is most accurate for interpolating values of $x$ that are near the beginning of the data set (i.e., close to $x_0$).
  • Newton's Backward Formula: Uses the bottom diagonal of the table ($y_n, \nabla y_n, \dots$). It is most accurate for interpolating values near the end of the data set (i.e., close to $x_n$).
  • Central Difference Formulas (e.g., Stirling's Formula): These more advanced formulas use a zig-zag path through the *center* of the difference table. They provide the highest accuracy for interpolating values near the middle of the data set.

Choosing the right formula for the job is a key skill. If your data is a time series and you want to predict the near future, you would use the most recent data points and the backward formula. If you are trying to fill in a missing data point at the beginning of a test, you would use the forward formula.


Solving Numerical Problems: A Step-by-Step Guide

Problem 1: Using Newton's Backward Formula

Question: The following data represents the velocity of a rocket. Use a third-order Newton's backward difference polynomial to estimate the velocity at $t=27$ seconds.

Time (t)10152025
Velocity (v)35486590

Step 1: Construct the Difference Table
The step size is $h=5$. We need to find the backward differences.

tv∇v∇²v∇³v
1035
154813
2065174
25902584

The coefficients from the bottom diagonal (starting from the last point, $t_n=25$) are $y_n=90$, $\nabla y_n=25$, $\nabla^2 y_n=8$, and $\nabla^3 y_n=4$.

Step 2: Calculate the parameter 's'
Our last point is $x_n=25$. We want to estimate at $x=27$. The step size is $h=5$. $$s = \frac{x - x_n}{h} = \frac{27 - 25}{5} = \frac{2}{5} = 0.4$$

Step 3: Substitute into the Backward Formula
We use the third-order formula: $$P(x) = y_n + s\nabla y_n + \frac{s(s+1)}{2!}\nabla^2 y_n + \frac{s(s+1)(s+2)}{3!}\nabla^3 y_n$$ Plugging in our values: $$P(27) = 90 + (0.4)(25) + \frac{(0.4)(0.4+1)}{2}(8) + \frac{(0.4)(0.4+1)(0.4+2)}{6}(4)$$ $$P(27) = 90 + 10 + \frac{(0.4)(1.4)}{2}(8) + \frac{(0.4)(1.4)(2.4)}{6}(4)$$ $$P(27) = 90 + 10 + (0.28)(8) + (0.224)(4)$$ $$P(27) = 90 + 10 + 2.24 + 0.896 = 103.136$$

Final Answer: The estimated velocity at $t=27$ seconds is 103.136.


Test Your Intuition!

Newton Interpolation Quiz

No comments

No comments yet. Be the first!

Post a Comment

Search This Blog

Explore More Topics

Loading topics…