Prove quadratic convergence and limit

Click For Summary
SUMMARY

The discussion focuses on proving quadratic convergence of a fixed point iteration defined by the function g, which is twice continuously differentiable near the fixed point xi. The iterative formula is given as z(n+1) = g(g(z(n))) - [g(g(z(n))) - g(z(n))]^2 / [g(g(z(n))) - 2*g(z(n)) + z(n)]. Participants utilize Taylor series expansions of g(z(n)) and g(g(z(n))) to derive expressions that demonstrate convergence. The final expression for z(n+1) incorporates derivatives of g at xi, indicating the complexity of the convergence proof.

PREREQUISITES
  • Understanding of fixed point theory and convergence concepts
  • Familiarity with Taylor series expansions
  • Knowledge of calculus, particularly derivatives and their applications
  • Experience with iterative methods in numerical analysis
NEXT STEPS
  • Study the properties of quadratic convergence in numerical methods
  • Learn about Taylor series and their applications in function approximation
  • Explore fixed point iteration methods and their convergence criteria
  • Investigate simplification techniques for complex derivative expressions
USEFUL FOR

Mathematicians, numerical analysts, and students studying fixed point iteration methods and convergence proofs will benefit from this discussion.

radutanasa
Messages
1
Reaction score
0
Hello,

I have this problem that I simply cannot nail down. Please help.
xi - fixed point of g, g is twice continuously differentiable in a vecinity of xi.

z(n+1) = g(g(z(n))) - [g(g(z(n))) - g(z(n))]^2 / [g(g(z(n))) - 2*g(z(n)) + z(n)]

Using taylor series expansion of g(z(n)) and g(g(z(n))) in a vecinity of xi I have to prove that xi is z's limit and that the convergence is quadratic.

Thank you for your help!

Radu
 
Physics news on Phys.org
Interesting problem. I don't know a smart way to do it, so I just dive right in and calculate the series up to two terms around x_i (I won't bother to indicate the higher order terms):

[tex] g(z) = \frac{1}{2} g''\left(x_i\right) \left(z-x_i\right){}^2+g'\left(x_i\right) \left(z-x_i\right)+g\left(x_i\right)[/tex]

[tex] g(g(z)) = \frac{1}{2} \left(g''\left(g\left(x_i\right)\right) g'\left(x_i\right){}^2+g'\left(g\left(x_i\right)\right) g''\left(x_i\right)\right)<br /> \left(z-x_i\right){}^2+g'\left(g\left(x_i\right)\right) g'\left(x_i\right) \left(z-x_i\right)+g\left(g\left(x_i\right)\right)<br /> [/tex]

Now if we go ahead and replace these into your expression and use the fact that g(g(x_i) = g(x_i) = x_i we get the following for z_{n + 1}:

[tex]\frac{g''\left(x_i\right){}^2 \left(x_i-z\right){}^3+4 x_i \left(g'\left(x_i\right)-1\right){}^2+2 \left(z-x_i\right) \left(g'\left(x_i\right)-1\right) \left(2 x_i+z<br /> g'\left(x_i\right)\right) g''\left(x_i\right)}{4 \left(g'\left(x_i\right)-1\right){}^2+2 \left(z-x_i\right)<br /> \left(g'\left(x_i\right){}^2+g'\left(x_i\right)-2\right) g''\left(x_i\right)}[/tex]

Note: I dropped the subscript z_n because of formatting problems with LaTeX in these forums.

I do not know how to simplify the expressions involving the derivatives, but maybe this is already to the point where simplifying into a particular form might make it easy to see the convergence.
 
Last edited:

Similar threads

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