DifferentialEquations
Chapter 7 · Section 7.3

7.3 Basic Theory of Linear Systems in Normal Form (2×2)

The operator method of Section 7.1 produces answers, but it says little about why the answers have the form they do, or when a solution exists at all. In Chapter 4 we answered such questions for a single linear equation: existence and uniqueness, superposition, the Wronskian, fundamental sets. This section does the same for systems. We work with two unknowns, where everything can be written out in full, but almost every statement carries over word for word to unknowns (Sections 7.7 and 7.8).

The theory is cleanest for systems in which each equation is solved for one derivative. That is the normal form, and it is the form used throughout the rest of the chapter, in numerical methods (Section 8.8) and in engineering software.

Systems in normal form

Definition 7.3.1 (Linear system in normal form)

A linear system of two first-order equations in normal form is

where the coefficients and the forcing terms are given functions on an interval . If and on the system is homogeneous; otherwise it is nonhomogeneous. A solution on is a pair of differentiable functions , satisfying both equations at every .

We will often write a solution as a column, , or as an ordered pair . Think of it as a point that moves in the -plane as time goes on. Section 7.4 makes much use of that picture.

For example, , is homogeneous with constant coefficients, while , is nonhomogeneous with variable coefficients.

Every higher-order linear equation is a system

Normal form is more general than it looks. Take a second-order linear equation written with leading coefficient ,

Introduce the position and the velocity as two new unknowns: and . Then by definition, and the equation itself says . So

a system in normal form. Conversely, the first component of any solution of this system solves the original equation. An th-order equation becomes a system of equations in the same way (, , …, ), and two second-order equations, such as those of the coupled springs in Section 7.2, become four first-order equations. This is how computers solve higher-order equations (Section 8.8).

Example 7.3.1

Write the initial-value problem as a system in normal form, and state the initial conditions for the new unknowns.

Solution
  1. Normalize. For divide by :
  2. New unknowns. Put , . Then Here , , , , , .
  3. Initial conditions. and .
  4. Remark. The coefficients are continuous everywhere except at . We will see in a moment that this guarantees a unique solution on , the largest interval containing on which they are continuous.

Existence and uniqueness

For a system we prescribe the value of each unknown at the initial time: two numbers for two unknowns.

Theorem 7.3.1 (Existence and uniqueness)

Suppose are continuous on an interval . Let and let be any two numbers. Then the system (7.3.1) has exactly one solution , with and this solution exists on the whole interval .

Three remarks.

  • As in Theorem 4.1.1, a linear system cannot "blow up" inside : the solution lives on the entire interval of continuity. Nonlinear systems do not have this property.
  • When all coefficients are constants, .
  • Uniqueness has a geometric meaning for homogeneous systems with constant coefficients: two different solution curves in the -plane can never cross, because at a crossing point both would satisfy the same initial condition. Figure 7.3.2 illustrates this.

A useful special case: the only solution of a homogeneous system with , is the trivial solution , , since the zero pair obviously solves it and there is only one solution.

Homogeneous systems

From now on until the last part of the section we consider the homogeneous system

with coefficients continuous on .

Theorem 7.3.2 (Superposition)

If and are solutions of (7.3.2) on , then for any constants the pair is also a solution on .

The proof is one line per equation. For the first one,

and the second is identical. Note that the same constant multiplies both components of a solution.

Linear independence of solutions

Superposition builds new solutions from old ones. To build all solutions we need two that are "really different".

Definition 7.3.2 (Linear dependence and independence)

Two pairs of functions and are linearly dependent on if there are constants , not both zero, such that Otherwise they are linearly independent on .

In words: two pairs are dependent exactly when one is a constant multiple of the other, with the same constant for both components. For instance, and are dependent (the second is times the first), while and are independent even though their first components coincide.

The Wronskian

Definition 7.3.3 (Wronskian of two solution pairs)

The Wronskian of the pairs and is the determinant whose columns are the two pairs:

Compare this with Chapter 4: there the second row contained derivatives . For the system obtained from by , , the second components are the derivatives, so the two Wronskians coincide.

How does change in time? Differentiate the determinant row by row (product rule):

Now substitute the system. In the first bracket use :

In the second bracket use ; the same cancellation leaves . Therefore

a first-order linear equation for (Section 2.3). Solving it gives the following formula.

Theorem 7.3.3 (Abel's formula)

If and are solutions of (7.3.2) on and , then In particular, is either zero everywhere on or zero nowhere on .

The sum of the diagonal coefficients is called the trace of the system. For constant coefficients Abel's formula reads simply .

Theorem 7.3.4 (Wronskian test)

Two solutions of the homogeneous system (7.3.2) are linearly independent on if and only if their Wronskian is different from zero at one (and then, by Abel's formula, at every) point of .

Why. If the solutions are dependent, one is a multiple of the other, so the columns of are proportional and . Conversely, suppose at some . Then the algebraic system , has a solution , because its determinant is zero. With these constants, , is a solution (superposition) that vanishes at . By uniqueness it is the trivial solution, so and : the solutions are dependent.

The test is for solutions only

For pairs of functions that are not solutions of a common system with continuous coefficients, a Wronskian that vanishes (even identically) does not imply dependence. Example 7.3.5 shows two independent pairs with .

Fundamental sets and the general solution

Definition 7.3.4 (Fundamental set)

Two linearly independent solutions of the homogeneous system (7.3.2) on form a fundamental set of solutions on .

Theorem 7.3.5 (General solution of the homogeneous system)

The homogeneous system (7.3.2) with coefficients continuous on has a fundamental set on . If , is a fundamental set, then every solution on can be written as for suitable constants . This expression is called the general solution.

Why. Existence: pick and let be the solution with initial values and the solution with initial values (Theorem 7.3.1). Their Wronskian at is . Every solution is a combination: given a solution , solve for ; this is possible because the determinant is . Then and are two solutions with the same initial values, so by uniqueness they are equal.

The argument is also the practical recipe for initial-value problems: once a fundamental set is known, an IVP is just a linear system for and .

Nonhomogeneous systems

Theorem 7.3.6 (General solution of the nonhomogeneous system)

Let be any one solution of the nonhomogeneous system (7.3.1) on , and let , be a fundamental set of the corresponding homogeneous system (7.3.2). Then every solution of (7.3.1) on has the form

The reason: if is any solution of (7.3.1), the difference solves the homogeneous system, because the forcing terms cancel when we subtract. By Theorem 7.3.5 that difference is . This is the same structure "complementary + particular" that you know from Chapter 4.

Method in a nutshell: using the theory
  1. Check that the coefficients and forcing are continuous on an interval containing ; the solution then exists and is unique on all of .
  2. Find (or verify) two solutions of the homogeneous system. Compute their Wronskian at one convenient point. If it is not zero, they form a fundamental set. Abel's formula is a free check of your algebra.
  3. For a nonhomogeneous system add one particular solution.
  4. For an IVP substitute and solve the linear system for .

How do we find the two homogeneous solutions? For constant coefficients this is the subject of Section 7.4. You can already experiment with the solver below: enter the coefficient matrix row by row () and an initial point, and it returns a fundamental set, the solution of the IVP and its phase-plane picture. The preset is the system of Example 7.3.2.

Example 7.3.2

Consider , .

  1. Verify that and are solutions.
  2. Show that they form a fundamental set and check Abel's formula.
  3. Write the general solution and solve the IVP , .
Solution
  1. First pair. give and ; and . Both equations hold.
  2. Second pair. , : and ; and . Both hold.
  3. Wronskian. So the solutions are independent on and form a fundamental set.
  4. Abel's check. The trace is , so Abel predicts . It matches.
  5. General solution.
  6. IVP. At : and . Subtracting, , so and : Figure 7.3.1 shows the two components; Figure 7.3.2 shows the same solution as a curve in the -plane.
Figure 7.3.1. The solution of Example 7.3.2 as two functions of time. Salt-in-two-tanks behaviour again: starts at zero, is "fed" by , peaks, and then both decay.
Figure 7.3.2. The same system drawn in the -plane. Each arrow shows the direction at that point; each curve is one solution, traced as increases. The highlighted curve starts at (Example 7.3.2). By uniqueness, no two curves cross.
Example 7.3.3 (Variable coefficients)

Show that and form a fundamental set of and solve the IVP , . On which interval is the solution guaranteed?

Solution
  1. First pair. , : and . ✓. and . ✓
  2. Second pair. , : and . ✓. and . ✓
  3. Wronskian. So the pairs form a fundamental set on .
  4. Abel's check. The trace is , so . ✓
  5. IVP. , . At : and , so , :
  6. Interval. The coefficients , , are continuous on , which contains . Theorem 7.3.1 guarantees exactly one solution there. (The formulas happen to make sense for all , but the theorem promises nothing across , where the coefficients are undefined.)
Example 7.3.4 (A nonhomogeneous system)

Show that , is a solution of and find the solution with , .

Solution
  1. Check the particular solution. and . ✓. and . ✓
  2. Homogeneous part. The homogeneous system is the one of Example 7.3.2, so by Theorem 7.3.6
  3. Initial conditions. and . Subtracting, : , .
  4. Answer. , . For large the solution follows the particular solution .
Example 7.3.5 (A vanishing Wronskian that proves nothing)

Let and . Show that and are linearly independent on although their Wronskian is identically zero. What does this say about systems they could solve?

Solution
  1. Wronskian. for every .
  2. Independence. Suppose and for all . Put in the second: . Then for all , so . Hence the pairs are independent.
  3. Conclusion. If and were solutions of one homogeneous system (7.3.2) with coefficients continuous on an interval , Theorem 7.3.4 would force on . So no such system exists on any interval. The Wronskian test works for solutions, not for arbitrary pairs of functions.

Exercises

  1. Consider , . (a) Verify that and are solutions. (b) Compute their Wronskian and check it with Abel's formula. (c) Solve the IVP , .
    Answer (b) , trace . (c) , .
  2. Write as a system in normal form.
    Answer , (with , ).
  3. Write as a system of three first-order equations in normal form.
    Answer , , .
  4. Find the largest interval on which Theorem 7.3.1 guarantees a unique solution of , , , .
    Answer : needs , and is discontinuous at .
  5. Without solving, find the Wronskian of any two solutions of , , up to a constant factor.
    Answer .
  6. Without solving, find the Wronskian of any two solutions of , on , up to a constant factor.
    Answer .
  7. Verify that and are solutions of , on . Compute , check Abel's formula, and solve the IVP , .
    Answer (trace ). , .
  8. Find a constant solution of , and then the general solution (use Exercise 1).
    Answer . , .
  9. Can two solutions of a homogeneous system with coefficients continuous on have Wronskian ? Explain.
    Answer No. By Abel's formula is either never zero or identically zero, but vanishes only at .
  10. Show that the pairs and have , and find constants, not both zero, that show they are dependent.
    Answer ; .
  11. Let and be solutions of the nonhomogeneous system (7.3.1). Show that solves the homogeneous system, but in general does not solve (7.3.1).
    Answer Subtracting, the forcing terms cancel. Adding, the forcing appears twice: the sum solves the system with forcing .
  12. For the system of Exercise 1, find the solutions with initial values and at , and check that their Wronskian at is .
    Answer and ; , so .