Real Eigenvalues — Question 1

PDF ↗

Question 1

Let X=(x,y)TX=(x,y)^T solve X′=(−211−2)X,X(0)=(20).X'=\begin{pmatrix}-2&1\\1&-2\end{pmatrix}X,\qquad X(0)=\binom 20. For geometric interpretation use modal coordinates u=(x+y)/2u=(x+y)/2 and v=(x−y)/2v=(x-y)/2. A stable node means that all nearby states tend to the origin forward in time, with two real negative eigenvalues.

Tasks

  1. Find both eigenpairs and solve the IVP. Explain why the two modes give every solution.

  2. Eliminate time for this IVP in the (u,v)(u,v) coordinates, retaining the domain. Translate the result to the original phase plane.

  3. Determine the time direction and the limiting tangent line at the origin. Find lim⁡t→∞y(t)/x(t)\lim_{t\to\infty}y(t)/x(t).

  4. For arbitrary X(0)=(p,q)TX(0)=(p,q)^T, identify precisely the initial states that approach along the faster eigendirection instead. Classify the origin.

Original worksheet page 1: question and worked solution for 5-7-001
Show solutionHide solution

Question 1 – Solution

Strategy. Decompose into independent eigendirections, then distinguish the dominant slow mode from initial data that cancel it.

Step 1: Resolve the modes. The characteristic polynomial is (λ+1)(λ+3)(\lambda+1)(\lambda+3). Eigenvectors are v1=(1,1)Tv_1=(1,1)^T for −1-1 and v2=(1,−1)Tv_2=(1,-1)^T for −3-3. They form a basis, so modal coordinates satisfy u′=−uu'=-u, v′=−3vv'=-3v. Initially u=v=1u=v=1, giving X=e−t(11)+e−3t(1−1).\boxed{X=e^{-t}\binom 11+e^{-3t}\binom 1{-1}.} This decoupling also proves completeness for arbitrary initial data.

Step 2: Retain the orbit domain. Here u=e−t>0u=e^{-t}>0 and v=e−3t=u3v=e^{-3t}=u^3 for all real tt. Thus the orbit is x−y2=(x+y2)3,x+y>0.\frac{x-y}{2}=\left(\frac{x+y}{2}\right)^3,\qquad x+y>0. The origin is a limiting point, not a point attained at finite time.

Step 3: Read the approach direction. As time increases, uu decreases toward zero. Moreover y/x=(1−e−2t)/(1+e−2t)→1y/x=(1-e^{-2t})/(1+e^{-2t})\to 1. Both the state direction and the velocity direction approach the line y=xy=x; indeed dy/dx=(1−3e−2t)/(1+3e−2t)→1dy/dx=(1-3e^{-2t})/(1+3e^{-2t})\to 1. At (2,0)(2,0) the velocity is (−4,2)(-4,2), consistent with the arrows.

Step 4: Keep the exceptional line. For initial (p,q)(p,q), the slow and fast coefficients are (p+q)/2(p+q)/2 and (p−q)/2(p-q)/2. Every solution tends to zero. A nonzero solution approaches along the fast line y=−xy=-x exactly when p+q=0p+q=0; otherwise its limiting line is y=xy=x. The zero state remains fixed and has no approach direction. Both eigenvalues are negative, so the origin is a stable node.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 5-7-001

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