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 - Fourier Series: Eigenvalues and Eigenfunctions - ii


Example 2 Find all the eigenvalues and eigenfunctions for the following BVP.y+λy=0y(0)=0y(2π)=0

Here we are going to work with derivative boundary conditions. The work is pretty much identical to the previous example however so we won’t put in quite as much detail here. We’ll need to go through all three cases just as the previous example so let’s get started on that.
λ>0_
The general solution to the differential equation is identical to the previous example and so we have,

y(x)=c1cos(λx)+c2sin(λx)Applying the first boundary condition gives us,
0=y(0)=λc2c2=0Recall that we are assuming that λ>0 here and so this will only be zero if c2=0. Now, the second boundary condition gives us,
0=y(2π)=λc1sin(2πλ)Recall that we don’t want trivial solutions and that λ>0 so we will only get non-trivial solution if we require that,
sin(2πλ)=02πλ=nπn=1,2,3,Solving for λ and we see that we get exactly the same positive eigenvalues for this BVP that we got in the previous example.
λn=(n2)2=n24n=1,2,3,The eigenfunctions that correspond to these eigenvalues however are,
yn(x)=cos(nx2)n=1,2,3,So, for this BVP we get cosines for eigenfunctions corresponding to positive eigenvalues.
Now the second case.
λ=0_
The general solution is,

y(x)=c1+c2xApplying the first boundary condition gives,
0=y(0)=c2Using this the general solution is then,
y(x)=c1and note that this will trivially satisfy the second boundary condition,
0=y(2π)=0Therefore, unlike the first example, λ=0 is an eigenvalue for this BVP and the eigenfunctions corresponding to this eigenvalue is,
y(x)=1Again, note that we dropped the arbitrary constant for the eigenfunctions.
Finally let’s take care of the third case.
λ<0_
The general solution here is,

y(x)=c1cosh(λx)+c2sinh(λx)Applying the first boundary condition gives,
0=y(0)=λc1sinh(0)+λc2cosh(0)=λc2c2=0Applying the second boundary condition gives,
0=y(2π)=λc1sinh(2πλ)As with the previous example we again know that 2πλ0 and so sinh(2πλ)0. Therefore, we must have c1=0.
So, for this BVP we again have no negative eigenvalues.
In summary then we will have the following eigenvalues/eigenfunctions for this BVP.
λn=n24yn(x)=cos(nx2)n=1,2,3,λ0=0y0(x)=1Notice as well that we can actually combine these if we allow the list of n’s for the first one to start at zero instead of one. This will often not happen, but when it does we’ll take advantage of it. So the “official” list of eigenvalues/eigenfunctions for this BVP is,
λn=n24yn(x)=cos(nx2)n=0,1,2,3,
So, in the previous two examples we saw that we generally need to consider different cases for λ as different values will often lead to different general solutions. Do not get too locked into the cases we did here. We will mostly be solving this particular differential equation and so it will be tempting to assume that these are always the cases that we’ll be looking at, but there are BVP’s that will require other/different cases.
Also, as we saw in the two examples sometimes one or more of the cases will not yield any eigenvalues. This will often happen, but again we shouldn’t read anything into the fact that we didn’t have negative eigenvalues for either of these two BVP’s. There are BVP’s that will have negative eigenvalues.

Let’s take a look at another example with a very different set of boundary conditions. These are not the traditional boundary conditions that we’ve been looking at to this point, but we’ll see in the next chapter how these can arise from certain physical problems.


Example 3 Find all the eigenvalues and eigenfunctions for the following BVP.y+λy=0y(π)=y(π)y(π)=y(π)

So, in this example we aren’t actually going to specify the solution or its derivative at the boundaries. Instead we’ll simply specify that the solution must be the same at the two boundaries and the derivative of the solution must also be the same at the two boundaries. Also, this type of boundary condition will typically be on an interval of the form [-L,L] instead of [0,L] as we’ve been working on to this point.
As mentioned above these kind of boundary conditions arise very naturally in certain physical problems and we’ll see that in the next chapter.
As with the previous two examples we still have the standard three cases to look at.
λ>0_
The general solution for this case is,

y(x)=c1cos(λx)+c2sin(λx)Applying the first boundary condition and using the fact that cosine is an even function (i.e.cos(x)=cos(x)) and that sine is an odd function (i.e. sin(x)=sin(x)). gives us,
c1cos(πλ)+c2sin(πλ)=c1cos(πλ)+c2sin(πλ)c1cos(πλ)c2sin(πλ)=c1cos(πλ)+c2sin(πλ)c2sin(πλ)=c2sin(πλ)0=2c2sin(πλ)This time, unlike the previous two examples this doesn’t really tell us anything. We could have sin(πλ)=0 but it is also completely possible, at this point in the problem anyway, for us to have c2=0 as well.
So, let’s go ahead and apply the second boundary condition and see if we get anything out of that.
λc1sin(πλ)+λc2cos(πλ)=λc1sin(πλ)+λc2cos(πλ)λc1sin(πλ)+λc2cos(πλ)=λc1sin(πλ)+λc2cos(πλ)λc1sin(πλ)=λc1sin(πλ)2λc1sin(πλ)=0So, we get something very similar to what we got after applying the first boundary condition. Since we are assuming that λ>0 this tells us that either sin(πλ)=0 or c1=0.
Note however that if sin(πλ)0 then we will have to have c1=c2=0 and we’ll get the trivial solution. We therefore need to require that sin(πλ)=0 and so just as we’ve done for the previous two examples we can now get the eigenvalues,
πλ=nπλ=n2n=1,2,3,Recalling that λ>0 and we can see that we do need to start the list of possible n’s at one instead of zero.
So, we now know the eigenvalues for this case, but what about the eigenfunctions. The solution for a given eigenvalue is,
y(x)=c1cos(nx)+c2sin(nx)and we’ve got no reason to believe that either of the two constants are zero or non-zero for that matter. In cases like these we get two sets of eigenfunctions, one corresponding to each constant. The two sets of eigenfunctions for this case are,
yn(x)=cos(nx)yn(x)=sin(nx)n=1,2,3,Now the second case.
λ=0_
The general solution is,

y(x)=c1+c2xApplying the first boundary condition gives,
c1+c2(π)=c1+c2(π)2πc2=0c2=0Using this the general solution is then,
y(x)=c1and note that this will trivially satisfy the second boundary condition just as we saw in the second example above. Therefore, we again have λ=0 as an eigenvalue for this BVP and the eigenfunctions corresponding to this eigenvalue is,
y(x)=1Finally let’s take care of the third case.
λ<0_
The general solution here is,

y(x)=c1cosh(λx)+c2sinh(λx)Applying the first boundary condition and using the fact that hyperbolic cosine is even and hyperbolic sine is odd gives,
c1cosh(πλ)+c2sinh(πλ)=c1cosh(πλ)+c2sinh(πλ)c2sinh(πλ)=c2sinh(πλ)2c2sinh(πλ)=0Now, in this case we are assuming that λ<0 and so we know that πλ0 which in turn tells us that sinh(πλ)0. We therefore must have c2=0.
Let’s now apply the second boundary condition to get,
λc1sinh(πλ)=λc1sinh(πλ)2λc1sinh(πλ)=0c1=0By our assumption on λ we again have no choice here but to have c1=0.
Therefore, in this case the only solution is the trivial solution and so, for this BVP we again have no negative eigenvalues.
In summary then we will have the following eigenvalues/eigenfunctions for this BVP.
λn=n2yn(x)=sin(nx)n=1,2,3,λn=n2yn(x)=cos(nx)n=1,2,3,λ0=0y0(x)=1Note that we’ve acknowledged that for λ>0 we had two sets of eigenfunctions by listing them each separately. Also, we can again combine the last two into one set of eigenvalues and eigenfunctions. Doing so gives the following set of eigenvalues and eigenfunctions.
λn=n2yn(x)=sin(nx)n=1,2,3,λn=n2yn(x)=cos(nx)n=0,1,2,3,
Once again, we’ve got an example with no negative eigenvalues. We can’t stress enough that this is more a function of the differential equation we’re working with than anything and there will be examples in which we may get negative eigenvalues.
Now, to this point we’ve only worked with one differential equation so let’s work an example with a different differential equation just to make sure that we don’t get too locked into this one differential equation.

Before working this example let’s note that we will still be working the vast majority of our examples with the one differential equation we’ve been using to this point. We’re working with this other differential equation just to make sure that we don’t get too locked into using one single differential equation.

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