Fixed Point Iteration, does it converge?

11.1k Views Asked by At

Find the fixed point(s) of $g(x) = x^2 + 3x - 3$. Does the fixed point iteration(s) converge(s) to the fixed points if you start with a close enough first approximation?

I set $g(x) = x$ and got $g(x) = x^2 + 2x - 3$. So either $x=1$ or $x=-3$. How do I find out if it converges or diverges?

4

There are 4 best solutions below

0
On

Hint: the answer depends on the derivative of $g$ at the fixed point

3
On

If we write

$$x_{n+1} = G(x_n)$$

The error term is given by:

$$e_n = [G'(r)]^ne_0$$

where $r$ is the root.

The error reduces if $|G'(r)| \lt 1$ and increases if $|G'(r)| \gt 1$.

Also, are you sure you wrote $G(x)$ correctly?

Update

  • Let's see, we have a root at $r = 1$ or $r = -3$.
  • $|G'(r)| = |2r+2|$.
  • At $r= 1, ~~~|G'(r)| = |4|>1 \rightarrow$ diverges.
  • At $r= -3, |G'(r)| = |4|>1 \rightarrow$ diverges.
1
On

One way to understand the qualitative behaviour of iteration of $g$ is a cobweb diagram. Read that article and try to apply it to your situation. After doing that, you should understand Robert Israels hint to look at the derivative.

0
On

First, the fixed points of $g(x)=x^2+2x-3$ are not $1$ and $-3$, those are the roots. The fixed points are the solutions to $x^2+2x-3=x$ or, by the quadratic formula, $(-1\pm\sqrt{13})/2$.

To show that the larger of these $(\sqrt{13}-1)/2$ is repulsive, we can re-write $g(x)$ as

$$g(x)=\left(x+\frac{1}{2} \left(1-\sqrt{13}\right)\right)^2+\left(1+\sqrt{13}\right) \left(x+\frac{1}{2} \left(1-\sqrt{13}\right)\right)+\frac{1}{2} \left(\sqrt{13}-1\right).$$

You can verify this equality by expanding it out; you can even have WolframAlpha do the algebra for you.

Now, consider a quadratic of the form

$$f(x) = (x-x_0)^2 + m(x-x_0) + x_0,$$

where $m>1$. This is exactly what we have for $x_0 = (\sqrt{13}-1)/2$ and $m=1+\sqrt{13}$.

Now, I claim that $$|f(x)-x_0|>|x-x_0|$$ for $x$ sufficiently close to $x_0$. This follows from the equality $$|f(x)-x_0| = |(x-x_0) + m| \cdot |x-x_0|.$$ All we need is $|(x-x_0) + m|>1$, which will certainly be true for $x$ close to $x_0$, since $m>1$.

Note that $m$ is exactly the value of the derivative of $f$ at $x_0$ and what we see near the fixed point is that the linear term in the expansion dominates the behavior of the function. This is exactly the hint that the others were trying to offer.