The [second-order hyperbolic theory](/pages/1075) treats a single scalar equation $u_{tt} + Lu = f$ for one unknown [function](/page/Function) $u$. Many important physical systems — the equations of gas dynamics, electromagnetism, elasticity — naturally involve **several coupled unknowns** evolving simultaneously. These are described by first-order hyperbolic systems, where the unknown is a vector $u = (u^1, \ldots, u^m)$ and the spatial [derivatives](/page/Derivative) are coupled through matrix coefficients.

{width=80%}

The passage from scalar to system is not merely a matter of generality: scalar second-order equations can always be **reduced** to first-order systems (by introducing auxiliary variables for derivatives), and conversely, first-order systems often arise directly from conservation laws rather than from higher-order equations. The first-order formulation makes the wave-propagation structure explicit: the eigenvalues of the coefficient matrices are the **characteristic speeds** at which information travels, and the eigenvectors determine which combinations of the unknowns propagate along each characteristic.

## The General First-Order System

[definition: First-Order Hyperbolic System] A **first-order system of PDEs** in $m$ unknowns is: \begin{align*} u_t + \sum_{k=1}^n A^k(x,t)\,u_{x_k} + B(x,t)\,u = f(x,t) \qquad \text{in } \mathbb{R}^n \times (0, T], \end{align*} with initial condition $u(\cdot, 0) = g$ in $\mathbb{R}^n$. Here $u: \mathbb{R}^n \times [0, T] \to \mathbb{R}^m$ is the unknown vector, $A^k(x,t)$ are $m \times m$ matrices for $k = 1, \ldots, n$, $B(x,t)$ is an $m \times m$ matrix of lower-order coefficients, and $f$ is a vector-valued source. [/definition]

The system is written more compactly as $u_t + \sum_k A^k u_{x_k} + Bu = f$. The term "first-order" refers to the fact that only first derivatives of $u$ appear (both in $t$ and in $x$), in contrast to the [second-order scalar equation](/pages/1075) $u_{tt} + Lu = f$ where second derivatives of $u$ occur.

The character of the system — whether it exhibits wave-like propagation or diffusive behaviour — is determined entirely by the coefficient matrices $A^k$. The key concept is:

[definition: Symmetric Hyperbolicity] The system $u_t + \sum_k A^k u_{x_k} + Bu = f$ is **symmetric hyperbolic** if the matrices $A^k(x,t)$ are all **symmetric** ($A^k = (A^k)^\top$) for every $(x,t)$. [/definition]

Symmetry of the coefficient matrices is the first-order analogue of the uniform ellipticity condition $\sum a_{ij}\xi_i\xi_j \ge \theta|\xi|^2$ that appears throughout the [elliptic](/pages/1033), [parabolic](/pages/1074), and [hyperbolic](/pages/1075) scalar theories. Just as ellipticity ensures coercivity of the bilinear form, symmetry of the $A^k$ ensures that the natural energy $\frac{1}{2}\|u(t)\|_{L^2}^2$ satisfies a useful energy identity — the antisymmetric part would destroy this.

[remark: Why Symmetry?] The reason symmetry is the right condition becomes clear when one derives the energy estimate (see below). Testing the equation with $u$ and integrating by parts produces terms of the form $\int u \cdot A^k u_{x_k}\,dx$. [Integration](/page/Integral) by parts gives: \begin{align*} \int_{\mathbb{R}^n} u \cdot A^k u_{x_k}\,dx = -\int_{\mathbb{R}^n} u_{x_k} \cdot A^k u\,dx - \int_{\mathbb{R}^n} u \cdot A^k_{x_k} u\,dx. \end{align*} If $A^k$ is symmetric, the first term on the right equals $-\int u \cdot A^k u_{x_k}\,dx$, so the two sides combine to give $2\int u \cdot A^k u_{x_k}\,dx = -\int u \cdot A^k_{x_k} u\,dx$, which is bounded by lower-order terms. Without symmetry, this cancellation fails. [/remark]

## Characteristics and Propagation Speeds

The physical content of a hyperbolic system is encoded in its **characteristics** — the directions along which information propagates. For a symmetric hyperbolic system, these are determined by the eigenvalues of the coefficient matrices.

[definition: Characteristic Speeds] For a direction $\xi \in \mathbb{R}^n \setminus \{0\}$, the **symbol** of the system at $(x,t)$ in direction $\xi$ is the $m \times m$ matrix: \begin{align*} A(x,t,\xi) := \sum_{k=1}^n \xi_k\,A^k(x,t). \end{align*} Since each $A^k$ is symmetric, $A(x,t,\xi)$ is symmetric for every $\xi$ and hence has $m$ real eigenvalues $\lambda_1(x,t,\xi) \le \cdots \le \lambda_m(x,t,\xi)$. These are the **characteristic speeds** of the system in direction $\xi$. [/definition]

{width=75%}

Each eigenvalue $\lambda_k$ represents a speed at which a particular "mode" of the solution propagates in direction $\xi$. The corresponding eigenvector determines which linear combination of the components $u^1, \ldots, u^m$ travels at that speed. The extreme eigenvalues $\lambda_1$ and $\lambda_m$ bound the speed of propagation: no information travels faster than $\max_{|\xi|=1}|\lambda_m(x,t,\xi)|$.

[remark: Reduction of Scalar Equations to Systems] The [second-order wave equation](/pages/1075) $u_{tt} - u_{xx} = 0$ (in one spatial dimension) can be written as a first-order symmetric hyperbolic system by introducing $v := u_t$ and $w := u_x$. The equation $u_{tt} = u_{xx}$ becomes $v_t = w_x$, and the compatibility condition $u_{tx} = u_{xt}$ gives $w_t = v_x$. In matrix form: \begin{align*} \begin{pmatrix} v \\ w \end{pmatrix}_t - \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} \begin{pmatrix} v \\ w \end{pmatrix}_x = \begin{pmatrix} 0 \\ 0 \end{pmatrix}. \end{align*} The coefficient matrix $A = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}$ is symmetric with eigenvalues $\lambda = \pm 1$ — the two wave speeds — and eigenvectors $(1, \pm 1)^\top$. The combinations $v + w = u_t + u_x$ and $v - w = u_t - u_x$ are the **Riemann invariants**, each propagating along one family of characteristics. {width=80%} More generally, any [second-order hyperbolic equation](/pages/1075) $u_{tt} + Lu = f$ in $n$ spatial dimensions can be converted to a first-order symmetric system in $m = n + 1$ unknowns $u = (u_t, u_{x_1}, \ldots, u_{x_n})$. This is one motivation for developing the first-order theory. [/remark]

## Energy Estimates

The fundamental tool for symmetric hyperbolic systems is, as in all the preceding theories, an energy estimate. The symmetry of the $A^k$ makes this particularly clean.

[theorem: Energy Estimate for Symmetric Hyperbolic Systems] Let $u$ be a smooth solution of the symmetric hyperbolic system $u_t + \sum_k A^k u_{x_k} + Bu = f$ in $\mathbb{R}^n \times (0, T]$ with $u(\cdot, 0) = g$. Assume $A^k, B \in C^1$ with bounded derivatives. Then there exists a constant $C > 0$, depending on $T$, $\|B\|_{L^\infty}$, and $\|\partial_{x_k} A^k\|_{L^\infty}$, such that: \begin{align*} \max_{0 \le t \le T} \|u(t)\|_{L^2(\mathbb{R}^n)}^2 \le C\left(\|g\|_{L^2(\mathbb{R}^n)}^2 + \int_0^T \|f(t)\|_{L^2(\mathbb{R}^n)}^2\,dt\right). \end{align*} [/theorem]

[proof] **Step 1: Energy computation.** Define $e(t) := \frac{1}{2}\int_{\mathbb{R}^n} |u(x,t)|^2\,d\mathcal{L}^n$. Differentiate: \begin{align*} e'(t) = \int_{\mathbb{R}^n} u \cdot u_t\,d\mathcal{L}^n = \int_{\mathbb{R}^n} u \cdot \left(f - \sum_k A^k u_{x_k} - Bu\right) d\mathcal{L}^n. \end{align*} **Step 2: [Integration by parts](/theorems/210).** The key term is $\int u \cdot A^k u_{x_k}\,d\mathcal{L}^n$. Since $A^k$ is symmetric, integration by parts (with the [boundary](/page/Boundary) terms vanishing since $u$ is defined on all of $\mathbb{R}^n$ and decays at infinity) gives: \begin{align*} \int_{\mathbb{R}^n} u \cdot A^k u_{x_k}\,d\mathcal{L}^n = -\frac{1}{2}\int_{\mathbb{R}^n} u \cdot (\partial_{x_k} A^k)\,u\,d\mathcal{L}^n. \end{align*} This is the step where symmetry is essential: it converts a first-order derivative term into a zeroth-order term. **Step 3: Gronwall.** Combining: \begin{align*} e'(t) \le \frac{1}{2}\sum_k \|\partial_{x_k} A^k\|_{L^\infty}\int |u|^2\,d\mathcal{L}^n + \|B\|_{L^\infty}\int |u|^2\,d\mathcal{L}^n + \int |u||f|\,d\mathcal{L}^n. \end{align*} By Cauchy-Schwarz and Young's inequality on the last term: $e'(t) \le C_0\,e(t) + \frac{1}{2}\|f(t)\|_{L^2}^2$, where $C_0$ depends on the coefficient bounds. Gronwall's inequality gives the result. [/proof]