【正文】 p。 以下都是針對(duì)相互獨(dú)立的事件，P(A,B)=P(B|A)*P(A)P(A,B,C)=P(C)*P(A,B|C)=P(A,C|B)*P(B)=P(B,C|A)*P(A)P(A,B,C,D)=P(D)*P(A,B|D)*P(C|A)=P(D)*P(A,B|D)*P(C|B)P(A,B|C)=P(D1,A,B|C)+P(D2,A,B|C) 預(yù)備知識(shí)首先你得對(duì)下面的公式表示認(rèn)同。 如此循環(huán)迭代，直到前后兩次參數(shù)的變化量小于某個(gè)值為止。 // 向后變量 double[][] gamma = new double[T][]。 A[1][0] = 。 t T 1。 } alpha[t][i] *= B[i][(out_seq[t])]。 t++) { double frac = 。 i++) { frac += alpha[t][i] * beta[t][i]。 i++) { double frac = 。 } A[i][j] = d離散數(shù)學(xué) / frac。 // 迭代次數(shù) while (iter 0) { // EStep updateAlpha()。 for (int i = 0。 } ()。 i 。 t++) { if (out_seq[t] == (j)) d離散數(shù)學(xué) += gamma[t][i]。 k++) // frac += xi[t][i][k]。 i M。 t T 1。 j M。 i++) { alpha[t][i] = 0。 t++) { alpha[t] = new double[M]。 i++) { A[i] = new double[M]。 // 測試用的隱藏狀態(tài)序列 int T = 32。除以 C)則會(huì)顯示該合適公式是 A^B(1amp。離散數(shù)學(xué)算法的基本思路是：隨機(jī)初始化一組參數(shù)θ(0)，根據(jù)后驗(yàn)概率Pr(Y|X。是從狀態(tài)i轉(zhuǎn)移出去的次數(shù)的期望值。 // 輸出活動(dòng)的種數(shù) double[] PI。 PI = new double[M]。 B[1][1] = 。 } } // 更新向前變量 public void updateAlpha() { for (int i = 0。 i M。 j++) { frac += alpha[t][i] * A[i][j] * B[j][(out_seq[t + 1])] * beta[t + 1][j]。 } // for(int i=0。 j M。 t T。 updateXi()。 } ()。 for (int i = 0。 updateB()。 for (int j