三硬币模型的变分贝叶斯EM解法

题目：假设有3枚硬币，分别记做A，B，C。这些硬币正面出现的概率分别是$\pi$,$p$和$q$。进行如下掷硬币实验：先掷硬币A，根据其结果选出硬币B或C，正面选硬币B，反面选硬币C；然后投掷选中的硬币，出现正面记作1，反面记作0；独立地重复$n$次$(n=10)$，结果为1111110000…
我们只能观察投掷硬币的结果，而不知其过程，估计这三个参数$\pi$,$p$和$q$。

我们假设A出现正面的概率$\pi$满足一个贝塔分布$Be(\alpha)$，其中$\alpha=(\alpha_1,\alpha_2)$，B和C出现正面的概率$p$和$q$分别满足$Be(\beta_i),i=1,2$，其中

$$\beta= \left [ \begin{matrix} \beta_{11}&\beta_{12}\\ \beta_{21}&\beta_{22} \end{matrix} \right ]$$

那么，对于任意一列抽样数据$X=(x_1,x_2,...,x_n)$，对应存在一列 Z=(