---
title: Reduced Bussgang Kalman Filter for 1-Bit Data
url: https://www.emergentmind.com/topics/reduced-bussgang-aided-kalman-filter
type: topic
---

# Reduced Bussgang Kalman Filter for 1-Bit Data

Searching arXiv for the specified papers to ground the article in the cited literature.
Reduced Bussgang-aided Kalman filter (RBKF) is a computationally efficient variant of the Bussgang-aided Kalman filter (BA-KF) for state estimation from 1-bit observations. It is designed for settings in which the observation dimension is large, often with many 1-bit analog-to-digital converters, so that the dominant computational cost of the Bussgang-based update is the inversion of a high-dimensional innovation covariance. RBKF preserves the Bussgang linearization and adaptive thresholding strategy of BA-KF, but reduces the measurement dimension by a fixed linear projection before the update. In the formulation reported in "State Estimation with 1-Bit Observations and Imperfect Models: Bussgang Meets Kalman in Neural Networks" [2507.17284], the method targets both linear and nonlinear state-space models, while its computational rationale aligns with the broader catalog of structure-exploiting Kalman filter optimizations surveyed in "On Computational Complexity Reduction Methods for Kalman Filter Extensions" [1512.03077].

## 1. Problem formulation and scope

RBKF is posed for a possibly nonlinear state-space model with state evolution
$$
x_k = f(x_{k-1}) + w_k,
$$
and pre-quantization observation
$$
u_k = y_k = h(x_k) + v_k \in \mathbb{R}^n,
$$
where \(w_k \sim N(0,Q_k)\) and \(v_k \sim N(0,R_k)\), mutually independent, and independent of \(x_0\) [2507.17284]. The filter operates after 1-bit quantization, with the analog-to-digital converter applying the elementwise sign function with threshold \(\tau_k\):
$$
r_k = Q(u_k-\tau_k), \qquad Q(\cdot)=\mathrm{sign}(\cdot)
$$
applied coordinatewise [2507.17284].

A defining feature of the method is adaptive thresholding, also described as dithering. The threshold is set to the one-step-ahead prediction of the measurement,
$$
\tau_k = \hat y_{k|k-1} = h(\hat x_{k|k-1}),
$$
so that
$$
z_k := u_k-\tau_k
$$
is zero-centered in the conditional sense \(E[z_k \mid r_{1:k-1}] = 0\) [2507.17284]. This step is central because it restores the simple zero-mean Bussgang coefficients for 1-bit quantization and avoids the more expensive non-zero-mean formulas involving Gaussian integrals [2507.17284].

RBKF itself is a model-based method. The formulation assumes that \(f(\cdot)\), \(h(\cdot)\), \(Q_k\), and \(R_k\) are known, possibly with EKF-style linearization for nonlinear models [2507.17284]. Robustness to model mismatch is limited in comparison with the accompanying Bussgang-aided KalmanNet, which was introduced in the same work for partial model knowledge [2507.17284].

## 2. Bussgang linearization under 1-bit quantization

The filter is built on the Bussgang decomposition for Gaussian inputs. If \(u \sim N(0,C_u)\) and \(y=\mathrm{sign}(u)\) is applied elementwise, then the cross-covariance is
$$
C_{yu}=E[yu^T]=\sqrt{\frac{2}{\pi}}\,D_u^{-1/2}C_u,
$$
with \(D_u:=\mathrm{diag}(C_u)\) [2507.17284]. The corresponding equivalent linearization is
$$
y \approx A u + e, \qquad A=C_{yu}C_u^{-1}=\sqrt{\frac{2}{\pi}}\,D_u^{-1/2},
$$
where \(e\) is uncorrelated with \(u\), so \(E[e\,u^T]=0\) [2507.17284].

For the 1-bit output covariance, the method uses the arcsine law. Writing
$$
\rho_u=D_u^{-1/2}C_uD_u^{-1/2},
$$
the correlation matrix of \(u\), one has
$$
C_{yy}=\frac{2}{\pi}\arcsin(\rho_u),
$$
with the \(\arcsin\) operator applied elementwise to entries in \([-1,1]\) [2507.17284].

Within the Kalman filtering context, these identities are applied to the dithered residual
$$
z_k=u_k-\hat y_{k|k-1}.
$$
For linear \(h\), or for a Jacobian-based residual in the nonlinear case,
$$
z_k = H_k(x_k-\hat x_{k|k-1}) + v_k,
$$
and the quantized observation is approximated as
$$
r_k = Q(z_k) \approx B_k z_k + \eta_k,
$$
with
$$
B_k=\sqrt{\frac{2}{\pi}}\,D_{z,k}^{-1/2}, \qquad D_{z,k}:=\mathrm{diag}(P_{k|k-1}),
$$
where \(\eta_k\) is uncorrelated with \(z_k\) [2507.17284]. The effective linear measurement model is then expressed on the state deviation \(\delta x_k:=x_k-\hat x_{k|k-1}\) as
$$
r_k \approx H_{b,k}\,\delta x_k + w_k^b,
$$
with
$$
H_{b,k}:=B_k H_k, \qquad w_k^b:=B_k v_k+\eta_k
$$
[2507.17284].

This construction is directly related to the statistical linearization viewpoint surveyed in [1512.03077], where a Jacobian-like matrix is written as
$$
J=P_{g(z)z}P_{zz}^{-1},
$$
and explicitly noted to be the Bussgang-type gain for memoryless nonlinearities when inputs are Gaussian [1512.03077]. This suggests that RBKF can be understood as a specialized, quantization-aware instance of the broader statistical-linearization family.

## 3. From BA-KF to the reduced formulation

The parent BA-KF first performs prediction in EKF style for nonlinear models:
$$
\hat x_{k|k-1}=f(\hat x_{k-1|k-1}),
$$
$$
F_k=\frac{\partial f}{\partial x}\ \text{evaluated at}\ \hat x_{k-1|k-1},
$$
$$
\Sigma_{k|k-1}=F_k\Sigma_{k-1|k-1}F_k^T+Q_k,
$$
$$
\hat y_{k|k-1}=h(\hat x_{k|k-1}),
$$
$$
H_k=\frac{\partial h}{\partial x}\ \text{evaluated at}\ \hat x_{k|k-1},
$$
and
$$
P_{k|k-1}=H_k\Sigma_{k|k-1}H_k^T+R_k
$$
[2507.17284].

Using the Bussgang approximation, BA-KF forms the output covariance
$$
S_k:=\mathrm{Cov}(r_k)=\frac{2}{\pi}\arcsin\!\left(D_{z,k}^{-1/2}P_{k|k-1}D_{z,k}^{-1/2}\right),
$$
or equivalently defines
$$
R_{b,k}:=S_k-H_{b,k}\Sigma_{k|k-1}H_{b,k}^T
$$
[2507.17284]. The Kalman-like gain is
$$
K_k:=\Sigma_{k|k-1}H_{b,k}^T S_k^{-1},
$$
and the update is
$$
\hat x_{k|k}=\hat x_{k|k-1}+K_k r_k,
$$
$$
\Sigma_{k|k}=(I-K_kH_{b,k})\Sigma_{k|k-1}
$$
[2507.17284]. Because adaptive dithering enforces \(E[r_k]=0\), no centering term is needed in the innovation [2507.17284].

RBKF modifies only the measurement side of this construction. It reduces the measurement dimension before the Bussgang update by a fixed projection
$$
A \in \mathbb{R}^{n_r \times n}, \qquad n_r = a n,\quad 0<a\le 1,
$$
for example by averaging \(L=1/a\) replicated ADCs per feature:
$$
r_k^* = A r_k
$$
[2507.17284]. The projected covariance and effective measurement matrix become
$$
S_k^* = A S_k A^T,
$$
$$
H_{b,k}^* = A H_{b,k} = (A B_k)H_k \equiv B_k^* H_k
$$
[2507.17284]. The reduced update is then
$$
K_k^*=\Sigma_{k|k-1}(H_{b,k}^*)^T(S_k^*)^{-1},
$$
$$
\hat x_{k|k}=\hat x_{k|k-1}+K_k^* r_k^*,
$$
$$
\Sigma_{k|k}=(I-K_k^*H_{b,k}^*)\Sigma_{k|k-1}
$$
[2507.17284].

The reduction is therefore not a modification of the Bussgang constants themselves, but a projection of the quantized observation, its covariance, and the effective linearized measurement operator.

## 4. Computational reduction and relation to generic KF optimization

The explicit motivation for RBKF is that the dominant cost in BA-KF is inverting \(S_k \in \mathbb{R}^{n\times n}\), which is \(O(n^3)\). This is particularly restrictive when many 1-bit ADCs are used and often \(n \gg m\) [2507.17284]. The reported per-step complexity is
- BA-KF: \(O(m^3 + n m^2 + n^3)\),
- RBKF: \(O(m^3 + (a n)m^2 + (a n)^3)\) plus the projection cost \(O(n n_r)\),
with dominant term \(O((a n)^3)=O(a^3 n^3)\) [2507.17284].

Memory is also reduced: storing \(S_k\) and its factors changes from \(O(n^2)\) to \(O((a n)^2)\) [2507.17284]. In applications with replicated measurements, the projection can be chosen as block-averaging,
$$
A=\frac{1}{L}(1_L^T\otimes I_p),
$$
when \(n=Lp\) and \(p\) is the number of physical features [2507.17284].

This reduction mechanism is closely aligned with the general principles in [1512.03077]. That tutorial emphasizes sequential measurement processing for diagonal or block-diagonal \(R_k\), grouped or partitioned measurement updates, low-rank updates via the Woodbury identity, square-root filters, and partially linear or conditionally linear reductions as standard routes for lowering Kalman-filter complexity [1512.03077]. In that framework, RBKF can be interpreted as a measurement-side dimension-reduction strategy specialized to Bussgang-effective observation models. A plausible implication is that RBKF belongs to the same structural-optimization class as grouped measurement updates and reduced-dimensional measurement transforms, although the 2025 formulation presents it directly as a projection of \(r_k\), \(S_k\), and \(H_{b,k}\) rather than through the generic notation \(g(z)=A g_n(Tz)+Hz\) used in [1512.03077].

The broader tutorial also notes that exploitation of the structure of the problem can lead to improved estimation accuracy while reducing the computational load [1512.03077]. For RBKF, the article’s own empirical claims are more specific: when projections average multiple i.i.d. replicas with identical \(R_k\), information loss can be negligible, whereas heterogeneous noise or spatially varying statistics cause mild performance loss [2507.17284].

## 5. Nonlinear integration, numerical conditioning, and implementation

RBKF supports nonlinear \(f(\cdot)\) and \(h(\cdot)\) through EKF linearizations. The state Jacobian is
$$
F_k=\frac{\partial f}{\partial x}\Big|_{x=\hat x_{k-1|k-1}},
$$
and the measurement Jacobian is
$$
H_k=\frac{\partial h}{\partial x}\Big|_{x=\hat x_{k|k-1}},
$$
after which the method proceeds through \(P_{k|k-1}=H_k\Sigma_{k|k-1}H_k^T+R_k\), the Bussgang mapping, and the RBKF update [2507.17284]. A UKF variant is stated to be possible by plugging sigma-point-based \(P_{k|k-1}\), but the reported implementation uses EKF [2507.17284].

The implementation guidance emphasizes regularization and conditioning. The recommended stabilized quantities are
$$
D_{z,k}=\mathrm{diag}(P_{k|k-1})+\epsilon I,
$$
followed by safe inverse square root evaluation; the normalized covariance
$$
\rho_k=D_{z,k}^{-1/2}P_{k|k-1}D_{z,k}^{-1/2}
$$
is clipped entrywise to \([-1+\delta,\,1-\delta]\) before applying \(\arcsin\); \(S_k\) is symmetrized and regularized by small diagonal loading; and \(S_k^*\) is regularized by \(+\epsilon I\) if needed [2507.17284]. The update covariance \(\Sigma_{k|k}\) is then symmetrized and, if necessary, adjusted by tiny diagonal loading to enforce positive semidefiniteness [2507.17284].

Square-root filtering is recommended for improved numerical robustness [2507.17284]. This recommendation is consistent with the general survey in [1512.03077], where square-root filters based on Cholesky or QR are described as numerically favorable, replacing inversions with triangular solves. The same tutorial also recommends the Joseph-stabilized covariance update
$$
P_{k|k}=(I-K_kH_k)P_{k|k-1}(I-K_kH_k)^T+K_kR_kK_k^T
$$
as preferred for stability [1512.03077]. The RBKF formulas reported in [2507.17284] use the simpler covariance update \((I-K_k^*H_{b,k}^*)\Sigma_{k|k-1}\). This suggests that, in implementations where numerical robustness is critical, the generic stabilization guidance from [1512.03077] is relevant even though the 2025 paper presents the reduced method in compact form.

The practical notes on the projection \(A\) are also explicit. Block-averaging is presented as a good default when there are \(L\) replicas per physical measurement; it preserves scale and keeps \(S_k^*\) well-conditioned. For heterogeneous noise, more sophisticated \(A\), such as weighted averaging, can improve performance at modest extra cost [2507.17284].

## 6. Empirical behavior, comparative position, and limitations

The reported experiments cover the Lorenz-Attractor model and the Michigan NCLT dataset [2507.17284]. On the Lorenz attractor with a single 1-bit ADC per feature and no reduction, the following mean-squared error results are reported:
- EKF (ideal measurements): \(-19.31\) dB,
- KalmanNet (ideal): \(-19.49\) dB,
- EKF with 1-bit: \(17.85\) dB,
- KalmanNet with 1-bit: \(12.95\) dB,
- BKF with 1-bit: \(-17.38\) dB,
- BKNet with 1-bit: \(-17.31\) dB
[2507.17284].

For multiple 1-bit ADCs per feature, the comparison between RBKF and BKF is more directly informative. With identical noise across ADCs, RBKF matches BKF MSE across noise levels, while inference time drops dramatically as the number of ADCs increases; an example given is 128 ADCs per feature, where BKF is approximately \(5.43\) s and RBKF approximately \(0.89\) s on the same setup [2507.17284]. With heterogeneous noise across ADCs, where variance is uniform in \([-20,-10]\) dB, the reported Lorenz results are:
- EKF (ideal): \(-22.41\) dB,
- BKF for \(1,8,64,128\) ADCs per feature: \(-20.53,-26.59,-31.80,-33.79\) dB,
- RBKF for the same cases: \(-20.53,-26.33,-31.45,-33.33\) dB
[2507.17284].

On Michigan NCLT with partial model knowledge and a single ADC per feature, the reported values are:
- EKF (ideal): \(33.41\) dB,
- KalmanNet (ideal): \(19.15\) dB,
- with 1-bit, EKF: \(37.79\) dB,
- KalmanNet: \(34.67\) dB,
- BKF: \(32.58\) dB,
- BKNet: \(18.62\) dB
[2507.17284]. These results are used in the source to distinguish the intended operating regimes: RBKF when models are known and efficiency is needed, BA-KF when measurement dimension is modest and runtime is not critical, and BKNet when \(f\), \(h\), \(Q\), and \(R\) are partially known or mismatched [2507.17284].

The limitations are stated in model and distributional terms. Degradation is expected under strongly heterogeneous noise among replicas, correlations not preserved by simple averaging, severe model mismatch in \(f\), \(h\), \(Q\), and \(R\), and non-Gaussian pre-quantization inputs, since Bussgang is second-order optimal for Gaussian inputs [2507.17284]. The method is also sensitive to the centering step: dithering that re-centers the residual is described as critical [2507.17284]. More generally, [1512.03077] cautions that partitioning and structural reductions must be balanced against numerical conditioning and linearization drift, especially in strongly nonlinear regimes.

In summary, RBKF is a reduced-dimensional Bussgang-linearized Kalman update for 1-bit observations. Its defining ingredients are adaptive thresholding \(\tau_k=\hat y_{k|k-1}\), the Bussgang coefficient \(B_k=\sqrt{2/\pi}\,D_{z,k}^{-1/2}\), the arcsine-law covariance \(S_k\), and a fixed projection \(A\) that compresses the quantized measurement prior to inversion [2507.17284]. Within the broader theory of Kalman filter complexity reduction, it is an instance of exploiting observation structure to replace a high-dimensional update by a reduced one, while retaining recursive state estimation under the same underlying model assumptions [1512.03077].

Source: https://www.emergentmind.com/topics/reduced-bussgang-aided-kalman-filter