A digression into optimization, needed to explain how each iteration refines its approximate solution (43:35, 44:23).
The setting
- Goal: iteratively compute
x* = argmin_{x ∈ ℝⁿ} f(x)forf : ℝⁿ → ℝ(44:43). - First-order method assumptions:
- Hessian at
x: then × nmatrix with(∇²f(x))_{ij} = ∂²f(x)/∂xᵢ∂xⱼ(46:56).
The Hessian condition
- Assume constants
0 < μ < Lwithμ·I ⪯ ∇²f(x) ⪯ L·Ifor allx ∈ ℝⁿ(48:38,49:13). - Loewner ordering:
A ⪯ BmeansB − Ais positive semidefinite, equivalentlyxᵀ(B − A)x ≥ 0for allx(49:13,49:56,50:56). - Interpretation: the Hessian's largest eigenvalue is at most
L, its smallest is at leastμ— the matrix is well-conditioned, with no eigenvalue too large or too close to zero (50:30). - Justification sketched via the Rayleigh quotient: the largest eigenvalue is
max_{‖x‖=1} xᵀAx, i.e.xᵀAx / xᵀxin general; comparing againstxᵀ(L·I)x = L‖x‖²gives the eigenvalue bound. Taken on faith for this course (51:37,52:07,52:29,52:50).
Orders of methods
- k-th order method = the update rule has access to k-th order partial derivatives; using the Hessian in the update rule makes a method second order (
57:26). - Interior point will use a second-order method (Newton's) and will not assume control over third-order partials — a different kind of control is used instead (
57:44,58:02). - Newton's method is not a pure first-order method here: its update rule genuinely needs the Hessian (
58:24).
What an iterative method is
- Start at
x₀ ∈ ℝⁿ, computex_{k+1}fromx_kvia an update rule, stop after some number of iterations (53:36). - Target guarantee:
f(x_{k+1}) − f(x*) ≤ α·( f(x_k) − f(x*) )withαsmall (55:08,56:38). - By induction the error after
ksteps is multiplied byαᵏ, soα < 1means exponentially decaying error (56:38,57:03).