Harmonic Mean Density Fusion

• Harmonic Mean Density (HMD) fusion is an interpolation-based probabilistic fusion rule that integrates multiple local density functions via a reciprocal formulation, ensuring consistency and robustness.
• The method leverages the power mean at exponent -1 to suppress spurious peaks and avoid overconfidence, providing mathematically principled error reductions in Gaussian and non-Gaussian contexts.
• HMD fusion is widely applied in distributed tracking and multisensor inference, yielding tighter covariances and improved RMSE performance compared to naive and covariance intersection methods.

Harmonic Mean Density (HMD) fusion is an interpolation-based probabilistic fusion rule designed for robust combination of multiple local probability density functions (pdfs)—notably in distributed and track fusion scenarios where unknown cross-correlations or mismatches between local densities make conventional fusion methods inconsistent or overly conservative. HMD fusion is grounded in the power mean (Hölder mean) formalism at exponent $\alpha = -1$, and features a reciprocal formulation that prioritizes suppression of spurious or over-confident peaks, avoids double-counting of priors, and produces tight and consistent fused posteriors in both Gaussian and non-Gaussian contexts. Its mathematically principled foundations, computational tractability for Gaussian and Gaussian mixture models, and empirically demonstrated error reductions have led to its adoption and further development in distributed estimation, consensus tracking, and multisensor inference architectures (Sharma et al., 2024, Dubey et al., 5 Nov 2025, Koliander et al., 2022, Sharma et al., 2024).

1. Mathematical Formulation

Given $N$ local densities $f_i(x)$ with weights $w_i$ summing to one, the harmonic mean density is defined as

$$f_{HMD}(x) = \left( \sum_{i=1}^N w_i\, f_i(x)^{-1} \right)^{-1}$$

with normalization constant $$Z_{HMD} = \int \left( \sum_{i=1}^N w_i\,f_i(x)^{-1} \right)^{-1} dx$$, yielding the normalized fused density. For the two-density, two-weight case (with $\omega \in [0,1]$),

$$f_{HMD}(x) = \frac{f_1(x) f_2(x)}{(1-\omega) f_1(x) + \omega f_2(x)} \Bigg/ \int \frac{f_1(u) f_2(u)}{(1-\omega) f_1(u) + \omega f_2(u)} du$$

For more than two densities, a recursive structure or generalized weighted sum is used; for three densities with weights $(\omega_1, \omega_2, \omega_3)$, recursion is realized as

$$\mathcal{M}^h_{\omega_1,\omega_2,\omega_3}(f_1, f_2, f_3) = \mathcal{M}^h_{\alpha,\omega_3}\left(\mathcal{M}^h_{\omega_1,\omega_2}(f_1,f_2),\,f_3\right)$$

with $N$0. This reciprocal structure imposes “veto” behavior: the fused density vanishes wherever any local density vanishes, strongly penalizing isolated high-likelihood spikes.

2. Theoretical Properties

HMD fusion arises from multiple foundational perspectives:

• Axiomatic Characterization: HMD is the unique member of the power-mean ($N$1-mean) family at $N$2 satisfying joint likelihood commutation (Likelihood Principle), weak setwise function property (WSFP), and unanimity preservation, but not strong linearity (Koliander et al., 2022).
• Optimization Principle: It uniquely minimizes the weighted sum of Pearson $N$3 divergences,

$N$4

over normalized, nonnegative densities $N$5, yielding $N$6 as a variational solution (Sharma et al., 2024, Koliander et al., 2022).

• Supra-Bayesian Interpretation: In supra-Bayesian pooling of posteriors for linear-Gaussian models under certain priors, HMD fusion emerges as the consensus estimate for uncorrelated or weakly correlated noise models (Koliander et al., 2022).
• Support and Shape: HMD satisfies the tight bound

$N$7

for all $N$8, preserving the support intersection of all input densities and producing robustly conservative interpolation.

• Divergence and Consistency: For all $N$9, $f_i(x)$0, i.e., the fused density is always closer (in KL sense) to each local density than any two local densities are to each other. The fused density is nondecreasing in any argument, monotonic in all inputs, and consistent (in the sense of Normalized Estimation Error Squared tracking) on empirical benchmarks (Sharma et al., 2024).

3. Extension to Gaussian and Mixture Densities

For univariate or multivariate Gaussian inputs, HMD fusion supports closed-form implementation with mean and covariance updates:

Given $f_i(x)$1, the denominator mixture $f_i(x)$2 is approximated by a single Gaussian matched in moments, with

$f_i(x)$3

$f_i(x)$4

The fused covariance and mean are given by

$f_i(x)$5

$f_i(x)$6

Generalization to $f_i(x)$7 densities yields analogous, recursively computed expressions, with explicit formulas available via leave-one-out Gaussian products and moment matching (Dubey et al., 5 Nov 2025).

For Gaussian mixture models (GMMs), pairwise products of mixture components generate a new mixture which is divided by a moment-matched denominator Gaussian, followed by pruning and merging steps to control mixture cardinality (Sharma et al., 2024, Sharma et al., 2024).

4. Algorithmic Implementation

Efficient and scalable implementations target both closed-form (Gaussian) and sampling-based (mixture/nonparametric) scenarios:

• Gaussian Case: Sequentially compute, for each node, moment-matched denominator and apply Sherman–Morrison or blockwise inversion to obtain fused covariance/mean as above. Recommended weights can be uniform or computed by optimizing a criterion over symmetric KL divergence between the fused and local densities (Dubey et al., 5 Nov 2025).
• Mixture Densities: For each fusion step, compute all pairwise products for the numerator, approximate the denominator as a single Gaussian, divide, and normalize. Merge/prune mixture components to handle complexity growth.
• Sampling-based (HMD-S algorithm): For arbitrary densities, perform self-normalized importance sampling using local densities as proposal distributions. For each sample, compute weights proportional to the ratio of the other densities to the denominator, then normalize. This approach requires no additional proposal and generalizes naturally to GMM inputs (Sharma et al., 2024).
• Consensus/Distributed Fusion: Each agent computes and broadcasts local estimates and covariance; global weights may be optimized at each time step via numerical minimization of a KL-divergence–based cost function (Dubey et al., 5 Nov 2025).

5. Comparative Analysis: HMD, Naive Fusion, and Covariance Intersection

Method	Consistency	Covariance Tightness	Computational Structure
Naive/Kalman Fusion	Over-confident	Often singular	Product of densities
Covariance Intersection	Consistent, conservative	Loose (deflated)	Geometric mean
Harmonic Mean Density	Consistent, robust	Tighter but valid	Reciprocal (harmonic) mean

• Naive Fusion applies direct multiplication of posteriors, leading to over-confident covariances due to double-counting of common prior information. For Gaussians: $f_i(x)$8, which often leads to divergence (Sharma et al., 2024).
• Covariance Intersection (CI)/Geometric Mean Density (GMD), $f_i(x)$9, provides a consistent fusion strategy under unknown correlations but is recognized to be conservative, with expanded covariances: $w_i$0.
• HMD explicitly subtracts an information-theoretic denominator term and produces $w_i$1, resulting in covariances that are always at least as tight as those of CI or the inverse covariance intersection (ICI), but with consistency validated in both NEES and RMSE metrics (Sharma et al., 2024, Sharma et al., 2024).

6. Applications and Empirical Performance

Distributed target tracking, sensor consensus, and multi-sensor scenario benchmarks demonstrate that HMD fusion offers superior trade-offs in real-world tasks:

• 3D and 2D Target Tracking: In both nearly constant velocity (NCV) and maneuvering multi-model (IMM) regimes, HMD outperforms CI, often showing 20–50% RMSE improvement and maintaining NEES tightly within 95% confidence bands (Sharma et al., 2024).
• Underwater Consensus Tracking: Weighted HMD fusion, with KL-divergence–optimized weights, achieves minimum average RMSE and drives track divergence to zero, sustaining ANEES within confidence bounds (Dubey et al., 5 Nov 2025).
• Real-world Multi-Target Fusion: Across multi-radar and multi-track fusion with clutter and non-Gaussian densities, HMD-GA yields consistently lower positional RMSE and NEES near theoretical means, dominating both CI and naive fusion for most targets (Sharma et al., 2024).

7. Limitations, Computational Considerations, and Future Directions

HMD fusion’s primary limitations include the need for accurate approximation of the denominator mixture in non-Gaussian settings, and the complex optimization of fusion weights when targeting further covariance control or entropy minimization. Normalization requires one-dimensional or higher-order integration, though sampling-based implementations relieve the need for explicit proposal distributions. A plausible implication is that advances in mixture reduction and moment-matching techniques, or adaptive and decentralized weight optimization, will further enhance HMD’s applicability to high-dimensional and resource-constrained systems (Sharma et al., 2024, Dubey et al., 5 Nov 2025, Sharma et al., 2024).

---

References:

• "On Pooling-Based Track Fusion Strategies: Harmonic Mean Density" (Sharma et al., 2024)
• "Consensus Tracking of an Underwater Vehicle Using Weighted Harmonic Mean Density" (Dubey et al., 5 Nov 2025)
• "Fusion of Probability Density Functions" (Koliander et al., 2022)
• "Harmonic Mean Density Fusion in Distributed Tracking: Performance and Comparison" (Sharma et al., 2024)