Solving this type of recurrence equation

  • Context: Graduate 
  • Thread starter Thread starter samuelandjw
  • Start date Start date
  • Tags Tags
    Recurrence Type
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
4 replies · 3K views
samuelandjw
Messages
22
Reaction score
0
Hi,

The problem is to solve
[tex](1-g_{i+1})P_{i+1}-P_{i}+g_{i-1}P_{i-1}=0[/tex]
for [tex]P_{i}[/tex]
with boundary condition
[tex]P_{i}=P_{i+L}, g_{i}=g_{i+L}[/tex]
Can anyone provide any guide of solving this type of recurence equation?
Thank you!
 
Last edited:
Physics news on Phys.org
I don't know of any general methods for solving recurrence relations in several unknown functions.

Since setting g and P to be constant sequences gives a solution to the problem, I suggest that you state additional conditions that g and P must satistfy if having them constant isn't what you are after.
 
Stephen Tashi said:
I don't know of any general methods for solving recurrence relations in several unknown functions.

Since setting g and P to be constant sequences gives a solution to the problem, I suggest that you state additional conditions that g and P must satistfy if having them constant isn't what you are after.

Thanks for your reply.

In my problem, [tex]g_i[/tex] is given and arbitrary, so in general [tex]g_i[/tex] is not a constant sequence.
 
Then, as I interpret the problem, it amounts to solving [itex]L[/itex] simultaneous linear equations with constant coefficients and unknowns [itex]P_1,P_2,..P_L[/itex].

Are you looking for a closed form symbolic answer instead of a numerical one?

( The wikipedia article on "recurrence relation" says that linear constant coefficient difference equations can be solved with z-transforms. I, myself, have never done that. )
 
Because both [itex]\lbrace P_{i}\rbrace[/itex] and [itex]\lbrace g_{i}\rbrace[/itex] are period with a common period of L, you should use the discrete Fourier transform:

[tex] P_{i} = \frac{1}{\sqrt{L}} \, \sum_{k = 0}^{L - 1}{\tilde{P}_{k} \, \exp{\left(\frac{2 \pi j k \, i}{L}\right)}}[/tex]
and similarly for [itex]g_{i}[/itex]. Then you will use the convolution theorem for the products.