Module for Fixed Point Iteration

Check out the new Numerical Analysis Projects page.

 

    A fundamental principle in computer science is iteration.  As the name suggests, a process is repeated until an answer is achieved. Iterative techniques are used to find roots of equations, solutions of linear and nonlinear systems of equations, and solutions of differential equations.

    A rule or function [Graphics:Images/FixPointIterationMod_gr_1.gif] for computing successive terms is needed, together with a starting value [Graphics:Images/FixPointIterationMod_gr_2.gif].  Then a sequence of values [Graphics:Images/FixPointIterationMod_gr_3.gif] is obtained using the iterative rule [Graphics:Images/FixPointIterationMod_gr_4.gif].   The sequence has the pattern
       [Graphics:Images/FixPointIterationMod_gr_5.gif]            (starting value)
       [Graphics:Images/FixPointIterationMod_gr_6.gif]
       [Graphics:Images/FixPointIterationMod_gr_7.gif]
    [Graphics:Images/FixPointIterationMod_gr_8.gif]  
       [Graphics:Images/FixPointIterationMod_gr_9.gif]
    [Graphics:Images/FixPointIterationMod_gr_10.gif]
    [Graphics:Images/FixPointIterationMod_gr_11.gif]  

    What can we learn from an unending sequence of numbers?  If the numberd tend to a limit, we suspect that it is the answer.  

 

Finding Fixed Points

Definition (FixedPoint). A fixed point of a function [Graphics:Images/FixPointIterationMod_gr_12.gif] is a number [Graphics:Images/FixPointIterationMod_gr_13.gif] such that [Graphics:Images/FixPointIterationMod_gr_14.gif].  

Caution.  A fixed point is not a root of the equation  [Graphics:Images/FixPointIterationMod_gr_15.gif],  it is a solution of the equation  [Graphics:Images/FixPointIterationMod_gr_16.gif].  

    Geometrically, the fixed points of a function  [Graphics:Images/FixPointIterationMod_gr_17.gif]  are the point(s) of intersection of the curve  [Graphics:Images/FixPointIterationMod_gr_18.gif]  and the line  [Graphics:Images/FixPointIterationMod_gr_19.gif].  

        

[Graphics:Images/FixPointIterationMod_gr_20.gif]

  

Definition (Fixed Point Iteration). The iteration [Graphics:Images/FixPointIterationMod_gr_21.gif] for [Graphics:Images/FixPointIterationMod_gr_22.gif] is called fixed point iteration.  

Theorem (For a converging sequence). Assume that [Graphics:Images/FixPointIterationMod_gr_23.gif] is a continuous function and that [Graphics:Images/FixPointIterationMod_gr_24.gif] is a sequence generated by fixed point iteration.

If  [Graphics:Images/FixPointIterationMod_gr_25.gif],  then [Graphics:Images/FixPointIterationMod_gr_26.gif] is a fixed point of [Graphics:Images/FixPointIterationMod_gr_27.gif].  

Proof.

 

    The following two theorems establish conditions for the existence of a fixed point and the convergence of the fixed-point iteration process to a fixed point.

Theorem (First Fixed Point Theorem). Assume that [Graphics:Images/FixPointIterationMod_gr_35.gif], i. e.  [Graphics:Images/FixPointIterationMod_gr_36.gif]  is continuous on  [Graphics:Images/FixPointIterationMod_gr_37.gif].  
Then we have the following conclusions.
(i).    If the range of the mapping [Graphics:Images/FixPointIterationMod_gr_38.gif] satisfies [Graphics:Images/FixPointIterationMod_gr_39.gif] for all [Graphics:Images/FixPointIterationMod_gr_40.gif], then  [Graphics:Images/FixPointIterationMod_gr_41.gif] has a fixed point in [Graphics:Images/FixPointIterationMod_gr_42.gif].
(ii).    Furthermore, suppose that [Graphics:Images/FixPointIterationMod_gr_43.gif] is defined over [Graphics:Images/FixPointIterationMod_gr_44.gif] and that a positive constant [Graphics:Images/FixPointIterationMod_gr_45.gif] exists with
    [Graphics:Images/FixPointIterationMod_gr_46.gif]  for all  [Graphics:Images/FixPointIterationMod_gr_47.gif],  then [Graphics:Images/FixPointIterationMod_gr_48.gif] has a unique fixed point [Graphics:Images/FixPointIterationMod_gr_49.gif] in [Graphics:Images/FixPointIterationMod_gr_50.gif].    

Proof of (i).

Proof of (ii).

 

Theorem (Second Fixed Point Fheorem). Assume that the following hypothesis hold true.  
(a)    [Graphics:Images/FixPointIterationMod_gr_80.gif] is a fixed point of a function [Graphics:Images/FixPointIterationMod_gr_81.gif],
(b)    [Graphics:Images/FixPointIterationMod_gr_82.gif],
(c)    [Graphics:Images/FixPointIterationMod_gr_83.gif] is a positive constant,
(d)    [Graphics:Images/FixPointIterationMod_gr_84.gif], and
(e)    [Graphics:Images/FixPointIterationMod_gr_85.gif]  for all  [Graphics:Images/FixPointIterationMod_gr_86.gif].  
Then we have the following conclusions.
(i).    If [Graphics:Images/FixPointIterationMod_gr_87.gif]  for all   [Graphics:Images/FixPointIterationMod_gr_88.gif],  then the iteration  [Graphics:Images/FixPointIterationMod_gr_89.gif]  will converge to the
    unique fixed point [Graphics:Images/FixPointIterationMod_gr_90.gif].  In this case, [Graphics:Images/FixPointIterationMod_gr_91.gif] is said to be an attractive fixed point.  
(ii).    If [Graphics:Images/FixPointIterationMod_gr_92.gif]  for all   [Graphics:Images/FixPointIterationMod_gr_93.gif],  then the iteration  [Graphics:Images/FixPointIterationMod_gr_94.gif]  will not converge to [Graphics:Images/FixPointIterationMod_gr_95.gif].  
    In this case, [Graphics:Images/FixPointIterationMod_gr_96.gif] is said to be a repelling fixed point and the iteration exhibits local divergence.  

 

Proof.

Remark 1.  It is assumed that [Graphics:Images/FixPointIterationMod_gr_129.gif] in statement (ii).

Remark 2.  Because [Graphics:Images/FixPointIterationMod_gr_130.gif]  is continuous on an interval containing [Graphics:Images/FixPointIterationMod_gr_131.gif], it is permissible to use the simpler criterion [Graphics:Images/FixPointIterationMod_gr_132.gif]  and  [Graphics:Images/FixPointIterationMod_gr_133.gif] in (i) and (ii), respectively.  

Corollary. Assume that [Graphics:Images/FixPointIterationMod_gr_134.gif]  satisfies hypothesis (a)-(e)  of the previous theorem.  Bounds for the error involved when using  [Graphics:Images/FixPointIterationMod_gr_135.gif]  to approximate  [Graphics:Images/FixPointIterationMod_gr_136.gif]  are given by         

    [Graphics:Images/FixPointIterationMod_gr_137.gif]              for  [Graphics:Images/FixPointIterationMod_gr_138.gif],  
and
    [Graphics:Images/FixPointIterationMod_gr_139.gif]      for  [Graphics:Images/FixPointIterationMod_gr_140.gif].  

Graphical Interpretation of Fixed-point Iteration

    Since we seek a fixed point [Graphics:Images/FixPointIterationMod_gr_141.gif] to [Graphics:Images/FixPointIterationMod_gr_142.gif],  it is necessary that the graph of the curve  [Graphics:Images/FixPointIterationMod_gr_143.gif]  and the line  [Graphics:Images/FixPointIterationMod_gr_144.gif]  intersect at the point [Graphics:Images/FixPointIterationMod_gr_145.gif].  
The following animations illustrate two types iteration: monotone and oscillating.      

Animations (Fixed Point Iteration  Fixed Point Iteration).  Internet hyperlinks to animations.

Algorithm (Fixed Point Iteration).  To find a solution to the equation  [Graphics:Images/FixPointIterationMod_gr_146.gif] by starting with  [Graphics:Images/FixPointIterationMod_gr_147.gif] and iterating  [Graphics:Images/FixPointIterationMod_gr_148.gif].

 

Pseudo-code

 

Computer Programs (Fixed Point Iteration).

Mathematica Subroutine (Fixed Point Iteration).

[Graphics:Images/FixPointIterationMod_gr_159.gif]

Example 1.  Use fixed point iteration to find the fixed point(s) for the function   [Graphics:Images/FixPointIterationMod_gr_160.gif].

Solution 1.

 

Example 2.  Use Mathematica's built in subroutine "FixedPointList" and experiment finding the fixed point(s) for the function   [Graphics:Images/FixPointIterationMod_gr_278.gif].

Solution 2.

 

Example 3.  Use Mathematica's built in subroutine "FixedPointList" and continue investigation of the "repulsive fixed point" for the function   [Graphics:Images/FixPointIterationMod_gr_294.gif].

Solution 3.

 

Research Experience for Undergraduates

Fixed Point Iteration  Fixed Point Iteration  Internet hyperlinks to web sites and a bibliography of articles.  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004