Fact-checked by Grok 2 weeks ago

Bounded variation

In mathematical analysis, a function of bounded variation (BV function) on a closed interval [a, b] is a real-valued function f: [a, b] \to \mathbb{R} whose total variation V_a^b(f)—defined as the supremum over all partitions P = \{a = x_0 < x_1 < \cdots < x_n = b\} of the sums \sum_{i=1}^n |f(x_i) - f(x_{i-1})|—is finite. This condition captures functions that do not oscillate excessively, generalizing monotone functions while excluding highly irregular ones like the Weierstrass nowhere-differentiable function. Functions of bounded variation were introduced by Camille Jordan in 1881 to analyze the pointwise convergence of Fourier series, forming the basis for Jordan's criterion on such convergence. A fundamental property is that every BV function can be expressed as the difference of two non-decreasing functions, and conversely, any such difference is of bounded variation; this decomposition separates the function into its absolutely continuous part, jump discontinuities (captured by a saltus function), and a singular continuous part. BV functions are bounded on [a, b], have at most countably many discontinuities, and are differentiable almost everywhere with respect to Lebesgue measure. The space of BV functions on [a, b], denoted BV[a, b], is a Banach space under the norm \|f\|_{BV} = |f(a)| + V_a^b(f), and it contains important subclasses such as Lipschitz and absolutely continuous functions, though the converse does not hold. Applications extend to integration theory, where BV functions serve as integrators in Riemann-Stieltjes integrals, and to geometry, as the graph of a BV function has finite arc length if and only if the function is BV. In higher dimensions and metric measure spaces, BV functions generalize to spaces with weak derivatives in L^1, playing a key role in calculus of variations and partial differential equations.

Historical Context

Early Foundations

The concept of functions of bounded variation emerged in the late 19th century as a tool to address convergence issues in Fourier series. In 1881, Camille Jordan introduced the notion while investigating the convergence of trigonometric series, motivated by the need for sufficient conditions ensuring pointwise convergence to the generating function. In his work on Fourier series, Jordan defined the total variation of a function f on a closed interval [a, b] as the least upper bound (supremum) of sums \sum |f(x_i) - f(x_{i-1})| taken over all finite partitions a = x_0 < x_1 < \cdots < x_n = b of the interval. He observed that functions with finite total variation possess at most countably many discontinuities and that their Fourier series converge to the function at points of continuity, a result that marked a significant advance over Dirichlet's earlier uniform continuity criterion. Jordan emphasized that bounded variation implies a certain regularity sufficient for convergence in the context of trigonometric representations. This definition immediately connected to monotone functions, which Jordan recognized as having bounded variation equal to their total net change over the interval, since increasing or decreasing functions yield variation precisely |f(b) - f(a)|. Monotone functions thus form a foundational subclass of bounded variation functions, with Jordan proving that any bounded variation function can be decomposed as the difference of two monotone functions—a result now known as the . This property underscored the utility of bounded variation in analysis, providing a bridge between continuous and discontinuous behaviors. Shortly thereafter, in 1885, Cesare Arzelà extended these ideas in his work on rectifiable curves and functions, exploring monotone increasing and decreasing functions as primitives for understanding variation in geometric contexts. Arzelà's contributions highlighted how bounded variation captures the "length" of a function's graph, linking it to arc length in the plane and paving the way for applications in differential geometry. His analysis reinforced Jordan's framework by showing that functions of bounded variation admit a canonical representation via monotone components, further solidifying their role in early real analysis. The framework laid by Jordan and Arzelà found natural extension to the , introduced by Thomas Joannes Stieltjes in 1894, where the integrator function requires bounded variation to ensure integrability for continuous integrands. This connection emphasized that monotone functions, as a subset of bounded variation functions, suffice for many integration purposes, with the total variation providing a measure of the integrator's "roughness." While initial developments focused on one dimension, brief considerations of multivariable generalizations appeared in subsequent works, adapting the supremum over partitions to higher-dimensional domains.

Key Developments in the 20th Century and Beyond

In the early 20th century, the theory of (BV) functions was extended from one dimension to several variables through integral representations of the total variation. Leonida Tonelli introduced this generalization in 1926, defining the total variation for functions of multiple variables via the integral of the absolute value of the gradient, which laid the groundwork for multidimensional . This approach was further refined by Lamberto Cesari in 1936, who formalized in higher dimensions using a perimeter-based measure, establishing key properties like closure under weak convergence and applications to . A significant advancement came in 1967 with A.I. Vol'pert's work, which developed a calculus for BV functions including the chain rule and introduced the concept of approximate limits to handle discontinuities precisely. Vol'pert also defined special BV (SBV) functions, decomposing the variation into absolutely continuous, jump, and Cantor parts, which became essential for analyzing . During the 1970s and 1980s, Ennio De Giorgi contributed to the integration of BV theory into partial differential equations (), particularly through techniques for relaxation of variational integrals involving discontinuous solutions. His frameworks enabled the study of free discontinuity problems and evolution of interfaces in PDEs, influencing the treatment of BV solutions in elliptic and parabolic settings. In the 1990s and early 2000s, Luigi Ambrosio advanced BV theory by addressing relaxation in the calculus of variations, proving integral representations for relaxed functionals on BV spaces and establishing compactness results for sequences with bounded variation. Ambrosio's work extended BV concepts to , providing tools for non-Euclidean geometries. A comprehensive synthesis appears in the 2000 monograph by Ambrosio, Fusco, and Pallara, which details the structure theorem for BV functions and applications to free discontinuity problems. Post-2020 developments have further generalized BV theory. In 2025, research on —where A is a constant-coefficient differential operator—established fine properties analogous to classical BV, including slicing decompositions and rectifiability of singular sets. Extensions to , introduced in the sense of Shiba in 2025, adapt the variation measure to exponents varying with position, enhancing flexibility for non-standard growth conditions in variational integrals. Additionally, 2024–2025 works have defined BV functions in 2-normed spaces, endowing them with Banach-like structures and proving completeness, thus broadening applications to multilinear analysis.

Formal Definition and Notation

One-Dimensional BV Functions

A partition of a closed interval [a, b] is a finite ordered set of points a = x_0 < x_1 < \cdots < x_n = b. For a function f: [a, b] \to \mathbb{R} and such a partition P, the associated variation is defined as V_P(f) = \sum_{i=0}^{n-1} |f(x_{i+1}) - f(x_i)|. For example, consider the interval [0, 1] and the partition P = \{0, 0.5, 1\}; then V_P(f) = |f(0.5) - f(0)| + |f(1) - f(0.5)|. The total variation of f over [a, b] is the supremum of V_P(f) taken over all possible partitions P of [a, b], V(f; [a, b]) = \sup_P V_P(f). The function f is said to be of bounded variation on [a, b], denoted f \in BV[a, b], if V(f; [a, b]) < \infty. This condition ensures that the "oscillations" of f across the interval remain controlled, distinguishing BV functions from more irregular ones like the Weierstrass nowhere-differentiable function. Every function of bounded variation admits a Jordan decomposition, expressing it as the difference of two non-decreasing functions: f = f^+ - f^-, where f^+ and f^- are increasing. Specifically, one can define f^+(x) = V(f; [a, x]) and f^-(x) = V(f; [a, x]) - f(x), both of which are increasing on [a, b]. This decomposition is unique up to the addition of a constant to one function and subtraction of the same constant from the other. The total variation satisfies V(f; [a, b]) = V(f^+; [a, b]) + V(f^-; [a, b]). In one dimension, functions of bounded variation can be equivalently characterized as f = g + j + s, where g is an absolutely continuous function, j is a jump function capturing the discontinuous part across countably many points, and s is a singular continuous function. The absolutely continuous component g integrates an L^1 density, while the jump function j sums the sizes of discontinuities, and s accounts for continuous but non-absolutely continuous singular behavior, such as the . This representation highlights how BV functions combine smooth variation, finite jumps, and controlled singular parts.

Multi-Dimensional BV Functions

In the context of Sobolev-like spaces on an open domain \Omega \subset \mathbb{R}^n, the distributional derivative of a function u \in L^1(\Omega) is defined in the sense of distributions by \langle Du, \phi \rangle = -\int_\Omega u \operatorname{div} \phi \, dx for all test functions \phi \in C_c^1(\Omega, \mathbb{R}^n). This framework extends the classical notion of derivatives to functions that may lack smoothness, capturing both absolutely continuous and singular behaviors through measure theory. A function u \in L^1(\Omega) belongs to the space BV(\Omega) of functions of bounded variation if its distributional gradient Du is a finite Radon vector measure on \Omega, meaning that the total variation |Du|(\Omega) < \infty. The total variation of Du is quantified by |Du|(\Omega) = \sup \left\{ \left| \int_\Omega u \operatorname{div} \phi \, dx \right| : \phi \in C_c^1(\Omega, \mathbb{R}^n), \|\phi\|_{L^\infty(\Omega)} \leq 1 \right\}, which represents the supremum over all admissible test fields and serves as a seminorm on BV(\Omega). This definition generalizes the one-dimensional case, where bounded variation corresponds to finite total change in the function, but requires vector measures to handle higher-dimensional geometry. The measure Du admits a Lebesgue decomposition Du = D^a u + D^s u, where D^a u is the absolutely continuous part with respect to Lebesgue measure, given by \nabla u \, dx for an L^1 density \nabla u, and D^s u is the singular part, which includes both jump and Cantor-like components. This decomposition is unique and follows from the general Lebesgue decomposition theorem for Radon measures relative to Lebesgue measure. A fundamental compactness result states that if (u_k) \subset BV(\Omega) satisfies u_k \to u in L^1(\Omega) and \sup_k |Du_k|(\Omega) < \infty, then there exists a subsequence (u_{k_j}) such that Du_{k_j} \rightharpoonup^* Du weakly* in the space of Radon measures, where u \in BV(\Omega). Additionally, up to a further subsequence, u_{k_j} \to u almost everywhere. Weak* convergence of Du_k to Du is tested against continuous functions with compact support, reflecting the duality between bounded variation measures and the space C_c(\Omega, \mathbb{R}^n).

Locally BV Functions

In mathematical analysis, a function u is said to be locally of bounded variation (locally BV) on an open set \Omega \subset \mathbb{R}^n if u \in L^1_{\mathrm{loc}}(\Omega) and, for every compact subset K \subset \Omega, the distributional derivative Du satisfies |Du|(K) < \infty, where |Du| denotes the . This condition ensures that the variation is controlled on bounded regions, allowing the study of functions with potentially unbounded global behavior. Equivalently, u is locally BV if Du is a locally finite vector-valued , meaning |Du| assigns finite mass to every compact set in \Omega. This local finiteness contrasts with the global BV condition, where |Du|(\Omega) < \infty must hold over the entire domain; thus, every globally BV function is locally BV, but the converse fails in general, particularly on unbounded domains. A classic example of a function that is locally BV but not globally BV is u(x) = \log |x| in \mathbb{R}^2, where the distributional derivative has finite variation on compact sets but accumulates infinite variation globally due to the slow decay at infinity. Another illustrative case in one dimension is u(x) = \sin(1/x) on (0,1), extended appropriately, which exhibits bounded variation on subintervals [\epsilon, 1] for \epsilon > 0 but infinite on (0,1). The space of locally BV functions, often denoted \mathrm{BV}_{\mathrm{loc}}(\Omega) or \mathrm{LBV}(\Omega), plays a crucial role in the analysis of partial differential equations and variational problems on unbounded domains such as \mathbb{R}^n, where global BV would impose overly restrictive conditions like sublinear growth at infinity. For instance, the identity function u(x) = |x| belongs to \mathrm{BV}_{\mathrm{loc}}(\mathbb{R}^n) since |\nabla u| = 1 almost everywhere yields finite variation on bounded sets, but |Du|(\mathbb{R}^n) = \infty, excluding it from global BV. These spaces facilitate extensions of compactness and approximation results from bounded to unbounded settings, with the distributional derivative remaining a locally finite Radon measure to ensure well-posedness in applications like image processing on infinite domains.

Standard Notation

In the context of functions of bounded variation, standard notation relies on the Lebesgue space L^1(\Omega) and the space of , where \Omega \subset \mathbb{R}^n is an . The space of functions of bounded variation, denoted BV(\Omega), consists of all functions u \in L^1(\Omega) such that the distributional Du is a vector-valued with finite , i.e., |Du|(\Omega) < \infty. The associated norm is given by \|u\|_{BV(\Omega)} = \|u\|_{L^1(\Omega)} + |Du|(\Omega), where |Du| denotes the measure of Du. The total variation measure |Du| is defined as the smallest Radon measure such that |Du|(E) = \sup\left\{\left|\int_E u \operatorname{div} \phi \, dx \right| : \phi \in C_c^\infty(\Omega; \mathbb{R}^n), \|\phi\|_{L^\infty} \leq 1\right\} for every Borel set E \subset \Omega. The distributional derivative Du admits a decomposition Du = D^a u + D^j u + D^c u, where D^a u is the absolutely continuous part, D^j u is the jump part across the jump set J_u, and D^c u is the Cantor part, with the total variation satisfying |Du| = |D^a u| + |D^j u| + |D^c u|. The subspace of special functions of bounded variation, denoted SBV(\Omega), is the set of all u \in BV(\Omega) such that the Cantor part vanishes, i.e., D^c u = 0, so that Du = D^a u + D^j u. The variation seminorm, often denoted V(u, \Omega), is defined as V(u, \Omega) = |Du|(\Omega), which measures the total variation of u over \Omega.

Basic Properties

Discontinuities and Differentiability

Functions of bounded variation (BV) on an open set \Omega \subseteq \mathbb{R}^n are approximately differentiable \mathcal{L}^n-almost everywhere, where \mathcal{L}^n denotes the Lebesgue measure, and the approximate gradient \nabla u belongs to L^1(\Omega; \mathbb{R}^n). In the one-dimensional case, for u: [a, b] \to \mathbb{R} of bounded variation, classical differentiability holds Lebesgue-almost everywhere, with the derivative u' \in L^1[a, b]. This differentiability property stems from the decomposition of the distributional derivative Du into an absolutely continuous part D^a u = \nabla u \, \mathcal{L}^n and a singular part, ensuring the existence of the approximate differential at almost every point. A BV function u admits a precise representative \tilde{u}, which equals u \mathcal{L}^n-almost everywhere and is defined pointwise via approximate limits: at points of approximate continuity, \tilde{u}(x) is the approximate limit of u. This precise representative captures the fine structure of u, distinguishing regions of approximate differentiability from singularities, and is essential for analyzing pointwise behavior without altering the total variation. The discontinuities of a BV function are precisely characterized by its jump set J_u, the set of points where the approximate limits from opposite sides of the approximate normal differ. In the one-dimensional setting, these discontinuities are either removable or of jump type, with no Cantor-type singularities (i.e., no continuous points of essential discontinuity), and there are at most countably many such points. In higher dimensions, the jump set J_u is (n-1)-rectifiable and satisfies \mathcal{H}^{n-1}(J_u) < \infty, where \mathcal{H}^{n-1} is the (n-1)-dimensional Hausdorff measure, ensuring that jumps contribute finitely to the total variation. This finiteness implies that discontinuities are "negligible" in a measure-theoretic sense, concentrating on a set of codimension one.

Semi-Continuity and Compactness

One fundamental property of the space of functions of , denoted BV(Ω), is the lower semi-continuity of the total variation functional with respect to convergence in . Specifically, if {u_k} is a sequence in BV(Ω) converging to u in , then liminf_{k→∞} |D u_k|(Ω) ≥ |D u|(Ω), where |D u|(Ω) denotes the total variation of the distributional derivative D u over the domain Ω. This result ensures that the total variation acts as a relaxed functional in variational problems, preserving minimizers under weak convergence. The proof relies on the decomposition of the derivative measures. Since u_k → u in L¹(Ω), the absolutely continuous parts converge appropriately, while the singular parts are controlled via the weak* compactness of bounded Radon measures. In particular, the measures D u_k converge weak* to D u as Radon measures on Ω, and by the lower semi-continuity of the total mass under weak* convergence for positive measures, the liminf of the total variations follows directly. More precisely, up to a subsequence, |D u_k| weak* converges to some positive Radon measure μ ≥ |D u|, implying liminf |D u_k|(Ω) ≥ μ(Ω) ≥ |D u|(Ω). This argument extends the classical lower semi-continuity for smooth functions to the BV setting. A key consequence of this semi-continuity is a compactness theorem for BV(Ω): any bounded set in BV(Ω), equipped with the norm ||u||{BV} = ||u||{L¹} + |D u|(Ω), is precompact in L¹(Ω). In one dimension, this is Helly's selection theorem, which states that from a sequence of functions of bounded variation on an interval, one can extract a subsequence converging pointwise almost everywhere to a function of bounded variation. The generalization to higher dimensions follows from the weak* compactness of the unit ball in the space of Radon measures with finite total mass, combined with the semi-continuity result. Specifically, for {u_k} ⊂ BV(Ω) with sup_k ||u_k||{BV} < ∞, a subsequence u{k_j} converges in L¹(Ω) to some u ∈ BV(Ω) with |D u|(Ω) ≤ liminf |D u_{k_j}|(Ω) ≤ sup |D u_k|(Ω). This embedding BV(Ω) ↪ L¹(Ω) is compact, making BV spaces suitable for approximation and existence theorems in calculus of variations.

Space Structure and Completeness

The space of functions of , denoted BV(Ω), on an open set Ω ⊂ ℝⁿ is a Banach space when equipped with the norm \|u\|_{\mathrm{BV}(\Omega)} = \|u\|_{L^1(\Omega)} + |\mathrm{D}u|(\Omega), where |\mathrm{D}u| denotes the total variation measure of the distributional derivative \mathrm{D}u. This norm topology induces strong convergence in BV(Ω), with respect to which the space is complete. BV(Ω) continuously embeds into L¹(Ω), but the space is not separable, even when Ω is bounded; for instance, the uncountable family of characteristic functions of disjoint subintervals with fixed positive length yields an uncountable set where pairwise distances are bounded below by a positive constant. The dual space (BV(Ω))^* contains L^∞(Ω), which acts via integration against functions in BV(Ω), as well as certain Radon measures that account for the variation component. More precisely, elements of the dual can be represented by divergence-measure vector fields whose divergence is a Radon measure. For locally bounded variation functions, BV_loc(Ω) consists of those u ∈ L¹_loc(Ω) such that |Du|(K) < ∞ for every compact K ⊂ Ω; this space carries a natural metric structure, defined via an exhaustion of Ω by compact sets {Ω_k}, as d(u,v) = \sum_{k=1}^\infty 2^{-k} \min\{ \|u-v\|_{\mathrm{BV}(\Omega_k)}, 1 \}, making BV_loc(Ω) a complete metric space.

Chain Rule and Composition

In the context of functions of (BV), the chain rule provides a formula for the distributional derivative of compositions, extending classical differentiation rules to nonsmooth settings. For a Lipschitz continuous function g: \mathbb{R} \to \mathbb{R} and a locally BV function u: \Omega \to \mathbb{R} on an open interval \Omega \subseteq \mathbb{R}, the composition g \circ u is also locally BV, and its distributional derivative satisfies D(g \circ u) = g'(u) \, Du in the sense of measures, where g' is extended constantly on sets where it is undefined (e.g., at points of nondifferentiability of g). This holds almost everywhere with respect to the total variation measure |Du|, accounting for the absolutely continuous, jump, and Cantor parts of Du. A key consequence is the control on the total variation of the composition: for any compact subinterval K \subset \Omega, |D(g \circ u)|(K) \leq \operatorname{Lip}(g) \, |Du|(K), where \operatorname{Lip}(g) denotes the Lipschitz constant of g. This inequality follows directly from the chain rule and the boundedness of g', ensuring that the BV seminorm is preserved up to the Lipschitz factor, which is crucial for stability in approximation schemes and compactness arguments in local BV spaces. Extensions of the chain rule to non-Lipschitz outer functions g require additional growth or regularity conditions to ensure g \circ u remains BV. For instance, if g is C^1 with at most linear growth (i.e., |g'(t)| \leq C(1 + |t|) for some constant C), and u is locally BV, then D(g \circ u) = g'(u) \, Du holds in the measure sense, provided the singular parts of Du are controlled by the growth of g'. Such conditions prevent explosion of the total variation due to unbounded derivatives, as established in generalized formulations that approximate g by Lipschitz truncations. These results find applications in the analysis of nonsmooth compositions within local BV spaces, particularly in solving quasilinear PDEs where the right-hand side involves composite nonlinearities. For example, in conservation laws with flux functions that are non-Lipschitz but satisfy growth bounds, the chain rule enables the derivation of entropy inequalities and stability estimates for entropy solutions in BV. This framework supports the study of shock formations and interface evolutions in one-dimensional hyperbolic systems, where local BV regularity preserves the structure of discontinuities under composition.

Examples and Illustrations

Canonical Examples

Monotone functions provide a fundamental class of functions of bounded variation. For an increasing function f on a closed interval [a, b], the total variation V(f; [a, b]) equals the net change f(b) - f(a), which is finite since f is bounded on the compact interval. Similarly, for a decreasing function, V(f; [a, b]) = |f(b) - f(a)|. Absolutely continuous functions form another core example within the space of bounded variation functions. If f is absolutely continuous on [a, b], then it is differentiable almost everywhere, and its total variation is given by V(f; [a, b]) = \int_a^b |f'(x)| \, dx, which is finite. Smooth functions with bounded derivatives also belong to this class, as they are continuously differentiable and thus absolutely continuous. For a C^1 function f on [a, b] with \|f'\|_\infty < \infty, the total variation satisfies V(f; [a, b]) = \int_a^b |f'(x)| \, dx \leq \|f'\|_\infty (b - a). Step functions with finitely many jumps exemplify functions of bounded variation that are discontinuous. For a step function constant on subintervals of a partition of [a, b] with jumps at finitely many points, the total variation equals the sum of the absolute sizes of the jumps. A canonical instance is the H(x), defined as H(x) = 0 for x < 0 and H(x) = 1 for x \geq 0, which has total variation V(H; \mathbb{R}) = 1 due to its single jump of size 1 at x = 0. The Cantor function (or devil's staircase) is a canonical example of a singular continuous function of bounded variation. It is constant on the intervals removed in the construction of the Cantor set, increases only on the Cantor set (of Lebesgue measure zero), and has total variation 1 on [0,1], yet its derivative is zero almost everywhere.

Counterexamples

A classic example of a continuous function that fails to be of bounded variation is f(x) = \sin(1/x) for x \in (0,1], extended by f(0) = 0 if desired, though the behavior near 0 dominates. This function exhibits infinite oscillations as x \to 0^+, causing its total variation V(f; (0,1]) = \infty. To see this, consider the partition points x_n = 1/(n\pi + \pi/2) for n = 0,1,2,\dots, where \sin(1/x_n) = (-1)^n. The successive differences satisfy |f(x_n) - f(x_{n+1})| = 2, and since there are infinitely many such points accumulating at 0, the variation over these points alone sums to infinity, yielding a lower bound V(f) \geq \sum 2 = \infty. In contrast, the function g(x) = x^2 \sin(1/x) for x \in (0,1] with g(0) = 0 is continuous and of on [0,1], despite similar infinite oscillations near 0. This serves as a counterexample to the intuition that any highly oscillatory function must have unbounded variation. To compute the total variation, note that g is absolutely continuous (hence BV) with V(g) = \int_0^1 |g'(x)| \, dx, where g'(x) = 2x \sin(1/x) - \cos(1/x). The term |2x \sin(1/x)| \leq 2x integrates to 1, while \int_0^1 |\cos(1/x)| \, dx = \int_1^\infty |\cos u| / u^2 \, du < \infty by integration by parts or comparison to \sum \int_{n\pi}^{(n+1)\pi} 1/u^2 \, du \sim \sum 1/n^2 < \infty. Alternatively, using the partition at extrema x_n = 1/(n\pi + \pi/2), the differences |g(x_n) - g(x_{n+1})| \approx 2 x_n^2 \sim 2/(n\pi)^2, and \sum_n 2/(n\pi)^2 = (2/\pi^2) \sum 1/n^2 = (2/\pi^2)(\pi^2/6) = 1/3 < \infty, providing a finite lower bound that confirms boundedness when combined with monotonicity on [x_0,1]. The Weierstrass function w(x) = \sum_{n=0}^\infty a^n \cos(b^n \pi x) with $0 < a < 1 and ab > 1 + 3\pi/2 is another prominent : it is continuous everywhere but nowhere differentiable on \mathbb{R}, and thus not of bounded variation on any interval, as BV functions are differentiable . The infinite "wiggles" at all scales lead to unbounded variation, analogous to the fractal nature preventing finite for its graph. In higher dimensions, consider the characteristic function \chi_E of a domain E \subset \mathbb{R}^2 bounded by a like the , which has \log 4 / \log 3 > 1 and infinite . Such \chi_E fails to be BV because the total variation of its distributional gradient |D \chi_E|(\mathbb{R}^2) equals the perimeter of E, which is infinite; BV functions in \mathbb{R}^d require the distributional derivative to be a finite .

Generalizations and Extensions

Weighted and Variable-Exponent BV

Weighted bounded variation spaces generalize the classical BV spaces by incorporating a positive weight function w to account for non-uniform metrics in the domain \Omega \subset \mathbb{R}^n. A function u \in L^1(\Omega) belongs to the weighted BV space \mathrm{BV}_w(\Omega) if the weighted total variation is finite, defined as |D u|_w(\Omega) = \int_\Omega w \, d|D u| < \infty, where |D u| is the total variation measure of the distributional derivative D u. This formulation, introduced by Baldi, allows the weight w > 0 to belong to a subclass of Muckenhoupt's A_1 weights to ensure desirable properties. The full norm on \mathrm{BV}_w(\Omega) is typically \|u\|_{L^1(\Omega; w)} + |D u|_w(\Omega), where \|u\|_{L^1(\Omega; w)} = \int_\Omega |u| w \, dx. Equivalently, under suitable conditions on w, the weighted variation can be expressed in dual form as |D u|_w(\Omega) = \sup \left\{ \int_\Omega u \operatorname{div} \phi \, dx : \phi \in C_c^1(\Omega; \mathbb{R}^n), \, |\phi| \leq w \right\}. This space recovers the standard BV space when w \equiv 1. Recent investigations have confirmed that \mathrm{BV}_w(\Omega) is a Banach space, complete under the given norm, provided w is lower semicontinuous and satisfies appropriate admissibility conditions like belonging to the A_1^* class. Variable-exponent bounded variation spaces further extend this framework by allowing the exponent in the variation to depend on the spatial variable x \in \Omega, denoted as \mathrm{BV}^{p(\cdot)}(\Omega) where p: \Omega \to (1, \infty) is a satisfying standard log-Hölder continuity for well-posedness. For functions on intervals, the p(\cdot)-variation in the sense is defined as the supremum over tagged partitions \pi = \{ (t_{j-1}, t_j, x_j) \}_{j=1}^n of [a, b] with tags x_j \in [t_{j-1}, t_j]: V_{p(\cdot)}(f; [a, b]) = \sup_\pi \sum_{j=1}^n \left| f(t_j) - f(t_{j-1}) \right|^{p(x_j)}, with the space consisting of functions where this supremum is finite, adjusted via the modular \rho_{p(\cdot)}(f) = \inf \{ \lambda > 0 : V_{p(\cdot)}(f / \lambda; [a, b]) \leq 1 \} < \infty to form a Banach space under the Luxemburg norm \|f\|_{p(\cdot)} + V_{p(\cdot)}(f). In higher dimensions, analogous definitions use rectifiable currents or approximate the variation with variable p(x)-norms on the derivative measure. These spaces are complete under the adapted modular norm, forming a Banach function space useful for non-homogeneous problems. A notable recent development incorporates variable exponents into higher-order variations, specifically the second bounded variation in the Shiba sense. Introduced in 2025, the space \Lambda \mathrm{BV}^{2, p(\cdot)}[a, b] consists of functions f: [a, b] \to \mathbb{R} where the second p(\cdot)-variation is finite, generalizing Shiba's original construction by replacing constant p with variable p(x). The variation is given by the supremum over partitions \pi \in \Pi_3[a, b] (partitions with at least three points) of V_{\Lambda, p(\cdot)}^2(f; [a, b]) = \sup_\pi \sum_i \left( |Q(f; x_{i-1}, x_i)|^ {p(x_i)} + |Q(f; x_i, x_{i+1})|^{p(x_{i+1})} \right)^{1/p(x_i)}, where Q(f; \alpha, \beta) = f(\beta) - f(\alpha) captures second-order differences via consecutive first differences, with the supremum ensuring boundedness. The norm is \|f\|_\infty + V_{\Lambda, p(\cdot)}^2(f; [a, b]). This space is complete, forming a Banach algebra under pointwise multiplication, extending classical results for constant exponents and addressing variable growth in applications requiring higher regularity.

Special BV Classes

Special functions of bounded variation, denoted SBV(Ω), form a subclass of BV(Ω) where the distributional derivative decomposes as D u = \nabla u \, \mathcal{L}^n + D^j u, with the Cantor part D^c u = 0, ensuring no diffuse singularities. This structure confines the singular contributions to jumps, making SBV particularly suitable for modeling problems with concentrated discontinuities, such as free discontinuity sets in variational models. The jump set J_u, defined as the set where u exhibits a clear jump discontinuity, is (n-1)-rectifiable with finite (n-1)-dimensional Hausdorff measure, \mathcal{H}^{n-1}(J_u) < \infty. On J_u, the jump \llbracket u \rrbracket is well-defined and \mathcal{H}^{n-1}-measurable, allowing precise description of the discontinuity locus. Vol'pert's 1967 theorem establishes that BV functions, and thus those in SBV, are approximately continuous \mathcal{L}^n-almost everywhere, with the approximate gradient coinciding with the precise representative where it exists. For vector-valued fields, the space of bounded deformation, BD(Ω), extends the BV framework to displacement fields u: \Omega \to \mathbb{R}^n where the symmetric gradient \varepsilon(u)_{ij} = \frac{1}{2} (\partial_j u_i + \partial_i u_j) has components that are Radon measures of bounded total variation. Unlike full BV on the gradient, BD focuses on the deformation tensor, capturing symmetric strains relevant to linearized elasticity and plasticity without requiring boundedness of the full distributional derivative. This space admits a decomposition analogous to BV, with absolutely continuous, jump, and Cantor parts, but applied to the strain measure.

Discrete BV and Sequences

In the discrete setting, a real sequence \{a_n\}_{n=1}^\infty is said to be of bounded variation, denoted a \in \mathrm{bv}, if its total variation V(a) = \sum_{n=1}^\infty |a_{n+1} - a_n| < \infty. This definition extends the continuous analog by considering the supremum of sums of absolute differences over finite consecutive segments, ensuring the series of first differences converges absolutely. Sequences in \mathrm{bv} are necessarily convergent and bounded, forming a Banach space under the norm \|a\|_{\mathrm{bv}} = |a_1| + V(a). The space \mathrm{bv} is isometrically isomorphic to \ell^1 via the forward difference operator \Delta a = (\Delta a_n)_{n=1}^\infty, where \Delta a_n = a_{n+1} - a_n, as any such sequence satisfies a_n = a_1 + \sum_{k=1}^{n-1} \Delta a_k and the norm aligns with the \ell^1-norm of the differences (up to the constant term). This equivalence highlights that bounded variation sequences are precisely those whose first differences belong to \ell^1, providing a direct link to summable increments. Subspaces like \mathrm{bv}_0 = \{a \in \mathrm{bv} : a_1 = 0\} correspond exactly to the image of \ell^1 under partial summation, preserving the isometric structure. In numerical schemes, discrete BV concepts arise in finite difference discretizations of total variation functionals, where the sum of absolute differences approximates the continuous semi-norm to ensure stability and convergence in optimization algorithms. For time series analysis, BV sequences model data with sparse changes, such as abrupt shifts in trends, enabling efficient online forecasting algorithms that achieve near-optimal cumulative squared error bounds of \tilde{O}(n^{1/3} C_n^{2/3}) for variation bounded by C_n. Recent extensions from 2024–2025 unify discrete BV variations, such as Waterman \Lambda-variation and Chanturia classes, within ideal spaces on \mathbb{N}, characterizing compact embeddings between these spaces via properties like summability ideals and exhaustion conditions. For instance, compactness in Waterman-type spaces holds when \sum b_k = o(\sum a_k) for generating sequences, while Chanturia classes require g(n) = o(h(n)), all framed through ideal inclusions. This ideal perspective provides a general framework for analyzing boundedness and convergence in generalized discrete BV spaces.

Vector-Valued and Measure-Theoretic BV

The concept of functions of bounded variation extends naturally from the scalar case to vector-valued functions. For a function u: \Omega \to \mathbb{R}^m, where \Omega \subset \mathbb{R}^n is an open set, the space BV(\Omega, \mathbb{R}^m) consists of those u \in L^1(\Omega, \mathbb{R}^m) whose distributional derivative Du is a matrix-valued Radon measure of finite total variation. When m=1, this reduces to the classical scalar BV space. The total variation measure |Du| is defined such that |Du|(\Omega) = \sup \left\{ \int_\Omega \operatorname{trace}(Du \, \phi) : \phi \in C_c^\infty(\Omega, \mathbb{R}^{n \times m}), \|\phi\|_\infty \leq 1 \right\}, where \|\phi\|_\infty denotes the essential supremum of the operator norm of \phi(x), and the trace is taken with respect to the standard inner product on matrices. The BV norm is then given by \|u\|_{BV} = \|u\|_{L^1} + |Du|(\Omega), making BV(\Omega, \mathbb{R}^m) a Banach space. In the measure-theoretic framework, functions of bounded variation are intimately connected to vector measures of bounded total variation. A signed Radon measure \mu on \Omega is of bounded variation if its total variation |\mu|(\Omega) < \infty, where |\mu| is the unique positive Radon measure satisfying |\mu|(E) = \sup \left\{ \sum |\mu(E_i)| : \{E_i\} \text{ partition of } E \right\} for Borel sets E \subset \Omega. For vector-valued measures \mu: \mathcal{B}(\Omega) \to \mathbb{R}^m, the total variation is similarly defined as |\mu|(E) = \sup \left\{ \sum_{i} \|\mu(E_i)\|_{\mathbb{R}^m} : \{E_i\} \text{ partition of } E \right\}, using the Euclidean norm on \mathbb{R}^m. The space M(\Omega, \mathbb{R}^m) of all such vector measures with finite total variation, equipped with the norm \|\mu\|_M = |\mu|(\Omega), forms a Banach space that serves as the dual of C_0(\Omega, \mathbb{R}^m). In this setting, the derivative Du of a BV function belongs to M(\Omega, \mathbb{R}^{m \times n}), and the total variation |Du| captures the "size" of the distributional gradient. Recent advancements have generalized vector-valued BV spaces to broader contexts, such as functions taking values in more abstract structures or defined on non-standard domains. In particular, a 2025 contribution introduces a new space of generalized vector-valued functions of bounded variation, extending classical definitions to handle discontinuities and singularities in a unified manner while preserving key compactness and approximation properties. These generalizations maintain the Banach space structure of the underlying measures, ensuring duality with continuous functions on the respective spaces.

Applications

Mathematical Analysis and PDEs

In mathematical analysis, functions of bounded variation (BV) play a crucial role in the study of variational problems and partial differential equations (PDEs), particularly where classical smoothness assumptions fail and discontinuities or jumps are inherent. BV spaces extend by allowing weak derivatives to be measures rather than functions, enabling the treatment of problems involving free discontinuities or interfaces. This framework is essential for establishing existence, regularity, and compactness results in the calculus of variations and PDE theory, where the total variation seminorm controls the "size" of singularities. A key application arises in the relaxation of non-convex integral functionals with linear growth at infinity, where the natural space for relaxed minimizers is . For a functional I(u) = \int_\Omega f(x, u, \nabla u) \, dx defined on , non-convexity in the gradient variable can lead to lack of lower semicontinuity under weak convergence; relaxation yields an equivalent lower semicontinuous functional on BV, often expressed via the of f. Specifically, if f grows linearly in |\nabla u|, the relaxed functional incorporates the |Du| as \int_\Omega f^{**}(x, u, \nabla u) \, dx + \int_\Omega g(x, u) \, d|Du|^s, where f^{**} is the convexification and g accounts for singular parts. This result ensures the existence of minimizers in BV for problems like or phase transitions, as established in foundational relaxation theorems. Trace theorems for BV functions provide boundary continuity essential for Dirichlet problems and elliptic PDEs. For a bounded open set \Omega \subset \mathbb{R}^n with Lipschitz boundary, the trace operator T: BV(\Omega) \to L^1(\partial \Omega) is well-defined, linear, and continuous with respect to the BV norm, satisfying \|T u\|_{L^1(\partial \Omega)} \leq C \|u\|_{BV(\Omega)} for some constant C depending on \Omega. Moreover, this operator is surjective: every v \in L^1(\partial \Omega) extends to a BV function in \Omega with controlled total variation. These properties, extending Gagliardo's traces for , facilitate the analysis of boundary value problems in BV, such as those arising in elliptic regularity theory. In the context of the Plateau problem, BV functions underpin the parametrization and regularity of minimal surfaces spanning a given boundary curve. Post-1967 developments, building on Almgren's existence theorems for area-minimizing currents, leverage BV to study graphical representations and the structure of singular sets in higher codimensions. For instance, minimal surfaces can be approximated by graphs of BV functions whose total variation bounds the area, enabling proofs of density of regular points and partial regularity via monotonicity formulas adapted to BV estimates. This connection highlights BV's role in resolving the Plateau problem beyond smooth cases, as detailed in treatments linking varifold theory to BV parametrizations. For hyperbolic PDEs, particularly scalar conservation laws \partial_t u + \div F(u) = 0 in \mathbb{R}^n \times (0,T), BV provides a framework for entropy solutions when initial data u_0 \in BV. Such solutions remain in BV_{\mathrm{loc}} and satisfy the Kružkov entropy condition: for any convex entropy-entropy flux pair (\eta, q), \partial_t \eta(u) + \div q(u) \leq 0 holds in the distributional sense, with the total variation decaying according to the one-sided Lipschitz condition on the flux. In multiple dimensions, existence is proved via DiPerna-Lions regularization, mollifying the equation and passing to the limit using renormalized solutions for transport equations with BV velocities, ensuring entropy dissipation concentrates on jump sets. This yields uniqueness and stability in L^1, with BV compactness implying convergence to entropy solutions.

Image Processing and Computer Vision

In image processing, the Rudin-Osher-Fatemi (ROF) model, introduced in 1992, serves as a foundational approach for denoising by minimizing the total variation of the image while preserving edges. The model solves the optimization problem \min_u \int_\Omega |\nabla u| \, dx + \frac{\lambda}{2} \int_\Omega (u - f)^2 \, dx, where u is the denoised image, f is the noisy input, \Omega is the image domain, and \lambda > 0 balances smoothness and data fidelity. This formulation leverages the space to penalize oscillations while allowing discontinuities, effectively removing without blurring sharp boundaries, as demonstrated in applications to images. Total variation minimization extends to image inpainting, where missing regions are filled by solving \min_u \| \nabla u \|_{TV} + \lambda \int_{\Omega \setminus D} (u - f)^2 \, dx subject to u = f on the known domain D^c, promoting smooth extensions across gaps while respecting BV regularity. This method excels in texture synthesis and scratch removal, with numerical solvability via primal-dual algorithms. For segmentation, total variation underpins models like the Chan-Vese active contour framework, which approximates the Mumford-Shah functional through level sets: \min_{\phi,c_1,c_2} \int_\Omega |\nabla H(\phi)| + \lambda \int_\Omega (f - c_1 H(\phi) - c_2 (1 - H(\phi)))^2 \, dx, where H is the Heaviside function and c_1, c_2 are regional means; this partitions images into piecewise constant regions, robust to intensity inhomogeneities. Vector-valued extensions of BV spaces enable color image restoration by generalizing the total variation norm to \| \nabla \mathbf{u} \|_{TV} = \int_\Omega \sqrt{ \sum_{j=1}^3 |\nabla u^j|^2 } \, dx for \mathbf{u} = (u^1, u^2, u^3) in RGB channels, coupling components to avoid color artifacts. Seminal work in 2007 analyzed such models for super-resolution and denoising, proving existence via lower in vector spaces. Recent extensions (2015–2025) incorporate anisotropic or structure-tensor-based norms, improving handling of chromatic edges in hyperspectral tasks. From 2020 onward, has integrated with BV principles, notably through neural networks approximating BV functions. A 2025 study demonstrates the approximation capabilities of single-hidden-layer ReLU networks for BV functions on the unit circle, enabling efficient representation of discontinuous signals in vision tasks like . For example, a 2018 study on electrical capacitance (ECT) for multiphase flows used deep autoencoders and reported relative image errors of 10.88% compared to 22.40% for the Landweber iteration method in reconstructing permittivity distributions. These hybrids leverage BV's edge-preserving properties within learned priors for robust inverse problems.

Physics and Engineering

In fracture mechanics, functions of bounded variation (BV) provide a mathematical for modeling in brittle materials, where paths are represented as discontinuities with finite perimeter to capture the energy associated with fracture surfaces. This approach is central to variational models of brittle fracture, such as the nonlinear Griffith model, which combines in the uncracked domain with a term proportional to the length, ensuring the existence of minimizers through Γ-convergence to linear elastic limits. By incorporating BV spaces, these models handle free discontinuity problems without prescribing paths a priori, allowing for and in arbitrary dimensions while maintaining bounded for the displacement field. In , BV functions are employed in rate-independent models to describe transitions and associated phenomena, such as those in shape memory alloys, where the evolution of microstructures involves jumps in state variables and energy measured by . Seminal formulations treat the system as a generalized gradient flow in metric spaces, with solutions of bounded variation in time capturing irreversible transformations and loops through energetic and variational inequalities. For scalar dissipation models, the BV norm quantifies the interfacial energy during changes, enabling the analysis of low- reversible martensitic transformations as limits of viscous approximations, akin to Stefan problems with kinetic undercooling. BV solutions arise in chemical kinetics through reaction-diffusion equations, particularly for systems exhibiting propagating interfaces or fronts, where the total variation bounds the number and strength of discontinuities in concentration profiles. In such models, BV regularity ensures the existence of entropy-admissible weak solutions for Neumann boundary conditions, facilitating the generation and propagation of phase boundaries in multi-species reactions without classical smoothness assumptions. This framework is crucial for analyzing traveling waves in autocatalytic reactions, where the BV structure controls the sharpness of transition layers and prevents unphysical oscillations. In applications, BV constraints regularize control problems by penalizing excessive switching in signals, as in where bounded variation controls mitigate chattering phenomena in variable structure systems, ensuring stable trajectories with finite jumps. For instance, in of hybrid systems, BV regularization promotes smooth approximations to bang-bang controls while preserving . In , post-1967 extensions by Vol'pert enable the definition of products in BV spaces for laws, allowing entropy solutions to the Euler equations with shocks, where the measures the strength of discontinuities across wave fronts. This is essential for modeling compressible flows with finite propagation speed and bounded .