Skip to content

Question

Math 6610 • Analysis of Numerical Methods I

Fall 1997 • Question 7

Broyden’s Method Consider the nonlinear system: \[ F(x)=0\qquad \text {where }F:\mathbb {R}^n\to \mathbb {R}^n \] and the iteration \[ x^{(k+1)}=x^{(k)}-(A^{(k)})^{-1}F(x^{(k)}) \] where \(x^{(0)}\in \mathbb {R}^n\) is given.

Newton’s Method is defined by the choice \(A^{(k)}=F'(x)\) where \(F'(x)\) is the Jacobean matrix of \(F\). An alternative choice for \(A^{(k+1)}\) is made in Broyden’s Method and is motivated as follows:

Let \(A^{(0)}\) be a given and suitably chosen matrix. The scalar Secant Method suggests imposing on \(A^{(k+1)}\) the condition \[ A^{(k+1)}(x^{(k+1)}-x^{(k)})=F(x^{(k+1)})-F(x^{(k)}). \] This however does not uniquely determine \(A^{k+1}\) as it does not give any information on how \(A^{k+1}\) acts on vectors perpendicular to \((x^{(k+1)}-x^{(k)})\). Thus require also that: \[ A^{(k+1)}t=A^{(k)}t \qquad \text {for all }t\in \mathbb {R}^n\text { such that } t^T(x^{(k+1)}-x^{(k)})=0 \] Use these conditions to derive an explicit expression for \(A^{(k+1)}\) in terms of \(A^{(k)}\), \(x^{(k)}\), and \(x^{(k+1)}\).

Raw LaTeX

{Broyden's Method}
Consider the nonlinear system:
\[
  F(x)=0\qquad \text{where }F:\R^n\to\R^n
\]
and the iteration
\[
  x^{(k+1)}=x^{(k)}-(A^{(k)})^{-1}F(x^{(k)})
\]
where \(x^{(0)}\in\R^n\) is given.

Newton's Method is defined by the choice \(A^{(k)}=F'(x)\) where \(F'(x)\) is the
Jacobean matrix of \(F\). An alternative choice for \(A^{(k+1)}\) is made in Broyden's
Method and is motivated as follows:

Let \(A^{(0)}\) be a given and suitably chosen matrix. The scalar Secant Method
suggests imposing on \(A^{(k+1)}\) the condition
\[
  A^{(k+1)}(x^{(k+1)}-x^{(k)})=F(x^{(k+1)})-F(x^{(k)}).
\]
This however does not uniquely determine \(A^{k+1}\) as it does not give any
information on how \(A^{k+1}\) acts on vectors perpendicular to
\((x^{(k+1)}-x^{(k)})\). Thus require also that:
\[
  A^{(k+1)}t=A^{(k)}t
  \qquad
  \text{for all }t\in\R^n\text{ such that }
  t^T(x^{(k+1)}-x^{(k)})=0
\]
Use these conditions to derive an explicit expression for \(A^{(k+1)}\) in terms of
\(A^{(k)}\), \(x^{(k)}\), and \(x^{(k+1)}\).