---
title: 'MAPS: Mode-Aware Probabilistic Scheduling'
url: https://www.emergentmind.com/topics/mode-aware-probabilistic-scheduling-maps
type: topic
---

# MAPS: Mode-Aware Probabilistic Scheduling

Mode-Aware Probabilistic Scheduling (MAPS) is an adaptive control framework that couples an Interacting Multiple Model (IMM) estimator with a Linear Parameter-Varying (LPV) control law and uses updated mode probabilities directly as interpolation weights for online gain synthesis. In the formulation introduced for a brushed DC motor with varying friction, MAPS treats friction-induced operating conditions as discrete modes, performs Bayesian mode inference in real time, and schedules controller gains probabilistically rather than through a measured deterministic scheduling variable. The resulting architecture is intended for friction-aware adaptive control in uncertain, time-varying environments and was validated in a Hardware-in-the-Loop Simulation (HILS) setting on a QUBE-Servo 2 platform [2509.12695].

## 1. Defining mechanism

The defining feature of MAPS is the replacement of traditional deterministic LPV scheduling weights with posterior mode probabilities generated online by an IMM estimator. In the DC motor formulation, the controller does not rely on explicit identification of the viscous friction coefficient at runtime. Instead, it infers the probability of each friction mode and forms a convex combination of precomputed vertex gains:
$$
K_k = \sum_{i=1}^M \mu_k^{(i)} K^{[i]}.
$$
For regulation, the control law is
$$
u_k = -K_k x_k,
$$
and for tracking, which was used in the experiments,
$$
e_k = x_{\mathrm{ref},k} - x_k,\qquad u_k = K_k e_k.
$$
The sign convention was chosen to match the hardware controller [2509.12695].

This construction tightly couples estimation and control. Classical LPV scheduling uses a measured or computed parameter $\rho_k$ and associated weights $\xi_i(\rho_k)$, whereas MAPS uses $\mu_k^{(i)}$ inferred from the IMM. The distinction is operationally important: the framework adapts to latent, unmeasured changes such as friction variation without explicit friction model identification. It also differs from hard switching among mode-specific controllers, because gains are blended continuously rather than selected by winner-take-all logic. This suggests a design emphasis on smooth adaptation during mode transitions rather than abrupt controller replacement.

## 2. System model and LPV representation

The reference application is a brushed DC motor with armature voltage input $v$ and states angular position $\theta$, angular velocity $\omega$, and armature current $i$. The coupled electrical and mechanical dynamics are
$$
L_m \frac{di}{dt} = -R_m i - K_e \omega + v,\qquad
J_{eq} \frac{d\omega}{dt} = K_t i - b\,\omega,
$$
where $L_m$, $R_m$, $K_e$, $K_t$, and $J_{eq}$ are motor parameters and $b$ is the viscous friction coefficient. With
$$
x := \begin{bmatrix}\theta\\ \omega\\ i\end{bmatrix},\qquad
u := v,\qquad
y := \theta,
$$
the continuous-time state-space model is
$$
\dot{x} = A(b)x + Bu,\qquad y = Cx,
$$
with
$$
A(b)=
\begin{bmatrix}
0 & 1 & 0\\
0 & -\dfrac{b}{J_{eq}} & \dfrac{K_t}{J_{eq}}\\
0 & -\dfrac{K_e}{L_m} & -\dfrac{R_m}{L_m}
\end{bmatrix},\quad
B=
\begin{bmatrix}
0\\0\\ \dfrac{1}{L_m}
\end{bmatrix},\quad
C=\begin{bmatrix}1&0&0\end{bmatrix}.
$$
Variations in $b$, for example due to temperature or load, induce distinct operating modes. The MAPS implementation uses two representative friction modes identified experimentally as $b_{\min}$ and $b_{\max}$ [2509.12695].

After forward-Euler discretization with sampling time $T$,
$$
x_{k+1} = \Phi(\rho_k)x_k + \Gamma u_k + w_k,\qquad
y_k = Hx_k + v_k,
$$
where
$$
\Phi(\rho_k)=I+T A(\rho_k),\qquad \Gamma = TB,\qquad H=C.
$$
The scheduling parameter $\rho_k$ represents the viscous friction coefficient at time $k$. In polytopic LPV form,
$$
\Phi(\rho_k)=\sum_{i=1}^{N_v}\xi_i(\rho_k)\Phi^{[i]},\qquad
\sum_{i=1}^{N_v}\xi_i(\rho_k)=1,\qquad
\xi_i(\rho_k)\ge 0,
$$
with vertices $\Phi^{[1]}=\Phi(b_{\min})$ and $\Phi^{[2]}=\Phi(b_{\max})$ for $N_v=2$ [2509.12695].

For each mode $i$, the process and measurement noises are zero-mean Gaussian:
$$
w_k^{(i)} \sim \mathcal{N}(0,Q_i),\qquad
v_k^{(i)} \sim \mathcal{N}(0,R_i).
$$
In the HILS validation, identical covariances were used for both modes:
$$
Q_1 = Q_2 = \mathrm{diag}(10^{-6},10^{-6},10^{-6}),\qquad
R_1 = R_2 = 10^{-5}.
$$
The two-mode construction is minimal rather than exhaustive. A plausible implication is that MAPS depends on whether the chosen vertex set captures the dominant friction variation relevant to control performance.

## 3. IMM estimation and probabilistic scheduling of gains

The estimator layer uses an IMM with $M=2$ modes corresponding to $b_{\min}$ and $b_{\max}$. Mode transitions are modeled by a Markov matrix
$$
\Pi=[\pi_{ij}] \in \mathbb{R}^{M\times M},\qquad
\pi_{ij}=\mathbb{P}\{\text{mode } j \text{ at } k \mid \text{mode } i \text{ at } k-1\}.
$$
In the experiments, the self-transition probabilities were set to $0.9$ and the switching probabilities to $0.1$, so $\pi_{11}=\pi_{22}=0.9$ and $\pi_{12}=\pi_{21}=0.1$ [2509.12695].

Given previous mode probabilities $\mu_{k-1}^{(i)}$, the IMM forms one-step-ahead predicted mode probabilities
$$
p_{k|k-1}^{(j)} = \sum_{i=1}^M \pi_{ij}\mu_{k-1}^{(i)},
$$
and conditional mixing weights
$$
\mu_k^{(i|j)} =
\frac{\pi_{ij}\mu_{k-1}^{(i)}}{p_{k|k-1}^{(j)}}.
$$
These define the mixed initial state and covariance for each mode-conditioned Kalman filter. For each mode $i$, prediction and update are carried out with the model $(\Phi^{[i]},\Gamma,H,Q_i,R_i)$, using innovation $e_i(k)$, innovation covariance $S_i(k)$, and gain $K_i(k)$. The corresponding likelihood is
$$
L_i(k)=\frac{1}{\sqrt{(2\pi)^m|S_i(k)|}}
\exp\left(-\frac{1}{2}e_i(k)^\top S_i(k)^{-1}e_i(k)\right),
$$
where $m=\dim(y_k)$ [2509.12695].

Posterior mode probabilities are then computed by Bayesian normalization:
$$
\mu_k^{(i)} =
\frac{L_i(k)\,p_{k|k-1}^{(i)}}
{\sum_{\ell=1}^M L_\ell(k)\,p_{k|k-1}^{(\ell)}}.
$$
The fused IMM state estimate is
$$
x_k = \sum_{i=1}^M \mu_k^{(i)} x_k^{(i)},
$$
with a corresponding covariance fusion formula that includes the within-mode covariances and the spread of mode-conditioned estimates. MAPS uses these updated $\mu_k^{(i)}$ immediately as convex weights for control gain interpolation [2509.12695].

A common misunderstanding is to treat MAPS as an estimator of a continuous friction parameter. The formulation in fact operates through discrete friction modes and posterior mode probabilities. Another common misunderstanding is to equate probabilistic scheduling with mode switching. MAPS explicitly contrasts its convex blending with hard mode switching and attributes to the latter the risk of chattering and transient performance degradation during mode transitions.

## 4. Gain synthesis, constraints, and stability guarantees

The controller gains at the LPV vertices are synthesized by discrete-time LQR on the vertex models $(\Phi^{[i]},\Gamma)$ with infinite-horizon cost
$$
J = \sum_{k=0}^{\infty} \left(x_k^\top Q x_k + u_k^\top R u_k\right),
$$
where $Q \succeq 0$ and $R \succ 0$. The discrete algebraic Riccati equation yields $P$, and the vertex gain is
$$
K^{[i]} = \left(R + \Gamma^\top P \Gamma\right)^{-1}\Gamma^\top P \Phi^{[i]}.
$$
In the HILS implementation,
$$
Q=\mathrm{diag}(100,1,1),\qquad R=10.
$$
Voltage is limited by the QUBE-Servo 2 hardware, reference amplitudes were chosen in $[-4,4]$ V, and saturation is handled by clamping the computed $u_k$. No integral action is present, so anti-windup is not needed [2509.12695].

The stability analysis is framed through a common Lyapunov argument. The lemma stated for the vertex-closed loops is: if there exists $P \succ 0$ such that for each vertex $i$,
$$
\left(\Phi^{[i]} + \Gamma K^{[i]}\right)^\top
P
\left(\Phi^{[i]} + \Gamma K^{[i]}\right) - P \prec 0,
$$
then any convex combination
$$
\Phi^{\mathrm{cl}}_k =
\sum_{i=1}^{N_v} \mu_k^{(i)}
\left(\Phi^{[i]} + \Gamma K^{[i]}\right)
$$
is quadratically stable for all admissible $\mu_k$ on the probability simplex. The associated theorem states that with $\mu_k$ generated by the IMM and $K(\rho_k)=\sum_i \mu_k^{(i)}K^{[i]}$, the closed loop
$$
x_{k+1} =
\left(\Phi(\rho_k) + \Gamma K(\rho_k)\right)x_k
$$
is quadratically stable under the lemma’s conditions. The paper further states that, under bounded parameter estimation error and Lipschitz continuity of $\Phi(\cdot)$ and $K(\cdot)$, MAPS preserves exponential stability when $\|\hat{\rho}_k-\rho_k\|$ and $\|\rho_{k+1}-\rho_k\|$ are sufficiently small [2509.12695].

These guarantees are structurally important because they formalize the status of probability-weighted interpolation as a closed-loop design mechanism rather than a heuristic. The stability claim is not attached to a specific friction identification procedure; it is attached to convexity of the closed loop and the existence of a common Lyapunov certificate across the vertices.

## 5. Algorithmic workflow, HILS results, and implementation

The per-sample workflow is specified as: measure $z_k$; perform IMM interaction to compute $p_{k|k-1}^{(j)}$, $\mu_k^{(i|j)}$, and mixed $(x_{k-1}^{0(j)},P_{k-1}^{0(j)})$; run mode-conditioned Kalman prediction and update; evaluate likelihoods $L_i(k)$ and posterior probabilities $\mu_k^{(i)}$; fuse the state estimate; interpolate the gain; compute $e_k$ and $u_k$; apply saturation; and send $u_k$ to the plant [2509.12695].

The HILS platform consists of a QUBE-Servo 2 DC motor with optical encoder and amplifier, operated at a sampling rate of $500$ Hz with $T=2$ ms. Signals recorded were $\theta$, $\omega$, $i$, and $u$. Friction was emulated through a minimum-friction condition with no load and a maximum-friction condition with manual external load. The computational platform was MATLAB/Simulink Real-Time, and the reported per-step latency at $500$ Hz was negligible [2509.12695].

State-estimation accuracy in HILS was reported through RMSE as follows.

| State | MAPS RMSE | KF RMSE |
|---|---:|---:|
| $\theta$ [rad] | 0.0014 | 0.0038 |
| $\omega$ [rad/s] | 0.5725 | 0.9825 |
| $i$ [A] | 0.0223 | 0.0523 |

These correspond to an approximately $63\%$ reduction for $\theta$, approximately $42\%$ for $\omega$, and approximately $57\%$ for $i$ relative to the KF baseline [2509.12695].

Tracking performance was compared between MAPS-gain LQR and fixed-gain LQR. For a step input with no load, the two methods showed similar performance, with fixed-gain slightly lower IAE by approximately $1.67\%$. For a sine input with no load, MAPS improved RMSE to $0.5110$ versus $0.5348$ and IAE to $14.1458$ versus $14.7622$. For a step input with load, MAPS reduced IAE to $17.6559$ versus $18.7269$ and MAE to $0.5886$ versus $0.6243$. For a sine input with load, MAPS reduced RMSE to $0.7009$ versus $0.8356$, MAE to $0.6344$ versus $0.7485$, and IAE to $19.0336$ versus $22.4560$. Transient responses showed reduced overshoot and faster settling under friction changes with MAPS-gain LQR [2509.12695].

Implementation details reported for tuning and feasibility include $M=2$, raw IMM probabilities without optional smoothing, and complexity $O(Mn^3)$ per step for KF factorizations. With $n=3$ and $M=2$, computations were described as trivial, and the framework ran comfortably at $500$ Hz with low latency. Practical conditions include observability and controllability of both vertex models, sensitivity of adaptation speed to the choice of $\Pi$, and the usual estimator trade-off that larger $Q$ increases adaptability but may amplify noise while larger $R$ prioritizes the model over measurements [2509.12695].

## 6. Generalization, limitations, and related uses of probabilistic scheduling

Within the LPV-adaptive-control setting, MAPS is presented as applicable to systems with mode-dependent dynamics, including varying damping in mechanical systems, aerodynamic regimes in flight control, tire-road friction in vehicles, and other hybrid or LPV plants. The stated prerequisites are a set of representative vertex models and stabilizing gains for each. The reported limitations are equally specific: limited vertex coverage may reduce performance if the true dynamics fall outside the polytopic hull; rapid, frequent switching can challenge the IMM and cause lag in $\mu_k$; poorly tuned $Q/R$ or $\Pi$ can degrade estimation; and weakly observable modes can produce unreliable probabilities, in which case additional sensing or richer models may be needed [2509.12695].

The term “probabilistic scheduling” also appears in distinct but related control and decision-making contexts. In networked control systems, a probabilistic scheduling algorithm was proposed for plants sharing a communication network of limited capacity. That framework models each plant as a two-mode Markovian jump linear system with connected/closed-loop and disconnected/open-loop modes, constructs disjoint subsets of plants of cardinality equal to the network capacity, assigns stationary probabilities to those subsets, and derives necessary and sufficient LMIs for mean-square stability of the resulting homogeneous Markov chains [2204.00788]. The paper explicitly notes that it does not use the term MAPS, but its design is described as intrinsically mode-aware through MJLS characterization and probabilistic scheduling of network access.

A further development appears in fault-aware communications scheduling for robotic fleets, where an IMM-MPC framework maintains a probabilistic belief over discrete fault modes with time-inhomogeneous dynamics and optimizes a two-term objective coupling acquisition value with information gain. In that setting, mode-aware probabilistic scheduling addresses observational aliasing, such as failed contacts that cannot distinguish a lethal hardware fault from a benign communications outage under standard actions. The framework uses belief updates over discrete modes,
$$
\bar{\mu}_{k+1}^i(m') = \sum_{m \in M} \Pi_{mm'}(\tau_k^i)\mu_k^i(m),
$$
followed by Bayesian updating under action-conditioned observation likelihoods, and was reported to recover $59.8\%$ of lethal-fault spacecraft versus $9.0\%$ for binary-MPC in a rideshare scenario and $24.1\%$ versus $3.0\%$ in a megaconstellation scenario, across $200$ randomized trials [2605.09158].

These related uses clarify a terminological point. In LPV-based adaptive control, MAPS denotes probabilistic gain scheduling from inferred operating modes; in networked control, probabilistic scheduling refers to stochastic allocation of communication access; in fleet communications, it denotes receding-horizon action scheduling under mode uncertainty. The shared structure is mode awareness, probabilistic inference or probabilistic allocation, and a control or scheduling law that depends explicitly on those probabilities.

Source: https://www.emergentmind.com/topics/mode-aware-probabilistic-scheduling-maps