Simple curve - not so simple function?

  • Context: Undergrad 
  • Thread starter Thread starter Miffymycat
  • Start date Start date
  • Tags Tags
    Curve Function
Click For Summary
SUMMARY

The discussion centers around identifying the function f(x) based on a given sequence of y-values corresponding to x-values. The sequence is defined by a second-order difference equation, leading to a characteristic equation with solutions that suggest a general solution involving exponential terms. Participants explore whether the function can be classified as a polynomial, hyperbola, or a variation of a power function, ultimately concluding that it resembles a power function with specific characteristics. The conversation highlights the complexities of curve fitting and the mathematical relationships involved.

PREREQUISITES
  • Understanding of second-order difference equations
  • Familiarity with characteristic equations and their solutions
  • Knowledge of polynomial functions and their properties
  • Basic linear algebra concepts, including matrix diagonalization
NEXT STEPS
  • Study the properties of second-order difference equations
  • Learn about characteristic equations and their applications in sequences
  • Explore polynomial fitting techniques in software like Excel or R
  • Investigate power functions and their variations in mathematical modeling
USEFUL FOR

Mathematicians, data analysts, chemists working with data modeling, and anyone interested in advanced mathematical functions and curve fitting techniques.

Miffymycat
Messages
47
Reaction score
0
What is f(x) when eg

x = 0, 1, 2, 3, 4, 5, 6, etc
y = 0, 8, 12, 14, 15, 15.5, 15.75 etc (the y value at x =1 is arbitrary)

ie each successive y value adds half the difference of the preceding 2 values.

(It is effectively the inverse of a first order exponential decay where the y value halves at constant x intervals).

But it doesn’t fit an exponential ...! Could it be a type of hyperbola? or polynomial? or both? or something else entirely?

I have no decent curve fitting software and struggling to find this.
Thank you!
 
Physics news on Phys.org
Try drawing it on a graph paper.
 
Thanks ... but how would that help, if I've already tried to fit it in Excel?
 
You can try finding the nth term of the sequence(That's what I do).
 
So you are saying that x_{n+2}= x_{n+1}+ (x_n+ x_{n+1})/2= (3/2)x_{n+1}+ (1/2)x_n

That's a "second order difference equation" which has associated characteristic equation s^2= (3/2)s+ 1/2 or s^2- (3/2)s- 1/2= 0. The solutions to that equation are s= \frac{3\pm \sqrt{17<br /> }}{4}. That means that the general solution to the difference equation is 8^{n/4}(C2^{\sqrt{17}n/4}+ D2^{-\sqrt{17}n/4})
 
Last edited by a moderator:
So you have a sequence defined by

x_0 = 0,~x_1 = 8,~x_{n+1} = x_n + \frac{1}{2}(x_n - x_{n-1}) = \frac{3}{2}x_n - \frac{1}{2}x_{n-1}

Let's find a general term. We can express this sequence as follows:

\left(<br /> \begin{array}{cc}<br /> x_{n+1}\\ x_n<br /> \end{array}<br /> \right)<br /> =<br /> \left(<br /> \begin{array}{cc}<br /> \frac{3}{2} &amp; -\frac{1}{2}\\ 1 &amp; 0<br /> \end{array}<br /> \right)<br /> \left(<br /> \begin{array}{cc}<br /> x_{n}\\ x_{n-1}<br /> \end{array}<br /> \right)<br />

So, we get

\left(<br /> \begin{array}{cc}<br /> x_{n+1}\\ x_n<br /> \end{array}<br /> \right)<br /> =<br /> \left(<br /> \begin{array}{cc}<br /> \frac{3}{2} &amp; -\frac{1}{2}\\ 1 &amp; 0<br /> \end{array}<br /> \right)^n<br /> \left(<br /> \begin{array}{cc}<br /> x_1\\ x_0<br /> \end{array}<br /> \right)<br />

By using the diagonalization theory of linear algebra, we can write this as

\left(<br /> \begin{array}{cc}<br /> x_{n+1}\\ x_n<br /> \end{array}<br /> \right)<br /> =<br /> \left(<br /> \begin{array}{cc}<br /> \frac{1}{2} &amp; 1\\ 1 &amp; 1<br /> \end{array}<br /> \right)<br /> \left(<br /> \begin{array}{cc}<br /> \frac{1}{2^n} &amp; 0\\ 0 &amp; 1<br /> \end{array}<br /> \right)<br /> \left(<br /> \begin{array}{cc}<br /> -2 &amp; 2\\ 2 &amp; -1<br /> \end{array}<br /> \right)<br /> \left(<br /> \begin{array}{cc}<br /> x_1\\ x_0<br /> \end{array}<br /> \right)<br />

and thus

\left(<br /> \begin{array}{cc}<br /> x_{n+1}\\ x_n<br /> \end{array}<br /> \right)<br /> =<br /> \left(<br /> \begin{array}{cc}<br /> 2 - \frac{1}{2^n} &amp; \frac{1}{2^n}- 1\\ 2 - \frac{1}{2^{n-1}} &amp; \frac{1}{2^{n-1}}-1<br /> \end{array}<br /> \right)<br /> \left(<br /> \begin{array}{cc}<br /> x_1\\ x_0<br /> \end{array}<br /> \right)<br />

Thus we get

x_{n+1} = 2x_1 - \frac{x_1}{2^n} + \frac{x_0}{2^n} - x_0

Since ##x_1 = 8## and ##x_0 = 0##, we get

x_{n+1} = 16 - \frac{8}{2^n}
 
Thanks both! I'm not a mathematician ... just a chemist struggling to make sense of some data!

So Hallsofivy, is a "second order difference equation" the same thing as a second order polynomial? I tried to fit the data to a polynomial and it required 6th order before R2= 1 in Excel! (And not sure how you got from the series to the associated characteristic equation)

I see you guys have both written it as a series, but can it be expressed in terms of y=f(x)? And described as eg a power function or polynomial etc?
 
You can manipulate Micromass's approach to get, for example:

f(x) = 16 - 2^{4-x}

And, if f(1) = a, then:

f(x) = a(2 - 2^{1-x}) = 2a(1 - \frac{1}{2^x})
 
Last edited:
... building on micromass's reply ... something like y = 2a -(a/2^x)? where a is an arbitrary constant. Is it a hyperbolic function??! Or a hyperbolic polynomial! does that exist?
 
  • #10
Thanks Perok ... I was sending mine while yours was arriving - I was close! So does the function fit one of those categories?
 
  • #11
I would say it's just a variation of a power function (powers of 1/2).
 
  • #12
OK - a bit weird though, that it fits precisely to polynomial order 6! Coincidence?
 
  • #13
Miffymycat said:
OK - a bit weird though, that it fits precisely to polynomial order 6! Coincidence?

Any set of n values will fit a polynomial of order n-1.
 
  • Like
Likes   Reactions: 1 person

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
915
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 8 ·
Replies
8
Views
4K
  • · Replies 1 ·
Replies
1
Views
913
  • · Replies 6 ·
Replies
6
Views
3K