Mellin transform

From formulasearchengine
Jump to navigation Jump to search

Template:Differential equations

In mathematics, separation of variables (also known as the Fourier method) is any of several methods for solving ordinary and partial differential equations, in which algebra allows one to rewrite an equation so that each of two variables occurs on a different side of the equation.

Ordinary differential equations (ODE)

Suppose a differential equation can be written in the form

ddxf(x)=g(x)h(f(x))

which we can write more simply by letting y=f(x):

dydx=g(x)h(y).

As long as h(y) ≠ 0, we can rearrange terms to obtain:

dyh(y)=g(x)dx,

so that the two variables x and y have been separated. dx (and dy) can be viewed, at a simple level, as just a convenient notation, which provides a handy mnemonic aid for assisting with manipulations. A formal definition of dx as a differential (infinitesimal) is somewhat advanced.

Alternative notation

Some who dislike Leibniz's notation may prefer to write this as

1h(y)dydx=g(x),

but that fails to make it quite as obvious why this is called "separation of variables". Integrating both sides of the equation with respect to x, we have

1h(y)dydxdx=g(x)dx,(1)

or equivalently,

1h(y)dy=g(x)dx

because of the substitution rule for integrals.

If one can evaluate the two integrals, one can find a solution to the differential equation. Observe that this process effectively allows us to treat the derivative dydx as a fraction which can be separated. This allows us to solve separable differential equations more conveniently, as demonstrated in the example below.

(Note that we do not need to use two constants of integration, in equation (1) as in

1h(y)dy+C1=g(x)dx+C2,

because a single constant C=C2C1 is equivalent.)

Example (I)

The ordinary differential equation

ddxf(x)=f(x)(1f(x))

may be written as

dydx=y(1y).

If we let g(x)=1 and h(y)=y(1y), we can write the differential equation in the form of equation (1) above. Thus, the differential equation is separable.

As shown above, we can treat dy and dx as separate values, so that both sides of the equation may be multiplied by dx. Subsequently dividing both sides by y(1y), we have

dyy(1y)=dx.

At this point we have separated the variables x and y from each other, since x appears only on the right side of the equation and y only on the left.

Integrating both sides, we get

dyy(1y)=dx,

which, via partial fractions, becomes

1ydy+11ydy=1dx,

and then

ln|y|ln|1y|=x+C

where C is the constant of integration. A bit of algebra gives a solution for y:

y=11+Bex.

One may check our solution by taking the derivative with respect to x of the function we found, where B is an arbitrary constant. The result should be equal to our original problem. (One must be careful with the absolute values when solving the equation above. It turns out that the different signs of the absolute value contribute the positive and negative values for B, respectively. And the B = 0 case is contributed by the case that y = 1, as discussed below.)

Note that since we divided by y and (1y) we must check to see whether the solutions y(x)=0 and y(x)=1 solve the differential equation (in this case they are both solutions). See also: singular solutions.

Example (II)

Population growth is often modeled by the differential equation

dPdt=kP(1PK)

where P is the population with respect to time t, k is the rate of growth, and K is the carrying capacity of the environment.

Separation of variables may be used to solve this differential equation.

dPdt=kP(1PK)
dPP(1PK)=kdt

To evaluate the integral on the left side, we simplify the fraction

1P(1PK)=KP(KP)

and then, we decompose the fraction into partial fractions

KP(KP)=1P+1KP

Thus we have

(1P+1KP)dP=kdt

ln|P|ln|KP|=kt+C

ln|KP|ln|P|=ktC

ln|KPP|=ktC

|KPP|=ektC

|KPP|=eCekt

KPP=±eCekt

Let A=±eC.

KPP=Aekt

KP1=Aekt

KP=1+Aekt

PK=11+Aekt

P=K1+Aekt

Therefore, the solution to the logistic equation is

P(t)=K1+Aekt

To find A, let t=0 and P(0)=P0. Then we have

P0=K1+Ae0

Noting that e0=1, and solving for A we get

A=KP0P0

Partial differential equations<pde>...</pde>

The method of separation of variables is also used to solve a wide range of linear partial differential equations with boundary and initial conditions, such as heat equation, wave equation, Laplace equation and Helmholtz equation.

Homogeneous case

Consider the one-dimensional heat equation.The equation is

Template:NumBlk

The boundary condition is homogeneous, that is Template:NumBlk

Let us attempt to find a solution which is not identically zero satisfying the boundary conditions but with the following property: u is a product in which the dependence of u on x, t is separated, that is:

Template:NumBlk

Substituting u back into equation and using the product rule,

Template:NumBlk

Since the right hand side depends only on x and the left hand side only on t, both sides are equal to some constant value − λ. Thus:

Template:NumBlk

and

Template:NumBlk

− λ here is the eigenvalue for both differential operators, and T(t) and X(x) are corresponding eigenfunctions.

We will now show that solutions for X(x) for values of λ ≤ 0 cannot occur:

Suppose that λ < 0. Then there exist real numbers B, C such that

X(x)=Beλx+Ceλx.

From Template:EqNote we get

Template:NumBlk

and therefore B = 0 = C which implies u is identically 0.

Suppose that λ = 0. Then there exist real numbers B, C such that

X(x)=Bx+C.

From Template:EqNote we conclude in the same manner as in 1 that u is identically 0.

Therefore, it must be the case that λ > 0. Then there exist real numbers A, B, C such that

T(t)=Aeλαt,

and

X(x)=Bsin(λx)+Ccos(λx).

From Template:EqNote we get C = 0 and that for some positive integer n,

λ=nπL.

This solves the heat equation in the special case that the dependence of u has the special form of Template:EqNote.

In general, the sum of solutions to Template:EqNote which satisfy the boundary conditions Template:EqNote also satisfies Template:EqNote and Template:EqNote. Hence a complete solution can be given as

u(x,t)=n=1DnsinnπxLexp(n2π2αtL2),

where Dn are coefficients determined by initial condition.

Given the initial condition

u|t=0=f(x),

we can get

f(x)=n=1DnsinnπxL.

This is the sine series expansion of f(x). Multiplying both sides with sinnπxL and integrating over [0,L] result in

Dn=2L0Lf(x)sinnπxLdx.

This method requires that the eigenfunctions of x, here {sinnπxL}n=1, are orthogonal and complete. In general this is guaranteed by Sturm-Liouville theory.

Nonhomogeneous case

Suppose the equation is nonhomogeneous,

Template:NumBlk

with the boundary condition the same as Template:EqNote.

Expand h(x,t), u(x,t) and f(x,t) into

Template:NumBlk

Template:NumBlk

Template:NumBlk

where hn(t) and bn can be calculated by integration, while un(t) is to be determined.

Substitute Template:EqNote and Template:EqNote back to Template:EqNote and considering the orthogonality of sine functions we get

u'n(t)+αn2π2L2un(t)=hn(t),

which are a sequence of linear differential equations that can be readily solved with, for instance, Laplace transform, or Integrating factor. Finally, we can get

un(t)=eαn2π2L2t(bn+0thn(s)eαn2π2L2sds).

If the boundary condition is nonhomogeneous, then the expansion of Template:EqNote and Template:EqNote is no longer valid. One has to find a function v that satisfies the boundary condition only, and subtract it from u. The function u-v then satisfies homogeneous boundary condition, and can be solved with the above method.

In orthogonal curvilinear coordinates, separation of variables can still be used, but in some details different from that in Cartesian coordinates. For instance, regularity or periodic condition may determine the eigenvalues in place of boundary conditions. See spherical harmonics for example.

Matrices

The matrix form of the separation of variables is the Kronecker sum.

As an example we consider the 2D discrete Laplacian on a regular grid:

L=DxxDyy=DxxI+IDyy,

where Dxx and Dyy are 1D discrete Laplacians in the x- and y-directions, correspondingly, and I are the identities of appropriate sizes. See the main article Kronecker sum of discrete Laplacians for details.

References

  • A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9.
  • 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.

    My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
  • 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.

    My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

External links