我的模型系统:一个各向同性扩散粒子,它在各种扩散系数(D1 <-> D2 <-> D3 <-> ...)之间进行随机切换。
由于可以根据高斯分布对沿该假设粒子轨迹的位移进行建模,因此使用高斯+模型选择的混合似乎很自然,以提取有关不同“状态”数量或存在的扩散系数的信息,这将表现为混合物中的不同成分。
似乎有很多代码可以在协方差矩阵不受约束的 GMM 上执行 EM。然而,在我的特定应用中,各向同性扩散意味着我的矩阵不仅是对角线,而且对角线的所有分量对于每个混合分量都是相等的,这意味着在 x、y、z 方向上的扩散速率是相同的。
任何人都可以就这种特殊情况下的期望和最大化步骤将如何变化提供指导吗?