21.04.2015 Views

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

390 LESSON 34. CRITICAL POINTS<br />

S<strong>in</strong>ce the eigenvalues are purely imag<strong>in</strong>ary, the orig<strong>in</strong> is a center.<br />

If we perturb either of the diagonal coefficients, the eigenvalues develop a<br />

real part. For example, the system<br />

}<br />

x ′ = 1.01x + 2y<br />

y ′ (34.95)<br />

= −3x − y<br />

has eigenvalues λ ≈ 0.005 ± 2.23383i, mak<strong>in</strong>g it an unstable spiral; and the<br />

system<br />

}<br />

x ′ = 0.99x + 2y<br />

y ′ (34.96)<br />

= −3x − y<br />

has eigenvalues λ ≈ −0.005 ± 2.2383i, for a stable spiral.<br />

The magnitude of the real part grows approximately l<strong>in</strong>early as the perturbation<br />

grows. In general, the perturbed system<br />

}<br />

x ′ = (1 + ε)x + 2y<br />

y ′ (34.97)<br />

= −3x − y<br />

will have eigenvalues<br />

λ = 1 2<br />

[ε ± √ −20 + 4ε + ε 2 ]<br />

= ε 2 ± i√ 5<br />

√<br />

1 − ε 5 − ε2<br />

20<br />

The results for perturbations of ε = ±0.05 are shown <strong>in</strong> figure 34.8.<br />

(34.98)<br />

(34.99)<br />

Non-isolated Critical Po<strong>in</strong>ts<br />

When the matrix A <strong>in</strong> equation (34.36) is s<strong>in</strong>gular, the critical po<strong>in</strong>ts will<br />

not, <strong>in</strong> general, be isolated, and they will not fall <strong>in</strong>to any of the categories<br />

that we have discussed so far.<br />

Set ∆ = ad−bc = 0, and suppose that all four of a, b, c, and d are nonzero.<br />

Then we can solve for any one of the four coefficients <strong>in</strong> terms of the others,<br />

e.g., d = bc/a, so that the trajectories are def<strong>in</strong>ed by<br />

dy cx + dy cx + (bc/a)y<br />

= = = c dx ax + by ax + by a = d b<br />

(34.100)<br />

The trajectories are parallel l<strong>in</strong>es with slope c/a. There is only one critical

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!