4

I have the following recurrence relation :

$$g(0) = c $$ $$g(i+1) = g(i) + (1-g(i))*g(i)^{2}$$

where $0 < c < 1$. Is there any closed form for this relation? If not can you give me an upper bound on $g(n)$?

Did
  • 279,727
Kostas
  • 143
  • 5

2 Answers2

7

For every $i\geqslant0$, $$ 1-g(i+1)=(1-g(i)^2)(1-g(i)). $$ This shows that the sequence $(g(i))_{i\geqslant0}$ is increasing and $g(0)=c$ in $(0,1)$ shows that $g(i)$ is in $[c,1)$ for every $i\geqslant0$. Thus the sequence $(g(i))$ converges and its limit $\gamma$ is such that $\gamma\gt c\gt0$ and $1-\gamma=(1-\gamma^2)(1-\gamma)$, that is, $$ \color{red}{\lim\limits_{i\to+\infty}g(i)=1}. $$ For every $i\geqslant0$, $$ 1-g(i+1)=(1+g(i))(1-g(i))^2. $$ This shows that $x(i)=2^{-i}\log(1-g(i))$ solves the recursion $$ x(i+1)=x(i)+2^{-i-1}\log(1+g(i)). $$ Iterating this and starting from $x(0)=\log(1-g(0))=\log(1-c)$, one gets, for every $i\geqslant0$,, $$ \ \qquad\qquad2^{-i}\log(1-g(i))=\log(1-c)+\sum_{k=0}^{i-1}2^{-k-1}\log(1+g(k)).\qquad\qquad(*) $$ Since $\log(1+c)\lt\log(1+g(k))\lt\log2$ for every $k\geqslant0$, the sum on the RHS of $(*)$ is a partial sum of a convergent series, hence $2^{-i}\log(1-g(i))$ converges to a finite limit $-K_c$ such that $0\leqslant K_c\leqslant-\log(1-c^2)$. That is, $$ \color{red}{g(i)=1-\exp(-2^iK_c+o(2^i))}. $$ Finally, the sum on the RHS of $(*)$ being nonnegative, $2^{-i}\log(1-g(i))\geqslant\log(1-c)$, that is, $$ \color{red}{g(i)\leqslant1-(1-c)^{2^i}}. $$ Note 1: The last estimation above can be refined since $\sum\limits_{k=0}^{i-1}2^{-k-1}=1-2^{-i}$ and $g(k)\geqslant c$ for every $k\geqslant0$, hence $$ 2^{-i}\log(1-g(i))\geqslant\log(1-c)+(1-2^{-i})\log(1+c), $$ that is, $$ g(i)\leqslant1-(1-c)^{2^i}(1+c)^{2^i-1}. $$ Note 2: All this does not determine whether $K_c=0$ or $K_c\gt0$. Nevertheless, the semi-explicit formula one gets from $(*)$, namely, $$ -K_c=\log(1-c)+\sum_{k=0}^{+\infty}2^{-k-1}\log(1+g(k)), $$ and the upper bound $g(k)\lt1$ for every $k\geqslant0$ yield $-K_c\lt\log(2(1-c))$ hence $K_c\gt0$ for every $c\geqslant1/2$. From here, using the function $u:c\mapsto c+c^2(1-c)$ and the notation $u^{k+1}=u^k\circ u$ for every $k\geqslant0$, one sees that $K_{u(c)}=2K_c$ for every $c$ in $(0,1)$. For every $c$ in $(0,1)$, there exists some finite nonnegative integer $k(c)$ such that $u^{k(c)}(c)\geqslant1/2$, hence $K_c\geqslant2^{-k(c)}K_{u^{k(c)}(c)}\gt0$.

Finally, $c\mapsto K_c$ is nondecreasing on $[0,1)$, $K_0=0$, $K_c\to+\infty$ when $c\to1$, and $\color{red}{K_c\gt0}$ for every $\color{red}{0\lt c\lt1}$.

Did
  • 279,727
  • Strange this answer only has one upvote (yes, mine). Perhaps because it is very terse. I am guessing we can also show that $K_c = |\log (2 - 2c)|$ – Aryabhata Feb 25 '12 at 21:09
  • Thanks for answering! I am still in the process of understanding the answer. Could you explain in more detail the second and last relation? – Kostas Feb 26 '12 at 00:31
  • Ok! I understood the second relation (i am blind)!!! – Kostas Feb 26 '12 at 01:07
  • I understood the answer!! Thanks a lot, it was extremely helpful!!!! I was actually trying to find an upper bound on $max_{i}(n-i)g(i-1)$ for $1 \leq i \leq n$. I will try maximizing using your bound on g(i)! Thanks again! :) – Kostas Feb 26 '12 at 01:18
  • @Aryabhata: Thanks for the kind words. Tried to untersify a bit. – Did Feb 26 '12 at 10:36
  • Kostas: please tell me if some steps are still unclear. – Did Feb 26 '12 at 10:42
  • @Didier : Thanks a lot for answering! I am not a mathematician, so a lot of the things you wrote might take me some time to realy understand. The second upper bound you gave, was really helpuful! I am currently having trouble finding an upper bound on $max_{i}(n-i)g(i-1)$ using your second upper bound on $g(i)$, that is $g(i-1) \leq 1 - (1-c)^{2^{i-1}}(1+c)^{2^{i-1}-1}$. If you could help me with that i would be most greatful! Thanks again! – Kostas Feb 27 '12 at 01:09
  • What is $n$ in $(n-i)g(i)$? 2. If this is the bound that interests you, you could have asked THAT...
  • – Did Feb 27 '12 at 01:13
  • Assume that $1 \leq i \leq n$....I mentioned it in a comment under my initial post, and as a comment under your first answer...but i thought that provided a bound on $g(i)$ it would then be easy to find the one that interests me (not the case though)!! – Kostas Feb 27 '12 at 01:29