next up previous
Next: Descent methods Up: No Title Previous: Newton's method

Broyden's Method

Broyden's method is a generalization of the secant method to nonlinear systems. Roughly speaking, the secant method replaces the derivative by a finite difference

eqnarray50

Broyden gives a systems generalization of this formula

displaymath144

Instead of using the Jacobian matrix, Broyden uses a matrix tex2html_wrap_inline146 satisfying the same formula. Note the system provides only N equations to determine the NXN matrix. The ``best possible'' choice for tex2html_wrap_inline146 is a minimal modification of tex2html_wrap_inline150 , tex2html_wrap_inline152 . This leads to a uniquely defined matrix. Set tex2html_wrap_inline154 Then

displaymath156

Broyden's method then iterates

displaymath158

Of course, once again we are faced with finding an inverse matrix, something we would like to avoid. rather, one can solve tex2html_wrap_inline160 and then update tex2html_wrap_inline162 .

Broyden also suggested using the Sherman-Morrison formula (sometimes called the `bad Broyden' method) which tells us how to compute the inverse in an efficient manner. The formula says that, for a nonsingular matrix A and vectors u, v, such that tex2html_wrap_inline164 , then

displaymath166

To implement the Sherman-Morrison formula in Broyden, then, let tex2html_wrap_inline168 . A count shows that this formula leads to an update of the inverse costing only tex2html_wrap_inline170 operations.

Broyden's method is not as fast as the quadratic convergence of Newton's method. But the smaller operation count per iteration is often worth the trade-off.



E. Bruce Pitman
Tue Nov 17 10:17:29 EST 1998