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 - First Order: Intervals of Validity



We’ve called this section Intervals of Validity because all of the examples will involve them. However, there is a lot more to this section. We will see a couple of theorems that will tell us when we can solve a differential equation. We will also see some of the differences between linear and nonlinear differential equations.
First let's take a look at a theorem about linear first order differential equations. This is a very important theorem although we’re not going to really use it for its most important aspect.

Theorem 1

Consider the following IVP.
y+p(t)y=g(t)y(t0)=y0If p(t) and g(t) are continuous functions on an open interval α<t<β and the interval contains to, then there is a unique solution to the IVP on that interval.

So, just what does this theorem tell us? First, it tells us that for nice enough linear first order differential equations solutions are guaranteed to exist and more importantly the solution will be unique. We may not be able to find the solution but do know that it exists and that there will only be one of them. This is the very important aspect of this theorem. Knowing that a differential equation has a unique solution is sometimes more important than actually having the solution itself!

Next, if the interval in the theorem is the largest possible interval on which p(t) and g(t) are continuous then the interval is the interval of validity for the solution. This means, that for linear first order differential equations, we won't need to actually solve the differential equation in order to find the interval of validity. Notice as well that the interval of validity will depend only partially on the initial condition. The interval must contain to, but the value of yo, has no effect on the interval of validity.
Let’s take a look at an example.


Example 1 Without solving, determine the interval of validity for the following initial value problem.(t29)y+2y=ln|204t|y(4)=3

First, in order to use the theorem to find the interval of validity we must write the differential equation in the proper form given in the theorem. So we will need to divide out by the coefficient of the derivative.
y+2t29y=ln|204t|t29Next, we need to identify where the two functions are not continuous. This will allow us to find all possible intervals of validity for the differential equation. So, p(t) will be discontinuous at t=±3 since these points will give a division by zero. Likewise, g(t) will also be discontinuous at t=±3 as well as t=5 since at this point we will have the natural logarithm of zero. Note that in this case we won't have to worry about natural log of negative numbers because of the absolute values.
Now, with these points in hand we can break up the real number line into four intervals where both p(t) and g(t) will be continuous. These four intervals are,
<t<33<t<33<t<55<t<The endpoints of each of the intervals are points where at least one of the two functions is discontinuous. This will guarantee that both functions are continuous everywhere in each interval.
Finally, let's identify the actual interval of validity for the initial value problem. The actual interval of validity is the interval that will contain to=4. So, the interval of validity for the initial value problem is.
3<t<5

In this last example we need to be careful to not jump to the conclusion that the other three intervals cannot be intervals of validity. By changing the initial condition, in particular the value of to, we can make any of the four intervals the interval of validity.
The first theorem required a linear differential equation. There is a similar theorem for non-linear first order differential equations. This theorem is not as useful for finding intervals of validity as the first theorem was so we won’t be doing all that much with it.
Here is the theorem.

Theorem 2

Consider the following IVP.
y=f(t,y)y(t0)=y0If f(t,y) and fy are continuous functions in some rectangle α<t<βγ<y<δ containing the point (toyo) then there is a unique solution to the IVP in some interval toh<t<to+h that is contained in α<t<β.

That’s it. Unlike the first theorem, this one cannot really be used to find an interval of validity. So, we will know that a unique solution exists if the conditions of the theorem are met, but we will actually need the solution in order to determine its interval of validity. Note as well that for non-linear differential equations it appears that the value of y0 may affect the interval of validity.
Here is an example of the problems that can arise when the conditions of this theorem aren’t met.


Example 2 Determine all possible solutions to the following IVP.y=y13y(0)=0

First, notice that this differential equation does NOT satisfy the conditions of the theorem.
f(y)=y13dfdy=13y23So, the function is continuous on any interval, but the derivative is not continuous at y=0 and so will not be continuous at any interval containing y=0. In order to use the theorem both must be continuous on an interval that contains yo=0 and this is problem for us since we do have yo=0.
Now, let’s actually work the problem. This differential equation is separable and is fairly simple to solve.
y13dy=dt32y23=t+cApplying the initial condition gives c=0 and so the solution is.
32y23=ty23=23ty2=(23t)3y(t)=±(23t)32So, we’ve got two possible solutions here, both of which satisfy the differential equation and the initial condition. There is also a third solution to the IVP. y(t)=0 is also a solution to the differential equation and satisfies the initial condition.

In this last example we had a very simple IVP and it only violated one of the conditions of the theorem, yet it had three different solutions. All the examples we’ve worked in the previous sections satisfied the conditions of this theorem and had a single unique solution to the IVP. This example is a useful reminder of the fact that, in the field of differential equations, things don’t always behave nicely. It’s easy to forget this as most of the problems that are worked in a differential equations class are nice and behave in a nice, predictable manner.

Let’s work one final example that will illustrate one of the differences between linear and non-linear differential equations.


Example 3 Determine the interval of validity for the initial value problem below and give its dependence on the value of yo
y=y2y(0)=y0

Before proceeding in this problem, we should note that the differential equation is non-linear and meets both conditions of the Theorem 2 and so there will be a unique solution to the IVP for each possible value of yo.
Also, note that the problem asks for any dependence of the interval of validity on the value of yo. This immediately illustrates a difference between linear and non-linear differential equations. Intervals of validity for linear differential equations do not depend on the value of yo. Intervals of validity for non-linear differential can depend on the value of yo as we pointed out after the second theorem.
So, let’s solve the IVP and get some intervals of validity.
First note that if yo=0 then y(t)=0 is the solution and this has an interval of validity of
<t<So, for the rest of the problem let's assume that y00. Now, the differential equation is separable so let's solve it and get a general solution.
y2dy=dt1y=t+cApplying the initial condition gives
c=1y0The solution is then.
1y=t1y0y(t)=11y0ty(t)=y01y0tNow that we have a solution to the initial value problem we can start finding intervals of validity. From the solution we can see that the only problem that we’ll have is division by zero at
t=1y0This leads to two possible intervals of validity.
<t<1y01y0<t<That actual interval of validity will be the interval that contains to=0. This however, depends on the value of yo. If yo<0 then 1y0<0 and so the second interval will contain to=0. Likewise, if yo>0 then 1y0>0 and in this case the first interval will contain to=0.
This leads to the following possible intervals of validity, depending on the value of yo.
If y0>0<t<1y0 is the interval of validity.If y0=0<t< is the interval of validity.If y0<01y0<t< is the interval of validity.On a side note, notice that the solution, in its final form, will also work if yo=0.

So, what did this example show us about the difference between linear and non-linear differential equations?
First, as pointed out in the solution to the example, intervals of validity for non-linear differential equations can depend on the value of yo, whereas intervals of validity for linear differential equations don’t.

Second, intervals of validity for linear differential equations can be found from the differential equation with no knowledge of the solution. This is definitely not the case with non-linear differential equations. It would be very difficult to see how any of these intervals in the last example could be found from the differential equation. Knowledge of the solution was required in order for us to find the interval of validity.

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