In this section we need to take a look at the equation of a line in R 3 R 3 . As we saw in the previous section the equation y = m x + b y = m x + b does not describe a line in R 3 R 3 , instead it describes a plane. This doesn’t mean however that we can’t write down an equation for a line in 3-D space. We’re just going to need a new way of writing down the equation of a curve. So, before we get into the equations of lines we first need to briefly look at vector functions. We’re going to take a more in depth look at vector functions later. At this point all that we need to worry about is notational issues and how they can be used to give the equation of a curve. The best way to get an idea of what a vector function is and what its graph looks like is to look at an example. So, consider the following vector function. → r ( t ) = ⟨ t , 1 ⟩ r → ( t ) = ⟨ t , 1 ⟩ A vector function is a function that takes one or more variables, one in this case, and returns a vector. Note as we
Example 3 Find the solution to the following system.
We first need to convert this into matrix form. This is easy enough. Here is the matrix form of the system.
This is just the system from the first example and so we’ve already got the solution to this system. Here it is.
Now, since we want the solution to the system not in matrix form let’s go one step farther here. Let’s multiply the constants and exponentials into the vectors and then add up the two vectors.
Now, recall,
So, the solution to the system is then,
This is just the system from the first example and so we’ve already got the solution to this system. Here it is.
Now, since we want the solution to the system not in matrix form let’s go one step farther here. Let’s multiply the constants and exponentials into the vectors and then add up the two vectors.
Now, recall,
So, the solution to the system is then,
Let’s work another example.
Example 4 Solve the following IVP.
So, the first thing that we need to do is find the eigenvalues for the matrix.
Now let’s find the eigenvectors for each of these.
:
We’ll need to solve,
The eigenvector in this case is,
:
We’ll need to solve,
The eigenvector in this case is,
Then general solution is then,
Now, we need to find the constants. To do this we simply need to apply the initial conditions.
Now solve the system for the constants.
The solution is then,
Now let’s find the eigenvectors for each of these.
:
We’ll need to solve,
The eigenvector in this case is,
:
We’ll need to solve,
The eigenvector in this case is,
Then general solution is then,
Now, we need to find the constants. To do this we simply need to apply the initial conditions.
Now solve the system for the constants.
The solution is then,
Now let’s find the phase portrait for this system.
Example 5 Sketch the phase portrait for the following system.
From the last example we know that the eigenvalues and eigenvectors for this system are,
This one is a little different from the first one. However, it starts in the same way. We’ll first sketch the trajectories corresponding to the eigenvectors. Notice as well that both of the eigenvalues are negative and so trajectories for these will move in towards the origin as increases. When we sketch the trajectories we’ll add in arrows to denote the direction they take as increases. Here is the sketch of these trajectories.
Now, here is where the slight difference from the first phase portrait comes up. All of the trajectories will move in towards the origin as increases since both of the eigenvalues are negative. The issue that we need to decide upon is just how they do this. This is actually easier than it might appear to be at first.
The second eigenvalue is larger than the first. For large and positive ’s this means that the solution for this eigenvalue will be smaller than the solution for the first eigenvalue. Therefore, as increases the trajectory will move in towards the origin and do so parallel to . Likewise, since the second eigenvalue is larger than the first this solution will dominate for large and negative ’s. Therefore, as we decrease the trajectory will move away from the origin and do so parallel to .
Adding in some trajectories gives the following sketch.
In these cases we call the equilibrium solution a node and it is asymptotically stable. Equilibrium solutions are asymptotically stable if all the trajectories move in towards it as increases.
This one is a little different from the first one. However, it starts in the same way. We’ll first sketch the trajectories corresponding to the eigenvectors. Notice as well that both of the eigenvalues are negative and so trajectories for these will move in towards the origin as increases. When we sketch the trajectories we’ll add in arrows to denote the direction they take as increases. Here is the sketch of these trajectories.
The second eigenvalue is larger than the first. For large and positive ’s this means that the solution for this eigenvalue will be smaller than the solution for the first eigenvalue. Therefore, as increases the trajectory will move in towards the origin and do so parallel to . Likewise, since the second eigenvalue is larger than the first this solution will dominate for large and negative ’s. Therefore, as we decrease the trajectory will move away from the origin and do so parallel to .
Adding in some trajectories gives the following sketch.
Note that nodes can also be unstable. In the last example if both of the eigenvalues had been positive all the trajectories would have moved away from the origin and in this case the equilibrium solution would have been unstable.
Before moving on to the next section we need to do one more example. When we first started talking about systems it was mentioned that we can convert a higher order differential equation into a system. We need to do an example like this so we can see how to solve higher order differential equations using systems.
Example 6 Convert the following differential equation into a system, solve the system and use this solution to get the solution to the original differential equation.
So, we first need to convert this into a system. Here’s the change of variables,
The system is then,
where,
Now we need to find the eigenvalues for the matrix.
Now let’s find the eigenvectors.
:
We’ll need to solve,
The eigenvector in this case is,
:
We’ll need to solve,
The eigenvector in this case is,
The general solution is then,
Apply the initial condition.
This gives the system of equations that we can solve for the constants.
The actual solution to the system is then,
Now recalling that,
we can see that the solution to the original differential equation is just the top row of the solution to the matrix system. The solution to the original differential equation is then,
Notice that as a check, in this case, the bottom row should be the derivative of the top row.
The system is then,
where,
Now we need to find the eigenvalues for the matrix.
Now let’s find the eigenvectors.
:
We’ll need to solve,
The eigenvector in this case is,
:
We’ll need to solve,
The eigenvector in this case is,
The general solution is then,
Apply the initial condition.
This gives the system of equations that we can solve for the constants.
The actual solution to the system is then,
Now recalling that,
we can see that the solution to the original differential equation is just the top row of the solution to the matrix system. The solution to the original differential equation is then,
Notice that as a check, in this case, the bottom row should be the derivative of the top row.
Comments
Post a Comment