Fact-checked by Grok 2 weeks ago

Phase space

In physics, phase space is a multidimensional space that represents all possible states of a dynamical system, with each axis corresponding to a coordinate such as position or momentum, allowing a complete description of the system's configuration at any instant.^[1] For a classical system of N particles in three-dimensional space, the phase space is typically 6N-dimensional, comprising 3N position coordinates and 3N momentum coordinates, forming a comprehensive framework for analyzing trajectories and evolutions under Hamiltonian mechanics.^[2] The concept originated in the late 19th century amid efforts to formalize statistical mechanics, with early contributions from Ludwig Boltzmann, who introduced the concept of "phase" in 1872 to describe the distribution of molecular states, and the term "phase space" coined by J. Willard Gibbs in 1902, with refinements by Henri Poincaré emphasizing its role in ensemble theory and ergodic behavior.^[3] In classical statistical mechanics, phase space underpins key principles such as Liouville's theorem, which states that the phase space volume occupied by an ensemble of systems remains constant over time due to the incompressible flow of trajectories, enabling predictions of thermodynamic properties from microscopic dynamics.^[4] This conservation property highlights phase space's utility in bridging deterministic mechanics with probabilistic descriptions of large systems. In quantum mechanics, phase space adapts to the uncertainty principle through quasi-probability distributions, such as the Wigner function, which provides a phase-space representation of the density operator while revealing interference effects and the classical limit. These formulations, developed prominently by Eugene Wigner in 1932, allow quantum states to be visualized in a continuous phase space despite the discrete nature of observables, facilitating comparisons between classical and quantum behaviors in areas like quantum optics and many-body physics.^[5] Beyond fundamental theory, phase space concepts extend to applications in chaos theory, where Poincaré sections map complex trajectories to lower dimensions, and in engineering fields like control systems, aiding the analysis of stability and attractors.^[6]

Fundamental Principles

Definition and Coordinates

In classical mechanics and statistical mechanics, the concept of phase space, developed in the late 19th century including contributions from Ludwig Boltzmann, was advanced by Josiah Willard Gibbs in his 1902 book Elementary Principles in Statistical Mechanics, where he explicitly termed it "phase space" while representing the states of a mechanical system in a multidimensional space.^[7] This framework was further formalized within Hamiltonian mechanics, building on William Rowan Hamilton's 1834 reformulation of dynamics, where phase space provides a complete description of a system's evolution. Phase space is mathematically defined as the cotangent bundle T^*Q of the configuration space Q, which parameterizes all possible positions of the system; for a system with n degrees of freedom, the phase space has dimension $2n, and each point (q, p) uniquely specifies the state of the system at a given time, encompassing both positional and momentum information.^[1] In this construction, the configuration space Q is typically a manifold of dimension n, such as \mathbb{R}^n for unconstrained particles, and the cotangent bundle equips it with fiber coordinates representing momenta, enabling a symplectic structure that preserves the geometry under dynamics.^[8] The standard coordinates of phase space consist of generalized position coordinates q_i (for i = 1, \dots, n) and their conjugate momenta p_i, defined via the Lagrangian as p_i = \frac{\partial L}{\partial \dot{q}_i}, where L is the system's Lagrangian.^[8] For a single particle in one dimension, phase space is two-dimensional with coordinates (q, p), where q is position and p = m \dot{q} is momentum; for N non-interacting particles in three dimensions, it extends to $6N-dimensional space with coordinates ( \mathbf{q}_1, \mathbf{p}_1, \dots, \mathbf{q}_N, \mathbf{p}_N ).^[1] The dynamics on this space are governed by the Hamiltonian H(q, p, t), the total energy expressed in terms of these coordinates, through Hamilton's equations:

\frac{dq_i}{dt} = \frac{\partial H}{\partial p_i}, \quad \frac{dp_i}{dt} = -\frac{\partial H}{\partial q_i}.

These equations describe trajectories in phase space as integral curves of the Hamiltonian vector field.^[8] Canonical transformations are coordinate changes (q, p) \to (Q, P) that preserve the form of Hamilton's equations and the symplectic structure of phase space, ensuring that the new coordinates Q_i, P_i also satisfy the conjugate pairing and the fundamental Poisson brackets \{q_i, p_j\} = \delta_{ij}.^[8] Such transformations maintain the geometric integrity of phase space, allowing equivalent descriptions of the same dynamics.

Conjugate Variables

In Hamiltonian mechanics, conjugate variables form pairs consisting of generalized coordinates q_i and their corresponding momenta p_i, where the momenta are defined through the Legendre transformation of the Lagrangian L(q, \dot{q}, t) to the Hamiltonian H(q, p, t), specifically p_i = \frac{\partial L}{\partial \dot{q}_i}.^[9]^[10] This pairing elevates the velocities \dot{q}_i to independent dynamical variables p_i, enabling a symmetric formulation of the equations of motion in phase space.^[11] The symplectic structure of phase space arises from these conjugate pairs, endowing the space with a Poisson bracket that encodes the fundamental algebraic relations: \{q_i, p_j\} = \delta_{ij}, \{q_i, q_j\} = 0, and \{p_i, p_j\} = 0, where \delta_{ij} is the Kronecker delta.^[12]^[13] More generally, for smooth functions f and g on the phase space, the Poisson bracket is defined as

\{f, g\} = \sum_i \left( \frac{\partial f}{\partial q_i} \frac{\partial g}{\partial p_i} - \frac{\partial f}{\partial p_i} \frac{\partial g}{\partial q_i} \right),

which satisfies bilinearity, antisymmetry, and the Jacobi identity, thereby defining a Poisson algebra on the manifold.^[12]^[14] This bracket governs the time evolution of any function f via \frac{df}{dt} = \{f, H\} + \frac{\partial f}{\partial t}, where H is the Hamiltonian, yielding Hamilton's equations as a special case for f = q_i or p_i.^[8]^[13] Canonical coordinates refer to any set of variables (Q_k, P_k) obtained from the original (q_i, p_i) via point transformations that preserve the fundamental Poisson brackets, ensuring the symplectic structure remains invariant.^[15]^[16] Such transformations maintain the form of Hamilton's equations and are essential for simplifying problems in classical mechanics.^[17] Darboux's theorem guarantees that any symplectic manifold admits a local coordinate system where the symplectic form takes the canonical Darboux form \sum_i dq_i \wedge dp_i, affirming the local existence of conjugate coordinates around any point.^[18]^[19] This result underscores the uniformity of symplectic geometry, implying no local invariants beyond the dimension in such spaces.^[20]

Liouville's Theorem

Liouville's theorem states that in Hamiltonian systems, the volume of any region in phase space remains invariant under time evolution, implying an incompressible flow where phase space volumes neither expand nor contract despite the deformation of their shapes.^[21]^[22] This invariance underscores the deterministic and reversible nature of classical Hamiltonian dynamics, ensuring that the measure of accessible states is preserved along trajectories.^[21] The theorem was originally proved by Joseph Liouville in 1838 in the context of differential equations, without explicit reference to phase space or mechanics.^[23] An independent proof was provided by Carl Gustav Jacob Jacobi around 1842, later published in 1866, where he applied it to mechanical systems using Hamilton's equations.^[24] To sketch the proof, consider the velocity field in phase space derived from Hamilton's equations for conjugate variables q_i and p_i: \dot{q}_i = \frac{\partial H}{\partial p_i} and \dot{p}_i = -\frac{\partial H}{\partial q_i}, so the field is \mathbf{v} = \left( \frac{\partial H}{\partial p_1}, \dots, -\frac{\partial H}{\partial q_1}, \dots \right).^[22] The divergence is

\nabla \cdot \mathbf{v} = \sum_i \left( \frac{\partial}{\partial q_i} \frac{\partial H}{\partial p_i} + \frac{\partial}{\partial p_i} \left( -\frac{\partial H}{\partial q_i} \right) \right) = 0,

since the mixed partial derivatives commute.^[22] By the divergence theorem, the rate of change of a phase space volume V is

\frac{dV}{dt} = \int_V (\nabla \cdot \mathbf{v}) \, dV = 0.

^[22] Geometrically, the time evolution corresponds to a canonical transformation that preserves the symplectic 2-form \omega = \sum_i dq_i \wedge dp_i, ensuring the flow maintains the volume element dV = \omega^n / n! for an n-degree-of-freedom system.^[21] This preservation highlights how Hamiltonian flows act like shear deformations in phase space, distorting shapes without altering volumes.^[21]

Low-Dimensional Representations

One-Dimensional Systems

In one-dimensional systems, the phase space is a two-dimensional plane parameterized by the position coordinate q and its conjugate momentum p, representing the state of a system with a single degree of freedom.^[25] Trajectories in this space trace the evolution of the system under the governing dynamics, with the Hamiltonian H(q, p) determining the flow.^[26] A canonical example is the simple harmonic oscillator, governed by the Hamiltonian H = \frac{p^2}{2m} + \frac{1}{2} k q^2, where m is the mass and k is the spring constant.^[27] The constant-energy contours form closed elliptical trajectories in the (q, p) plane, as the total energy E = H remains fixed, yielding p = \pm \sqrt{2m \left( E - \frac{1}{2} k q^2 \right)}.^[26] The parametric equations for the motion are q(t) = A \cos(\omega t + \phi) and p(t) = -m \omega A \sin(\omega t + \phi), where A is the amplitude, \omega = \sqrt{k/m}, and \phi is the initial phase; plotting these parametrically produces the ellipse.^[27] These closed orbits reflect periodic motion, with the area enclosed by the trajectory proportional to the action variable, preserved under Hamiltonian dynamics as per Liouville's theorem.^[25] For a free particle, the Hamiltonian simplifies to H = \frac{p^2}{2m}, with no potential term.^[25] Since \dot{p} = -\frac{\partial H}{\partial q} = 0, momentum p is constant, resulting in straight-line trajectories parallel to the q-axis at fixed p = \sqrt{2m E}.^[28] This illustrates unbounded motion without oscillation, where position evolves linearly as q(t) = q_0 + \frac{p}{m} t. In contrast, a damped harmonic oscillator introduces dissipation, modifying the equations to \ddot{q} + 2\beta \dot{q} + \omega^2 q = 0, where \beta > 0 is the damping coefficient.^[25] This non-Hamiltonian system produces spiraling inward trajectories in phase space, with amplitude decaying exponentially and energy loss causing the orbit to contract toward the origin.^[29] The phase space volume contracts at a rate \nabla \cdot \mathbf{v} = -2\beta < 0, violating the area preservation of Liouville's theorem applicable to conservative systems.^[25]

Two-Dimensional Systems

In classical mechanics, systems with two degrees of freedom possess a four-dimensional phase space spanned by the generalized coordinates q_1, q_2 and their conjugate momenta p_1, p_2.^[30] This structure arises because each degree of freedom contributes two dimensions to the phase space, allowing trajectories to evolve on three-dimensional energy hypersurfaces within the full four-dimensional manifold.^[30] Due to the high dimensionality, direct visualization is challenging, so common projections onto two-dimensional planes, such as the (q_1, p_1) or (q_2, p_2) subspace, are employed to reveal qualitative behaviors while marginalizing over the other variables.^[31] For uncoupled harmonic oscillators, the Hamiltonian takes the separable form

H = \sum_{i=1}^2 \left( \frac{p_i^2}{2m} + \frac{1}{2} k q_i^2 \right),

where m is the mass and k the spring constant, assuming identical parameters for simplicity.^[32] The trajectories in the full phase space are then the Cartesian product of independent elliptical orbits in each (q_i, p_i) plane, with constant areas enclosed by each ellipse determined by the action variables.^[33] Projections onto a single (q_i, p_i) plane thus show isolated ellipses, while the full dynamics exhibit Lissajous figures when projecting onto coordinate-momentum mixed planes, reflecting the incommensurate frequencies if \omega_1 \neq \omega_2.^[34] In coupled oscillators, such as two masses connected by springs, the coupling introduces off-diagonal terms in the potential, leading to normal modes where the system decouples into independent oscillators after a linear transformation to normal coordinates.^[32] These modes manifest as in-phase (beating with lower frequency) and out-of-phase (higher frequency) oscillations, with the phase space trajectories becoming separable in action-angle variables after diagonalization.^[34] The action variables J_i quantify the energy in each mode, while angle variables \theta_i track the phases, enabling a toroidal representation of the invariant tori on which motion occurs.^[33] Projections of these trajectories often reveal quasi-periodic filling of annular regions, hinting at the underlying separability.^[32] A canonical example of non-separable dynamics is the double pendulum, with two degrees of freedom corresponding to the angles \theta_1, \theta_2 and angular momenta p_{\theta_1}, p_{\theta_2}. In this four-dimensional phase space, trajectories for small amplitudes approximate tori from linearized normal modes, but larger energies produce chaotic orbits that appear as tangled, space-filling curves in projections like (\theta_1, p_{\theta_1}), indicating exponential sensitivity to initial conditions. These projections capture the transition from regular to chaotic motion, where nearby trajectories diverge rapidly, though the full chaos analysis resides in higher-dimensional applications. To manage the four-dimensional complexity, especially in periodically driven two-degree-of-freedom systems, Poincaré sections provide a reduction to two-dimensional stroboscopic maps by recording intersections of trajectories with a chosen hypersurface, such as when a phase variable equals a constant.^[31] For instance, in a driven double oscillator, sampling at driving period intervals yields points in a (q_1, p_1) plane that reveal invariant curves for integrable cases or scattered points for chaotic ones, effectively lowering the dimensionality for analysis.^[31] This technique preserves the symplectic structure and highlights periodic orbits as fixed points in the map.^[31]

Phase Plot

A phase plot is a graphical representation of a dynamical system's trajectory in phase space, depicting the evolution of state variables such as position and velocity over time to form a curve that traces the system's path.^[29] This visualization eliminates the explicit dependence on time, focusing instead on the relationship between the variables themselves.^[35] Phase plots are constructed through parametric plotting of the solutions to the governing differential equations, where one state variable is plotted against another as derived from the system's dynamics.^[35] For instance, in a one-dimensional system like simple harmonic motion, the plot of velocity versus position yields an ellipse for undamped cases, illustrating the closed orbit corresponding to periodic behavior.^[29] In higher dimensions, however, direct visualization of phase plots becomes limited, as projecting multi-variable trajectories onto two or three dimensions can obscure the full structure. One key advantage of phase plots is their ability to reveal qualitative aspects of the system's dynamics—such as periodic, quasi-periodic, or aperiodic motion—without requiring explicit analytical or numerical solutions to the equations.^[35] This approach allows for intuitive insight into the trajectory's shape and stability, aiding in the analysis of oscillatory or rotational behaviors. The early development of phase plots is attributed to Henri Poincaré's work in celestial mechanics around 1890, where he used such graphical methods to explore trajectories in the three-body problem, laying foundational techniques for visualizing complex orbital dynamics.

Phase Portrait

A phase portrait provides a qualitative representation of the global structure of trajectories in the phase space of a dynamical system, illustrating the flow defined by the vector field across the entire space. It consists of a collection of representative trajectories that reveal the overall behavior, including how solutions evolve from various initial conditions without solving the equations explicitly. This visualization is particularly useful in low-dimensional systems, such as two-dimensional phase planes, where the portrait captures the topological features of the flow.^[36]^[37] Key elements of a phase portrait include equilibrium points, where the vector field vanishes, classifying the local dynamics as sinks (stable nodes or foci), sources (unstable nodes or foci), or saddles based on the nature of nearby trajectories. Stable and unstable manifolds emanate from these points, with stable manifolds attracting trajectories and unstable manifolds repelling them, often forming separatrices that divide the phase space into regions of distinct behavior. Additionally, homoclinic orbits connect a single equilibrium point to itself, while heteroclinic orbits link distinct equilibria, both contributing to the portrait's global topology by highlighting invariant sets and barriers to flow.^[38]^[39] To analyze these elements, linearization is employed around each equilibrium point by computing the Jacobian matrix of the vector field, whose eigenvalues determine the local phase portrait type. For a two-dimensional autonomous system given by

\dot{x} = f(x, y), \quad \dot{y} = g(x, y),

equilibrium points satisfy f(x^*, y^*) = 0 and g(x^*, y^*) = 0. The Jacobian J at (x^*, y^*) is

J = \begin{pmatrix} \frac{\partial f}{\partial x} & \frac{\partial f}{\partial y} \\ \frac{\partial g}{\partial x} & \frac{\partial g}{\partial y} \end{pmatrix}_{(x^*, y^*)},

with trace \tau = \operatorname{tr}(J) and determinant \Delta = \det(J). The eigenvalues are roots of \lambda^2 - \tau \lambda + \Delta = 0; the equilibrium is a stable node if \tau < 0 and \Delta > 0 with real eigenvalues of the same sign, a saddle if \Delta < 0, and a focus if the eigenvalues are complex with negative real part for stability. Basins of attraction, regions converging to a particular attractor, are delineated by these manifolds and orbits in the full portrait.^[40]^[37]^[41] A canonical example is the Van der Pol oscillator, modeled by \ddot{x} - \mu (1 - x^2) \dot{x} + x = 0 for \mu > 0, rewritten in phase space as \dot{x} = y, \dot{y} = \mu (1 - x^2) y - x. The origin is the sole equilibrium, an unstable focus, surrounded by a stable limit cycle that appears as a closed orbit in the phase portrait, attracting all nearby trajectories and demonstrating self-sustained oscillations central to nonlinear dynamics.^[42]

Phase Integral

The phase integral, also known as the action integral, is defined as the line integral \oint p \, dq taken over a closed path in phase space, where p is the momentum and q is the position coordinate.^[43] This integral quantifies the action associated with a periodic orbit in the system.^[43] In classical mechanics, the phase integral plays a central role in the formulation of action-angle variables, which are canonical coordinates suited for integrable Hamiltonian systems with periodic motion.^[43] The action variable J is given by J = \frac{1}{2\pi} \oint p(q) \, dq, where the integral is evaluated over one complete cycle of the motion.^[44] In these coordinates, the Hamiltonian depends only on the actions H = H(J), and the conjugate angle variables \phi evolve linearly with time.^[43] The frequency \omega of the motion is then \omega = \frac{\partial H}{\partial J}.^[45] Moreover, the action J serves as an adiabatic invariant, remaining constant under slow variations of system parameters, such as in gradually changing potentials.^[46] Historically, the phase integral was instrumental in the old quantum theory developed by Niels Bohr in 1913 and extended by Arnold Sommerfeld in 1916, where it formed the basis for the quantization condition \oint p \, dq = n h (with n an integer and h Planck's constant) to predict atomic energy levels.^[47] In two-dimensional phase space, the phase integral J corresponds directly to the area A enclosed by the orbit, via the relation J = A / 2\pi.^[44] This geometric interpretation underscores its utility in visualizing conserved quantities for bounded motion.

Applications in Physics

Dynamical Systems and Chaos Theory

In dynamical systems, phase space serves as the arena for visualizing and analyzing the time evolution of states governed by ordinary differential equations (ODEs), where each trajectory corresponds to a unique solution curve parametrized by time. These trajectories encapsulate the deterministic unfolding of the system, revealing fixed points, limit cycles, or more complex behaviors depending on the nonlinearity and dimensionality. In dissipative systems, characterized by contraction in certain phase space directions, trajectories from diverse initial conditions converge toward attractors—invariant sets of finite measure that capture the long-term dynamics, such as stable equilibria or periodic orbits.^[48]^[49] Chaos emerges in nonlinear dynamical systems when trajectories exhibit extreme sensitivity to initial conditions, a phenomenon quantified by Lyapunov exponents that measure the average exponential rates of divergence or convergence in phase space. The largest Lyapunov exponent \lambda is given by

\lambda = \lim_{t \to \infty} \frac{1}{t} \ln \left| \frac{\delta \mathbf{x}(t)}{\delta \mathbf{x}(0)} \right|,

where \delta \mathbf{x}(t) denotes the infinitesimal separation between nearby trajectories at time t; a positive \lambda signifies local instability and chaotic motion, as small perturbations grow exponentially. Strange attractors, which underpin chaotic regimes in dissipative systems, possess fractal dimensions—typically non-integer values between the topological and embedding dimensions—arising from the intricate folding and stretching of phase space volumes, leading to self-similar structures with infinite detail at finer scales. Illustrative examples highlight phase space's role in chaos. The Lorenz attractor, derived from a three-dimensional truncation of Navier-Stokes equations for atmospheric convection, manifests as a butterfly-shaped strange attractor in the projection of its phase space (variables x, y, z), where trajectories weave chaotically between lobes, never repeating yet bounded by dissipation.^[50] Similarly, the Hénon map, a discrete two-dimensional system defined by iterations x_{n+1} = 1 - a x_n^2 + y_n and y_{n+1} = b x_n (with parameters a=1.4, b=0.3), produces a fractal strange attractor in the (x, y) plane, demonstrating how simple quadratic maps can yield dense, non-periodic orbits sensitive to initial points. Ergodicity in chaotic dynamical systems implies that, for almost all initial conditions on an invariant measure, the time average of an observable along a single trajectory equals the spatial average over the relevant phase space subset, such as an energy surface. This equivalence facilitates the use of long-time simulations to approximate ensemble properties and connects to measure-preserving flows in Hamiltonian chaos via Liouville's theorem, ensuring incompressible evolution in phase space.

Statistical Mechanics and Thermodynamics

In statistical mechanics, phase space serves as the arena for describing ensembles of systems, where the positions and momenta of particles define the microscopic states. The microcanonical ensemble corresponds to an isolated system with fixed energy, volume, and particle number, represented by a uniform probability distribution over the hypersurface of constant energy in phase space. The entropy S is then given by S = k \ln \Omega, where k is Boltzmann's constant and \Omega is the volume of the accessible phase space divided by h^{3N} N!, with h being Planck's constant and N the number of particles, ensuring dimensional consistency and accounting for particle indistinguishability.^[51]^[52] The canonical ensemble, applicable to systems in thermal contact with a heat reservoir at temperature T, employs a Boltzmann distribution for the probability density \rho in phase space, given by \rho \propto \exp(-\beta H), where \beta = 1/(kT) and H(q,p) is the Hamiltonian. The partition function Z, which normalizes this distribution and encodes thermodynamic properties, is computed as the phase space integral Z = \frac{1}{h^{3N} N!} \int \exp(-\beta H(q,p)) \, dq \, dp. This integral over the entire phase space yields quantities like the Helmholtz free energy via F = -kT \ln Z, linking microscopic dynamics to macroscopic thermodynamics.^[51]^[52]^[53] Liouville's theorem plays a crucial role by guaranteeing the conservation of phase space volume under Hamiltonian evolution, which underpins the stationarity of equilibrium distributions in both ensembles. This incompressibility of phase space flow ensures that probability densities remain constant along trajectories, allowing time-independent ensemble averages for equilibrium states.^[51]^[54] The ergodic hypothesis further justifies the use of ensemble averages by positing that, for sufficiently chaotic systems, the time average of an observable equals the phase space average over the ensemble, enabling the equivalence of temporal and statistical descriptions in thermodynamics. This assumption, essential for deriving equilibrium properties from dynamical evolution, holds for most physical systems of interest.^[55]

Quantum Mechanics

In quantum mechanics, phase space formulations provide a bridge between classical and quantum descriptions by representing quantum states as quasi-probability distributions over position and momentum. The Wigner-Weyl transform, introduced by Eugene Wigner in 1932, maps quantum operators and states onto functions in phase space, enabling a reformulation of quantum mechanics in terms of these distributions.^[56] The central object is the Wigner function W(q, p), defined for a wave function \psi(q) as

W(q, p) = \frac{1}{\pi \hbar} \int_{-\infty}^{\infty} \psi^*(q + y) \psi(q - y) e^{2 i p y / \hbar} \, dy,

which transforms the density operator into a phase-space representation.^[56] This function generalizes the classical phase-space probability density, recovering classical statistical mechanics in the semiclassical limit as \hbar \to 0. Key properties of the Wigner function include its marginal distributions: integrating over momentum yields the position probability density \int W(q, p) \, dp = |\psi(q)|^2, while integrating over position gives the momentum space probability \int W(q, p) \, dq = |\tilde{\psi}(p)|^2 / (2\pi \hbar).^[56] Unlike classical probabilities, the Wigner function can take negative values, a hallmark of non-classical quantum behavior that signals interference effects without a direct classical analog; the volume of negativity quantifies this quantumness.^[57] The algebraic structure of quantum phase space relies on the Moyal bracket, the quantum counterpart to the classical Poisson bracket, which governs dynamics and operator products. Defined as

\{f, g\}_M = \frac{2}{\hbar} f \sin\left( \frac{\hbar}{2} (\partial_q \partial_{p'} - \partial_p \partial_{q'}) \right) g \bigg|_{p'=p, q'=q},

it arises from the antisymmetric part of the Moyal star product, deforming classical commutators into a non-commutative phase-space multiplication. Time evolution of the Wigner function follows the von Neumann equation in this framework, yielding \partial W / \partial t = \{ H, W \}_M, where H(q, p) is the Weyl-transformed Hamiltonian; this equation preserves the quasi-probability structure under unitary dynamics. Wigner's original work laid the foundation for applications in quantum optics, where the function characterizes non-classical light states like squeezed vacuum through negativity and higher-order correlations, and in quantum state tomography, enabling full reconstruction of quantum systems from phase-space measurements such as homodyne detection.^[56]^[58]^[59]

Optics

In geometric optics, phase space provides a framework for representing light rays as points in a coordinate system where the horizontal axis denotes position q and the vertical axis denotes optical momentum p = n \sin \theta, with n as the refractive index and \theta as the ray angle relative to the optical axis.^[60] This representation allows the evolution of ray bundles through optical elements to be tracked as transformations in phase space, preserving the structure of the bundle in lossless systems. The conservation of etendue, defined as the phase space volume occupied by the rays, follows from Liouville's theorem and holds for reversible optical transformations without absorption or scattering. For beam characterization in paraxial optics, the Wigner distribution function serves as a phase space representation of the beam's intensity and phase, particularly for partially coherent light. It is given by

W(q, p) = \frac{1}{\lambda^2} \iint \psi^*\left(q - \frac{r}{2}\right) \psi\left(q + \frac{r}{2}\right) e^{-i k p \cdot r} \, d^2 r,

where \psi is the complex field amplitude, [\lambda](/page/Lambda) is the wavelength, and k = 2\pi / \lambda.^[60] The area in phase space enclosed by the Wigner distribution relates directly to beam quality, quantified by the factor M^2, which measures deviation from an ideal Gaussian beam; for a diffraction-limited beam, M^2 = 1, and the phase space area is \lambda / 4 per transverse dimension.^[61] The ABCD matrix formalism describes linear optical transformations, such as those induced by lenses and mirrors, as symplectic mappings in phase space that preserve the determinant of the transformation matrix (equal to 1 for lossless systems). A ray vector \begin{pmatrix} q \\ p \end{pmatrix} at the output is obtained via \begin{pmatrix} q' \\ p' \end{pmatrix} = \begin{pmatrix} A & B \\ C & D \end{pmatrix} \begin{pmatrix} q \\ p \end{pmatrix}, enabling efficient simulation of beam propagation through complex systems.^[62] An illustrative example is the propagation of a Gaussian beam, where the invariant emittance \varepsilon = \frac{1}{2\pi} \iint dq \, dp remains constant and equals \lambda / 4\pi for a fundamental mode, capturing the minimum phase space volume due to diffraction.^[63] In lossless optical systems, the etendue

G = n^2 \iint \cos \theta \, dq \, d\theta

is conserved, linking the spatial extent and angular spread of the beam across transformations.^[64]

Applications in Other Fields

Engineering and Control Theory

In engineering and control theory, phase space provides a foundational framework for representing and analyzing dynamic systems through state-space models. These models describe the evolution of a system's state vector \mathbf{x}(t) \in \mathbb{R}^n, which captures all relevant variables in the phase space, via the linear differential equation \dot{\mathbf{x}} = A \mathbf{x} + B \mathbf{u}, where A is the state matrix, B is the input matrix, and \mathbf{u} is the control input.^[65] The output \mathbf{y} = C \mathbf{x} + D \mathbf{u} relates the phase space states to measurable quantities, with C and D as output matrices.^[65] This representation, introduced by Kalman, unifies time-domain analysis and enables the design of controllers that manipulate trajectories within the phase space to achieve desired behaviors.^[65] Stability analysis in phase space often employs Lyapunov's direct method, which assesses asymptotic stability without solving the differential equations explicitly. A Lyapunov function V(\mathbf{x}) is selected as positive definite in the phase space (i.e., V(\mathbf{0}) = 0 and V(\mathbf{x}) > 0 for \mathbf{x} \neq \mathbf{0}), and its time derivative along system trajectories must satisfy \dot{V}(\mathbf{x}) = \frac{\partial V}{\partial \mathbf{x}} (A \mathbf{x} + B \mathbf{u}) < 0 for \mathbf{x} \neq \mathbf{0} under appropriate control \mathbf{u}.^[66] For linear systems, quadratic forms V(\mathbf{x}) = \mathbf{x}^T P \mathbf{x} (with P > 0) lead to the Lyapunov equation A^T P + P A = -Q for some positive definite Q, ensuring phase space trajectories converge to the origin.^[66] This approach extends to nonlinear systems by constructing V to bound phase space regions, guaranteeing global stability when \dot{V} is negative definite.^[66] Reachability and observability concepts leverage phase space to evaluate control system capabilities. Reachability determines if any state in the phase space can be attained from the origin using admissible inputs, quantified by the controllability Gramian W_c(T) = \int_0^T e^{A \tau} B B^T e^{A^T \tau} d\tau, whose rank equals the dimension of the reachable subspace if full.^[67] Observability assesses whether initial phase space states can be reconstructed from outputs, via the observability Gramian W_o(T) = \int_0^T e^{A^T \tau} C^T C e^{A \tau} d\tau, ensuring full rank for complete state inference.^[67] These Gramians, integrated over phase space evolutions, guide controller design by identifying controllable and observable subspaces, as formalized in linear system theory.^[67] For nonlinear systems, phase space partitioning facilitates hybrid control, where the state space is divided into regions with distinct dynamics or control laws. This partitioning, often using hyperplanes or polyhedra, models interfaces between continuous and discrete modes in hybrid systems, enabling stability guarantees via piecewise Lyapunov functions across partitions.^[68] Covering halfspaces approximate these boundaries, extending linear methods to handle switching in phase space while preserving reachability properties.^[68] A practical example is PID control for second-order systems, visualized in the phase plane of error e and error rate \dot{e}. The proportional term acts radially toward the origin, the derivative term along switching lines to damp oscillations, and the integral term shifts the effective equilibrium, resulting in spiral trajectories converging to zero error; tuning ensures encirclements avoid limit cycles.^[69]

Medicine and Biology

In medicine and biology, phase space provides a framework for analyzing the nonlinear dynamics of physiological systems, enabling the reconstruction of hidden attractors from time series data to distinguish healthy from pathological states. For cardiac dynamics, electrocardiogram (ECG) signals are reconstructed into phase space using Takens' embedding theorem, which guarantees that a time-delay embedding of sufficient dimension preserves the topological structure of the underlying attractor for deterministic systems.^[70] This reconstruction typically involves selecting an embedding dimension m and time delay \tau, transforming the one-dimensional ECG into a higher-dimensional space where trajectories reveal periodic or chaotic behaviors.^[71] In normal sinus rhythm, the phase space attractor is compact and low-dimensional, reflecting stable oscillatory dynamics, whereas arrhythmic states like ventricular tachycardia or fibrillation exhibit strange attractors with higher fractal dimensions and disrupted trajectories, indicating loss of synchronization.^[72]^[73] To determine the appropriate embedding dimension for heart rate variability (HRV) analysis in reconstructed phase space, the false nearest neighbors (FNN) method is employed, which identifies the minimal m where the attractor unfolds without false crossings. The FNN algorithm proceeds as follows: for each point in the m-dimensional embedding, find its nearest neighbor; then project to m+1 dimensions and compute the ratio r_i(m) = \frac{|\mathbf{x}_{i}(m+1) - \mathbf{y}_{i}(m+1)|}{|\mathbf{x}_{i}(m) - \mathbf{y}_{i}(m)|}; a neighbor is deemed false if r_i(m) > R_t (threshold, often 10–50) or if boundary effects dominate. The percentage of false neighbors drops to near zero at the correct m, typically 3–5 for HRV signals in healthy hearts but increasing in pathological conditions like congestive heart failure, signifying higher dynamical complexity.^[74]^[75]^[76] In population biology, the Lotka-Volterra predator-prey model exemplifies phase space analysis of ecological interactions, with the two-dimensional space spanned by prey density x and predator density y. The system is governed by \frac{dx}{dt} = \alpha x - \beta x y and \frac{dy}{dt} = \delta x y - \gamma y, where parameters represent growth, predation, and death rates, yielding closed periodic orbits around the equilibrium point (x^*, y^*) = (\gamma/\delta, \alpha/\beta). These cycles in phase space illustrate oscillatory population fluctuations without damping, contrasting stable equilibria in non-interacting models and providing insights into biodiversity maintenance.^[77]^[78] Neurodynamics leverages phase space to model neuron firing, as in the Hodgkin-Huxley equations, which describe action potentials in a four-dimensional space of membrane potential V and gating variables m, h, n for sodium and potassium channels. The dynamics, captured by C \frac{dV}{dt} = -g_{Na} m^3 h (V - V_{Na}) - g_K n^4 (V - V_K) - g_L (V - V_L) + I, along with cubic activation and exponential inactivation terms, produce limit cycle attractors corresponding to repetitive spiking for suprathreshold currents. This 4D representation reveals bifurcations from resting states to oscillatory firing, foundational for understanding neural excitability in biological circuits. Clinical applications include detecting chaotic dynamics in electroencephalogram (EEG) signals for epilepsy diagnosis, with 1990s research establishing nonlinear measures like Lyapunov exponents to quantify pre-ictal transitions. In epileptic seizures, phase space reconstruction of EEG reveals low-dimensional attractors with positive exponents indicating chaos, differing from the higher-dimensional noise-like activity in healthy brains; this approach, pioneered in intracranial EEG analysis, enables seizure prediction by tracking dynamical instability.^[79]^[80]

References

[1]
[PDF] 1. Phase space - UCLA Mathematics
In classical mechanics, the phase space is the space of all possible states of a physical system; by “state” we do not simply mean the positions q of all ...
[2]
[PDF] Phase space in classical physics Quantum mechanically, we can ...
First, we define the notion of phase space—each momentum and position coordinate counts as one dimension of phase space. Thus, for a single particle, there are ...
[3]
[PDF] The tangled tale of phase space - Purdue Physics department
The origins of both the concept of phase space and its name are historically obscure—which is surprising in view of the central role it plays in practically ...
[4]
[PDF] 1 Classical statistical physics
For a single particle, phase space in the 6-dimensional space with coordinates. (q,p) where q is the particle's position and p is its momentum. As the particle ...
[5]
Overview of the phase space formulation of quantum mechanics ...
Feb 22, 2021 · Here, an overview of quantum mechanics in phase space is presented. The formulation to generate a generalized phase-space function for any arbitrary quantum ...
[6]
[PDF] Physics 5153 Classical Mechanics The Hamiltonian and Phase Space
In the Hamiltonian approach, the Hamiltonian is defined in terms of the spatial coordinates and momenta, and the equations of motion are first order ...
[7]
[PDF] The Physics of J. Willard Gibbs in his Time
Who was the "Mr. Josiah Willard Gibbs of New Haven" who was appointed professor of mathematical physics by the Yale Corporation on 13 July 1871?'
[8]
[PDF] Elementary Principles in Statistical Mechanics
Josiah Willard Gibbs (1839–1903) was the greatest American mathematician and ... * If we regard a phase as represented by a point in space of 2 n dimen-.
[9]
[PDF] 4. The Hamiltonian Formalism - DAMTP
The phase space of the pendulum is a cylinder R⇥S1, with the R factor corresponding to the momentum. We draw this by flattening out the cylinder. The two ...<|control11|><|separator|>
[10]
7.13: Hamiltonian in Classical Mechanics - Physics LibreTexts
Jun 30, 2021 · Chapter 8 uses the Legendre transformation between the conjugate variables ( q , q ˙ , t ) and ( q , p , t ) where the generalized coordinate q ...
[11]
[PDF] Hamiltonian Mechanics
. Hamiltonian formalism uses qi and pi as dynamical variables, where pi are generalized momenta defined by pi = ∂L. ∂ ˙qi . (0.1). The resulting 2N ...
[12]
[PDF] Review of Hamiltonian Mechanics - The Cockcroft Institute
In Hamiltonian mechanics, the canonical variables consist of a set of coordinates and their conjugate momenta defined by equation (18). In general: The word ...
[13]
15.2: Poisson bracket Representation of Hamiltonian Mechanics
Jun 28, 2021 · The Poisson bracket representation of Hamiltonian mechanics provides a direct link between classical mechanics and quantum mechanics.
[14]
[PDF] Generalized Phase Space 1 Symplectic Structure - Hitoshi Murayama
3 Poisson Bracket. The Poisson brackets are defined using the inverse of the symplectic matrix as. {A, B} = ωij ∂A. ∂xi. ∂B. ∂xj . (10). For the conventional ...<|separator|>
[15]
Poisson Bracket - an overview | ScienceDirect Topics
"A Poisson bracket is defined as a binary operation on a set of variables that satisfies specific algebraic properties, such as {x_i, x_j} = a_ij x_i x_j, ...
[16]
15.3: Canonical Transformations in Hamiltonian Mechanics
Mar 21, 2021 · Canonical transformations are the foundation of Hamiltonian mechanics; they underlie Hamilton-Jacobi theory and action-angle variable theory.Generating functions · Type 1: F = F 1 ⁡ ( q , Q , t ) · Type 2: F = F 2 ⁡ ( q , P , t ) − Q...
[17]
[PDF] Chapter 4 Canonical Transformations, Hamilton-Jacobi Equations ...
The simplest case is for a single dimension of canonical coordinates (q, p). If H(q, p) = α1 then p = p(q, α1). There are two types of periodic motion to ...
[18]
[PDF] Hamiltonian Dynamics - Lecture 1 - CERN Indico
Nov 10, 2020 · In Hamiltonian mechanics, the canonical momenta pi = δL are promoted to coordinates on equal footing with the generalized coordinates qi .
[19]
Darboux's theorem in nLab
Nov 17, 2017 · For the case of symplectic manifolds the Darboux theorem may also be read as saying that a G-structure for G = Sp ( 2 n ) the symplectic group ...
[20]
[PDF] LECTURE 24: DARBOUX THEOREM 1. Moser's trick Recall that ...
Theorem 2.1 (Darboux theorem). Let (M,ω) be a symplectic manifold of dimension 2n. Then for any p ∈ M, there exists a neighborhood U of p in M and a ...
[21]
proof of Darboux's theorem (symplectic geometry) - PlanetMath.org
Mar 22, 2013 · Indeed, if we have a symplectic manifold (M,η) ( M , η ) , and a point x0 x 0 , we can take a (smooth) coordinate chart about x0 x 0 .
[22]
4 The Hamiltonian Formalism - DAMTP
The central idea of Liouville's theorem – that volume of phase space is constant – is somewhat reminiscent of quantum mechanics. Indeed, this is the first of ...
[23]
[PDF] 3 Conservation of volume in phase space - DSpace@MIT
In conservative systems, the conservation of volumes in phase space is known as Liouville's theorem. 4 Damped oscillators and dissipative systems. 4.1 ...
[24]
[PDF] Note sur la Théorie de la Variation des constantes arbitraires
NOTE. Sur la Théorie de la Variation des constantes arbitraires : PAR J LIOUVILLE. Soient η un nombre entier positif, x une fonction de t dont nous ...Missing: Joseph | Show results with:Joseph
[25]
C. G. J. Jacobi's Vorlesungen über Dynamik - Internet Archive
Jun 9, 2008 · Jacobi's Vorlesungen über Dynamik: Gehalten an der Universität zu Königsberg im ... ... PDF download · download 1 file · SINGLE PAGE PROCESSED ...
[26]
[PDF] Contents - Physics Courses
Consider next the case of the damped harmonic oscillator. In this case, phase space volumes contract. For a one- dimensional oscillator obeying x+2β ˙x+Ω2. 0 ...
[27]
[PDF] PHY411 Lecture notes -Introduction to Hamiltonian and Lagrangian ...
Sep 18, 2023 · The phase curves for the harmonic oscillator are shown in Figure 1. ... Phase space is 4 dimensional and so one conserved quantity is not ...<|control11|><|separator|>
[28]
9. The Simple Harmonic Oscillator - Galileo
The best we can do is to place the system initially in a small cell in phase space, of size Δx⋅Δp=ℏ/2. In fact, we shall find that in quantum mechanics phase ...
[29]
[PDF] Problem 1. Phase-space and its characteristic flow
(i) Show by direct substitution that for a free particle H = P2/2m the solution to. ∂f(t, Q, P). ∂t. + {f,H}. P,Q. = 0. (7) is f(t, Q, P) = f0(Q −. P m t, P) ...<|control11|><|separator|>
[30]
Phase Space Diagrams for an Oscillator (undamped and damped)
May 28, 2008 · A phase-space plot is a parametric graph of the velocity v(t) plotted as a function of the displacement x(t), with the changing variable being time.Missing: definition | Show results with:definition
[31]
[PDF] Notes on Classical Mechanics: 2nd module Phase space ...
A system with 2 degrees of freedom has a 4d phase space. In the absence of any conserved quantity, the trajectory could explore the whole of phase space. If ...
[32]
[PDF] Poincare Sections - MIT OpenCourseWare
Oct 17, 2022 · 1 Poincaré sections. The dynamical systems we study are of the form d dt x(t) = F(x, t). Systems of such equations describe a flow in phase ...
[33]
6.1: Two Coupled Oscillators - Physics LibreTexts
Jan 26, 2022 · Let us discuss oscillations in systems with several degrees of freedom, starting from the simplest case of two linear (harmonic), dissipation-free, 1D ...Missing: dimensional | Show results with:dimensional
[34]
15.5: Action-angle Variables - Physics LibreTexts
Mar 14, 2021 · The action-angle variable approach uses a canonical transformation to action and angle variables which provide a powerful, and elegant method to exploit ...Missing: coupled | Show results with:coupled
[35]
[PDF] Modern stellar dynamics, lecture 5: action–angle variables
Motion is separable in x,y – two uncoupled simple harmonic oscillators, two integrals of motion Ex ,Ey , actions are Jx = Ex /ωx ,Jy = Ey /ωy . -1. -0.5. 0.
[36]
Differential Equations - Phase Plane - Pauls Online Math Notes
Nov 16, 2022 · In this section we will give a brief introduction to the phase plane and phase portraits. We define the equilibrium solution/point for a ...
[37]
[PDF] 26. Phase portraits in two dimensions - MIT OpenCourseWare
The x, y plane is called the phase y plane (because a point in it represents the state or phase of a system). The phase portrait is a representative sampling of ...Missing: key elements
[38]
[PDF] Dynamic Behavior
A phase portrait is constructed by plotting the flow of the vector field corresponding to the planar dynamical system. That is, for a set of initial conditions, ...
[39]
[PDF] Towards a Global Theory of Singularly Perturbed Dynamical Systems
The surface of periodic orbits and the stable and unstable manifolds of the saddles are shown. The homoclinic orbit is contained in the closure of all three of ...
[40]
[PDF] Untitled
2 Computing stable manifolds of equilibria. The stable and unstable manifolds of equilibrium points are important geometric objects in the phase portraits of ...
[41]
[PDF] Math 312 Lecture Notes Linearization - Department of Mathematics
Mar 23, 2005 · These notes discuss linearization, in which a linear system is used to approximate the behavior of a nonlinear system.Missing: dynamical | Show results with:dynamical
[42]
[PDF] Discrete and Continuous Dynamical Systems - MIT Mathematics
May 18, 2014 · Behavior of vector fields near a fixed point in a linearized system is determined by the determinant ∆ and the trace τ of the Jacobian matrix.<|separator|>
[43]
[PDF] The Van der Pol Oscillator
r = 2 is a stable limit cycle. A limit r = 2 cycle is a periodic orbit ... The phase portrait in u and y, shown above right, is more transparent. The y ...
[44]
[PDF] 6.8. ACTION-ANGLE VARIABLES 199 200 CHAPTER 6 ...
Action-angle variables are canonical coordinates on phase space for integrable systems, where each point is labeled by ~φ, and each φi is a coordinate on a ...
[45]
Phys410 Notes, Demos and Supplements - Fall 2011 - Physics - UMD
... 1/2π ∮p dq. The integral is around the closed orbit, and is equal to the phase space area enclosed by the orbit. If k(t) changes slowly enough, the system ...
[46]
[PDF] arXiv:2011.11684v1 [astro-ph.GA] 23 Nov 2020
Nov 23, 2020 · ˙Θ = ∂H/∂J is known as the frequency Ω(J). Therefore. Θ(t) = Θ0 + Ω t,. (9) and hence the angles are linearly dependent on time. Finally ...
[47]
[PDF] 13. Adiabatic Invariants and Action-Angle Variables ) ( )
An adiabatic invariant stays constant when system parameters change gradually, even though the system's energy changes.
[48]
[PDF] {How Sommerfeld extended Bohr's model of the atom (1913–1916)}
Jan 30, 2014 · Sommerfeld's early reaction to the Bohr atom, 1913–1914 If Bohr's paper in the Philosophical Magazine with the pretentious title “On the ...
[49]
[PDF] arXiv:2003.07174v1 [nlin.AO] 28 Feb 2020
Feb 28, 2020 · dynamical systems is the phase space, a space containing every possible unique state the system could be found in. The evolution of a system ...
[50]
[PDF] arXiv:2208.12101v1 [math.DS] 25 Aug 2022
Aug 25, 2022 · It discusses the classical and modern results of the theory of attractors for dissipative PDEs including attractors for autonomous and non- ...
[51]
Deterministic Nonperiodic Flow in - AMS Journals
Deterministic nonperiodic flow involves nonlinear equations where solutions are unstable, and slightly different initial states can lead to different states. ...
[52]
[PDF] Elementray Principles in Statistical Mechanics. - Project Gutenberg
Jan 22, 2016 · ∗If we regard a phase as represented by a point in space of 2n dimensions, the ... Mechanics, by Josiah Willard Gibbs. *** END OF THIS PROJECT ...
[53]
[PDF] 2. Classical Gases - DAMTP
Including the factor of N! in the partition function ensures that the entropy does not increase when identical atoms are mixed5. 2.2.4 The Ideal Gas in the ...
[54]
[PDF] Lecture 7: Ensembles
Ensembles are the possible microstates a system could be in, sharing a macroscopic property. Examples include microcanonical, canonical, Gibbs, and grand ...
[55]
[PDF] Chapter 4: Statistical Mechanics [version 1004.1.K] - Caltech PMA
the system's location in phase space is equal to the statistical-equilibrium ensemble average: ¯. A ≡ lim. T→0 & +T/2. −T/2A(q(t), p(t)) = &A% ≡ #n. Anρn ...<|control11|><|separator|>
[56]
Ergodic theorem, ergodic theory, and statistical mechanics - PNAS
Feb 17, 2015 · This hypothesis states that some orbit of the flow will pass arbitrarily close to every point of phase space, or in other words this orbit is ...
[57]
On the Quantum Correction For Thermodynamic Equilibrium
The formula is developed for this correction by means of a probability function and the result discussed. References (11). J. von Neumann, Gött. Nachr.
[58]
Wigner Function Negativity and Contextuality in Quantum ...
Apr 2, 2015 · The key difference is that Wigner functions can assume negative values, and this negativity is taken as an indication of quantumness.Missing: marginal | Show results with:marginal
[59]
Wigner function tomography via optical parametric amplification
Aug 22, 2023 · Wigner function tomography is indispensable for characterizing quantum states, but its commonly used version, balanced homodyne detection, ...
[60]
Generalized overlap quantum state tomography1
Oct 7, 2020 · In Section I, we provide the formalism for measuring the Wigner function overlap of quantum optical states. This formalism is extended in ...
[61]
[PDF] Numerical phase space optics methods and applications to the ...
Each point in phase space corresponds to an optical ray. Conservation of étendue. Given a second plane Σ2 parallel to Σ1, the phase space volume element 𝑑𝐺 = 𝑑 ...
[62]
M^2 Factor – M squared, laser beam, quality factor ... - RP Photonics
The M² factor is a parameter for quantifying the beam quality of laser beams. It is strongly related to the beam parameter product.
[63]
[PDF] ABCD Matrices as Similarity Transformations of Wigner ... - arXiv
May 16, 2008 · Abstract. The beam transfer matrix, often called the ABCD matrix, is a two-by-two matrix with unit determinant, and with three independent.
[64]
[PDF] Rays in phase space & Gaussian beams - CLASSE (Cornell)
M2 beam quality factor. M2 = 1 only for a pure Gaussian beam! Page 18. Gaussian beams. P3330 Exp Optics. 18. Links/References.Missing: M² | Show results with:M²
[65]
Application of the conservation of etendue theorem for 2-D ...
The etendue of the 2-D bundle is obtained by adding (integrating) each differential contribution dxdp + dydq (or dx′dp′ + dy′dq′) of the bundle. Considering ...
[66]
[PDF] Mathematical Description of Linear Dynamical Systems - Duke People
Control theory is supposed to deal with physical systems, and not merely with mathematical objects such as a differential equation or a trans- fer function. We ...
[67]
[PDF] 4 Lyapunov Stability Theory
In this section we review the tools of Lyapunov stability theory. These tools will be used in the next section to analyze the stability properties.<|separator|>
[68]
[PDF] Notes on Linear Systems Theory
This matrix is called the control- lability Gramian. (Note that FBS2e uses Wr to represent the reachability matrix [B AB A2B ...]. This is different than the ...Missing: phase | Show results with:phase
[69]
State Space Partitioning for Hybrid Control Systems - IEEE Xplore
In this paper, the state space partitioning in the interface of the hybrid control system is modeled using covering halfspaces; this formulation extends the ...Missing: phase | Show results with:phase
[70]
[PDF] PID Control
In this chapter we discuss the basic ideas of PID control and the methods for choosing the parameters of the controllers. Many aspects of control can be ...
[71]
OptRPC: A novel and optimized recurrence plot-based system for ...
In this work, we treated cardiac movement observed through ECG as a dynamical system. Takens' embedding theorem was used to reconstruct phase space containing ...
[72]
Experimental analysis of heart rate variability of long-recording ...
As mentioned above, the Takens theorem dictates that proper phase space reconstruction is achieved when m is greater than [2mc+1]. This criterion is difficult ...<|separator|>
[73]
Phase space reconstruction of an experimental model of cardiac ...
The shape of attractors remains similar but disturbed in case of arrhythmia compared to the normal conditions.Missing: states | Show results with:states
[74]
Analysis of an experimental model of in vitro cardiac tissue using ...
Phase space reconstructions in normal and arrhythmic signals highlight strange attractors. •. The time lag, the fractal and correlation dimensions could be ...
[75]
Increase in the embedding dimension in the heart rate variability ...
Oct 6, 2006 · The false nearest neighbor (FNN) method was proposed by Kennel et al. to determine m 0 ⁠.15 There are different ways to implement the FNN method ...
[76]
Time Series Analysis using Embedding Dimension on Heart Rate ...
False Nearest Neighbour (FNN) method of estimating ED has been adapted for analysing and predicting variables responsible for HRV time series.Missing: neighbors | Show results with:neighbors
[77]
New algorithm of mortality risk prediction for cardiovascular patients ...
In this paper, mutual information (MI) method for choosing the time delay and the false nearest neighbor method is used for selecting the embedding dimension.
[78]
The Lotka‐Volterra Predator‐Prey Model with Foraging–Predation ...
In the first case, the population equilibrium is in the part of the prey‐predator density phase space where predator activity is at its maximum ( v = 1 v = 1 ), ...
[79]
10.2: Modelling Predator-Prey Interactions - Biology LibreTexts
Nov 12, 2021 · In this chapter, we will explore the classic Lotka-Volterra predator-prey model (Rosenzweig and MacArthur 1963), which treats each population as if it were ...
[80]
REVIEW : Chaos Theory and Epilepsy - Leonidas D. Iasemidis, J ...
It is hypothesized that epilepsy may be an example of chaos. In this review, some of the basic concepts of nonlinear dynamics and chaos are illustrated.<|separator|>
[81]
Can Chaotic Analysis of Electroencephalogram Aid the Diagnosis of ...
The chaotic features of EEG have been shown in previous studies to be an indicator of the complexity of brain dynamics. The smaller values of chaotic features ...