从隐式方程中找到概率向量

计算科学 矩阵 可能性 逆问题 矩阵方程
2021-12-21 13:07:29

我有q n维向量yi和一个矩阵B^形状的n×m. 我在找q m维向量xi这样:

  1. yi=B^xi
  2. 每个向量xi表示概率分布(每个条目都是非负的,并且L1标准是 1)

相关问题是:

  • 矢量图yi很吵
  • 矩阵B^可以更准确地估计,但可能相当稀疏

如何解决这样的问题?

1个回答

我有一些想法,你可以试试。我会修一个i, 简化符号。如果(xi)i=1q是独立的,这应该没有问题。

  1. 最小二乘:解决最小化问题
    minBxy22,
    s.t.: x0,x1=1.
    这可以用任何标准的非线性优化库来完成。
  2. 贝叶斯方法:在贝叶斯方法中,您将尝试计算向量的后验度量。你假设你有一些关于向量的先验信息,然后将给出的信息合并到这个先验信息中。您拥有的先验信息是是具有个条目的概率向量。幸运的是,有代表这种概率向量的概率度量,例如Dirichlet 分布假设我们上有一个表示先验知识对于贝叶斯过程,您需要构建一个可能性。如果我们假设噪声数据是由 xyxmp0Rm
    yB(xtrue)+η,
    其中是高斯噪声。可能性是 您现在可以使用贝叶斯规则,例如重要性采样或马尔可夫链蒙特卡罗来计算后验测度。但为此,您可能应该查看维基百科页面:贝叶斯推理马尔可夫链蒙特卡洛后验平均值将为您最优估计量,测量本身将量化向量周围的不确定性。ηN(0,Γ)
    L(y|x)=exp((1/2)Γ1/2(yBx)22).
    L2xtruex