Md. Asif Uddin
    Problem I.9.B03

    Small eigenvalues can hide a large transient

    counterexample▲▲▲

    Exact values where possible; six decimal places for non-exact numerical results.

    STATEMENT

    Disprove the claim that recurrent eigenvalues below one prevent every finite-time amplification.

    GIVEN

    Consider the linear state update ht=Wht−1\vec h_t=\mat W\vec h_{t-1} with

    W=[1/2401/2].\mat W=\begin{bmatrix}1/2&4\\0&1/2\end{bmatrix}.

    Both eigenvalues are 1/21/2. Use Euclidean vector norms and let e2=(0,1)T\vec e_2=(0,1)^{\mathsf T}.

    FIND

    A formula for Wn\mat W^n, its effect on e2\vec e_2, and why this does not contradict eventual decay or the norm bound (I.9.6).

    STRATEGY

    Write W=(1/2)I+N\mat W=(1/2)\mat I+\mat N. Since N2=0\mat N^2=0, only two terms survive in the binomial expansion.

    SOLUTION

    For integer n≥1n\geq1,

    Wn=2−nI+n 2−(n−1)N=[2−n4n 2−(n−1)02−n].\begin{aligned} \mat W^n &=2^{-n}\mat I+n\,2^{-(n-1)}\mat N\\ &=\begin{bmatrix} 2^{-n}&4n\,2^{-(n-1)}\\ 0&2^{-n} \end{bmatrix}. \end{aligned}

    It follows that

    We2=(4,1/2)T,W2e2=(4,1/4)T.\begin{aligned} \mat W\vec e_2&=(4,1/2)^{\mathsf T},\\ \mat W^2\vec e_2&=(4,1/4)^{\mathsf T}. \end{aligned}

    Their norms are approximately 4.0311294.031129 and 4.0078054.007805, although the initial vector has norm one. At ten steps the first coordinate is 0.0781250.078125 and the second is 0.00097656250.0009765625.

    The matrix is non-normal: it does not commute with its transpose. Its eigenvalue magnitudes describe the eventual exponential rate, not a uniform bound on every finite power. The off-diagonal term introduces a factor of nn before the exponential decay wins. Both 2−n2^{-n} and n2−(n−1)n2^{-(n-1)} tend to zero, so Wn\mat W^n eventually tends to the zero matrix.

    Backward propagation uses (Wn)T(\mat W^n)^{\mathsf T}. Its spectral norm equals that of Wn\mat W^n, so a backward direction can be amplified too. For example, WT(1,0)T=(1/2,4)T\mat W^{\mathsf T}(1,0)^{\mathsf T}=(1/2,4)^{\mathsf T}.

    There is no conflict with (I.9.6). Indeed ∥W∥2≥∥We2∥2>4\|\mat W\|_2\geq\|\mat W\vec e_2\|_2>4. The eigenvalue 1/21/2 is not a valid replacement for the factor norm in that inequality.

    Answer

    Eigenvalues inside the unit circle permit substantial transient amplification. Here a unit perturbation grows by more than four in one step, then ultimately decays.

    Check — sanity

    At n=1n=1 the power formula returns the supplied matrix. At n=2n=2 its off-diagonal entry is four, matching direct multiplication.

    Where this breaks

    For a nonlinear cell the temporal Jacobians can change at every step. A formula for powers of one fixed matrix is then insufficient. Even a large product norm only guarantees an amplified direction, not that the actual loss gradient points along it.

    Variation

    Set the off-diagonal entry to zero. The matrix becomes (1/2)I(1/2)\mat I, and every vector contracts by 2−n2^{-n} after nn steps. Identify the term responsible for the transient in the original example.

    Draws on