Sturm–Liouville theory

Sturm–Liouville theory is a part of mathematics. It studies special kinds of equations called Sturm–Liouville problems. These problems involve solving a type of equation with certain conditions at the ends of the interval being studied. The main goals are to find special numbers, called eigenvalues, for which there are non-trivial solutions, and then to find the matching solutions, called eigenfunctions.

This theory is very important in applied mathematics because Sturm–Liouville problems appear often. They are especially useful when solving separable linear partial differential equations. For example, in quantum mechanics, the one-dimensional time-independent Schrödinger equation is a Sturm–Liouville problem.

The theory is named after Jacques Charles François Sturm and Joseph Liouville, who developed it. For regular Sturm–Liouville problems, there are infinitely many eigenvalues, each with a unique eigenfunction. These eigenfunctions can be used to build a special space of functions called a Hilbert space.

Main results

Sturm–Liouville theory looks at special kinds of math equations. These equations help us understand how certain functions change over a range of values.

The main goal is to solve these equations under specific conditions. We want to find special numbers, called eigenvalues, and the functions that go with them, called eigenfunctions. These solutions help explain many natural patterns and phenomena.

For these problems to work well, certain conditions must be met. The functions need to be smooth and positive over the range being studied. There are also rules about how the solutions behave at the ends of the range.

The theory tells us important facts about these solutions, such as whether they cross zero and how often they do so. These properties are useful in many areas of science and engineering.

d d x [ p ( x ) d y d x ] + q ( x ) y = − λ w ( x ) y {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}\left[p(x){\frac {\mathrm {d} y}{\mathrm {d} x}}\right]+q(x)y=-\lambda \,w(x)y}

{ α 1 y ( a ) + α 2 y ′ ( a ) = 0 , α 1 , α 2 not both 0 , β 1 y ( b ) + β 2 y ′ ( b ) = 0 , β 1 , β 2 not both 0. {\displaystyle {\begin{cases}\alpha _{1}y(a)+\alpha _{2}y'(a)&=0,\qquad \alpha _{1},\alpha _{2}{\text{ not both }}0,\\\beta _{1}y(b)+\beta _{2}y'(b)&=0,\qquad \beta _{1},\beta _{2}{\text{ not both }}0.\end{cases}}}

Reduction to Sturm–Liouville form

The equation in this special shape is called Sturm–Liouville form or self-adjoint form. All second-order linear homogenous ordinary differential equations can be changed to look like this. We do this by multiplying both sides by a special kind of number called an integrating factor. But we cannot do this for second-order partial differential equations or when y is a vector.

Some examples are shown below.

Bessel equation

This equation can be changed to Sturm–Liouville form by dividing and grouping its parts.

Legendre equation

We can also change this equation to Sturm–Liouville form by spotting a pattern in its parts.

Example using an integrating factor

This example shows how multiplying by an integrating factor helps us change the equation.

Integrating factor for general second-order homogeneous equation

There is a general way to use an integrating factor to change any second-order homogenous equation into Sturm–Liouville form.

Sturm–Liouville equations as self-adjoint differential operators

This section talks about a special kind of math problem. It shows how a certain math rule, called L, works with functions. We can think of L as a tool that changes one function into another.

In math, we often look for special numbers called eigenvalues and special functions called eigenfunctions. For this problem, we want to find numbers λ and functions u that follow the rule L u = λ u. This means when we use the tool L on the function u, we get the same function back, but multiplied by the number λ.

This math rule L has a special property: it is self-adjoint. This means that when we use L on two functions and then compare the results in a certain way, it’s the same as comparing the original functions using L in the opposite order. This helps mathematicians understand the problem better.

The study of these problems is important in areas like quantum mechanics, where similar math rules describe how particles behave.

Application to inhomogeneous second-order boundary value problems

When solving some math problems, we need to find a special answer that works at two different points. This can be hard, but Sturm–Liouville theory helps by using special functions called eigenfunctions. These functions can be added together in a series to match the needed shape at the two points.

For example, if we want a solution that looks like a straight line between two points, we can use a series of sine waves. Even though these waves might look messy when added together, they still give the right answer when we need it. This method works for many different shapes and conditions, making it a useful tool in solving tricky math problems.

Application to partial differential equations

Some big math ideas can be solved using Sturm–Liouville theory. For example, think about a thin piece of material stretched tight like a drumhead, held in a square shape. We can study how it moves using special math steps.

One way to understand the movement is by breaking it into simple patterns that repeat over time, called "normal modes". Each pattern moves at its own special pace, and together they make up all possible ways the material can move.

Representation of solutions and numerical calculation

The Sturm–Liouville equation with boundary conditions can be solved in different ways. One way is to use special math methods like the Rayleigh–Ritz method or the matrix-variational method.

There are also numerical methods to solve these problems. Some of these methods include shooting methods, the finite difference method, and the spectral parameter power series method. These methods help find the special numbers, called eigenvalues, that solve the equation correctly.

X ( n ) ( t ) = { − ∫ a x X ( n − 1 ) ( t ) p ( t ) − 1 y 0 ( t ) − 2 d t n odd , ∫ a x X ( n − 1 ) ( t ) y 0 ( t ) 2 w ( t ) d t n even {\displaystyle X^{(n)}(t)={\begin{cases}\displaystyle -\int _{a}^{x}X^{(n-1)}(t)p(t)^{-1}y_{0}(t)^{-2}\,dt&n{\text{ odd}},\\[6pt]\displaystyle \quad \int _{a}^{x}X^{(n-1)}(t)y_{0}(t)^{2}w(t)\,dt&n{\text{ even}}\end{cases}}}

X ~ ( n ) ( t ) = { ∫ a x X ~ ( n − 1 ) ( t ) y 0 ( t ) 2 w ( t ) d t n odd , − ∫ a x X ~ ( n − 1 ) ( t ) p ( t ) − 1 y 0 ( t ) − 2 d t n even. {\displaystyle {\tilde {X}}^{(n)}(t)={\begin{cases}\displaystyle \quad \int _{a}^{x}{\tilde {X}}^{(n-1)}(t)y_{0}(t)^{2}w(t)\,dt&n{\text{ odd}},\\[6pt]\displaystyle -\int _{a}^{x}{\tilde {X}}^{(n-1)}(t)p(t)^{-1}y_{0}(t)^{-2}\,dt&n{\text{ even.}}\end{cases}}}