Math 4 B Systems of Differential Equations Phase

Math 4 B Systems of Differential Equations Phase Plane Analysis Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

A 2 x 2 system of 1 st order linear differential equations will have the form The solution vector x(t) can be visualized by a phase plane diagram. At each point on a solution curve, the tangent is given by x’. The phase plane diagram looks like a vector field, where at each point the slope of the vector is just x 2/x 1 (i. e. rise/run). Let’s look at an example first, then get a more general description. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Unstable Node Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Here are the possibilities for equilibrium solutions to a 2 x 2 system. The outcome depends on the values of the eigenvalues. Real Eigenvalues r 1 and r 2: If both are negative the result is an asymptotically stable node. If both are positive the result is an unstable node. If they have opposite signs, there is an unstable saddle point. Complex Eigenvalues λ±iμ: If λ=0, a stable center point. If λ<0 the result is a stable spiral. If λ>0 the result is an unstable spiral. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Unstable Saddle Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Unstable Spiral Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Let’s work through the solution procedure for this one as well. Find eigenvectors for the eigenvalues: Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Let’s work through the solution procedure for this one as well. Find eigenvectors for the eigenvalues: Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Let’s work through the solution procedure for this one as well. Find eigenvectors for the eigenvalues: Here is the corresponding fundamental matrix: Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. The general solution can now be written in terms of this fundamental matrix: Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. The general solution can now be written in terms of this fundamental matrix: To solve the initial value problem, we can incorporate the initial value by writing the solution this way: Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. The general solution can now be written in terms of this fundamental matrix: To solve the initial value problem, we can incorporate the initial value by writing the solution this way: By coincidence this matrix is already the identity, so we don’t need to find the inverse. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Notice that with this formulation, we can quickly change the initial value to obtain a new solution, without recalculating the fundamental matrix. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Draw the Phase Plane diagram, and sketch the solution curve. Stable Center Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Find the solution to this initial value problem. Find eigenvectors for the eigenvalues: Eigenvector is split into real and imaginary parts. Solutions can be written in terms of complex exponentials, or sine/cosines. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Find the solution to this initial value problem. Here is the eigenvector for r=2 i. The other one is the conjugate. Eigenvector is split into real and imaginary parts. Solutions can be written in terms of complex exponentials, or sine/cosines. After some rearranging two independent solutions are: Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Find the solution to this initial value problem. Solution can be written as a Fundamental Matrix: The columns of this fundamental matrix X(t) are the solutions x(1) and x(2). This is the general solution. The next step is to incorporate the initial value. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Find the solution to this initial value problem. Here is the format for the solution with the initial value: The matrix Y(t) is another version of the fundamental matrix that also has the property that Y(0)=I. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Example: Find the solution to this initial value problem. Multiplying matrices, we find that Y(t) is: Now we can multiply this by the initial value vector to get the solution. this is x 1(t) this is x 2(t) Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB