Skip to main content

Calculus III - 3-Dimensional Space: Equations of Lines

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

Differential Equations - Second Order: Basic Concepts


In this chapter we will be looking exclusively at linear second order differential equations. The most general linear second order differential equation is in the form.
(1)p(t)y+q(t)y+r(t)y=g(t)
In fact, we will rarely look at non-constant coefficient linear second order differential equations. In the case where we assume constant coefficients we will use the following differential equation.
(2)ay+by+cy=g(t)
Where possible we will use (1) just to make the point that certain facts, theorems, properties, and/or techniques can be used with the non-constant form. However, most of the time we will be using (2) as it can be fairly difficult to solve second order non-constant coefficient differential equations.

Initially we will make our life easier by looking at differential equations with g(t)=0. When g(t)=0 we call the differential equation homogeneous and when g(t)0 we call the differential equation nonhomogeneous.

So, let’s start thinking about how to go about solving a constant coefficient, homogeneous, linear, second order differential equation. Here is the general constant coefficient, homogeneous, linear, second order differential equation.
ay+by+cy=0
It’s probably best to start off with an example. This example will lead us to a very important fact that we will use in every problem from this point on. The example will also give us clues into how to go about solving these in general.


Example 1 Determine some solutions toy9y=0

We can get some solutions here simply by inspection. We need functions whose second derivative is 9 times the original function. One of the first functions that I can think of that comes back to itself after two derivatives is an exponential function and with proper exponents the 9 will get taken care of as well.
So, it looks like the following two functions are solutions.
y(t)=e3tandy(t)=e3tWe’ll leave it to you to verify that these are in fact solutions.
These two functions are not the only solutions to the differential equation however. Any of the following are also solutions to the differential equation.
y(t)=9e3ty(t)=123e3ty(t)=56e3ty(t)=149e3ty(t)=7e3t6e3ty(t)=92e3t16e3tIn fact if you think about it any function that is in the form
y(t)=c1e3t+c2e3twill be a solution to the differential equation.

This example leads us to a very important fact that we will use in practically every problem in this chapter.


Principle of Superposition


If y1(t) and y2(t) are two solutions to a linear, homogeneous differential equation then so is
(3)y(t)=c1y1(t)+c2y2(t)

Note that we didn’t include the restriction of constant coefficient or second order in this. This will work for any linear homogeneous differential equation.
If we further assume second order and one other condition (which we’ll give in a second) we can go a step further.
If y1(t) and y2(t) are two solutions to a linear, second order homogeneous differential equation and they are “nice enough” then the general solution to the linear, second order homogeneous differential equation is given by (3).
So, just what do we mean by “nice enough”? We’ll hold off on that until a later section. At this point you’ll hopefully believe it when we say that specific functions are “nice enough”.
So, if we now make the assumption that we are dealing with a linear, second order homogeneous differential equation, we now know that (3) will be its general solution. The next question that we can ask is how to find the constants c1 and c2. Since we have two constants it makes sense, hopefully, that we will need two equations, or conditions, to find them.
One way to do this is to specify the value of the solution at two distinct points, or,
y(t0)=y0y(t1)=y1
These are typically called boundary values and are not really the focus of this course so we won’t be working with them here. We do give a brief introduction to boundary values in a later chapter if you are interested in seeing how they work and some of the issues that arise when working with boundary values.
Another way to find the constants would be to specify the value of the solution and its derivative at a particular point. Or,
y(t0)=y0y(t0)=y0
These are the two conditions that we’ll be using here. As with the first order differential equations these will be called initial conditions.


Example 2 Solve the following IVP.y9y=0y(0)=2y(0)=1

First, the two functions
y(t)=e3tandy(t)=e3tare “nice enough” for us to form the general solution to the differential equation. At this point, please just believe this. You will be able to verify this for yourself in a couple of sections.
The general solution to our differential equation is then
y(t)=c1e3t+c2e3tNow all we need to do is apply the initial conditions. This means that we need the derivative of the solution.
y(t)=3c1e3t+3c2e3tPlug in the initial conditions
2=y(0)=c1+c21=y(0)=3c1+3c2This gives us a system of two equations and two unknowns that can be solved. Doing this yields
c1=76c2=56The solution to the IVP is then,
y(t)=76e3t+56e3t

Up to this point we’ve only looked at a single differential equation and we got its solution by inspection. For a rare few differential equations we can do this. However, for the vast majority of the second order differential equations out there we will be unable to do this.
So, we would like a method for arriving at the two solutions we will need in order to form a general solution that will work for any linear, constant coefficient, second order homogeneous differential equation. This is easier than it might initially look.
We will use the solutions we found in the first example as a guide. All of the solutions in this example were in the form
y(t)=ert
Note, that we didn’t include a constant in front of it since we can literally include any constant that we want and still get a solution. The important idea here is to get the exponential function. Once we have that we can add on constants to our hearts content.
So, let’s assume that all solutions to
(4)ay+by+cy=0(5)y(t)=ert
To see if we are correct all we need to do is plug this into the differential equation and see what happens. So, let’s get some derivatives and then plug in.
y(t)=rerty(t)=r2erta(r2ert)+b(rert)+c(ert)=0ert(ar2+br+c)=0
So, if (5) is to be a solution to (4) then the following must be true
ert(ar2+br+c)=0
This can be reduced further by noting that exponentials are never zero. Therefore, (5) will be a solution to (4)provided r is a solution to
(6)ar2+br+c=0
This equation is typically called the characteristic equation for (4).
Okay, so how do we use this to find solutions to a linear, constant coefficient, second order homogeneous differential equation? First write down the characteristic equation, (6), for the differential equation, (4). This will be a quadratic equation and so we should expect two roots, r1 and r2. Once we have these two roots we have two solutions to the differential equation.
(7)y1(t)=er1tandy2(t)=er2t
Let’s take a look at a quick example.


Example 3 Find two solutions toy9y=0

This is the same differential equation that we looked at in the first example. This time however, let’s not just guess. Let’s go through the process as outlined above to see the functions that we guess above are the same as the functions the process gives us.
First write down the characteristic equation for this differential equation and solve it.
r29=0r=±3The two roots are 3 and -3. Therefore, two solutions are
y1(t)=e3tandy2(t)=e3tThese match up with the first guesses that we made in the first example.

You’ll notice that we neglected to mention whether or not the two solutions listed in (7) are in fact “nice enough” to form the general solution to (4). This was intentional. We have three cases that we need to look at and this will be addressed differently in each of these cases.
So, what are the cases? As we previously noted the characteristic equation is quadratic and so will have two roots, r1 and r2. The roots will have three possible forms. These are
  1. Real, distinct roots,r1r2.
  2. Complex root,r1,2=λ±μi.
  3. Double roots,r1=r2=r.
The next three sections will look at each of these in some more depth, including giving forms for the solution that will be “nice enough” to get a general solution.

Comments

Popular posts from this blog

Digital Signal Processing - Basic Continuous Time Signals

To test a system, generally, standard or basic signals are used. These signals are the basic building blocks for many complex signals. Hence, they play a very important role in the study of signals and systems. Unit Impulse or Delta Function A signal, which satisfies the condition,   δ ( t ) = lim ϵ → ∞ x ( t ) δ ( t ) = lim ϵ → ∞ x ( t )   is known as unit impulse signal. This signal tends to infinity when t = 0 and tends to zero when t ≠ 0 such that the area under its curve is always equals to one. The delta function has zero amplitude everywhere except at t = 0. Properties of Unit Impulse Signal δ(t) is an even signal. δ(t) is an example of neither energy nor power (NENP) signal. Area of unit impulse signal can be written as; A = ∫ ∞ − ∞ δ ( t ) d t = ∫ ∞ − ∞ lim ϵ → 0 x ( t ) d t = lim ϵ → 0 ∫ ∞ − ∞ [ x ( t ) d t ] = 1 Weight or strength of the signal can be written as; y ( t ) = A δ ( t ) y ( t ) = A δ ( t ) Area of the weighted impulse signal can

Differential Equations - First Order: Bernoulli

In this section we are going to take a look at differential equations in the form, y ′ + p ( x ) y = q ( x ) y n y ′ + p ( x ) y = q ( x ) y n where  p ( x ) p ( x )  and  q ( x ) q ( x )  are continuous functions on the interval we’re working on and  n n  is a real number. Differential equations in this form are called  Bernoulli Equations . First notice that if  n = 0 n = 0  or  n = 1 n = 1  then the equation is linear and we already know how to solve it in these cases. Therefore, in this section we’re going to be looking at solutions for values of  n n  other than these two. In order to solve these we’ll first divide the differential equation by  y n y n  to get, y − n y ′ + p ( x ) y 1 − n = q ( x ) y − n y ′ + p ( x ) y 1 − n = q ( x ) We are now going to use the substitution  v = y 1 − n v = y 1 − n  to convert this into a differential equation in terms of  v v . As we’ll see this will lead to a differential equation that we can solve. We are going to have to be c

Differential Equations - Systems: Solutions

Now that we’ve got some of the basics out of the way for systems of differential equations it’s time to start thinking about how to solve a system of differential equations. We will start with the homogeneous system written in matrix form, → x ′ = A → x (1) (1) x → ′ = A x → where,  A A  is an  n × n n × n  matrix and  → x x →  is a vector whose components are the unknown functions in the system. Now, if we start with  n = 1 n = 1 then the system reduces to a fairly simple linear (or separable) first order differential equation. x ′ = a x x ′ = a x and this has the following solution, x ( t ) = c e a t x ( t ) = c e a t So, let’s use this as a guide and for a general  n n  let’s see if → x ( t ) = → η e r t (2) (2) x → ( t ) = η → e r t will be a solution. Note that the only real difference here is that we let the constant in front of the exponential be a vector. All we need to do then is plug this into the differential equation and see what we get. First notice that

Calculus III - 3-Dimensional Space: Equations of Lines

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

Digital Signal Processing - Miscellaneous Signals

There are other signals, which are a result of operation performed on them. Some common type of signals are discussed below. Conjugate Signals Signals, which satisfies the condition  x ( t ) = x ∗ ( − t ) are called conjugate signals. Let  x ( t ) = a ( t ) + j b ( t ) So,  x ( − t ) = a ( − t ) + j b ( − t ) And  x ∗ ( − t ) = a ( − t ) − j b ( − t ) By Condition,  x ( t ) = x ∗ ( − t ) If we compare both the derived equations 1 and 2, we can see that the real part is even, whereas the imaginary part is odd. This is the condition for a signal to be a conjugate type. Conjugate Anti-Symmetric Signals Signals, which satisfy the condition  x ( t ) = − x ∗ ( − t ) are called conjugate anti-symmetric signal Let  x ( t ) = a ( t ) + j b ( t ) So  x ( − t ) = a ( − t ) + j b ( − t ) And  x ∗ ( − t ) = a ( − t ) − j b ( − t ) − x ∗ ( − t ) = − a ( − t ) + j b ( − t ) By Condition  x ( t ) = − x ∗ ( − t ) Now, again compare, both the equations just as w

Differential Equations - First Order: Modeling - i

We now move into one of the main applications of differential equations both in this class and in general. Modeling is the process of writing a differential equation to describe a physical situation. Almost all of the differential equations that you will use in your job (for the engineers out there in the audience) are there because somebody, at some time, modeled a situation to come up with the differential equation that you are using. This section is not intended to completely teach you how to go about modeling all physical situations. A whole course could be devoted to the subject of modeling and still not cover everything! This section is designed to introduce you to the process of modeling and show you what is involved in modeling. We will look at three different situations in this section : Mixing Problems, Population Problems, and Falling Objects. In all of these situations we will be forced to make assumptions that do not accurately depict reality in most cases, but wi

Differential Equations - Basic Concepts: Definitions

Differential Equation The first definition that we should cover should be that of  differential equation . A differential equation is any equation which contains derivatives, either ordinary derivatives or partial derivatives. There is one differential equation that everybody probably knows, that is Newton’s Second Law of Motion. If an object of mass  m m  is moving with acceleration  a a  and being acted on with force  F F  then Newton’s Second Law tells us. F = m a (1) (1) F = m a To see that this is in fact a differential equation we need to rewrite it a little. First, remember that we can rewrite the acceleration,  a a , in one of two ways. a = d v d t OR a = d 2 u d t 2 (2) (2) a = d v d t OR a = d 2 u d t 2 Where  v v  is the velocity of the object and  u u  is the position function of the object at any time  t t . We should also remember at this point that the force,  F F  may also be a function of time, velocity, and/or position. So, with all these things in

Differential Equations - Partial: Summary of Separation of Variables

Throughout this chapter we’ve been talking about and solving partial differential equations using the method of separation of variables. However, the one thing that we’ve not really done is completely work an example from start to finish showing each and every step. Each partial differential equation that we solved made use somewhere of the fact that we’d done at least part of the problem in another section and so it makes some sense to have a quick summary of the method here. Also note that each of the partial differential equations only involved two variables. The method can often be extended out to more than two variables, but the work in those problems can be quite involved and so we didn’t cover any of that here. So with all of that out of the way here is a quick summary of the method of separation of variables for partial differential equations in two variables. Verify that the partial differential equation is linear and homogeneous. Verify that the boundary condi

Differential Equations - Systems: Repeated Eigenvalues - i

This is the final case that we need to take a look at. In this section we are going to look at solutions to the system, → x ′ = A → x x → ′ = A x → where the eigenvalues are repeated eigenvalues. Since we are going to be working with systems in which  A A  is a  2 × 2 2 × 2  matrix we will make that assumption from the start. So, the system will have a double eigenvalue,  λ λ . This presents us with a problem. We want two linearly independent solutions so that we can form a general solution. However, with a double eigenvalue we will have only one, → x 1 = → η e λ t x → 1 = η → e λ t So, we need to come up with a second solution. Recall that when we looked at the double root case with the second order differential equations we ran into a similar problem. In that section we simply added a  t t  to the solution and were able to get a second solution. Let’s see if the same thing will work in this case as well. We’ll see if → x = t e λ t → η x → = t e λ t η → will also be a

Differential Equations - First Order: Modeling - ii

Example 4  A 50 kg object is shot from a cannon straight up with an initial velocity of 10m/s off a bridge that is 100 meters above the ground. If air resistance is given by 5 v v  determine the velocity of the mass when it hits the ground. First, notice that when we say straight up, we really mean straight up, but in such a way that it will miss the bridge on the way back down. Here is a sketch of the situation. Notice the conventions that we set up for this problem. Since the vast majority of the motion will be in the downward direction we decided to assume that everything acting in the downward direction should be positive. Note that we also defined the “zero position” as the bridge, which makes the ground have a “position” of 100. Okay, if you think about it we actually have two situations here. The initial phase in which the mass is rising in the air and the second phase when the mass is on its way down. We will need to examine both situations and set up an IVP for