Looking for General solution for a difference equation

Click For Summary

Discussion Overview

The discussion revolves around finding a general solution for a specific difference equation of the form At+1=(At+r)/(At+r+1), with participants exploring various approaches and methods to derive solutions, including special cases and asymptotic behavior. The scope includes theoretical exploration and mathematical reasoning.

Discussion Character

  • Exploratory
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant suggests setting At+1=At=A to find a special solution, but seeks a general solution.
  • Another participant presents a modified equation (1+r)At+1 - At + At+1*At = r, expressing skepticism about its simplicity.
  • There is a clarification regarding the formulation of the second equation, with participants confirming its intended structure.
  • A participant rewrites the original equation in a different form, suggesting a method to express A2 in terms of A1 and iteratively derive A_t.
  • One participant proposes that while a general solution is difficult to find, asymptotic behavior can be analyzed by assuming a form A_k = A + p_k, where p_k approaches zero as k increases.
  • Another participant elaborates on the Taylor series expansion to analyze the asymptotic behavior, leading to a recurrence relation for p_k.
  • A participant mentions that the series will only be a real approximation if the absolute value of (1-A) is less than 1, indicating dependence on the values of A and r.
  • One participant introduces the concept of a geometric series and discusses its implications for linear approximations.

Areas of Agreement / Disagreement

Participants express varying degrees of confidence regarding the difficulty of finding a general solution, with some focusing on specific cases and others exploring asymptotic behavior. No consensus is reached on a definitive general solution.

Contextual Notes

The discussion highlights limitations in deriving a general solution and emphasizes the dependence on specific values of A and r for stability in the asymptotic analysis.

ivyhawk
Messages
12
Reaction score
0
At+1=(At+r)/(At+r+1)

A1=constant

I know I can set At+1=At=A and solve for a special solution.

What would be a general solution?

I am not taking a course in Difference Equation, and this is not my homework but I encounter a similar question and I reduce it to this form.

Thanks
 
Last edited:
Physics news on Phys.org
Or One step further, how to solve this one?

(1+r)A_t+1-A_t+A_t+1*A_t=r

A_t means A with t subsciptI don't think it is as easy as it looks.
 
Last edited:
ivyhawk said:
I don't think it is as easy as it looks.

It doesn't even look like anything. :-p

Is your second one meant to be

(1+r)At+1-At+At+1*At=r ?
 
epenguin said:
It doesn't even look like anything. :-p

Is your second one meant to be

(1+r)At+1-At+At+1*At=r ?

Yes That is what I mean. And it is just algebra manipulation of the original equation
 
Your starting form is

[tex]A_{t+1} = \frac{A_t + r}{A_t + r + 1},[/tex]
which is perhaps more convenient to write as

[tex]A_{t+1} = 1 - \frac{1}{A_t+r+1}[/tex]

What you could to do try to solve it is find A_2 in terms of A_1, plug that into A_3 to get it in terms of A_1, and so on, to see if you can guess what A_t for general t will look like. (I don't know if this one will turn out nice. I did the first couple of steps myself, and it looks like you may want to look into http://en.wikipedia.org/wiki/Continued_fraction (go to the page on periodic continued fractions too). There may be nice solutions for specific values of r or A_1. (Of course, A_1 will be irrelevant as t -> infinity)
 
Last edited by a moderator:
Specific solution for A_t is easy to solve

Set A_t=A=A_t+1

We can solve that r=A/(1-A) or the other way, A=(1/2)(-r+(-)(r^2+4r)^1/2

But I still find it very hard to find a general solution.
 
A general solution is hard to find.
But you can get some idea about the asymptotic behaviour when the sequence is approaching equilibrium A.

So assume your solution is

[tex]A_{k}=A+p_{k}[/tex]

Where the sequence p goes to zero as k increases. At the limit of large k, you can expand the RHS to a taylor series. Because we want a linear easy solution, we will only care about the first order. So:

[tex]A+p_{k+1}=1-\frac{1}{A+r+1}\frac{1}{1+\frac{p_{k}}{A+r+1}}=1-\frac{1}{A+r+1}(1-\frac{p_{k}}{A+r+1})[/tex]

Let A+r+1=b to make it look nicer. Also 1-1/b=A. So:

[tex]A+p_{k+1}=A+(1-A)^{2}p_{k}[/tex]

Or

[tex]p_{k+1}-(1-A)^{2}p_{k}=0[/tex]

So at large k's, the series will be:

[tex]A_{k}=A+(1-A)^{2k}[/tex]

This will be a real approximation only if abs(1-A)<1 (the solution is really stable), and that depends on what value of A you have chose, and what is r.
 
How did u expand it?


elibj123 said:
A general solution is hard to find.
But you can get some idea about the asymptotic behaviour when the sequence is approaching equilibrium A.

So assume your solution is

[tex]A_{k}=A+p_{k}[/tex]

Where the sequence p goes to zero as k increases. At the limit of large k, you can expand the RHS to a taylor series. Because we want a linear easy solution, we will only care about the first order. So:

[tex]A+p_{k+1}=1-\frac{1}{A+r+1}\frac{1}{1+\frac{p_{k}}{A+r+1}}=1-\frac{1}{A+r+1}(1-\frac{p_{k}}{A+r+1})[/tex]

Let A+r+1=b to make it look nicer. Also 1-1/b=A. So:

[tex]A+p_{k+1}=A+(1-A)^{2}p_{k}[/tex]

Or

[tex]p_{k+1}-(1-A)^{2}p_{k}=0[/tex]

So at large k's, the series will be:

[tex]A_{k}=A+(1-A)^{2k}[/tex]

This will be a real approximation only if abs(1-A)<1 (the solution is really stable), and that depends on what value of A you have chose, and what is r.
 
This is a geometrical series:

[tex]\frac{1}{1-q}=1+q+q^{2}+q^{3}+...[/tex]

setting q-> (-q) the terms will occur with alternating signs so the linear approximation will be

[tex]\frac{1}{1+q}=1-q[/tex]
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 52 ·
2
Replies
52
Views
9K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 0 ·
Replies
0
Views
4K
  • · Replies 36 ·
2
Replies
36
Views
6K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K