function?hmm=baum_welch（hmmobs）

mix=hmm.mix;?%高斯混合模型
N=hmm.N;?%HMM的狀態數
K=length（obs）;?%訓練數據樣本數
SIZE=size（obs（1）.fea2）;?%特征矢量的個數

for?loop?=?1:40
???%?----計算前向?后向概率矩陣
???for?k=1:K
?????param（k）=getparam（hmmobs（k）.fea）;
???end

???%----重估轉移概率矩陣A
???for?i=1:N-1
?????demon=0;
?????for?k=1:K??
????????tmp=param（k）.ksai（:i:）;
????????demon=demon+sum（tmp（:））;?%對時間t，j求和
?????end
?????for?j=i:i+1??
????????nom=0;
????????for?k=1:K??
????????????tmp=param（k）.ksai（:ij）;
????????????nom=nom+sum（tmp（:））;??%對時間t求和
????????end
????????hmm.trans（ij）=nom/demon;
?????end
???end

???%----重估輸出觀測值概率B
???for?j=1:N?%狀態循環
?????for?l=1:hmm.M（j）?%混合高斯的數目
????????%計算各混合成分的均值和協方差矩陣
????????nommean=zeros（1SIZE）;
????????nomvar=zeros（1SIZE）;
????????denom=0;
????????for?k=1:K??%訓練數目的循環
???????????T=size（obs（k）.fea1）;??%幀數
???????????for?t=1:T???%幀數（時間）的遍歷
?????????????x=obs（k）.fea（t:）;
?????????????nommean=nommean+param（k）.gama（tjl）*x;
?????????????nomvar=nomvar+param（k）.gama（tjl）*（x-mix（j）.mean（l:））.^2;
?????????????denom=denom+param（k）.gama（tjl）;
???????????end
????????end
????????hmm.mix（j）.mean（l:）=nommean/denom;
????????hmm.mix（j）.var（l:）=nomvar/denom;
???
????????%計算各混合成分的權重
????????nom=0;
????????denom=0;
????????for?k=1:K
??????????tmp=param（k）.gama（:jl）;
??????????nom=nom+sum（tmp（:））;
??????????tmp=param（k）.gama（:j:）;
??????????denom=denom+sum（tmp（:））;
????????end
????????hmm.mix（j）.weight（l）=nom/denom;
?????end
???end???
??????
end
??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????185??2014-06-29?11:26??CMN.m
?????文件?????????670??2018-06-05?13:28??DTWScores.m
?????文件?????1083706??2018-06-05?13:26??Vectors1.mat
?????文件?????1084076??2018-06-05?13:26??Vectors2.mat
?????文件?????1083809??2018-06-05?13:26??Vectors3.mat
?????文件????????1655??2014-08-03?16:51??baum_welch.m
?????文件?????????234??2014-06-28?22:25??deltacoeff.m
?????文件????????2236??2016-03-17?14:46??disteusq.m
?????文件?????????375??2016-03-17?14:45??enframe.m
?????文件?????????132??2014-06-27?17:19??freq2mel.m
?????文件????????2170??2009-10-28?20:47??getparam.m
?????文件???????14854??2018-06-01?13:46??iconcumt.jpg
?????文件????????1290??2018-06-02?23:56??inithmm.m
?????目錄???????????0??2018-06-02?19:34??kk\
?????文件???????64044??2018-06-02?19:33??kk\1.wav
?????文件???????64044??2018-06-02?19:33??kk\10.wav
?????文件???????64044??2018-06-02?19:33??kk\11.wav
?????文件???????64044??2018-06-02?19:33??kk\12.wav
?????文件???????64044??2018-06-02?19:33??kk\13.wav
?????文件???????64044??2018-06-02?19:33??kk\14.wav
?????文件???????64044??2018-06-02?19:33??kk\15.wav
?????文件???????64044??2018-06-02?19:33??kk\16.wav
?????文件???????64044??2018-06-02?19:33??kk\17.wav
?????文件???????64044??2018-06-02?19:34??kk\18.wav
?????文件???????64044??2018-06-02?19:34??kk\19.wav
?????文件???????64044??2018-06-02?19:33??kk\2.wav
?????文件???????64044??2018-06-02?19:34??kk\20.wav
?????文件???????64044??2018-06-02?19:33??kk\3.wav
?????文件???????64044??2018-06-02?19:33??kk\4.wav
?????文件???????64044??2018-06-02?19:33??kk\5.wav
?????文件???????64044??2018-06-02?19:33??kk\6.wav
............此處省略58個文件信息