Repeated Eigenvalues — Question 6

PDF ↗

Question 6

Let A=(−1100−1100−1),X′=AX,X(0)=(pqr).A=\begin{pmatrix}-1&1&0\\0&-1&1\\0&0&-1\end{pmatrix},\qquad X'=AX,\qquad X(0)=\begin{pmatrix}p\\q\\r\end{pmatrix}. Use the Euclidean vector norm. For N=A+IN=A+I, the operator bounds ∥NX∥≤∥X∥\|NX\|\le\|X\| and ∥N2X∥≤∥X∥\|N^2X\|\le\|X\| may be used.

Tasks

  1. Compute the powers of NN needed to find the normalized evolution and the complete solution.

  2. For initial (0,0,1)T(0,0,1)^T, find the maximum of each of the first two components on t≥0t\ge 0 and its time.

  3. For arbitrary nonzero initial data, identify when the largest polynomial factor in the solution is quadratic, linear, or constant. Do all solutions decay?

  4. Can one constant CC satisfy ∥X(t)∥≤Ce−t∥X(0)∥\|X(t)\|\le C e^{-t}\|X(0)\| for all t≥0t\ge 0 and all initial data? Prove instead a bound with decay e−t/2e^{-t/2} and an explicit constant.

Original worksheet page 1: question and worked solution for 5-9-006
Show solutionHide solution

Question 6 – Solution

Strategy. A longer nilpotent chain truncates the exponential after a quadratic term; its degree matters for sharp decay estimates.

Step 1: Truncate the nilpotent exponential. Here N2N^2 has only entry (1,3)(1,3) equal to one; N3=0N^3=0. Direct differentiation verifies F(t)=e−t(I+tN+t2N2/2)F(t)=e^{-t}(I+tN+t^2N^2/2) and F(0)=IF(0)=I. Consequently X=e−t(p+qt+rt2/2q+rtr).\boxed{X=e^{-t}\begin{pmatrix}p+qt+rt^2/2\\q+rt\\r\end{pmatrix}.} This gives all initial states and shows why a length-three chain needs t2/2t^2/2, rather than another unrelated eigenvector.

Step 2: Locate the component peaks. For (p,q,r)=(0,0,1)(p,q,r)=(0,0,1), the first components are t2e−t/2t^2e^{-t}/2 and te−tte^{-t}. Their derivatives are e−tt(2−t)/2e^{-t}t(2-t)/2 and e−t(1−t)e^{-t}(1-t). Thus their maxima are 2e−2 at t=2\boxed{2e^{-2}\text{ at }t=2} and e−1 at t=1\boxed{e^{-1}\text{ at }t=1}, respectively. These peaks are not simultaneous.

Step 3: Classify polynomial degrees. If r≠0r\ne 0, the leading vector factor is (r/2)t2(1,0,0)T(r/2)t^2(1,0,0)^T. If r=0r=0 but q≠0q\ne 0, it is qt(1,0,0)Tqt(1,0,0)^T. If q=r=0q=r=0 and p≠0p\ne 0, it is the constant p(1,0,0)Tp(1,0,0)^T. Every polynomial factor times e−te^{-t} tends to zero, so all solutions decay; algebraic multiplicity three does not force degree two for every state.

Step 4: Prove an honest exponential estimate. For initial (0,0,1)(0,0,1), et∥X(t)∥≥t2/2e^t\|X(t)\|\ge t^2/2, ruling out the proposed constant CC. The given operator bounds instead imply ∥X(t)∥≤e−t(1+t+t2/2)∥X(0)∥\|X(t)\|\le e^{-t}(1+t+t^2/2)\|X(0)\|. Using e−t/2≤1e^{-t/2}\le 1, te−t/2≤2/ete^{-t/2}\le 2/e and (t2/2)e−t/2≤8/e2(t^2/2)e^{-t/2}\le 8/e^2 for t≥0t\ge 0 gives ∥X(t)∥≤(1+2/e+8/e2)e−t/2∥X(0)∥.\boxed{\|X(t)\|\le(1+2/e+8/e^2)e^{-t/2}\|X(0)\|.} A slightly slower exponential absorbs the polynomial transient.

Original worksheet page 2: question and worked solution for 5-9-006

Original worksheet layout. Use Enlarge or open the PDF for a closer view.