Floquet theory explained

Floquet theory is a branch of the theory of ordinary differential equations relating to the class of solutions to periodic linear differential equations of the form

x

=A(t)x,

with

x\in{Rn

} and

\displaystyleA(t)\in{R{n x

}} being a piecewise continuous periodic function with period

T

and defines the state of the stability of solutions.

\displaystyley=Q-1(t)x

with

\displaystyleQ(t+2T)=Q(t)

that transforms the periodic system to a traditional linear system with constant, real coefficients.

When applied to physical systems with periodic potentials, such as crystals in condensed matter physics, the result is known as Bloch's theorem.

Note that the solutions of the linear differential equation form a vector space. A matrix

\phi(t)

is called a fundamental matrix solution if the columns form a basis of the solution set. A matrix

\Phi(t)

is called a principal fundamental matrix solution if all columns are linearly independent solutions and there exists

t0

such that

\Phi(t0)

is the identity. A principal fundamental matrix can be constructed from a fundamental matrix using

\Phi(t)=\phi(t){\phi}-1(t0)

. The solution of the linear differential equation with the initial condition

x(0)=x0

is

x(t)=\phi(t){\phi}-1(0)x0

where

\phi(t)

is any fundamental matrix solution.

Floquet's theorem

Let

x

=A(t)x

be a linear first order differential equation,where

x(t)

is a column vector of length

n

and

A(t)

an

n x n

periodic matrix with period

T

(that is

A(t+T)=A(t)

for all real values of

t

). Let

\phi(t)

be a fundamental matrix solution of this differential equation. Then, for all

t\inR

,

\phi(t+T)=\phi(t)\phi-1(0)\phi(T).

Here

\phi-1(0)\phi(T)

is known as the monodromy matrix.In addition, for each matrix

B

(possibly complex) such that

eTB=\phi-1(0)\phi(T),

there is a periodic (period

T

) matrix function

t\mapstoP(t)

such that

\phi(t)=P(t)etBforallt\inR.

Also, there is a real matrix

R

and a real periodic (period-

2T

) matrix function

t\mapstoQ(t)

such that

\phi(t)=Q(t)etRforallt\inR.

In the above

B

,

P

,

Q

and

R

are

n x n

matrices.

Consequences and applications

This mapping

\phi(t)=Q(t)etR

gives rise to a time-dependent change of coordinates (

y=Q-1(t)x

), under which our original system becomes a linear system with real constant coefficients
y

=Ry

. Since

Q(t)

is continuous and periodic it must be bounded. Thus the stability of the zero solution for

y(t)

and

x(t)

is determined by the eigenvalues of

R

.

The representation

\phi(t)=P(t)etB

is called a Floquet normal form for the fundamental matrix

\phi(t)

.

The eigenvalues of

eTB

are called the characteristic multipliers of the system. They are also the eigenvalues of the (linear) Poincaré maps

x(t)\tox(t+T)

. A Floquet exponent (sometimes called a characteristic exponent), is a complex

\mu

such that

e\mu

is a characteristic multiplier of the system. Notice that Floquet exponents are not unique, since
(\mu+
2\piik
T
)T
e

=e\mu

, where

k

is an integer. The real parts of the Floquet exponents are called Lyapunov exponents. The zero solution is asymptotically stable if all Lyapunov exponents are negative, Lyapunov stable if the Lyapunov exponents are nonpositive and unstable otherwise.

References

. Ivar Ekeland. One. Convexity methods in Hamiltonian mechanics. Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)]. 19. Springer-Verlag. Berlin. 1990. x+247. 3-540-50613-6. 1051888.

. Teschl . Gerald . Gerald Teschl . Ordinary Differential Equations and Dynamical Systems . . . 2012 . 978-0-8218-8328-0 .