In this work, we propose a method for tracking multiple extended targets or unresolvable group targets in a clutter environment. First, based on the random matrix model (RMM), each target's joint kinematic–extent state is modelled as a gamma Gaussian inverse Wishart (GGIW) distribution. Considering the uncertainty of measurement origin caused by the clutters, we adopt the idea of probabilistic data association and describe the joint association event as an unknown parameter in the joint prior distribution. Then, variational Bayesian inference (VBI) is used to approximate the intractable posterior distribution. To improve practicality, we propose two lightweight schemes to reduce computational complexity. The first is clustering-based and effectively prunes joint association events. The second simplifies the variational posterior by using marginal association probabilities. Finally, we demonstrate effectiveness on simulations and real-data experiments and show that the method outperforms state-of-the-art baselines in accuracy and adaptability.