---
title: MF Stochastic Bounded Real Lemma
url: https://www.emergentmind.com/topics/mean-field-stochastic-bounded-real-lemma-mf-sbrl
type: topic
---

# MF Stochastic Bounded Real Lemma

Searching arXiv for the cited MF-SBRL papers and closely related formulations.
Mean-Field Stochastic Bounded Real Lemma (MF-SBRL) denotes a bounded-real characterization for stochastic systems whose dynamics depend not only on the state and disturbance, but also on their expectations. In the finite-horizon discrete-time setting, the MF-SBRL gives sufficient conditions under which the induced $\ell_2$ gain from an exogenous disturbance to a controlled output is strictly less than a prescribed level $\gamma$, using coupled backward Riccati-type recursions and positivity constraints [1607.00451]. In a later continuous-time formulation with affine terms, the same bounded-real question is expressed through coupled differential Riccati equations, together with BSDE and ODE components needed to accommodate affine drift and diffusion terms, and the bounded-real condition is stated as an equivalence [2507.19809].

## 1. Mean-field structure and system-theoretic setting

In the discrete-time formulation, the state $x_k \in \mathbb{R}^n$, control input $u_k \in \mathbb{R}^q$, exogenous disturbance $v_k \in \mathbb{R}^\ell$, and performance output $z_k \in \mathbb{R}^m$ evolve over a finite horizon $k \in \{0,1,\dots,K\}$ according to
\[
\begin{aligned}
x_{k+1} &= A(k)\,x_k + \bar{A}(k)\,\mathbb{E}[x_k] + B(k)\,v_k + \bar{B}(k)\,\mathbb{E}[v_k] \\
&\quad + \big(C(k)\,x_k + \bar{C}(k)\,\mathbb{E}[x_k] + D(k)\,v_k + \bar{D}(k)\,\mathbb{E}[v_k]\big)\,w_k + F_1(k)\,u_k, \\
z_k &= \Phi(k)\,x_k,
\end{aligned}
\]
with deterministic $x_0$ and scalar stochastic process $w_k$ satisfying
\[
\mathbb{E}[w_k]=0,\qquad \mathbb{E}[w_s w_t]=\delta_{st}.
\]
The disturbance $v_k$ belongs to $\ell_2$ over the horizon, is $\mathcal{F}_{k-1}$-measurable, and is uncorrelated with $w_k$ [1607.00451].

A standard notational device is the introduction of aggregated matrices
\[
\mathcal{A}(k):=A(k)+\bar{A}(k),\quad
\mathcal{B}(k):=B(k)+\bar{B}(k),\quad
\mathcal{C}(k):=C(k)+\bar{C}(k),\quad
\mathcal{D}(k):=D(k)+\bar{D}(k).
\]
These aggregated quantities encode the mean-field coupling and appear systematically in the bounded-real recursions.

The defining structural feature is the decomposition into centered and mean components. Taking expectations yields
\[
\mathbb{E}[x_{k+1}] = \mathcal{A}(k)\,\mathbb{E}[x_k] + \mathcal{B}(k)\,\mathbb{E}[v_k],
\]
while the fluctuation dynamics are written in terms of $x_k-\mathbb{E}[x_k]$ and $v_k-\mathbb{E}[v_k]$. This split is not cosmetic: it is the mechanism through which the mean-field coupling enters the bounded-real analysis and produces an additional matrix recursion absent from the classical non-mean-field stochastic BRL [1607.00451].

## 2. Bounded-real performance and induced-gain interpretation

For the discrete-time MF-SBRL, the disturbance-to-output map is defined with zero initial state:
\[
\mathcal{L}_K : \ell_2(N_K,\mathbb{R}^\ell)\to \ell_2(N_K,\mathbb{R}^m),\qquad
\mathcal{L}_K(v(\cdot)) := \Phi(k)\,x_k(0,v(\cdot)).
\]
Its induced $\ell_2$ gain is
\[
\|\mathcal{L}_K\| := \sup_{v\neq 0}\frac{\|z(\cdot)\|_{\ell_2(N_K,\mathbb{R}^m)}}{\|v(\cdot)\|_{\ell_2(N_K,\mathbb{R}^\ell)}}.
\]
The bounded real property with level $\gamma>0$ is the strict inequality
\[
\|\mathcal{L}_K\|<\gamma,
\]
equivalently expressed, when $x_0=0$, as
\[
\sum_{k=0}^K \mathbb{E}\big[\|z_k\|^2\big]
\le
\gamma^2 \sum_{k=0}^K \mathbb{E}\big[\|v_k\|^2\big].
\]
Here the mean-field terms enter through the state trajectory $x_k$, hence indirectly through the output $z_k=\Phi(k)x_k$ [1607.00451].

The induced-gain interpretation is the stochastic mean-field analogue of the deterministic and stochastic bounded-real viewpoint: the lemma certifies a dissipative inequality with attenuation level $\gamma$. The nontriviality lies in the fact that the energy balance must be handled separately on the centered and mean subspaces.

In the same discrete-time framework, the $H_2$ component used later for synthesis is built from the generic output
\[
z_k=\Phi_1(k)\,x_k+\Upsilon_1(k)\,u_k,
\]
with finite-horizon cost
\[
J_K(x_0,u) := \sum_{k=0}^K \mathbb{E}\big[\|z_k\|^2\big].
\]
This places the MF-SBRL inside a mixed $H_2/H_\infty$ architecture, where the $H_\infty$ side controls the disturbance attenuation and the $H_2$ side determines quadratic optimality [1607.00451].

## 3. Discrete-time MF-SBRL and its Riccati-type recursions

The discrete-time MF-SBRL is formulated through two coupled backward difference equations for symmetric matrices $P(k),Q(k)\in\mathcal{H}_n(\mathbb{R})$. At each time $k$, the paper defines
\[
\begin{aligned}
L\big(P(k+1)\big) &:= \mathcal{A}(k)^\top P(k+1)\,\mathcal{A}(k) + \mathcal{C}(k)^\top P(k+1)\,\mathcal{C}(k) - \Phi(k)^\top\Phi(k), \\
G\big(P(k+1)\big) &:= \mathcal{A}(k)^\top P(k+1)\,\mathcal{B}(k) + \mathcal{C}(k)^\top P(k+1)\,\mathcal{D}(k), \\
H\big(P(k+1)\big) &:= \gamma^2 I_\ell + \mathcal{B}(k)^\top P(k+1)\,\mathcal{B}(k) + \mathcal{D}(k)^\top P(k+1)\,\mathcal{D}(k),
\end{aligned}
\]
and similarly
\[
\begin{aligned}
\tilde{L}\big(P(k+1),Q(k+1)\big) &:= \mathcal{A}(k)^\top Q(k+1)\,\mathcal{A}(k) + \mathcal{C}(k)^\top P(k+1)\,\mathcal{C}(k) - \Phi(k)^\top\Phi(k), \\
\tilde{G}\big(P(k+1),Q(k+1)\big) &:= \mathcal{A}(k)^\top Q(k+1)\,\mathcal{B}(k) + \mathcal{C}(k)^\top P(k+1)\,\mathcal{D}(k), \\
\tilde{H}\big(P(k+1),Q(k+1)\big) &:= \gamma^2 I_\ell + \mathcal{B}(k)^\top Q(k+1)\,\mathcal{B}(k) + \mathcal{D}(k)^\top P(k+1)\,\mathcal{D}(k).
\end{aligned}
\]

The lemma states that if the constrained backward equations
\[
\begin{aligned}
P(k) &= L\big(P(k+1)\big) - G\big(P(k+1)\big)\,H\big(P(k+1)\big)^{-1}\,G\big(P(k+1)\big)^\top, \\
Q(k) &= \tilde{L}\big(P(k+1),Q(k+1)\big) - \tilde{G}\big(P(k+1),Q(k+1)\big)\,\tilde{H}\big(P(k+1),Q(k+1)\big)^{-1}\,\tilde{G}\big(P(k+1),Q(k+1)\big)^\top,
\end{aligned}
\]
with terminal conditions
\[
P(K+1)=0,\qquad Q(K+1)=0,
\]
satisfy
\[
H\big(P(k+1)\big)\succ 0,\qquad
\tilde{H}\big(P(k+1),Q(k+1)\big)\succ 0
\quad \text{for all } k,
\]
and admit a unique solution $(P_1(k),Q_1(k))$, then
\[
\|\mathcal{L}_K\|<\gamma
\qquad\text{and}\qquad
Q_1(k)\preceq 0,\ \forall k
\]
[1607.00451].

The structural novelty relative to the classical stochastic bounded real lemma is the second recursion $Q(k)$. It tracks the expected-state channel and is coupled to $P(k)$ through $\tilde{L}$, $\tilde{G}$, and $\tilde{H}$. A common misconception is that MF-SBRL is obtained by a straightforward replacement of state matrices by their mean-field aggregates. The discrete-time result shows otherwise: the mean dynamics require an additional backward equation and additional positivity constraints.

Another important point is the status of the result. In this discrete-time formulation, the MF-SBRL is a sufficient condition for $\|\mathcal{L}_K\|<\gamma$, not a necessary one. The reason given is that solvability of the backward recursions depends on positivity constraints that must hold at every step, and these constraints need not be implied solely by the induced-gain inequality [1607.00451].

## 4. Role in mean-field $H_2/H_\infty$ synthesis

The discrete-time MF-SBRL is used as the disturbance-attenuation component of a finite-horizon mixed $H_2/H_\infty$ design. On the $H_2$ side, the mean-field stochastic linear-quadratic problem is posed for
\[
\begin{aligned}
x_{k+1} &= A_1(k)\,x_k + \bar{A}_1(k)\,\mathbb{E}[x_k] + F_1(k)\,u_k \\
&\quad + \big(B_1(k)\,x_k + \bar{B}_1(k)\,\mathbb{E}[x_k]\big)\,w_k, \\
z_k &= \Phi_1(k)\,x_k + \Upsilon_1(k)\,u_k,
\end{aligned}
\]
with cost
\[
J_K(x_0,u) = \sum_{k=0}^K \mathbb{E}\big[\|z_k\|^2\big].
\]
A backward recursion obtained by completion of squares yields solvability conditions and the optimal control $u_k^*$, together with $Q_1(k)\succeq 0$ in that LQ setting [1607.00451].

For the mixed mean-field $H_2/H_\infty$ problem, the controller and the worst-case disturbance are taken in state-feedback form:
\[
u^*(k)=U(k)\,x_k+\bar{U}(k)\,\mathbb{E}[x_k],\qquad
v^*(k)=V(k)\,x_k+\bar{V}(k)\,\mathbb{E}[x_k].
\]
The existence result is expressed through four coupled matrix-valued backward equations: a controller-side Riccati pair $(P_1,Q_1)$ and a disturbance-side Riccati pair $(\tilde{P}_1,\tilde{Q}_1)$, all with zero terminal conditions and positivity conditions on
\[
H,\ \tilde{H},\ H_1,\ \tilde{H}_1.
\]
The worst-case disturbance gains are
\[
V(k) = -\,H\big(P_1(k+1)\big)^{-1}\,G_u\big(P_1(k+1)\big)^\top,\qquad
\bar{V}(k) = \tilde{H}\big(P_1(k+1),Q_1(k+1)\big)^{-1}\,\tilde{G}_u\big(P_1(k+1),Q_1(k+1)\big)^\top,
\]
whereas the controller gains are
\[
U(k) = -\,H_1\big(\tilde{P}_1(k+1)\big)^{-1}\,G_v\big(\tilde{P}_1(k+1)\big)^\top,\qquad
\bar{U}(k) = \tilde{H}_1\big(\tilde{P}_1(k+1),\tilde{Q}_1(k+1)\big)^{-1}\,\tilde{G}_v\big(\tilde{P}_1(k+1),\tilde{Q}_1(k+1)\big)^\top.
\]
Under solvability of these equations, the closed loop satisfies the $H_\infty$ bound $\|\mathcal{L}_K\|<\gamma$ and the $H_2$ optimality requirement against the worst-case disturbance $v^*$ [1607.00451].

This synthesis mechanism shows how MF-SBRL functions as more than a verification device. It is a constructive ingredient in a two-player finite-horizon design: first the disturbance channel is shaped through a bounded-real argument, and then the controller channel is optimized through an LQ argument, with both steps coupled by the mean-field terms.

## 5. Proof methodology and numerical realization

The proof strategy in the discrete-time setting combines mean-field decomposition, telescoping identities, and completion of squares. The first step rewrites the system in terms of the centered and mean quantities
\[
(x_k-\mathbb{E}[x_k],\ \mathbb{E}[x_k]),\qquad
(v_k-\mathbb{E}[v_k],\ \mathbb{E}[v_k]).
\]
The second step uses key telescoping identities for quadratic forms in $\tilde{x}_k:=x_k-\mathbb{E}[x_k]$ and $\mathbb{E}[x_k]$, exploiting the independence and zero-mean property of $w_k$ to eliminate cross terms. The resulting summed identity expresses
\[
\sum_{k=0}^K \mathbb{E}\big[\gamma^2\|v_k\|^2-\|z_k\|^2\big]
\]
as a quadratic form in the centered and mean variables, plus boundary terms depending on $P(K+1)$ and $Q(K+1)$. Completion of squares then produces the Riccati-type recursions and the positivity conditions ensuring nonnegativity of the residual terms [1607.00451].

The sign property $Q(k)\preceq 0$ in the MF-SBRL is not incidental. It reflects the fact that the expected-state channel enters the dissipation argument with its own storage contribution. This suggests that the mean dynamics are not a perturbative correction to the fluctuation dynamics, but a parallel channel requiring a separate energy accounting.

The synthesis proof repeats the same logic twice: once for the $H_\infty$ disturbance channel and once for the $H_2$ controller channel. That duplication is what yields four coupled backward equations rather than a single Riccati recursion.

A backward recursive algorithm is given for implementation. Starting from
\[
P_1(K+1)=Q_1(K+1)=0,\qquad \tilde{P}_1(K+1)=\tilde{Q}_1(K+1)=0,
\]
one computes
\[
H,\ \tilde{H},\ H_1,\ \tilde{H}_1,
\]
checks positive definiteness, inverts the matrices if possible, computes the terminal gains
\[
U(K),\ \bar{U}(K),\ V(K),\ \bar{V}(K),
\]
updates the Riccati variables at time $K$, and repeats backward for $k=K-1,K-2,\dots,0$. The numerical remarks recommend symmetric matrix factorizations such as Cholesky for positivity checks and inversion, exploitation of block structure, and describe the overall cost as $O(K)$ Riccati-like steps involving $n\times n$, $\ell\times \ell$, and $q\times q$ matrix operations [1607.00451].

The reported illustrative case is a two-step example with $K=2$ and $\gamma=0.8$, where the recursive algorithm computes $U(k)$, $\bar{U}(k)$, $V(k)$, $\bar{V}(k)$ and the matrices $P_1(k)$, $Q_1(k)$, $\tilde{P}_1(k)$, $\tilde{Q}_1(k)$, verifies the positivity conditions, and yields the mixed $H_2/H_\infty$ controller. Detailed numerical values are tabulated in the source [1607.00451].

## 6. Relation to classical stochastic BRL and intrinsic limitations

When the mean-field terms vanish,
\[
\bar{A}=\bar{B}=\bar{C}=\bar{D}=0,
\]
the discrete-time MF-SBRL reduces to the standard stochastic bounded real lemma for discrete-time systems with state- and disturbance-dependent noise. Under appropriate assumptions, the non-mean-field problem admits equivalence-type results, whereas the mean-field version in the discrete-time formulation remains sufficient only [1607.00451].

The main structural difference is therefore not merely the presence of expectation terms in the state equation, but the appearance of the additional recursion $Q(k)$ and the associated matrices $\tilde{L},\tilde{G},\tilde{H}$. These encode the coupling between the expected dynamics and the fluctuation dynamics. A plausible implication is that mean-field bounded-real analysis should be viewed as a genuinely two-layer Riccati problem, even when the original state equation appears close to a classical stochastic model.

Several limitations are explicit in the discrete-time result. The horizon is finite, the system is discrete-time, the multiplicative noise is scalar, the disturbance is exogenous and uncorrelated with the noise, and no explicit detectability or stabilizability assumptions are formulated beyond solvability and positivity of the backward recursions. Moreover, the Riccati-like equations may fail the intermediate positivity tests even when $\|\mathcal{L}_K\|<\gamma$ holds, which is precisely why the theorem is framed as sufficient rather than necessary [1607.00451].

A further misconception is that open-loop and closed-loop solvability should coincide automatically in mean-field stochastic $H_2/H_\infty$ problems. The later continuous-time affine theory explicitly distinguishes these notions and states that, in that setting, open-loop solvability does not imply closed-loop solvability because mean dynamics and affine terms require additional structural equations beyond the bounded-real Riccati pair [2507.19809].

## 7. Continuous-time affine extension and later MF-SBRL formulations

A continuous-time mean-field stochastic system with affine terms is given on $[t,T]$ by
\[
\begin{cases}
\displaystyle dX(s) = \big\{A_{1}(s)X(s) + \bar{A}_{1}(s)\mathbb{E}[X(s)] + B_{1}(s)u(s) + \bar{B}_{1}(s)\mathbb{E}[u(s)] \\
\qquad\qquad\quad + C_{1}(s)v(s) + \bar{C}_{1}(s)\mathbb{E}[v(s)] + b(s)\big\}\,ds \\
\quad\ \ + \big\{A_{2}(s)X(s) + \bar{A}_{2}(s)\mathbb{E}[X(s)] + B_{2}(s)u(s) + \bar{B}_{2}(s)\mathbb{E}[u(s)] \\
\qquad\qquad\quad + C_{2}(s)v(s) + \bar{C}_{2}(s)\mathbb{E}[v(s)] + \sigma(s)\big\}\,dW(s), \\
\displaystyle z(s) = \begin{pmatrix} Q(s)X(s) \\ N_1(s)u(s) \end{pmatrix},\qquad X(t)=\xi,
\end{cases}
\]
with deterministic bounded coefficients, affine terms $b(\cdot),\sigma(\cdot)\in L^2_{\mathbb F}(t,T;\mathbb R^n)$, and $N_1(s)^\top N_1(s)=I$ [2507.19809].

For the $H_\infty$ analysis, one sets $u\equiv 0$, $X(t)=0$, $b\equiv 0$, $\sigma\equiv 0$, and defines
\[
(Lv)(s):=z(s)=Q(s)X(s),
\qquad
\|L\|:=\sup_{v\neq 0}\frac{\|Lv\|}{\|v\|},
\]
where
\[
\|Lv\|=\Big(\mathbb{E}\int_t^T X(s)^\top Q(s)^\top Q(s)X(s)\,ds\Big)^{1/2},
\qquad
\|v\|=\Big(\mathbb{E}\int_t^T |v(s)|^2\,ds\Big)^{1/2}.
\]
The bounded-real inequality is
\[
\mathbb{E}\!\int_t^T |z(s)|^2\,ds \le \gamma^2\,\mathbb{E}\!\int_t^T |v(s)|^2\,ds
\iff \|L\|\le \gamma
\]
[2507.19809].

The continuous-time MF-SBRL introduces the aggregated matrices
\[
\widetilde{A}_i:=A_i+\bar{A}_i,\qquad \widetilde{C}_i:=C_i+\bar{C}_i,
\]
and the functions
\[
\Lambda^\gamma(P):=\gamma^2 I + C_2^\top P C_2,\qquad
\bar{\Lambda}^\gamma(P):=\gamma^2 I + \widetilde{C}_2^\top P \widetilde{C}_2,
\]
\[
\Phi(P):=PC_1 + A_2^\top P C_2,\qquad
\bar{\Phi}(P,\Pi):=\Pi\,\widetilde{C}_1 + \widetilde{A}_2^\top P \widetilde{C}_2.
\]
It then states that, for a given $\gamma>0$, the following are equivalent:

1. $\|L\|<\gamma$.
2. There exist continuous symmetric matrix solutions $P(\cdot),\Pi(\cdot)\in C([t,T];\mathbb S^n)$ to the coupled differential Riccati equations with $P(T)=0$, $\Pi(T)=0$,
\[
\begin{cases}
\dot{P} + PA_1 + A_1^\top P + A_2^\top P A_2 - Q^\top Q \\
\qquad - \big(PC_1 + A_2^\top PC_2\big)\big(\gamma^2 I + C_2^\top PC_2\big)^{-1}\big(C_1^\top P + C_2^\top P A_2\big) = 0,\\[0.3em]
\dot{\Pi} + \Pi\widetilde{A}_1 + \widetilde{A}_1^\top \Pi + \widetilde{A}_2^\top P \widetilde{A}_2 - Q^\top Q \\
\qquad - \big(\Pi\widetilde{C}_1 + \widetilde{A}_2^\top P \widetilde{C}_2\big)\big(\gamma^2 I + \widetilde{C}_2^\top P \widetilde{C}_2\big)^{-1}\big(\widetilde{C}_1^\top \Pi + \widetilde{C}_2^\top P \widetilde{A}_2\big) = 0,\\[0.2em]
\Lambda^\gamma(P)\ge \delta I,\qquad \bar{\Lambda}^\gamma(P)\ge \delta I\quad\text{for some }\delta>0.
\end{cases}
\]
This is a substantive contrast with the discrete-time 2016 formulation: the continuous-time affine result is stated as an equivalence rather than only a sufficiency condition [2507.19809].

Affine terms require an additional BSDE/ODE layer. If $(\eta,\zeta)$ and $\bar\eta$ solve the specified backward stochastic and deterministic equations, the minimizing disturbance is
\[
\begin{aligned}
v^*(s) &= - \Lambda^\gamma(P)^{-1} \Phi(P)^\top \big(X(s)-\mathbb{E}[X(s)]\big)
         - \bar{\Lambda}^\gamma(P)^{-1} \bar{\Phi}(P,\Pi)^\top \mathbb{E}[X(s)] \\
&\quad - \Lambda^\gamma(P)^{-1} \varphi(s)
         - \bar{\Lambda}^\gamma(P)^{-1} \bar{\varphi}(s),
\end{aligned}
\]
with
\[
\varphi := C_1^\top \eta + C_2^\top (P\sigma + \zeta),\qquad
\bar{\varphi} := \widetilde{C}_1^\top \mathbb{E}[X] + \widetilde{C}_2^\top (P\mathbb{E}[\sigma] + \mathbb{E}[\zeta]).
\]
The associated optimal cost is given explicitly in terms of $P,\Pi,\eta,\zeta,\bar\eta$ [2507.19809].

For the full joint $H_2/H_\infty$ problem, open-loop equilibrium is characterized by an MF-FBSDE together with stationarity conditions
\[
\begin{cases}
C_1^\top Y^*_1 + C_2^\top Z^*_1 + \bar{C}_1^\top \mathbb{E}[Y^*_1] + \bar{C}_2^\top \mathbb{E}[Z^*_1] + \gamma^2 v^* = 0,\\
B_1^\top Y^*_2 + B_2^\top Z^*_2 + \bar{B}_1^\top \mathbb{E}[Y^*_2] + \bar{B}_2^\top \mathbb{E}[Z^*_2] + u^* = 0,
\end{cases}
\]
while closed-loop synthesis is derived from four coupled CDREs $(P_1,\Pi_1,P_2,\Pi_2)$, two BSDEs, and two ODEs. The resulting feedback laws are
\[
u^*(s)=U(s)X^*(s)+\bar U(s)\mathbb E[X^*(s)]+U_0(s),\qquad
v^*(s)=V(s)X^*(s)+\bar V(s)\mathbb E[X^*(s)]+V_0(s),
\]
with
\[
V = -\varLambda^{-1}\Upsilon^\top,\qquad
\bar{V} = -\bar{\varLambda}^{-1}\bar{\Upsilon}^\top + \varLambda^{-1}\Upsilon^\top,
\]
\[
U = -\Theta^{-1}\Sigma^\top,\qquad
\bar{U} = -\bar{\Theta}^{-1}\bar{\Sigma}^\top + \Theta^{-1}\Sigma^\top.
\]
Under solvability of the four CDREs, the BSDEs, and the ODEs, the closed loop achieves the bounded-real property and $H_2$ optimality [2507.19809].

Taken together, the discrete-time and continuous-time formulations show a consistent pattern: mean-field coupling forces a separation between centered and mean dynamics; bounded-real certification is expressed through coupled Riccati objects rather than a single equation; and once affine or mixed $H_2/H_\infty$ features are introduced, additional backward equations become necessary. The precise status of the lemma, however, depends on the model class: the discrete-time finite-horizon version is sufficient, whereas the continuous-time affine version is formulated as necessary and sufficient under its stated conditions.

Source: https://www.emergentmind.com/topics/mean-field-stochastic-bounded-real-lemma-mf-sbrl